Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,908,512 A→C 6.1% intergenic (‑146/+37) 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,5120AC 6.1% 298.7 / 4.0 114intergenic (‑146/+37)pstB/pstAphosphate ABC transporter ATPase/phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (64/40);  new base C (4/4);  total (69/45)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 116 total aligned reads displayed.

CATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAA  >  NC_000913/3908364‑3908659
                                                                                                                                                    |                                                                                                                                                   
cATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGgag                                                                                                                                                    <  2:1158053/150‑1 (MQ=255)
 aTTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGCTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                                                   <  1:1350063/150‑1 (MQ=255)
   tGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCTCCCCAGAGGGGCGAgg                                                                                                                                                 <  1:613386/150‑1 (MQ=39)
   tGCAATCTTTTCCTCATCGTGCCTGATGCACTGCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCTCCCCAGAGGGGCGAgg                                                                                                                                                 <  2:538045/150‑1 (MQ=39)
    gCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAggg                                                                                                                                                <  2:1542776/150‑1 (MQ=255)
    gCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAggg                                                                                                                                                <  1:1240430/150‑1 (MQ=255)
     cAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCGTCCGGGCAAGTGGTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCGTCTGGGGAGAgggg                                                                                                                                               >  2:385621/1‑149 (MQ=255)
     cAATCTTTTCCTCATCCTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGGCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGt                                                                                                                                               <  2:818481/150‑1 (MQ=255)
      aaTCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGtt                                                                                                                                              >  1:369716/1‑150 (MQ=255)
       aTCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTa                                                                                                                                             >  1:1742147/1‑150 (MQ=255)
       aTCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTa                                                                                                                                             >  2:1190072/1‑150 (MQ=255)
        tCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAg                                                                                                                                            <  1:615726/150‑1 (MQ=255)
         cTTTTCCTCATCGTGCCTGCTGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCACAGGGGAGAGGGTTAgg                                                                                                                                           >  1:220926/1‑150 (MQ=255)
          ttttCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAggg                                                                                                                                          <  2:1848145/150‑1 (MQ=255)
          ttttCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCAGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAggg                                                                                                                                          <  1:715055/150‑1 (MQ=255)
              cgtCATCGTGCCTGCGGCACTTCGCTAATGAGGCTTACCGTAGGCCGGATCAGGGGTTCAGGCCGCATGCGGCAAAGTGGTACCCGTGTATAGTTTCTAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAg                                                                                                                                      <  1:264920/148‑1 (MQ=37)
               cTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg                                                                                                                                     <  1:1688750/150‑1 (MQ=255)
                tCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAACAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAggg                                                                                                                                    <  1:134347/150‑1 (MQ=255)
                tCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAggg                                                                                                                                    <  1:860928/150‑1 (MQ=255)
                tCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAggg                                                                                                                                    <  1:1995262/150‑1 (MQ=255)
                 cATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgggg                                                                                                                                   >  2:8459/1‑150 (MQ=255)
                   tCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaa                                                                                                                                 >  2:1466429/1‑150 (MQ=255)
                      tGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGGGTTAACCGGGTACAGAGGCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAc                                                                                                                              <  1:1153785/150‑1 (MQ=255)
                        ccTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAATCCg                                                                                                                            >  2:627030/1‑150 (MQ=255)
                            ttgccgttCGCTCATCTGCCTTACCGTAGGCCGGAGCAGGCGTTCACGCCGCTTCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAATGGCGCTGGGGCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTc                                                                                                                        <  1:1845945/144‑1 (MQ=37)
                               cACTTCGCTAATCAGGCTTACCGTAGGCCGGCTCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGCGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGc                                                                                                                     >  1:1460317/1‑150 (MQ=255)
                               cACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGTTTAGGGTGAGGGGAAAACCGTGTCAGc                                                                                                                     >  2:1561360/1‑150 (MQ=255)
                                   tCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCACATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCTCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGGGTCAGCAata                                                                                                                 >  2:1934797/1‑150 (MQ=255)
                                     gCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACCCCGCATCCGGCAAAGTTTTAACCGTGTACAGAGTCAAAAAGGCGCGTGGTCGGTCCCGTCGCCCCTCTGGGGAGAGGGTTAGGGGGAGGGGAAAACCGTGTCAGCAATATc                                                                                                               <  2:522962/150‑1 (MQ=39)
                                      cTAATCAGGATTACCGTAGGCCGTATCAGGCGTTCAGGCCGCATCCGGAAAAGTGTTAACCGTGTACTGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGTGGGTTTGGGTGAGGGGAAAACCGTGTCAGCAATATCa                                                                                                              <  2:1788917/150‑1 (MQ=38)
                                        aaTCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAACAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAAc                                                                                                            >  1:1702506/1‑150 (MQ=255)
                                           cAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCgt                                                                                                         >  2:1919134/1‑150 (MQ=255)
                                            aGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCgtg                                                                                                        <  1:8459/150‑1 (MQ=255)
                                             ggCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGGACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGGTAGGGTGAGGGGATAACCGGGTCAGCAATATCTACCgtgt                                                                                                       >  2:384764/1‑150 (MQ=39)
                                              gCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGTGAGGGTTAGGGTGAGGGGAAAACCGAGTCAGCAATATCAACCGTGtt                                                                                                      >  2:812610/1‑150 (MQ=255)
                                               cTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCTACCGTGgtt                                                                                                     >  1:822158/1‑150 (MQ=255)
                                                ttACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTTGGGTGAGGGGAAAACAGTGTCAGCATTATCAAGCGGGGTTa                                                                                                    >  2:685456/1‑150 (MQ=38)
                                                        ggCCGGATCAGGGGTTCACGCCGCATCCGGCAAAGTGTTAACCTTGTACAGAGGCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGc                                                                                            <  2:357011/150‑1 (MQ=255)
                                                         gCCGGATGCGGCGTGAACGCCTGATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGACAGCAATAACAACCGTGTTTATTCTTCGCa                                                                                           >  2:1273518/1‑150 (MQ=38)
                                                         gCCGGATCCGGCGGTCACGCCGCATCCGGCAAAGTGTTAACCGTGAACAGAGGCAAACAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCa                                                                                           <  1:1169777/150‑1 (MQ=38)
                                                          ccGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTGAGGGTGAGGGGAAACCCGTGTCAGCAATATCAACCGTGTTTATTCTTCGcca                                                                                          >  2:1856650/1‑148 (MQ=255)
                                                             gATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAACCCGTGTCAGCAATATCACCCGTGTTTATTCTTCGCaaaaa                                                                                       >  1:581532/1‑150 (MQ=255)
                                                                  ggCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcg                                                                                  >  2:628037/1‑150 (MQ=255)
                                                                  ggCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCGGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCACCCGTGTTTATTCTTCGCAAAACCACcg                                                                                  >  2:894484/1‑150 (MQ=255)
                                                                    cGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCACTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgcg                                                                                >  1:1591120/1‑150 (MQ=255)
                                                                           gccgcATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGa                                                                         >  2:468911/1‑150 (MQ=255)
                                                                               cATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGt                                                                     >  1:20029/1‑150 (MQ=255)
                                                                                 tCCGGCAAAGTGTTAAGCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTTGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTc                                                                   <  2:786629/150‑1 (MQ=255)
                                                                                 tCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTc                                                                   <  2:1068931/150‑1 (MQ=255)
                                                                                       aaaGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagcagc                                                             <  2:179495/150‑1 (MQ=255)
                                                                                              tAACGGTGTACTGAGTCAAAAAGGGGGGCGGTCGTTCGCCTCGCACCTCTGGGGCTTGGTTTAGGTTGAGGGGAAAACCGCGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAATCAACGCGCGCCAGAATGTTCAGCAGCAGTACGc                                                      <  1:1273518/150‑1 (MQ=255)
                                                                                              tAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGc                                                      <  1:222109/150‑1 (MQ=255)
                                                                                                 ccGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGcaca                                                   >  2:943546/1‑150 (MQ=255)
                                                                                                 ccGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGCTTAGGGTCAGGGGAAAACCGAGTAAGCAATATCAACCGTGTTTATTCTACGCAAAAACAACGCGCGCCAGAATGTTAAGCAGCAGTACGcaca                                                   >  1:867386/1‑150 (MQ=39)
                                                                                                   gtttACATAGTCAAAAAGGCGCGCGGGCGGTCCCCTCGCCCCTGTGGGTAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAgg                                                 <  1:1466429/147‑1 (MQ=39)
                                                                                                    gttACAGCGTCAAAAAGGCGCTCGGGCGGGCCCCGCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAggg                                                <  1:1571832/148‑1 (MQ=255)
                                                                                                       aCAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTTTTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTaa                                             >  1:379727/1‑150 (MQ=255)
                                                                                                       aCAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTaa                                             >  1:456741/1‑150 (MQ=255)
                                                                                                        cAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAAt                                            >  1:1659342/1‑150 (MQ=255)
                                                                                                         agagTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCGCTTGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATg                                           <  2:1525190/150‑1 (MQ=255)
                                                                                                          gagTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGa                                          >  2:194117/1‑150 (MQ=255)
                                                                                                           agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGAt                                         <  2:447542/150‑1 (MQ=255)
                                                                                                            gccAAAAAGGCGCGCGGTCGGCCCCCTCGCCCCTGTGGGGTGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATc                                        <  1:841226/148‑1 (MQ=255)
                                                                                                             tCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCa                                       <  1:1190072/150‑1 (MQ=255)
                                                                                                              cAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCTAAACCACCGCGCGCCAGAATGTTCCGCAGCAGTACGCCCAGGGTAATGATCaa                                      >  1:733148/1‑150 (MQ=38)
                                                                                                              cAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCaa                                      >  2:1952609/1‑150 (MQ=255)
                                                                                                              cAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGGAATGCTCaa                                      >  1:1110432/1‑150 (MQ=255)
                                                                                                               aaaaaGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCACAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACGTGGTAATGATCAAt                                     >  2:1215472/1‑150 (MQ=255)
                                                                                                                    ggCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTCCGCACAGGGTAATGATCAATAcccc                                >  2:569648/1‑150 (MQ=255)
                                                                                                                    ggCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATAcccc                                <  1:348225/150‑1 (MQ=255)
                                                                                                                     gcgcTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCg                               <  1:468911/150‑1 (MQ=255)
                                                                                                                      cgcTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCgg                              <  1:581537/150‑1 (MQ=255)
                                                                                                                      cgcTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCgg                              >  1:878850/1‑150 (MQ=255)
                                                                                                                          cggtcggtCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCa                          >  2:480949/1‑150 (MQ=255)
                                                                                                                          cggtcggtCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGGGTTTAATCTTCGCAAAAACACAGCGCGCCAGAATGGTCAGCAGCAGTACGCACAGGGAAATGATCAATACCCCGGCCCa                          >  2:322554/1‑150 (MQ=38)
                                                                                                                             tcggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGc                       >  1:1481916/1‑150 (MQ=255)
                                                                                                                             gcggccCCCTCGCCCCTCTGGGGCGCGGGTTAGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGc                       <  1:1260711/149‑1 (MQ=255)
                                                                                                                              cggtcCCCTCGCCCCTCTGGGTAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGcc                      <  2:369716/150‑1 (MQ=255)
                                                                                                                              cggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGcc                      <  2:1725767/150‑1 (MQ=255)
                                                                                                                              cggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGcc                      >  1:943070/1‑150 (MQ=255)
                                                                                                                               ggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCa                     >  2:720460/1‑150 (MQ=255)
                                                                                                                               ggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCa                     >  1:1756886/1‑150 (MQ=255)
                                                                                                                               ctaacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGGGTTTATTCTTCGCAAAATCAACGCGCGCCAGACTGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGcga                     >  1:1914016/5‑148 (MQ=255)
                                                                                                                                gtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGccca                    >  1:157531/1‑148 (MQ=255)
                                                                                                                                 tcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAAt                   >  2:708073/1‑150 (MQ=255)
                                                                                                                                   cccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGAAATAATGTTCAGCATCAGTACGCAAAGGGGAAGGATCAATACCCCGGGCGAGGCCAATtg                 >  2:321939/1‑150 (MQ=255)
                                                                                                                                   cccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAACCGTGTCACCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTCCGCACAGGGTAATGCTCAATACCCCGGCCCAGGCCAATtg                 >  1:1471556/1‑150 (MQ=255)
                                                                                                                                   cccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtg                 >  1:418450/1‑150 (MQ=255)
                                                                                                                                    ccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCCGCAATATCCCCCGTGTTTATTCTTCGCAAAAACAACGCTCGCCACCATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGTCCCAGGCCACTtgc                >  1:978511/1‑150 (MQ=38)
                                                                                                                                    ccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgc                >  2:1659197/1‑150 (MQ=255)
                                                                                                                                    ccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgc                >  1:286645/1‑150 (MQ=255)
                                                                                                                                    ccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgc                >  1:94328/1‑150 (MQ=255)
                                                                                                                                       cGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgctgc             >  1:661391/1‑150 (MQ=255)
                                                                                                                                       cGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGGCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgctgc             <  2:70182/150‑1 (MQ=255)
                                                                                                                                          cccTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCAt          >  2:1864190/1‑150 (MQ=255)
                                                                                                                                           ccTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCAtt         >  2:1692025/1‑150 (MQ=255)
                                                                                                                                           ccTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCAtt         >  1:1245010/1‑150 (MQ=255)
                                                                                                                                               tGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGc     >  1:1047388/1‑150 (MQ=255)
                                                                                                                                               tGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAAACGTGTTTATTCATCGCAAAAACAACGCGCGCCAGAAAGTTCAGCAGCAGAACGCACAGGGTAATGATCAATACCCCGGCCCATGCCAATTGCTGCCATTCCGc     >  1:1335651/1‑150 (MQ=255)
                                                                                                                                                  ggAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCCGCAGCAGTACGCACCGGGTAATGATCAAACCCCCGGCCCAGGCCAATTGCTGCCATTCCGCaaa  >  2:3131/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAA  >  NC_000913/3908364‑3908659

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: