Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 794,943 T→G 6.1% 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 6.1% 367.2 / 11.4 130intergenic (+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 (59/63);  new base G (4/4);  total (63/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.62e-01

Only 100 of 131 total aligned reads displayed.

TCGCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGTTACCAATG  >  NC_000913/794801‑795091
                                                                                                                                               |                                                                                                                                                    
tCGCCGTAATCATGGTACCTGTCGTGCTTGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAt                                                                                                                                               <  2:182281/150‑1 (MQ=255)
  gCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAttt                                                                                                                                             <  1:813898/150‑1 (MQ=255)
   ccGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGACAtttt                                                                                                                                            >  2:1211285/1‑150 (MQ=255)
   ccGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGAATTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAACGCCCGACTAAGCGGGCAttca                                                                                                                                            <  2:468665/150‑3 (MQ=255)
    cGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTTAAAAAGACCCGCCCGGACAAAATCATTCATTCCCTGAATGCCAGCTTAGTCGGGCAttttc                                                                                                                                           >  1:277916/1‑150 (MQ=255)
    cGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAttttc                                                                                                                                           >  1:1512615/1‑150 (MQ=255)
       aaTCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAttttcttt                                                                                                                                        <  1:982319/150‑1 (MQ=255)
        ttCATGGTAGGTGTCGTGATGGCCATCATCCTGGTTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAATATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAttttctttt                                                                                                                                       <  2:1128112/149‑1 (MQ=255)
        aTCATGGTACCTGACTTGATGGCCATCCTCCTGGGTCTGAGTTTCGGTCTTGGTGACGTATTAAACATCTTTTCTGGTGGTGGTAAAAAAGAACAGCCCGGCCTAAATCCTTGATTCCCTGAATGCCCGCTTAGTCGGGCAttttctttt                                                                                                                                       <  1:345801/150‑1 (MQ=255)
          cATGGTACCTGGCGTGATGGATATCATCCTGGGTCTGATTTACTGTCTTGGTGAATTATTCAAATTCTTTTCTGGTGTTGGTAAAAAAGACCAGACCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAttttctttttc                                                                                                                                     <  2:691030/150‑1 (MQ=255)
               tACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAAc                                                                                                                                <  1:1566244/150‑1 (MQ=255)
                aCCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACt                                                                                                                               <  1:62119/150‑1 (MQ=255)
                 ccTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACtt                                                                                                                              <  2:140430/150‑1 (MQ=255)
                  cTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTc                                                                                                                             >  2:1619517/1‑150 (MQ=255)
                   tGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTcc                                                                                                                            >  1:754297/1‑150 (MQ=255)
                     tCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTg                                                                                                                          >  2:1445164/1‑150 (MQ=255)
                     gcGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTg                                                                                                                          <  2:6448/149‑1 (MQ=255)
                      cGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGc                                                                                                                         >  1:763640/1‑150 (MQ=255)
                       gTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCt                                                                                                                        >  2:1687813/1‑150 (MQ=255)
                       gTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCt                                                                                                                        >  2:1373863/1‑150 (MQ=255)
                       gTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCt                                                                                                                        >  2:690190/1‑150 (MQ=255)
                        tGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTACCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTACTCAACTTACTGCtt                                                                                                                       >  1:248112/1‑150 (MQ=255)
                             gCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCt                                                                                                                  <  1:1654026/150‑1 (MQ=255)
                                atcatcCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGcc                                                                                                               <  2:710043/150‑1 (MQ=255)
                                atcatcCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGcc                                                                                                               <  2:831607/150‑1 (MQ=255)
                                atcatcCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGcc                                                                                                               <  1:1208728/150‑1 (MQ=255)
                                           tCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTc                                                                                                    <  2:1036871/150‑1 (MQ=255)
                                            cTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGTTGTTGGTAAAAAAGACCAGCCCTGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCATCTTCCTGCTTTTCCTGCCGATATTTTTTCt                                                                                                   <  1:462548/150‑1 (MQ=255)
                                             tGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCtt                                                                                                  >  1:548234/1‑150 (MQ=255)
                                             ggATTTACGGTGTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAATACCAGCCAGGACTAAATCTTTGATTCCCTGAATGCCCGCTTAGTCGTGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCtt                                                                                                  <  2:273743/149‑1 (MQ=255)
                                              gATTTACTGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCACGGACAAAATCATTGATCAACTGAATGCCCGCATAGTAGGGCATTTTCTTTTTCTCAGCTTCCGGCTGTTACTGCCGATATTTTGTCtat                                                                                                 >  2:1252266/1‑148 (MQ=255)
                                              gATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCTGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCtt                                                                                                  >  2:1399706/1‑150 (MQ=255)
                                                     ggTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCt                                                                                          >  2:1205656/1‑150 (MQ=255)
                                                      gTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGAGCATTTTCTTTTTCTCAACTTCCTGATTTACATGCCGATATTTTTTCTTATCTACCta                                                                                         >  1:13143/1‑149 (MQ=255)
                                                        cTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTcac                                                                                       <  1:1057534/150‑1 (MQ=255)
                                                        cTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTACCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTcac                                                                                       >  2:1442146/1‑150 (MQ=255)
                                                          tGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACaa                                                                                     >  1:556666/1‑150 (MQ=255)
                                                          tGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACaaa                                                                                    >  2:1556716/1‑150 (MQ=255)
                                                             tGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAgg                                                                                  <  1:557296/150‑1 (MQ=255)
                                                             tGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGt                                                                                 >  2:278624/1‑150 (MQ=255)
                                                             tGAAGTATTCAAAATCTTTTCTGGTGTTGGTAATAAAGACCAGCCGGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAgg                                                                                  <  1:246422/150‑1 (MQ=255)
                                                              gAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGt                                                                                 >  2:1133722/1‑150 (MQ=255)
                                                               aaGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGtt                                                                                <  1:1110926/150‑1 (MQ=255)
                                                                  tATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGc                                                                             <  2:475978/150‑1 (MQ=255)
                                                                  tATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGGTAGc                                                                             <  2:418940/150‑1 (MQ=255)
                                                                   aTTCAACATCTTTTCTGGTTTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCa                                                                            <  2:1284265/150‑1 (MQ=255)
                                                                      cAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATa                                                                         >  2:299709/1‑150 (MQ=255)
                                                                      cAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATa                                                                         >  2:1258613/1‑150 (MQ=255)
                                                                      cAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATa                                                                         >  2:462357/1‑150 (MQ=255)
                                                                       aaCATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATaa                                                                        >  2:648294/1‑150 (MQ=255)
                                                                         cATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGCTTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCACCTTCCGGCGTTTCCTGCCGAGATTTTTTCTGATCTCCCTCACATAGGTTAGCAATACct                                                                      >  2:435719/1‑150 (MQ=255)
                                                                         cATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTTAATGCCCGCTTAGTCGTGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAAct                                                                      <  1:1140250/150‑1 (MQ=255)
                                                                         cATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAAct                                                                      >  1:509045/1‑150 (MQ=255)
                                                                           tCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTTGCAATAActgc                                                                    >  1:1484476/1‑150 (MQ=255)
                                                                            cTTTTCTCGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGAATCCCTGAATGCCCGCTTAGGCGGGCATTTTCTTTTTCTCTACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAActgct                                                                   <  2:16170/150‑1 (MQ=255)
                                                                             ttttCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAActgctg                                                                  <  1:1484264/150‑1 (MQ=255)
                                                                              tttCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTgg                                                                 <  1:1514031/150‑1 (MQ=255)
                                                                                tCTGGTGTTTGTAAAAAAGACCAGCCCGTACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGa                                                               <  1:1137903/150‑1 (MQ=255)
                                                                                tCTGGTGTTGGTAAAAAAGACCAGCCCGGGCAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGa                                                               >  2:1610632/1‑150 (MQ=255)
                                                                                 cTGGTGTTGGTAAAAAAGCCCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGaa                                                              >  2:40285/1‑150 (MQ=255)
                                                                                    gtgtTGGTTAAAAAGACCAGCCCGGACTAAATGATTGATTCCCTGAATGCCCGCTTAGTCTGGCATTTTCTTTATCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACATAGGTTTGCAATAACTGCTGGGAAAAt                                                           <  1:1157975/150‑1 (MQ=255)
                                                                                    gtgtTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAAt                                                           >  2:645985/1‑150 (MQ=255)
                                                                                    gtgtTGGTAAAAAAGACCAGACCGGACAAAATCATTGATTCCCTGAATGCCCGCTTATTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAAt                                                           <  2:1079994/150‑1 (MQ=255)
                                                                                     tgtTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAAtt                                                          >  1:1346949/1‑150 (MQ=255)
                                                                                      gtTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATGCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTc                                                         <  1:1564794/150‑1 (MQ=255)
                                                                                       ttGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTcc                                                        <  2:1512615/150‑1 (MQ=255)
                                                                                        tGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCg                                                       <  2:1442579/150‑1 (MQ=255)
                                                                                        tGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCg                                                       >  2:584151/1‑150 (MQ=255)
                                                                                         ggTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGa                                                      >  2:1638614/1‑150 (MQ=255)
                                                                                          gTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCGGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAg                                                     >  2:15223/1‑150 (MQ=255)
                                                                                           tAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTTCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGt                                                    >  2:1058842/1‑150 (MQ=255)
                                                                                             aaaaaGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTa                                                  <  1:1435598/150‑1 (MQ=255)
                                                                                              aaaaGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAg                                                 >  2:1402618/1‑150 (MQ=255)
                                                                                                 aGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCg                                              <  2:413976/150‑1 (MQ=255)
                                                                                                 aGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCg                                              <  2:1261596/150‑1 (MQ=255)
                                                                                                  gACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGt                                             >  2:1562390/1‑150 (MQ=255)
                                                                                                    ccAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTta                                           <  2:763640/150‑1 (MQ=255)
                                                                                                         ccGGACAATATCTTTGAATCACTGAATGCCCGATTAGTCGGGCATTTTATTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCAACCTCACAAAGGTTAGCATTAACTGCTGGGAAAATTCCGAGGTAGTCGTTATATTg                                      <  1:933113/150‑1 (MQ=255)
                                                                                                             gACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCCCAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCg                                   >  1:1350300/1‑150 (MQ=255)
                                                                                                              aCAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGc                                  <  2:550560/150‑1 (MQ=255)
                                                                                                              aCAAAATCATAGATTCCCTGAATGCTCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGc                                  <  1:759167/150‑1 (MQ=255)
                                                                                                                        tGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACg                        <  1:1699362/150‑1 (MQ=255)
                                                                                                                         gATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGt                       >  1:1032821/1‑150 (MQ=255)
                                                                                                                             aagaaaaTGCCCGACTAAGCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGCGTTAGTCGTTATATTGTCGCCTACATAACGTTACa                   >  2:1172350/6‑150 (MQ=255)
                                                                                                                              ccTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACAt                  >  2:819334/1‑150 (MQ=255)
                                                                                                                              ccTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACAt                  >  2:1013400/1‑150 (MQ=255)
                                                                                                                              agaaaaTGCCCGACTAAGCGGGCATTTTCTTTTTCTCTACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTACGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACAt                  >  1:1039925/5‑150 (MQ=255)
                                                                                                                              agaaaaTGCCCGACTAAGCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACAt                  >  1:407966/5‑150 (MQ=255)
                                                                                                                              agaaaaTGCCCGACTAAGCGGGCATTCTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCGCAAAGGTTAGCAATTACTGCTGGGATATGTGCGAGTTAGTGGTTATATTGTCGCCTACATAACGTTGCAt                  >  1:251138/5‑150 (MQ=255)
                                                                                                                               cTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACAtt                 <  2:1466646/150‑1 (MQ=255)
                                                                                                                                aaaaTGCCCGCTTAGTCGGGCATATTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTa                >  1:322692/3‑150 (MQ=255)
                                                                                                                                   aTGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAgg             <  1:1533586/150‑1 (MQ=255)
                                                                                                                                   aTGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAgg             <  1:1580900/150‑1 (MQ=255)
                                                                                                                                      cccGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGt          >  1:1034351/1‑150 (MQ=255)
                                                                                                                                       gcGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGtt         >  2:1237512/2‑150 (MQ=255)
                                                                                                                                        cGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGTTa        <  2:543506/150‑1 (MQ=255)
                                                                                                                                             attCGGGCTTTTTCTTTTTCTCAACTTCCTCCTTTTCCTGCCGATATATTTTCTTATCTACCTCACAAAGGGTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGTTACCAAt   <  2:287953/148‑1 (MQ=255)
                                                                                                                                             aGTCGGACTTTTTCTTTTTATCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCAAAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGTTACCAAt   <  1:1211285/150‑1 (MQ=255)
                                                                                                                                              gTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGTTACCAATg  >  1:1670098/1‑150 (MQ=255)
                                                                                                                                              gTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGTTACCAATg  <  1:1450903/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
TCGCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGTTACCAATG  >  NC_000913/794801‑795091

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: