Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 915,332 C→T 5.8% intergenic (‑475/+20) ybjE ← / ← aqpZ putative transporter/aquaporin Z

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913915,3320CT 5.8% 276.1 / 3.1 86intergenic (‑475/+20)ybjE/aqpZputative transporter/aquaporin Z
Reads supporting (aligned to +/- strand):  ref base C (46/35);  new base T (3/2);  total (49/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

CGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACCGCGGTGCTG  >  NC_000913/915190‑915481
                                                                                                                                               |                                                                                                                                                     
cGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAATGCCTGACATTACGCCAGGc                                                                                                                                                >  2:658590/1‑150 (MQ=255)
  ggggTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCt                                                                                                                                              >  1:1341208/1‑150 (MQ=255)
      ggCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGGCGTAATGTCAGGCCttttt                                                                                                                                          <  2:980120/149‑4 (MQ=37)
      ggCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTg                                                                                                                                          <  2:195648/149‑1 (MQ=255)
       gCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGc                                                                                                                                         <  1:1197653/150‑1 (MQ=255)
         aGGGGGCGGCGGCGATTGGCCGCCCCGTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATATCAACATCTT‑AAAAAAAGGCCTGGCGTAATGTCAGGCCttttttta                                                                                                                                       <  2:22594/150‑7 (MQ=37)
         aGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGt                                                                                                                                       >  1:1687008/1‑150 (MQ=255)
               cggcggcgATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaatta                                                                                                                                 >  1:1197012/1‑150 (MQ=255)
               cggcggcgATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCAAGAAAGGAGCGCAACAAGAAGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTCACATTACGCCAGGCCTTCTGGGttaatta                                                                                                                                 >  1:193785/1‑150 (MQ=255)
                ggcggcgATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGC‑TAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTtaattaat                                                                                                                               <  2:294370/150‑1 (MQ=255)
                ggcggcgATTGGCCGCCCCCTGCGAGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGAGGAATCAACATATAACCACATCTG‑AAAAAAAGTCATGCCATTACGCCAGGCCTTCTGCGttaattaa                                                                                                                                >  1:257803/1‑150 (MQ=38)
                 gcggcgATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTtaattaat                                                                                                                               >  2:852238/1‑150 (MQ=255)
                 gcggcgATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTtaattaat                                                                                                                               >  2:641158/1‑150 (MQ=255)
                  cggcgATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATc                                                                                                                              <  2:54558/150‑1 (MQ=255)
                  cggcgATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑‑AAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAAATAATCa                                                                                                                             >  2:254622/1‑150 (MQ=255)
                     cgATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTA‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACg                                                                                                                           >  1:490066/1‑150 (MQ=255)
                          ggCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCttt                                                                                                                       >  2:1616097/1‑150 (MQ=255)
                              cccccgggcgcgcTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcag                                                                                                                  <  2:432830/143‑1 (MQ=255)
                              ccccccGCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcag                                                                                                                  <  2:179293/149‑1 (MQ=255)
                                ccccTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcagca                                                                                                                >  2:1358775/1‑150 (MQ=255)
                                      cgcgcTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTcc                                                                                                          >  1:1110463/1‑150 (MQ=255)
                                        cgcTCCTTGCGCCAGTGGCAATATGTAGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCgg                                                                                                        <  1:283965/150‑1 (MQ=255)
                                              ttGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCCCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCAATCAGGGTCCGGTAAAta                                                                                                  >  2:865364/1‑149 (MQ=255)
                                               tGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCa                                                                                                 >  1:1171154/1‑150 (MQ=255)
                                                  gcCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGac                                                                                              >  1:446220/1‑150 (MQ=255)
                                                  gcCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGac                                                                                              <  2:361230/150‑1 (MQ=255)
                                                  gcCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGac                                                                                              <  2:1637674/150‑1 (MQ=255)
                                                  gcCAGTGGAAATATGTTGCTTAGCTAATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGATATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGac                                                                                              >  1:1531077/1‑150 (MQ=255)
                                                    cAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCCCAACAAGATGATGAATCTACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGATAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGacca                                                                                            >  1:1029041/1‑150 (MQ=255)
                                                       tGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCg                                                                                         >  2:624782/1‑150 (MQ=255)
                                                           aaTATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATaa                                                                                     >  2:308555/1‑150 (MQ=255)
                                                              atGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATgc                                                                                  <  1:85699/150‑1 (MQ=255)
                                                               tGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATgcc                                                                                 <  1:1651053/150‑1 (MQ=255)
                                                                gTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATgcc                                                                                 <  1:1616097/150‑1 (MQ=255)
                                                                gTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATgccg                                                                                <  1:1101816/150‑1 (MQ=255)
                                                                 ttGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATgccgc                                                                               <  1:658590/150‑1 (MQ=255)
                                                                    cTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGa                                                                            >  1:951925/1‑150 (MQ=255)
                                                                        gCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAAt                                                                        >  2:1535063/1‑150 (MQ=255)
                                                                        gCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAAt                                                                        >  2:227202/1‑150 (MQ=255)
                                                                          tCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTg                                                                      >  2:11607/1‑150 (MQ=255)
                                                                          tCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTg                                                                      >  2:498657/1‑150 (MQ=255)
                                                                            aTGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGc                                                                    >  1:497053/1‑150 (MQ=255)
                                                                             tGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGca                                                                   >  2:67054/1‑150 (MQ=255)
                                                                                   gAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCAccc                                                             >  1:158816/1‑150 (MQ=255)
                                                                                   gAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCAccc                                                             <  1:1358775/150‑1 (MQ=255)
                                                                                    aGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCAcccc                                                            >  1:111599/1‑149 (MQ=255)
                                                                                    aGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCa                                                            <  2:1135962/150‑1 (MQ=255)
                                                                                     gcgcAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCag                                                           >  2:887799/1‑150 (MQ=255)
                                                                                       gcAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCAGACAATACGCCAGGCATTCTGCGTTAATTAATCACGCTTTTCAAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCagaa                                                         >  1:238200/1‑150 (MQ=255)
                                                                                          acaaGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAgaagaa                                                      <  1:1676079/150‑1 (MQ=255)
                                                                                            aaGATGATGAATCTACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAAcc                                                    <  2:1029041/150‑1 (MQ=255)
                                                                                                tgatgaATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGTCACCACCCAGAAGAACCACAg                                                >  2:329622/1‑150 (MQ=255)
                                                                                                 gatgaATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAgt                                               >  1:834424/1‑150 (MQ=255)
                                                                                                   tgaATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATgccgcc                                                                              <  2:678617/117‑1 (MQ=255)
                                                                                                   tgaATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATgccgcc                                                                              >  1:678617/1‑117 (MQ=255)
                                                                                                   tgaATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAgttg                                             >  1:1646502/1‑150 (MQ=255)
                                                                                                   tgaATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAgttg                                             >  1:1384210/1‑150 (MQ=255)
                                                                                                    gaATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCTGCAGGGTCCTGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAgttgt                                            <  2:362428/150‑1 (MQ=255)
                                                                                                    gaATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAgttgt                                            <  1:887799/150‑1 (MQ=255)
                                                                                                     aaTCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAgttgtt                                           <  1:641158/150‑1 (MQ=255)
                                                                                                      aTCAACATTTAACAACATCTT‑AAAAATACGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTc                                          <  2:126460/150‑1 (MQ=255)
                                                                                                      aTCAACATATAACAACATCTT‑AAAAAAAGGCCTGATATTACGCCAGGCCTTCTGCGTTAATTACTCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTc                                          <  2:1531077/150‑1 (MQ=255)
                                                                                                      aTCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTc                                          <  2:1564034/150‑1 (MQ=255)
                                                                                                           cATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATg                                     <  1:49636/150‑1 (MQ=255)
                                                                                                                aacaacATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAg                                <  1:1535063/150‑1 (MQ=255)
                                                                                                                aacaacATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAg                                >  2:1328790/1‑150 (MQ=255)
                                                                                                                 acaacaTCTT‑AAAAAAAGGCCTGACATTATGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAgc                               >  2:111155/1‑150 (MQ=255)
                                                                                                                 acaacaTCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAgc                               <  1:852238/150‑1 (MQ=255)
                                                                                                                 acaacaTCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAgc                               <  1:269549/150‑1 (MQ=255)
                                                                                                                 acaacaTCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAgc                               <  1:259205/150‑1 (MQ=255)
                                                                                                                  caacaTCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAgcc                              <  2:1341208/150‑1 (MQ=255)
                                                                                                                  aaacaTCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATAGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAgcc                              <  2:951925/149‑1 (MQ=255)
                                                                                                                   aacaTCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAgccg                             >  1:1611429/1‑150 (MQ=255)
                                                                                                                    acaTCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAgccgc                            >  1:1149605/1‑150 (MQ=255)
                                                                                                                    acaTCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAgccgc                            >  2:86004/1‑150 (MQ=255)
                                                                                                                     caTCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAgccgcc                           >  2:932527/1‑150 (MQ=255)
                                                                                                                       tCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCt                         >  1:1680108/1‑150 (MQ=255)
                                                                                                                       tCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAAATAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACAACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGcccg                         >  1:206383/1‑149 (MQ=255)
                                                                                                                        cTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTg                        >  2:1216500/1‑150 (MQ=255)
                                                                                                                        cTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTg                        >  2:1393534/1‑150 (MQ=255)
                                                                                                                            aaaaaaaGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGTTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGaa                     <  1:1089004/150‑1 (MQ=255)
                                                                                                                            aaaaaaaGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGaa                     <  1:308555/150‑1 (MQ=255)
                                                                                                                              cagaaGGCCTGGCGTAATGTCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGACCCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGa                   >  2:487784/4‑150 (MQ=255)
                                                                                                                                gaaGGCCTGGCGTAATGTCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATa                 >  1:726268/2‑150 (MQ=255)
                                                                                                                                aaaGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATa                 <  2:1149605/150‑1 (MQ=255)
                                                                                                                                 aaGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCCCCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATac                >  2:754246/1‑149 (MQ=255)
                                                                                                                                               cGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACCGCGGtgctg  <  2:446220/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
CGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTT‑AAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACCGCGGTGCTG  >  NC_000913/915190‑915481

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑