Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,324,952 A→G 5.4% intergenic (‑41/+49) folP ← / ← ftsH 7,8‑dihydropteroate synthase/protease, ATP‑dependent zinc‑metallo

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,324,9520AG 5.4% 435.4 / 10.8 147intergenic (‑41/+49)folP/ftsH7,8‑dihydropteroate synthase/protease, ATP‑dependent zinc‑metallo
Reads supporting (aligned to +/- strand):  ref base A (72/67);  new base G (3/5);  total (75/72)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.35e-01

Only 100 of 147 total aligned reads displayed.

CGAGTTATGCGTGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAG  >  NC_000913/3324804‑3325098
                                                                                                                                                    |                                                                                                                                                  
cGAGTTATGCGTGCCACCATCCGGAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAg                                                                                                                                                   <  2:1938072/150‑1 (MQ=255)
   gTTATGCGTGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACGCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCAc                                                                                                                                                >  1:188107/1‑150 (MQ=255)
   gTTATGCGTGACACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGGTGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGGTGCATGATGTTATCCCTGGTATGAATTGATTAGAAAAAACCCCGGAGCAc                                                                                                                                                >  1:115204/1‑150 (MQ=255)
         cGTGCCACCATCCGAATAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGACGTCGCCTGGGGAAAGAGTTTCATGAGGTTATCGCTGGTATGAAGTGATAAGAAAAAACCCCGGAGCACGCCCCg                                                                                                                                          >  1:301006/1‑150 (MQ=255)
         cGTGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCg                                                                                                                                          >  1:1824946/1‑150 (MQ=255)
         cGTGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCg                                                                                                                                          >  1:210119/1‑150 (MQ=255)
            gCCACCATCCTAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCgggg                                                                                                                                       <  2:1831485/150‑1 (MQ=255)
            gCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCgggg                                                                                                                                       <  1:1335107/150‑1 (MQ=255)
            gCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCgggg                                                                                                                                       <  1:1790667/150‑1 (MQ=255)
            gCCACCATCCGAAAACGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCgggg                                                                                                                                       <  1:607885/150‑1 (MQ=255)
             ccaccaTCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGGGCGTGCTCCGGGGt                                                                                                                                      <  1:861210/150‑1 (MQ=255)
              caccaTCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGtt                                                                                                                                     >  2:1557759/1‑150 (MQ=255)
               accaTCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGttt                                                                                                                                    >  2:848718/1‑150 (MQ=255)
                  aTCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCg                                                                                                                                 <  1:1442061/150‑1 (MQ=255)
                  aTCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCg                                                                                                                                 <  1:1942665/150‑1 (MQ=255)
                   tCCGAAAAGGAATCAGGCTTGACGTTGTGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGATGTTCCCTGGGCAAAGAGTTGCATGATGTTATCCCTGTTATGAAGTGATAAGAAAAAACCCCGGAGCTCGCCCCGGGGTTTTCgg                                                                                                                                <  2:365002/150‑1 (MQ=255)
                   accGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGGGCGTGCTCCGGGGttttttc                                                                                                                                <  1:324748/149‑4 (MQ=255)
                    ccGAAAATGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAATGTCCAGTGAAGTTCCCTGGGCAAAGAGTTTGATGATGTTATCCCTGGTATGAATTTATAAGAAAAAACCCCGGGGCGTGCTCCGGGGttttttct                                                                                                                               <  2:1582908/150‑5 (MQ=255)
                    ccGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGt                                                                                                                               >  2:1687821/1‑150 (MQ=255)
                     cGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTa                                                                                                                              >  2:1326198/1‑150 (MQ=255)
                     cGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCACTGGTATGAATTGATAAGAAAAAACCCCGGAGCCCGCCCCGGGGTTTTCGGTa                                                                                                                              >  1:176281/1‑150 (MQ=255)
                      gAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGGGCGTGCTCCGGGGttttttctta                                                                                                                             <  1:1801246/150‑7 (MQ=255)
                      gAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTAc                                                                                                                             >  2:1725729/1‑150 (MQ=255)
                      gAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTAc                                                                                                                             <  2:699187/150‑1 (MQ=255)
                      gAAAAGGAATCACGAGGGACGGTGAGGCTCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTGCAAGATGGTATCCCTGGTATGAATCGATAACAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTAc                                                                                                                             <  2:1330878/150‑1 (MQ=255)
                       aaaaGGAATCAGTCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGAAGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACa                                                                                                                            <  1:1302263/150‑1 (MQ=255)
                       aaaaGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACa                                                                                                                            <  1:518465/150‑1 (MQ=255)
                           ggAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAACGGCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATa                                                                                                                        <  1:595073/150‑1 (MQ=255)
                            gAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATAc                                                                                                                       >  2:1024540/1‑150 (MQ=255)
                            gAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATAc                                                                                                                       <  2:1644163/150‑1 (MQ=255)
                              aTCAGGCGTTACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAACGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAg                                                                                                                     <  1:1110980/150‑1 (MQ=255)
                              aTCAGGCGGGACGTTGAGGATCCCCCTTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCCTGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAg                                                                                                                     <  1:1117979/150‑1 (MQ=255)
                                 aGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTAGGTACAAATACAGTCa                                                                                                                  >  1:1652772/1‑150 (MQ=255)
                                 aGGCGTGACGTTGAGGATCCCCATTACCTGAGGCTGGCTACTGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAATACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCa                                                                                                                  <  2:120951/150‑1 (MQ=255)
                                    cGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGGTTTCGGTACAAATACAGTCATCt                                                                                                               >  2:22497/1‑150 (MQ=255)
                                     gTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTg                                                                                                              <  2:756829/150‑1 (MQ=255)
                                          gTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCg                                                                                                         <  1:1476891/150‑1 (MQ=255)
                                           ttGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATATGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCgg                                                                                                        >  2:1125614/1‑150 (MQ=255)
                                           ttGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCgg                                                                                                        >  2:1589190/1‑150 (MQ=255)
                                            tGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCggg                                                                                                       >  1:20582/1‑150 (MQ=255)
                                              aGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGaa                                                                                                     >  2:1024237/1‑150 (MQ=255)
                                              aGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGaa                                                                                                     >  1:838467/1‑150 (MQ=255)
                                               ggATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAac                                                                                                    >  1:1676059/1‑150 (MQ=255)
                                               ggATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAac                                                                                                    <  1:1217858/150‑1 (MQ=255)
                                                gATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAact                                                                                                   >  2:1777717/1‑150 (MQ=255)
                                                gATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGACATCTGATGCGGGAact                                                                                                   >  2:725633/1‑150 (MQ=255)
                                                   ccccATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAacttac                                                                                                >  2:1599266/1‑150 (MQ=255)
                                                     ccATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGTGGTTTTCGGTACAAATACAGTCATCTGATGCTGGAacttactt                                                                                              >  2:482245/1‑150 (MQ=255)
                                                      cATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTg                                                                                             >  1:1654760/1‑150 (MQ=255)
                                                      cATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTg                                                                                             <  2:221366/150‑1 (MQ=255)
                                                         tACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATTATGTTATCCCTGGTATGAATTCTTAAGAAAAAACCCCGCAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCg                                                                                          <  2:1810121/150‑1 (MQ=255)
                                                             tGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTa                                                                                      >  1:1906350/1‑150 (MQ=255)
                                                               aGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAAc                                                                                    <  2:1153734/150‑1 (MQ=255)
                                                                 gATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTg                                                                                  >  1:1960881/1‑150 (MQ=255)
                                                                 gATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGc                                                                                 <  2:1878003/150‑1 (MQ=255)
                                                                    ggCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGctc                                                                               <  2:613187/150‑1 (MQ=255)
                                                                       tAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGa                                                                            <  1:1620334/150‑1 (MQ=255)
                                                                         aGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACa                                                                          <  2:1906349/150‑1 (MQ=255)
                                                                          ggTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACAt                                                                         >  2:917592/1‑150 (MQ=255)
                                                                          ggTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACAt                                                                         >  2:1429059/1‑150 (MQ=255)
                                                                             ccATTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGgt                                                                      <  2:1794526/150‑1 (MQ=255)
                                                                             ccAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCATAACTGCTCTGACATGgt                                                                      >  2:1022507/1‑150 (MQ=255)
                                                                               tgTGAAGTACGCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAATCCCCGGAGCACGCCCGGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGgtgt                                                                    <  2:917988/149‑1 (MQ=255)
                                                                               aGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGgtgt                                                                    <  1:1594411/150‑1 (MQ=255)
                                                                               aGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGAATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGgtgt                                                                    >  1:18077/1‑150 (MQ=255)
                                                                                gTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGtt                                                                   >  2:1803716/1‑150 (MQ=255)
                                                                                gTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGtt                                                                   >  2:428720/1‑150 (MQ=255)
                                                                                  gAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTGTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTAc                                                                 >  1:366251/1‑150 (MQ=255)
                                                                                     gTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCg                                                              <  1:928326/150‑1 (MQ=255)
                                                                                           tGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCg                                                        <  2:1388655/150‑1 (MQ=255)
                                                                                            gggCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTTGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCgg                                                       <  1:1750400/150‑1 (MQ=255)
                                                                                              gCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCg                                                     <  2:1798150/150‑1 (MQ=255)
                                                                                              gCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCg                                                     >  2:373559/1‑150 (MQ=255)
                                                                                              gCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCg                                                     <  2:764197/150‑1 (MQ=255)
                                                                                              gCAAAGAGGTTCATGATGGTATCCCTGGTATGAATAGATAAGAAAACACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCg                                                     <  2:292813/150‑1 (MQ=255)
                                                                                                   gagTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTAcgc                                                >  2:1668160/1‑150 (MQ=255)
                                                                                                      tttCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGt                                             <  1:587500/150‑1 (MQ=255)
                                                                                                           tgatgTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTtcatc                                        >  1:1914872/1‑150 (MQ=255)
                                                                                                           tgatgTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGAGATACCCGGGTTCGGCGTACGCGTatcatc                                        >  1:1369213/1‑150 (MQ=255)
                                                                                                            gatgTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTtcatca                                       <  2:210119/150‑1 (MQ=255)
                                                                                                            gatgTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTtcatca                                       <  1:701541/150‑1 (MQ=255)
                                                                                                             atgTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCaa                                      >  1:509523/1‑150 (MQ=255)
                                                                                                              tgTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAAc                                     >  2:336983/1‑150 (MQ=255)
                                                                                                                 tATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCgg                                  >  1:1178082/1‑150 (MQ=255)
                                                                                                                 tATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCgg                                  >  1:270705/1‑150 (MQ=255)
                                                                                                                 tATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCgg                                  >  2:573180/1‑150 (MQ=255)
                                                                                                                    cccTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACg                               >  2:1521448/1‑150 (MQ=255)
                                                                                                                     ccTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACga                              >  1:621948/1‑150 (MQ=255)
                                                                                                                      cTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACgag                             <  2:1206789/150‑1 (MQ=255)
                                                                                                                        ggTATGAATTGATAAGAAAAAACCCCGGAGCACGCACCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACgagga                           >  1:129661/1‑150 (MQ=255)
                                                                                                                             gAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTATCGGTACAAATACAGTCATCTCATGCGGGAACTTACTTGTCGCCTAACTGCTCTCACAAGGTGTTACCCGGCTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCtt                      >  2:2023242/1‑150 (MQ=255)
                                                                                                                              aaTTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTATCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCttt                     <  1:364241/150‑1 (MQ=255)
                                                                                                                                ttGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTAAAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTgg                   <  2:828041/150‑1 (MQ=255)
                                                                                                                                  gATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGCGCCTTTGGAc                 >  2:184859/1‑150 (MQ=255)
                                                                                                                                  gATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGAGGACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGAc                 >  2:1599674/1‑150 (MQ=255)
                                                                                                                                     aaGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCCTCTGATGCGGGCACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTAc              <  1:1112797/150‑1 (MQ=255)
                                                                                                                                      aGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTAcc             >  1:1232662/1‑150 (MQ=255)
                                                                                                                                      aGAAAAAACCCCGGAGCACGCCCAGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTACCTGCTCTGACATTGTGGTACCCGGGTTCGGCGTACGCGGTTCATCACCCGGACGAGGAGCCTTTGGACTaac             >  1:744042/1‑148 (MQ=255)
                                                                                                                                             aCCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCg      >  1:893289/1‑150 (MQ=255)
                                                                                                                                               cccGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGcc    >  1:1790479/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CGAGTTATGCGTGCCACCATCCGAAAAGGAATCAGGCGTGACGTTGAGGATCCCCATTACGTGAGGATGGCTAAGGTCCAGTGAAGTACCCTGGGCAAAGAGTTTCATGATGTTATCCCTGGTATGAATTGATAAGAAAAAACCCCGGAGCACGCCCCGGGGTTTTCGGTACAAATACAGTCATCTGATGCGGGAACTTACTTGTCGCCTAACTGCTCTGACATGGTGTTACCCGGGTTCGGCGTACGCGGTTCATCAACCGGACGAGGAGCCTTTGGACTACCATTGTCGCCAG  >  NC_000913/3324804‑3325098

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: