Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,908,500 G→T 7.1% intergenic (‑134/+49) 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,5000GT 7.1% 259.1 / 9.3 99intergenic (‑134/+49)pstB/pstAphosphate ABC transporter ATPase/phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (41/49);  new base T (3/4);  total (44/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.93e-01

TCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGC  >  NC_000913/3908353‑3908648
                                                                                                                                                   |                                                                                                                                                    
tCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGcc                                                                                                                                                    <  2:1244547/150‑1 (MQ=255)
 cAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAACGCCGCTCGGTCGGTCCCCGCGccc                                                                                                                                                   >  2:1056887/1‑150 (MQ=255)
  aaCCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGcccc                                                                                                                                                  >  2:1337682/1‑150 (MQ=255)
  aaCCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGcccc                                                                                                                                                  >  1:256818/1‑150 (MQ=255)
       tACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTgg                                                                                                                                             <  2:1478747/150‑1 (MQ=255)
       tACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTgg                                                                                                                                             <  1:1328513/150‑1 (MQ=255)
        aCTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTggg                                                                                                                                            <  2:1410223/150‑1 (MQ=255)
        aCTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTggg                                                                                                                                            <  2:127869/150‑1 (MQ=255)
         cTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAAACAGGCTTACCGTAGGCCGGATCAGGCGATCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCACGGTCGGTCCCCTCGCCCCACTgggg                                                                                                                                           >  1:388772/1‑150 (MQ=255)
         cTCATTGCAATCTTTTCCGCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGAGCAGGCGTGCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTgggg                                                                                                                                           <  1:154471/150‑1 (MQ=255)
           cATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTAGGTCCACTCGCCCCTCTGGGgag                                                                                                                                         >  2:412568/1‑150 (MQ=255)
              tGCATTCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCTCCCCAGAGGGGCGAgg                                                                                                                                      <  2:1184828/150‑1 (MQ=39)
              tGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTACTCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAACGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                      <  2:111049/150‑1 (MQ=255)
              tGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                      <  2:1349317/150‑1 (MQ=255)
              tGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCACTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                      <  2:385218/150‑1 (MQ=255)
               gCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAggg                                                                                                                                     >  2:650074/1‑150 (MQ=255)
                cAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCTCCCCAGAGGGGCGAgggg                                                                                                                                    <  1:571514/150‑2 (MQ=39)
                cAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGGCGGTCCCCTCTCCCCAGAGGGGCGAgggg                                                                                                                                    <  1:1400068/150‑2 (MQ=39)
                 aaTCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGTATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGtt                                                                                                                                   <  2:871123/150‑1 (MQ=255)
                 aaTCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGtt                                                                                                                                   <  2:404832/150‑1 (MQ=255)
                   tCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACCGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAg                                                                                                                                 <  2:775096/150‑1 (MQ=255)
                    cTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAgg                                                                                                                                <  1:1275190/150‑1 (MQ=255)
                         cctcctCGTGCCTGATGCACTTCGCTAATCAGGCTTACGGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTCACCGTGTTCAGAGTCTAATAGGGGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAg                                                                                                                           <  1:14994/150‑1 (MQ=38)
                         ccTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAg                                                                                                                           >  1:1021506/1‑150 (MQ=255)
                            cATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgggg                                                                                                                        >  1:34237/1‑150 (MQ=255)
                            cATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgggg                                                                                                                        >  1:884872/1‑150 (MQ=255)
                             aTCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTTGGGTGAGGGGa                                                                                                                       >  2:792360/1‑150 (MQ=255)
                              tCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaa                                                                                                                      >  1:405238/1‑150 (MQ=255)
                                       aTGCACTTCGCTAATCAGGCTTACCGTAGGCCGGAGCCGGCGTTCACGCCGCATCGGGCAATGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGTTCCCCTCGCCCCTCTGGGGTGAGGGTTAGGGTGAGGGGAAAACCGTGTc                                                                                                             <  1:1692992/150‑1 (MQ=39)
                                        tGCACTTCGCTAATCAGGCTTACCGTCGGCCGGATCCGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCGCGGTCGGGCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGTGGGGAAAACCGTGTCa                                                                                                            <  1:437793/150‑1 (MQ=39)
                                         gCACTTCGCTAATCAGGCGTACCTTAGGCGGGATCTGGCGTTCAGGCCGCATCCGGCAAAGTGTTAACCGTGTTAAGTGTCAAAAAGGCGCGCGGGCGGCCCCCTCGCCCCTCTGGGGTTTGGGTTAGGGTGAGGGTAAAACCGTGTCAg                                                                                                           <  1:1068570/150‑1 (MQ=255)
                                             ctcccTAATCCGGCTAACCGTAGGCCGGCTCCGGCGTTCACGCCGCACCCGGCAATGTGTTAACCGTGTACAGATTCAAAAAGGCGCTCGGTCGGGCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAat                                                                                                       <  1:174862/146‑1 (MQ=37)
                                              tCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGGTAGGGGGAGGGAAAACCGTGTCAGCAatat                                                                                                     >  2:399112/1‑150 (MQ=255)
                                              tCGCAAATCAGGCTTCCGGTAGGCCGGCTCTGGCGTTCACGCCGCATCCGGAAAAGTGTTAACCGTGTCCAGCGGCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGTGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAata                                                                                                      <  2:1253523/150‑1 (MQ=37)
                                               cGCAAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTACCCGTGTACAGAGTCAAAAAGCCGCTCGGGCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAACAatat                                                                                                     <  1:916789/150‑1 (MQ=38)
                                                    aTCTGGCTTAGCGTAGGAGGGTTCAGGCGTTCACGCGGCTTCCGGCAAAGTGTTAACCGTGTACAGATTAAAAAAGGCGGTCGGTGGGTCCCCTCCCCCATCTGGGGAGTGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAAcc                                                                                                <  2:1423436/150‑1 (MQ=255)
                                                    aTCAGGCTTACCTTAGGCGGGATCAGGCGTTCACGCCGCATCCGGCAATGTGTTAACCGTGTACAGACTCAAAAAGGCGCTGGGTCGTTCCCCTCGCCCCTCTGGGGAGTGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAAcc                                                                                                <  1:1614529/150‑1 (MQ=38)
                                                     tccggCTTACCGTTGGCCGGATCTGCCGTTCCCGCCGCATCCGGCAAAGTGTTAACCGTGAACAGAGTCAAAATGGCGCTCGGGCGGTCCCCTCGCCCCTCGGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCg                                                                                               <  2:83845/147‑1 (MQ=38)
                                                        ggCTTACCGTAGGCAGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCCAGTAAAGAGTCAAAAAGGCGCTCGGTCAGTCCCCTCGCCCCTCTCAAGAGAGGTTTAGTGTGAGGGGAAAACAGTGTCATCAATATCAACCttgt                                                                                            >  2:931532/1‑150 (MQ=255)
                                                         gCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGtt                                                                                           >  1:889258/1‑150 (MQ=255)
                                                           ttACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGTTATGTCGGTCCCCTCGCCACTCTGTGGAGATGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTa                                                                                         <  1:276093/150‑1 (MQ=39)
                                                             aCCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAtt                                                                                       >  1:83766/1‑150 (MQ=255)
                                                                  cggCGGGACCAGGCGTGCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGTGTCAAAACGCCGCGCGGTCGGCCCCCTCGCCCCGCGGGGGTGAGGGTTAGGGGGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCg                                                                                  <  1:553146/149‑1 (MQ=255)
                                                                    gCCGGATGCGGCGTGAACGCCTGATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                      <  2:151654/96‑1 (MQ=37)
                                                                    gCCGGATGCGGCGTGAACGCCTGATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                      >  1:151654/1‑96 (MQ=37)
                                                                    gCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCa                                                                                >  1:1222189/1‑150 (MQ=255)
                                                                    gCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCa                                                                                >  1:1603291/1‑150 (MQ=255)
                                                                     ccGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaa                                                                               <  2:256818/150‑1 (MQ=255)
                                                                      gggATCAGGCGTTAACGCCGCATCCGGAAAAGTGTTAACCGTGTAAAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTGGCCCCTCTGGGGTGAGGGTTTGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaaa                                                                              <  2:1404345/149‑1 (MQ=255)
                                                                      cGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAATGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaaa                                                                              <  2:167212/150‑1 (MQ=255)
                                                                      cGGATCAGGCGGTCTCGCCGAATCCGTCAAAGTGTTATCCGTGTACAGAGTCAAAAACGCGCGCGGTCGGTCCCCTCGCCCCTCTGGTGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaaa                                                                              <  2:274334/150‑1 (MQ=38)
                                                                       ggATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaaaa                                                                             >  1:475495/1‑150 (MQ=255)
                                                                       ggATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaaaa                                                                             >  1:852703/1‑150 (MQ=255)
                                                                         ttaaGGCGTTCACGCGGCATGCGGCAAAGTGTTAACGGTGTACAGATTCAAAAAGGCGCTCGGGCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAaac                                                                           <  2:388772/147‑1 (MQ=38)
                                                                         aTCAGGGGTTCACGCCTCATCCGGCAAAGTGTTAACCGTGTAATGCGTCAAATAGCCGCGAGGGCGGCCCGCGCCCCCCTGTTGGTAGAGGGTTAGGGGGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAaac                                                                           <  2:556376/150‑1 (MQ=255)
                                                                             ggCGTTCAGGCCGCATCCGGCAAAGTATTAACCGTGTACAGAGGCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcg                                                                       <  2:251689/150‑1 (MQ=255)
                                                                               cGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTGCGCAACAACAAcgcg                                                                     >  1:84253/1‑150 (MQ=255)
                                                                               cGTGCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGGCAAAAAGGCTCTCGGTCGGTCCCCTCGCCCCTCTGGGGTGAGGGTTTGGGTGAGGGGAAAACCGTGGCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgcg                                                                     <  2:319566/150‑1 (MQ=39)
                                                                               agTTCACGCCGCATCCGGCAAAGTGTTAACCGTTTACAGAGTCAAACAGGCGCGCGGGCGGTCCCCGCGCCCCTCGGGGGAGAGGGTTTGGGGGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgcg                                                                     <  2:314888/149‑1 (MQ=255)
                                                                                 gtCACGCCGCATCCGGCAAAGTGGTAACCGGTTACAGAGTCAAAAAGGCGCGCGGGCGGTACCCTCGCCCCTCTGGGGAGAGGGTGAGGGGGAGGGGAAACCCGTGTCCGCACTATCAACCGTGGTTATTCTTCGCAAAAACAAcgcgcg                                                                   <  2:262636/149‑1 (MQ=37)
                                                                                   cACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGcc                                                                 >  2:357581/1‑150 (MQ=255)
                                                                                   cACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGcc                                                                 >  1:656737/1‑150 (MQ=255)
                                                                                     cgccgcATCCGGCAAAGTGTTAACCGTGTACAGATTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGAGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAg                                                               >  1:311218/1‑150 (MQ=255)
                                                                                     cgccgcATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGGGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAg                                                               <  1:650074/150‑1 (MQ=255)
                                                                                           aTCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGtt                                                         <  2:389333/150‑1 (MQ=255)
                                                                                            tCCTGCTAAGTGTTAACCGTGTACTGAGTCAAAAAGGCGCTGGGGCGGGCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGATAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTc                                                        <  2:1064135/150‑1 (MQ=38)
                                                                                            tCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTc                                                        >  2:1718228/1‑150 (MQ=255)
                                                                                              cGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcag                                                      >  2:379638/1‑150 (MQ=255)
                                                                                                gCAAAGTGATAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagca                                                    <  1:53380/150‑1 (MQ=255)
                                                                                                 cAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagcag                                                   >  2:711/1‑150 (MQ=255)
                                                                                                  aaaGTTTTAACCGTGAACAGCGTCACAACGGCGCTCGGGGGGTCCCCGCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTATATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagcagc                                                  <  2:1067292/150‑1 (MQ=38)
                                                                                                    aGTGTTAACCGTGTACAGAGTCAAAACGGCGCTCGGGCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTTAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagcagcag                                                <  2:1375506/150‑1 (MQ=255)
                                                                                                           aCCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTAGGcac                                         >  1:821534/1‑150 (MQ=255)
                                                                                                             cGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCGTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGGGTTGATTCTTCGCAAACCAACGCGCGCCAGAATGTTCAGCCGCAGTACGCACAgg                                      >  1:1150199/1‑150 (MQ=38)
                                                                                                                gtACAGAGTCAAACAGGCGCTCGGGCGGTCCCGTCGCCCCTCTGGGGCGAGGGTTAGGTTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGt                                    <  2:1713777/150‑1 (MQ=255)
                                                                                                                  aCAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCACAATGTTCAGCAGCAGTACGCACAGGGaaa                                  >  2:1536262/1‑150 (MQ=255)
                                                                                                                      agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGAt                              >  1:171572/1‑150 (MQ=255)
                                                                                                                       gTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACACCGCGCGCCAGAATGTTCAGCAGCAGTCCGCACAGGGTAATGATc                             >  2:1454675/1‑150 (MQ=255)
                                                                                                                       gTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATc                             >  1:1313132/1‑150 (MQ=255)
                                                                                                                         cAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCaa                           >  1:376704/1‑150 (MQ=255)
                                                                                                                          aaaaaGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGCTccag                          >  1:304422/1‑147 (MQ=255)
                                                                                                                              cggCGCTCGGTCGGGCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATAccc                      <  1:1238284/149‑1 (MQ=255)
                                                                                                                              aGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATAccc                      >  2:226382/1‑150 (MQ=255)
                                                                                                                                   cTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGcc                 <  2:405238/150‑1 (MQ=255)
                                                                                                                                    tcggtcggTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGccc                <  2:1148405/150‑1 (MQ=255)
                                                                                                                                      ggtcggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAACAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTTATGTGCAATACCCCGGCCCAg              >  2:263737/1‑150 (MQ=255)
                                                                                                                                      ggtcggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAAGGATCAATACCCCGGCCCAg              >  2:1440278/1‑150 (MQ=255)
                                                                                                                                       ggcggccCCCTCGCCCCTCTGGGGTGAGGGTTAGGGGGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAgg             <  2:753678/148‑1 (MQ=255)
                                                                                                                                        ccctaacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCa                            >  2:1283406/7‑134 (MQ=255)
                                                                                                                                        ccctaacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCa                            <  1:1283406/128‑1 (MQ=255)
                                                                                                                                        ccctaacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGc            >  2:1291207/7‑150 (MQ=255)
                                                                                                                                           taacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCaa         >  1:879468/4‑150 (MQ=255)
                                                                                                                                             ccccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTCCGCAAAAACAACGCGCGCCAGAATGTTCCGCAGCAGTACGCACAGGGTAAAGAGCACTACCCCGGCCCACGCCACtt       >  2:1342511/1‑150 (MQ=38)
                                                                                                                                             acccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAACCCGAGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACCGGGTAATGATCACTACCCCGGCCCAGGCCCAtt       >  1:237555/2‑150 (MQ=39)
                                                                                                                                                cTCGCCCCTCTGGGGTGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgct    >  1:361479/1‑150 (MQ=255)
                                                                                                                                                 tCGCCCCTCTGGGGAGAGGGTTAGGGTTAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgctg   <  1:1337682/150‑1 (MQ=255)
                                                                                                                                                 tCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgctg   <  2:34237/150‑1 (MQ=255)
                                                                                                                                                 gcGCCCCTCGGGGGTGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgctg   <  2:884872/149‑1 (MQ=255)
                                                                                                                                                  cGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgctgc  <  2:83766/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGC  >  NC_000913/3908353‑3908648

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: