Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,908,507 T→A 6.7% intergenic (‑141/+42) pstB ← / ← pstA phosphate ABC transporter ATPase/phosphate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,908,5070TA 6.7% 297.8 / 7.5 116intergenic (‑141/+42)pstB/pstAphosphate ABC transporter ATPase/phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (58/49);  new base A (4/4);  total (63/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.60e-01

Only 100 of 119 total aligned reads displayed.

ATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCG  >  NC_000913/3908359‑3908655
                                                                                                                                                    |                                                                                                                                                    
aTACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTg                                                                                                                                                     >  2:1421592/1‑150 (MQ=255)
  aCTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTggg                                                                                                                                                   >  1:1356801/1‑150 (MQ=255)
   cTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTgggg                                                                                                                                                  >  1:228680/1‑150 (MQ=255)
     cATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGgag                                                                                                                                                >  2:1174027/1‑150 (MQ=255)
        tGCACTCTTTTCCTCATCGTGCCTGATGCACTTCGCTCATCACGCTAACCGTAGGCCGGATCAGGCGATCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCTCCCCAGAGGGGCGAgg                                                                                                                                             <  2:698723/150‑1 (MQ=255)
        tGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCTCCCCAGAGGGGCGAgg                                                                                                                                             <  1:1628789/150‑1 (MQ=39)
         ggattCTTTTCCTCATCGTGCCCGTTGCACCTCGCTTATCACGCCTACCGTAGGCCGGAGCAGGCGTGCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGGCAAAAAGGCGCTCGGTCGGTCCCCTCTCCCCAGAGGGGCGAggg                                                                                                                                            <  1:658442/146‑1 (MQ=255)
         gCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAggg                                                                                                                                            <  1:365056/150‑1 (MQ=255)
         gCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAggg                                                                                                                                            <  1:1098380/150‑1 (MQ=255)
            aTCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCTCCCCAGAGGGGCGAggggac                                                                                                                                         <  1:1219623/150‑4 (MQ=39)
            aTCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTa                                                                                                                                         <  2:1335307/150‑1 (MQ=255)
             tCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCTCCCCTCTGGGGAGAGGGTTAg                                                                                                                                        >  2:1395423/1‑150 (MQ=255)
               ttttCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTgggg                                                                                                                                                  >  1:184204/1‑138 (MQ=255)
               ttttCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTgggg                                                                                                                                                  <  2:184204/138‑1 (MQ=255)
               ttttCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAggg                                                                                                                                      >  1:1445780/1‑150 (MQ=255)
               ttttCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAggg                                                                                                                                      <  2:712462/150‑1 (MQ=255)
                 ttCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTg                                                                                                                                    <  2:1683471/150‑1 (MQ=255)
                  tCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGa                                                                                                                                   >  1:673759/1‑150 (MQ=255)
                   ccTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGTGTTAGGGTGAg                                                                                                                                  <  1:1335290/150‑1 (MQ=255)
                    cTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg                                                                                                                                 <  1:625067/150‑1 (MQ=255)
                     tCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAggg                                                                                                                                <  2:1263621/150‑1 (MQ=255)
                      cATCGTGCCTGATGCACTTCGCTACTCCGGCTTACCGTCGGCCGGATCCGGCGGGCACGCCGCAGCCGGCAAAGTGGTCACCGTGTACCGAGGCAACACGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgggg                                                                                                                               <  1:920179/150‑1 (MQ=25)
                       aTCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTACCCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGCGGGTTAGGGTGAGGGGa                                                                                                                              >  2:1054542/1‑150 (MQ=255)
                        tCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaa                                                                                                                             >  2:745105/1‑150 (MQ=255)
                         cGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaa                                                                                                                            >  1:1382884/1‑150 (MQ=255)
                           ggcgtGATGCACTTCGCTAATCAGGCTTACCGTAGGCTGGATCAGGCGTTCACGCCGAATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGGCCCCTCGCCCCTCTGGGGAGAGGGTTAGTGTGTGGGGAAAAc                                                                                                                          <  2:1549894/146‑1 (MQ=39)
                            gCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTTTACAGAGTCAAAAAGGCGCTCGGTCGGGCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAcc                                                                                                                         <  1:826687/150‑1 (MQ=255)
                             ccTGATTCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGGGTTAACCGGGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCGCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCg                                                                                                                        <  1:1606675/150‑1 (MQ=255)
                               tGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTAAACCGTGTACAGAGTCAAAACGGCGCGCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCgtg                                                                                                                      <  2:552224/150‑1 (MQ=255)
                                gATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAACGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCgtgt                                                                                                                     <  1:722977/150‑1 (MQ=255)
                                                cAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCgt                                                                                                     >  1:1380723/1‑150 (MQ=255)
                                                    cTTACCGTAGGCCGGATCAGGCGTTCACCCCGCATCCGGAAATGTGTTAACCGTGTACAGAGTCAAAACGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGTGTGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGttt                                                                                                 <  2:1445780/150‑1 (MQ=39)
                                                      tACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAACCCGTGTCAGCAAGATCAACCGTGTTTAt                                                                                               >  2:897344/1‑150 (MQ=255)
                                                           tAGGCCGGAGCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttcttc                                                                                          <  1:25798/150‑1 (MQ=255)
                                                            ttgCCGGATGCGGCGTGAACGCCTGATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCg                                                                                         >  2:893265/3‑150 (MQ=38)
                                                            aGGCGGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTAAAGATTCAAAATGGGGCGCGGGCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTTGGGTGTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCg                                                                                         <  2:1341880/150‑1 (MQ=38)
                                                            aGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCg                                                                                         <  1:632373/150‑1 (MQ=255)
                                                            aGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCg                                                                                         >  2:922094/1‑150 (MQ=255)
                                                             ggCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGTAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGc                                                                                        <  2:1382884/150‑1 (MQ=255)
                                                             ggCCGGAGCAGGGGTGCACGCCGCATCGGGAAAAGTGTTAACCGTGTACAGAGTCAAAACGGCGCTGGGGCGGGCCCCTCGCCCCTCTGGGGTTAGGGTTAGGGGGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGc                                                                                        <  2:1253564/150‑1 (MQ=255)
                                                              gCCGGATGCGGCGTGAACGCCTGATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGTGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATAACAACCGTGTTTATTCTTCGCa                                                                                       >  2:319320/1‑150 (MQ=38)
                                                               ccGGATGCGGCGTGAACGCCTGATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTTGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaa                                                                                      >  1:747257/1‑150 (MQ=38)
                                                                cGGATCAGGCGTTCACGCCGCATCCGGAAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaaa                                                                                     <  2:388198/150‑1 (MQ=255)
                                                                  gATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGATGTGAAGACCGTGTCAGCAATATCAACAGGGTATATGCTTCGcctaac                                                                                   >  1:94778/1‑145 (MQ=38)
                                                                    tCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAaaca                                                                                 <  2:956856/150‑1 (MQ=255)
                                                                     cAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAaacaa                                                                                >  2:740336/1‑150 (MQ=255)
                                                                     cAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAaacaa                                                                                >  1:1294678/1‑150 (MQ=255)
                                                                     cAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGATTATTCTTCGCAAAaacac                                                                                >  2:934093/1‑149 (MQ=255)
                                                                       ggCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCATCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcg                                                                              >  2:251087/1‑150 (MQ=255)
                                                                       ggCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCATCGCCCCTCTGGGGAGAGGGTGAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAacaacacag                                                                              >  2:282038/1‑147 (MQ=255)
                                                                         cGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGATGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgcg                                                                            <  1:113794/150‑1 (MQ=255)
                                                                             cACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAACACCGTGTCTGCACTATCAACCGTGTGTATTCTTCGCAAACACAACGCGCGcc                                                                        >  2:842978/1‑150 (MQ=38)
                                                                             cACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGcc                                                                        >  2:191483/1‑150 (MQ=255)
                                                                               cgccgcATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAACCAACGCGCGCCAg                                                                      >  1:842014/1‑150 (MQ=255)
                                                                                    cATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGa                                                                     >  2:268823/1‑146 (MQ=255)
                                                                                      tCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCGGGGGAGAGGGTTAGGGTGAGGGGAAATCCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTc                                                               >  1:320662/1‑150 (MQ=255)
                                                                                       ccGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTca                                                              <  1:1031219/150‑1 (MQ=255)
                                                                                        cGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcag                                                             >  2:259544/1‑150 (MQ=255)
                                                                                        cGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcag                                                             >  1:1041895/1‑150 (MQ=255)
                                                                                          gCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATAAACAGTGTTTATTCTTAGCAAAAACAACGCTCGCCAGAATGTTcacca                                                           >  2:1460921/1‑150 (MQ=38)
                                                                                            aaaGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagcagc                                                         <  1:1174027/150‑1 (MQ=255)
                                                                                             aaGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagcagca                                                        <  2:228680/150‑1 (MQ=255)
                                                                                                       cGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAg                                              >  1:1151560/1‑150 (MQ=255)
                                                                                                           tACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTa                                          >  1:419125/1‑150 (MQ=255)
                                                                                                            aCAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTaa                                         >  1:1288208/1‑150 (MQ=255)
                                                                                                                agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGAt                                     >  2:393387/1‑150 (MQ=255)
                                                                                                                 gTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATc                                    >  2:1496969/1‑150 (MQ=255)
                                                                                                                 gTCAAAAAGGCGCTCGGTCGGGCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATc                                    <  2:427403/150‑1 (MQ=255)
                                                                                                                      aaaGGCGCTCGGTCGGTCCCCTCGCCCCGCTGGGGAGAGGGTAAGGGTGAGGGGACAACCGTGTCAGCAATATCAACCGTGTTTAGTATTCGCAAAACAAACGCGCGCCAGAGGGTGCAGCGGCAGTCCGCACAGGGTAATGATCaagtc                               >  1:292415/1‑147 (MQ=255)
                                                                                                                           cgcTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAACACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACACCGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCgg                          >  2:340108/1‑150 (MQ=255)
                                                                                                                           cgcTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCgg                          <  1:740336/150‑1 (MQ=255)
                                                                                                                            gcTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGc                         >  1:1638702/1‑150 (MQ=255)
                                                                                                                             cgcggtcggcCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGcc                        <  2:320662/148‑1 (MQ=255)
                                                                                                                             cTCGGTCGGTCCCCTCGCCCCTCTGGGGAGTGGGTTATGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTATTGATCATTAccctgtgc                        <  2:457414/150‑6 (MQ=39)
                                                                                                                             cTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGcc                        >  2:1413185/1‑150 (MQ=255)
                                                                                                                                 accctaacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAAGATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAgg                    >  1:491117/8‑150 (MQ=255)
                                                                                                                                  ccctaacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGACAACCGGGTCAGCAAGATCAACCGTGTTTATTCTTCGCAAACGCAGCGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTACTGATCTTTCCCCCGGCCCAGGc                   >  1:857631/7‑150 (MQ=255)
                                                                                                                                   cggtcCCCTCGCCCCTCTGGGGTGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCCACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGAACGCACAGGGTAATGATCAATACCCCGGCCCAGGcc                  >  1:1335820/1‑150 (MQ=255)
                                                                                                                                   cggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATAcccc                            >  1:171720/1‑140 (MQ=255)
                                                                                                                                   cggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATAcccc                            >  1:166129/1‑140 (MQ=255)
                                                                                                                                   cggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATAcccc                            <  2:171720/140‑1 (MQ=255)
                                                                                                                                   cggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGcc                  <  2:1638702/150‑1 (MQ=255)
                                                                                                                                   cctaacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTCCGCACAGGGTAATGATCACTACCCCGGCCCAGGcc                  >  1:537807/6‑150 (MQ=255)
                                                                                                                                    ggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAAt                                 >  1:30849/1‑134 (MQ=255)
                                                                                                                                    ggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAAt                                 <  2:30849/134‑1 (MQ=255)
                                                                                                                                     gtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCaa                >  2:580214/1‑150 (MQ=255)
                                                                                                                                       ccccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAAtt              >  2:1090917/1‑150 (MQ=255)
                                                                                                                                        cccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAATCTACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtg             >  1:266578/1‑150 (MQ=255)
                                                                                                                                        cccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGTGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtg             >  1:1574229/1‑150 (MQ=255)
                                                                                                                                          cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgct           >  2:1484454/1‑150 (MQ=255)
                                                                                                                                           tCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgctg          >  2:1571898/1‑150 (MQ=255)
                                                                                                                                             gCCCCTGTGGGGAGAGGGTTAGGGTTAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGcc        <  1:1054542/150‑1 (MQ=255)
                                                                                                                                             gCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGcc        <  2:1294678/150‑1 (MQ=255)
                                                                                                                                             gCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGGCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGcc        <  1:893265/150‑1 (MQ=255)
                                                                                                                                              ccccTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCa       >  1:907671/1‑150 (MQ=255)
                                                                                                                                              ccccTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCa       >  2:24930/1‑150 (MQ=255)
                                                                                                                                               cccTCTGGGGAGAGGTTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACATGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCAt      <  1:970235/150‑1 (MQ=255)
                                                                                                                                                 ctctGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTc    >  1:1262192/1‑150 (MQ=255)
                                                                                                                                                  tctGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTcc   <  1:611491/150‑1 (MQ=255)
                                                                                                                                                   ctGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCg  <  2:57949/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCG  >  NC_000913/3908359‑3908655

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