Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 794,939 T→C 7.0% intergenic (+17/‑150) 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,9390TC 7.0% 374.0 / 13.5 129intergenic (+17/‑150)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 (61/57);  new base C (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 = 9.96e-01

Only 100 of 131 total aligned reads displayed.

AAGTCTGGTATTCGCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGTTACC  >  NC_000913/794790‑795087
                                                                                                                                                     |                                                                                                                                                    
aaGTCTGGTATTCGCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCt                                                                                                                                                      <  1:1140473/150‑1 (MQ=255)
 aGTCTGGTATTCGCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCtt                                                                                                                                                     <  2:729435/150‑1 (MQ=255)
     tGGTATTCGCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTc                                                                                                                                                 <  2:524419/150‑1 (MQ=255)
      ggTATTCGCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCg                                                                                                                                                <  2:625252/150‑1 (MQ=255)
        tATTCTCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCGGGGGCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCggg                                                                                                                                              <  1:1833280/150‑1 (MQ=255)
           tCGCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAt                                                                                                                                           <  1:219713/150‑1 (MQ=255)
            cGCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTATTCGGGCAtt                                                                                                                                          >  1:573961/1‑150 (MQ=255)
            cGCCGTAATAATGGGACCTGTCGTGATGGCCATCATCCTGTGTCTGATTTACGGTCTAGGTGAAGTATTCAACATCTTTTGTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGACTAAGCGGGCAtt                                                                                                                                          <  2:259517/150‑1 (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:1385349/1‑149 (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)
                          tACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAAc                                                                                                                            >  1:1629158/1‑150 (MQ=255)
                             cTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGACCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTc                                                                                                                         >  2:2005016/1‑150 (MQ=255)
                              tGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGACTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCCGCCAGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTAGGGAATTTTCTATTTATCAACATcc                                                                                                                        >  1:75164/1‑150 (MQ=255)
                               gTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCt                                                                                                                       >  2:1786346/1‑150 (MQ=255)
                                 cGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCt                                                                                                                    >  2:796188/1‑150 (MQ=255)
                                       ggCCATCATCCTTGGGCTGATTTACGTTCTTGGTTAAGTATTCATCATGTTTTCTGGTGTTTGTAAAAAAGTCCAGCCCGGACAAAATCATTGATTCCCTGAATTCCAGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGGTTTTcc                                                                                                               <  2:573326/150‑1 (MQ=255)
                                         ccATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTg                                                                                                             >  2:267079/1‑150 (MQ=255)
                                          catcatCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGc                                                                                                            <  1:68845/150‑1 (MQ=255)
                                          catcatCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTACCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTATCCTGc                                                                                                            >  1:540022/1‑150 (MQ=255)
                                           atcatcCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTTCTTTTCCTGcc                                                                                                           <  1:2038799/150‑1 (MQ=255)
                                           atcatcCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGcc                                                                                                           <  2:1570940/150‑1 (MQ=255)
                                             catcCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGa                                                                                                         >  2:1597798/1‑150 (MQ=255)
                                               tcCTGCGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGATGGTAAAAAAGACCAGACAGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGTCATTTACTTTTTCTCAGCTTCCTGCTATTCCTGCCGatt                                                                                                       >  2:1242063/1‑149 (MQ=255)
                                                  tGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATAttt                                                                                                    >  1:778521/1‑150 (MQ=255)
                                                  tGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGGAAAAAAGCCCAGCCCGGACAAAATAATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCACCTTCCTGCTTTTCCTGCCGATAttt                                                                                                    >  1:1110748/1‑150 (MQ=255)
                                                    ggTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTACCTGAATGAACGCTTAGTCGGGCATTTTCTATTTCTCAACTTCCTGCTTTTCCTGCCGATAatttt                                                                                                  >  2:889134/1‑150 (MQ=255)
                                                     gTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATAtttttt                                                                                                 >  2:1382441/1‑150 (MQ=255)
                                                      tCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACGTCCTGCTTTTCCTGCCGATATTTTTTc                                                                                                >  2:1021526/1‑150 (MQ=255)
                                                      tCTGATTTACGGTCTTGGTGAAGTATTCAACATATTTTCTGGTGTTGGTAAAAAAGTCCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTc                                                                                                <  1:162472/150‑1 (MQ=255)
                                                       cTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCt                                                                                               <  2:1757998/150‑1 (MQ=255)
                                                       cTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCt                                                                                               <  1:1190440/150‑1 (MQ=255)
                                                        tGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCtt                                                                                              >  1:284404/1‑150 (MQ=255)
                                                         gATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTa                                                                                             >  2:170702/1‑150 (MQ=255)
                                                         gATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTa                                                                                             <  1:837472/150‑1 (MQ=255)
                                                           tttACGGTCTTGGTGAAGTATTCAACATCTTTTGTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAGTCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATc                                                                                           <  1:1709347/150‑1 (MQ=255)
                                                           tttACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATc                                                                                           <  2:1815041/150‑1 (MQ=255)
                                                               cGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGCCCAGCCCGGACAAAATCATTGATTACCTGAATGCCCGCTTATTCGGGCATTTTCTTTTTATCAACTTAATGATTTTCATTCAGATATTTTTTCTTATCTAcc                                                                                       >  1:895008/1‑150 (MQ=255)
                                                                 gTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTc                                                                                     <  2:1797442/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)
                                                                      ggTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACaaa                                                                                >  2:607311/1‑150 (MQ=255)
                                                                         gAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGt                                                                             >  1:1514675/1‑150 (MQ=255)
                                                                         gAAGTATTCAACATCTTTTCTGGAGTTGGTAAAAAAGACCAGCCCTGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGt                                                                             <  2:175333/150‑1 (MQ=255)
                                                                             tATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGc                                                                         <  2:1451494/150‑1 (MQ=255)
                                                                             tATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGc                                                                         <  2:1710151/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)
                                                                                            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)
                                                                                                   tGGTAAAAAAGACCAGCCCGGACAAAATAATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCg                                                   <  1:879944/150‑1 (MQ=255)
                                                                                                    ggTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGa                                                  <  2:1575506/150‑1 (MQ=255)
                                                                                                    ggTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGa                                                  >  1:1926943/1‑150 (MQ=255)
                                                                                                     ttAAAAAAGCCCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAg                                                 <  1:92891/149‑1 (MQ=255)
                                                                                                     gTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCAGATATTTTTTCGTATCTACCTCACAACGGTTAGCAATTCCTGAGGGTAAAATTCCCAg                                                 >  1:387376/1‑150 (MQ=255)
                                                                                                     gTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACATCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAg                                                 >  1:530611/1‑150 (MQ=255)
                                                                                                      tAAAAAAGACCAGCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGtt                                               <  1:796189/150‑1 (MQ=255)
                                                                                                       aaaaaaGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGATTAGTCGGGCATTTTATTTTTCTCAACTTCCTGCTTTTCCTGCAGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGtt                                               <  2:471117/150‑1 (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                                       >  2:1349725/1‑150 (MQ=255)
                                                                                                               ccAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTACGAGTTAGTCGTta                                       >  2:1909576/1‑150 (MQ=255)
                                                                                                               ccAGCCCGGACAAAATCAGTGATTCCCTGAATGCCCGCTTAGTCGGGCAGGTTCTTTTTCTCCCCTTCCGGCTTGTCCTGCCGATATTTTTTCTTATCGACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTta                                       >  1:36443/1‑150 (MQ=255)
                                                                                                                    ccGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCATATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCTAGTTAGTCGTTATATTg                                  >  2:1562179/1‑150 (MQ=255)
                                                                                                                       gACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCg                               <  2:906194/150‑1 (MQ=255)
                                                                                                                           aaaTCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTa                           >  2:1022637/1‑150 (MQ=255)
                                                                                                                           aaaTCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCATGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTTCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTa                           >  2:609292/1‑150 (MQ=255)
                                                                                                                              tCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACAt                        >  2:1066465/1‑150 (MQ=255)
                                                                                                                                   tcttCCCTTTATGCACGGGTATTCGGGGATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGt                   <  1:809351/148‑1 (MQ=255)
                                                                                                                                    aTTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGtt                  <  2:540022/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)
                                                                                                                                       cccTGAATGCCCGCTTTTTCGGTCTTTTTATTTTTCTCAACTTCCTGCTTTTCCTGGCGTTATTTTTTCTAATGTTCCTCTCAAAGGTTATCAATATCTGCTGTGAAATTTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACa               <  1:1276375/150‑1 (MQ=255)
                                                                                                                                       aagaaaaTGCCCGACTAAGCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCAGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACa               >  1:1246245/6‑150 (MQ=255)
                                                                                                                                        agaaaaTGCCCGACTAAGCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGAAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTAGCCTACATAACGTTACAt              >  1:483981/5‑150 (MQ=255)
                                                                                                                                         cTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACAtt             >  2:84999/1‑150 (MQ=255)
                                                                                                                                             aTGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTAACTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAgg         <  1:971330/150‑1 (MQ=255)
                                                                                                                                               gCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTATCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAgggg       <  1:1334207/150‑1 (MQ=255)
                                                                                                                                               gCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTATCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGACGTTATATTGTCGCCTACATAACGTTACATTAAgggg       <  1:1696759/150‑1 (MQ=255)
                                                                                                                                                 ccGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGtt     <  2:1655870/150‑1 (MQ=255)
                                                                                                                                                 ccGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGtt     <  2:79059/150‑1 (MQ=255)
                                                                                                                                                    cTTATTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGTTAcc  <  2:573961/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AAGTCTGGTATTCGCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGTTACC  >  NC_000913/794790‑795087

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: