Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 794,943 T→G 7.0% intergenic (+21/‑146) acrZ → / → modA AcrAB‑TolC efflux pump accessory protein, membrane‑associated/molybdate ABC transporter periplasmic binding protein; chlorate resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913794,9430TG 7.0% 360.1 / 14.1 129intergenic (+21/‑146)acrZ/modAAcrAB‑TolC efflux pump accessory protein, membrane‑associated/molybdate ABC transporter periplasmic binding protein; chlorate resistance protein
Reads supporting (aligned to +/- strand):  ref base T (62/58);  new base G (4/5);  total (66/63)
Fisher's exact test for biased strand distribution p-value = 7.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.12e-01

Only 100 of 130 total aligned reads displayed.

TGGTATTCGCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGTTACCAATGG  >  NC_000913/794795‑795092
                                                                                                                                                    |                                                                                                                                                     
tGGTATTCGCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTc                                                                                                                                                      <  2:524419/150‑1 (MQ=255)
 ggTATTCGCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCg                                                                                                                                                     <  2:625252/150‑1 (MQ=255)
      tCGCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAt                                                                                                                                                <  1:219713/150‑1 (MQ=255)
       cGCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTATTCGGGCAtt                                                                                                                                               >  1:573961/1‑150 (MQ=255)
         ccGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGACTAAGCGGGCAttca                                                                                                                                             <  2:1765926/150‑3 (MQ=255)
         ccGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCACGCTTAGTCGGGCAtttt                                                                                                                                             >  1:1180311/1‑150 (MQ=255)
         ccGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTACCTGAATGCCCGCTTAGTCGGGCCtttt                                                                                                                                             >  2:239556/1‑150 (MQ=255)
          cGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAttttc                                                                                                                                            <  2:1661881/150‑1 (MQ=255)
          cGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAtttta                                                                                                                                            >  1:79059/1‑149 (MQ=255)
           gTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAATAAGACCAGCCCGGACAAATTCATTGATTCCCTGAATGCCCGACTAAGCGGGCAttcagg                                                                                                                                           <  2:4036/150‑5 (MQ=255)
           gTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGACTAAGCGGGCAttcatg                                                                                                                                           <  1:1781597/150‑5 (MQ=255)
           gTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGACTAAGCGGGCAttcagg                                                                                                                                           <  2:1646696/150‑5 (MQ=255)
            tAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAttttctt                                                                                                                                          <  1:1942207/150‑1 (MQ=255)
              aTCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAttttctttt                                                                                                                                        >  1:882650/1‑150 (MQ=255)
               tCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAttttcttttt                                                                                                                                       >  2:1334206/1‑150 (MQ=255)
                cATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAttttctttttc                                                                                                                                      <  2:1764284/150‑1 (MQ=255)
                   ggTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCa                                                                                                                                   >  1:884340/1‑150 (MQ=255)
                   ggTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCATTTTCTCa                                                                                                                                   >  2:251187/1‑150 (MQ=255)
                     tACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAAc                                                                                                                                 >  1:1629158/1‑150 (MQ=255)
                         tGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGACTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCCGCCAGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTAGGGAATTTTCTATTTATCAACATcc                                                                                                                             >  1:75164/1‑150 (MQ=255)
                          gTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCt                                                                                                                            >  2:1786346/1‑150 (MQ=255)
                          gTCGTGATGGCCATCATCCTGGGTCTGATTAACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCt                                                                                                                            >  2:175581/1‑150 (MQ=255)
                            cGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGc                                                                                                                          >  2:1672645/1‑150 (MQ=255)
                            cGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCt                                                                                                                         >  2:796188/1‑150 (MQ=255)
                                  ggCCATCATCCTTGGGCTGATTTACGTTCTTGGTTAAGTATTCATCATGTTTTCTGGTGTTTGTAAAAAAGTCCAGCCCGGACAAAATCATTGATTCCCTGAATTCCAGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGGTTTTcc                                                                                                                    <  2:573326/150‑1 (MQ=255)
                                    ccATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTg                                                                                                                  <  1:161361/150‑1 (MQ=255)
                                    ccATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTg                                                                                                                  >  2:267079/1‑150 (MQ=255)
                                    ccATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTAGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTg                                                                                                                  >  1:1212389/1‑150 (MQ=255)
                                     catcatCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGc                                                                                                                 <  1:68845/150‑1 (MQ=255)
                                     catcatCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTACCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTATCCTGc                                                                                                                 >  1:540022/1‑150 (MQ=255)
                                        catcCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGa                                                                                                              >  2:1597798/1‑150 (MQ=255)
                                          tcCTGCGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGATGGTAAAAAAGACCAGACAGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGTCATTTACTTTTTCTCAGCTTCCTGCTATTCCTGCCGatt                                                                                                            >  2:1242063/1‑149 (MQ=255)
                                             tGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGGAAAAAAGCCCAGCCCGGACAAAATAATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCACCTTCCTGCTTTTCCTGCCGATAttt                                                                                                         >  1:1110748/1‑150 (MQ=255)
                                               ggTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTACCTGAATGAACGCTTAGTCGGGCATTTTCTATTTCTCAACTTCCTGCTTTTCCTGCCGATAatttt                                                                                                       >  2:889134/1‑150 (MQ=255)
                                                gTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATAtttttt                                                                                                      >  2:1382441/1‑150 (MQ=255)
                                                  cTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCt                                                                                                    >  1:76694/1‑150 (MQ=255)
                                                  cTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCt                                                                                                    <  1:1190440/150‑1 (MQ=255)
                                                   tGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCtt                                                                                                   >  1:284404/1‑150 (MQ=255)
                                                    gATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTa                                                                                                  <  1:837472/150‑1 (MQ=255)
                                                    gATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTa                                                                                                  >  2:170702/1‑150 (MQ=255)
                                                      tttACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATc                                                                                                <  2:1815041/150‑1 (MQ=255)
                                                       ttACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCt                                                                                               >  1:1695770/1‑150 (MQ=255)
                                                          cGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGCCCAGCCCGGACAAAATCATTGATTACCTGAATGCCCGCTTATTCGGGCATTTTCTTTTTATCAACTTAATGATTTTCATTCAGATATTTTTTCTTATCTAcc                                                                                            >  1:895008/1‑150 (MQ=255)
                                                          cGGGCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTAcc                                                                                            <  2:1565869/150‑1 (MQ=255)
                                                            gTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCCTTTCCTGCCGATATTTTTTCTTATCTACCTc                                                                                          <  2:804268/150‑1 (MQ=255)
                                                                tGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACaa                                                                                      >  2:884916/1‑150 (MQ=255)
                                                                 ggTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACaaa                                                                                     >  1:38732/1‑150 (MQ=255)
                                                                   tGAAGTATTCAACATGTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGATTGCCCGCGTAGTCGGGCATTTTCTTTTTCTCCACTTCCTGCTTTTCCTGCCGTTATTTTTTCTTATATAGCTGACAAAgg                                                                                   >  2:104144/1‑150 (MQ=255)
                                                                    gAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGt                                                                                  >  1:1514675/1‑150 (MQ=255)
                                                                    gAAGTATTCAACATCTTTTCTGGAGTTGGTAAAAAAGACCAGCCCTGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGt                                                                                  <  2:175333/150‑1 (MQ=255)
                                                                            cAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAGAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATa                                                                          >  1:782454/1‑150 (MQ=255)
                                                                            cAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATa                                                                          <  2:239665/150‑1 (MQ=255)
                                                                            cAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGACCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAAATTACTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCCCAAAGGTTAGAAATa                                                                          >  2:141373/1‑150 (MQ=255)
                                                                               cATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGTGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCt                                                                                               <  2:583706/126‑1 (MQ=255)
                                                                                     ttCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTggg                                                                 <  1:1346825/150‑1 (MQ=255)
                                                                                       cTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTACTGCTTTTCCTGCCGATATTTTTTCTTATATACCTCACAAAGGTTAGCAATAACTGCTGGGaa                                                               >  2:1461813/1‑150 (MQ=255)
                                                                                        tGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGaaa                                                              >  1:754410/1‑150 (MQ=255)
                                                                                          gtgtTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAAt                                                            <  2:284404/150‑1 (MQ=255)
                                                                                           tgtTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAAtt                                                           >  1:1396784/1‑150 (MQ=255)
                                                                                             ttGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTcc                                                         <  1:447120/150‑1 (MQ=255)
                                                                                              tGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCg                                                        <  2:1780031/150‑1 (MQ=255)
                                                                                              tGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCg                                                        <  2:1281167/150‑1 (MQ=255)
                                                                                               ggTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTTGCAATAACTGCTGGGAAAATTCCGa                                                       >  1:888383/1‑150 (MQ=255)
                                                                                               ggTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGa                                                       <  2:1575506/150‑1 (MQ=255)
                                                                                               ggTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGa                                                       >  1:1926943/1‑150 (MQ=255)
                                                                                                gTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCAGATATTTTTTCGTATCTACCTCACAACGGTTAGCAATTCCTGAGGGTAAAATTCCCAg                                                      >  1:387376/1‑150 (MQ=255)
                                                                                                gTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACATCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAg                                                      >  1:530611/1‑150 (MQ=255)
                                                                                                      aaGACCAGCCCTGACAAAATCTTTGATTCCCTGAATGCACGATTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGACGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTc                                                <  2:776480/150‑1 (MQ=255)
                                                                                                      aaGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTc                                                >  2:1084426/1‑150 (MQ=255)
                                                                                                       aGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCg                                               <  1:175581/150‑1 (MQ=255)
                                                                                                         aCCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGtt                                             <  2:1038841/150‑1 (MQ=255)
                                                                                                          ccAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTta                                            <  1:991162/150‑1 (MQ=255)
                                                                                                          ccAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTACGAGTTAGTCGTta                                            >  2:1909576/1‑150 (MQ=255)
                                                                                                               ccGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCATATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCTAGTTAGTCGTTATATTg                                       >  2:1562179/1‑150 (MQ=255)
                                                                                                                cGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGATGGGAAAATTCCGAGTTAGTCGTTATAtttt                                      >  1:1961108/1‑148 (MQ=255)
                                                                                                                  gACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCg                                    <  2:906194/150‑1 (MQ=255)
                                                                                                                    cAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGcc                                  <  2:1955919/150‑1 (MQ=255)
                                                                                                                      aaaTCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTa                                >  2:1022637/1‑150 (MQ=255)
                                                                                                                      aaaTCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCATGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTTCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTa                                >  2:609292/1‑150 (MQ=255)
                                                                                                                         tCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCACAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAAATCCGAGTTAGTCGTTATATTGTCGCCTACAt                             <  2:291162/150‑1 (MQ=255)
                                                                                                                              tcttCCCTTTATGCACGGGTATTCGGGGATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGt                        <  1:809351/148‑1 (MQ=255)
                                                                                                                              gATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGt                        >  2:1386178/1‑150 (MQ=255)
                                                                                                                               aTTCCCTGAATGACCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGtt                       <  2:1629157/150‑1 (MQ=255)
                                                                                                                                ttCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTa                      <  1:1672646/150‑1 (MQ=255)
                                                                                                                                 tCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTAc                     <  1:969504/150‑1 (MQ=255)
                                                                                                                                 aaagaaaaTGCCCGACTAAGCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTAc                     >  2:484746/7‑150 (MQ=255)
                                                                                                                                 aaagaaaaTGCCCGACTAAGCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTAc                     >  2:490144/7‑150 (MQ=255)
                                                                                                                                  cccTGAATGCCCGCTTTTTCGGTCTTTTTATTTTTCTCAACTTCCTGCTTTTCCTGGCGTTATTTTTTCTAATGTTCCTCTCAAAGGTTATCAATATCTGCTGTGAAATTTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACa                    <  1:1276375/150‑1 (MQ=255)
                                                                                                                                    cTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACAtt                  >  2:84999/1‑150 (MQ=255)
                                                                                                                                      gAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTaa                >  1:869141/1‑150 (MQ=255)
                                                                                                                                        aTGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACATCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAgg              <  1:1786347/150‑1 (MQ=255)
                                                                                                                                        aTGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTAACTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAgg              <  1:971330/150‑1 (MQ=255)
                                                                                                                                         tGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAggg             <  1:791082/150‑1 (MQ=255)
                                                                                                                                          gCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTATCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAgggg            <  1:1334207/150‑1 (MQ=255)
                                                                                                                                          gCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTATCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGACGTTATATTGTCGCCTACATAACGTTACATTAAgggg            <  1:1696759/150‑1 (MQ=255)
                                                                                                                                            ccGCTTAGTCGGGCTTTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACATCTCAAAGGTTAGCAATAACTGCTGGGAAAATTAATAGTTAGTCTATATATTGTCGCCTACATAACTTTACATTAAGTGGtt          >  1:333433/1‑150 (MQ=255)
                                                                                                                                            ccGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGtt          <  2:79059/150‑1 (MQ=255)
                                                                                                                                               cTTATTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGTTAcc       <  2:573961/150‑1 (MQ=255)
                                                                                                                                                  aGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGTGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGTTACCAAt    <  1:214339/150‑1 (MQ=255)
                                                                                                                                                    tCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGTTACCAATgg  <  2:1727562/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGGTATTCGCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGTTACCAATGG  >  NC_000913/794795‑795092

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