Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 794,943 T→G 5.5% 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 5.5% 370.6 / 3.2 126intergenic (+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/55);  new base G (3/6);  total (65/61)
Fisher's exact test for biased strand distribution p-value = 3.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 128 total aligned reads displayed.

GGTATTCGCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGTTACCAATG  >  NC_000913/794796‑795091
                                                                                                                                                   |                                                                                                                                                    
  tATTCGCCGTAATCATGGTTCCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAACTCATTGATTCCCTGAATGCCCGCTTAGTCggg                                                                                                                                                  <  1:1284547/150‑1 (MQ=255)
      cGCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAtt                                                                                                                                              >  2:1051897/1‑150 (MQ=255)
         cGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGACTAAGCGGGCAttcag                                                                                                                                           <  1:1951459/150‑4 (MQ=255)
          cgtaTCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGACTAAGCGGGCAttcagg                                                                                                                                          <  2:650907/147‑5 (MQ=255)
              tCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATTCCCGCTTAGTCGGGCAttttcttttt                                                                                                                                      >  1:1824324/1‑150 (MQ=255)
               cATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAttttctttttc                                                                                                                                     <  1:418360/150‑1 (MQ=255)
               cATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAttttctttttc                                                                                                                                     <  2:775631/150‑1 (MQ=255)
               cATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAttttctttttc                                                                                                                                     <  1:931689/150‑1 (MQ=255)
               cATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAttttcttttta                                                                                                                                     >  1:1069480/1‑149 (MQ=255)
                aTGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAttttctttttct                                                                                                                                    >  2:733197/1‑150 (MQ=255)
                aTGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCACTGAATGCCCGCTTAGTCGGGCAttttctttttct                                                                                                                                    >  2:608749/1‑150 (MQ=255)
                  ggTACCTTTCGGGATGGCAATCATCCTGCGTCTTATTTACTGTCTTTGTGAAGTATTCAACATCATTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCAAAGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCa                                                                                                                                  <  1:298359/150‑1 (MQ=255)
                      ccTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACtt                                                                                                                              >  1:1998144/1‑150 (MQ=255)
                        tGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCAATGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTcc                                                                                                                            <  1:436595/150‑1 (MQ=255)
                         gTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCt                                                                                                                           <  2:1715962/150‑1 (MQ=255)
                         gTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCt                                                                                                                           <  1:858958/150‑1 (MQ=255)
                                 ggCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCATGCTTTTcc                                                                                                                   >  2:908555/1‑150 (MQ=255)
                                   ccATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTg                                                                                                                 >  2:162492/1‑150 (MQ=255)
                                     atcatcCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGcc                                                                                                               <  2:505544/150‑1 (MQ=255)
                                       catcCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAACTTCATTGATTCCCTGAATGCCCGGTTAGTCGGGCTTTTTCGTTTTCTCAACTTCCTGCTTTTCCTGCCGa                                                                                                             >  2:87022/1‑150 (MQ=255)
                                        atcCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCCGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCACCTTCCTGCGTTTCCTGCCGat                                                                                                            >  2:1107747/1‑150 (MQ=255)
                                             gggTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATAtttt                                                                                                       >  2:1632526/1‑150 (MQ=255)
                                              ggTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATAttttt                                                                                                      >  1:1865411/1‑150 (MQ=255)
                                               gTCTTATTTACGGTCTTGTTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAATAGACCAGCCCGGACAAAATCATTGATTCCCTGTATGCCCGCTTATGCGGGCTTTTTCTTTTTCTGAACTTCCTGCTTTTCCTGCCGATAtttttt                                                                                                     <  2:1005519/150‑1 (MQ=255)
                                               gTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATAtttttt                                                                                                     >  2:1824212/1‑150 (MQ=255)
                                                tCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTc                                                                                                    <  2:1357130/150‑1 (MQ=255)
                                                 cTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCt                                                                                                   <  2:365856/150‑1 (MQ=255)
                                                   gATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTa                                                                                                 >  2:352204/1‑150 (MQ=255)
                                                   gATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTa                                                                                                 <  1:1348647/150‑1 (MQ=255)
                                                     tttACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATc                                                                                               >  2:1406455/1‑150 (MQ=255)
                                                         cGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTAcc                                                                                           >  2:916054/1‑150 (MQ=255)
                                                          ggTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCt                                                                                          >  1:392967/1‑150 (MQ=255)
                                                           gTCTTGGTGATGTATTCAACATCTTTTCTGGTGTTGGTTAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTATCCTGCCGATATTTTTTCTTATCTACCTc                                                                                         <  1:659085/150‑1 (MQ=255)
                                                                 gTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAg                                                                                   >  2:1079962/1‑150 (MQ=255)
                                                                   gAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCTTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGt                                                                                 >  2:41086/1‑150 (MQ=255)
                                                                    aaGTATTCAACATGTTTTCTGGTGTTGGGAAAAAAGCCCAGCCCGGACAAAATCATTGATTCACTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGtt                                                                                <  2:1793119/150‑1 (MQ=255)
                                                                    aaGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGtt                                                                                <  1:859347/150‑1 (MQ=255)
                                                                     aGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAAGTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTa                                                                               <  2:802897/150‑1 (MQ=255)
                                                                     aGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTa                                                                               >  2:260552/1‑150 (MQ=255)
                                                                       tATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGc                                                                             <  1:935717/150‑1 (MQ=255)
                                                                        aTTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCa                                                                            <  1:1588635/150‑1 (MQ=255)
                                                                        aTTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCAGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCa                                                                            <  1:397287/150‑1 (MQ=255)
                                                                         ttCAACATCTTTTCTGGTGTTGGTAAAAAAGACAAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCaa                                                                           >  1:1311010/1‑150 (MQ=255)
                                                                           cAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATa                                                                         >  1:1868519/1‑150 (MQ=255)
                                                                           cAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATa                                                                         >  2:31623/1‑150 (MQ=255)
                                                                             aCATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTTATTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAAc                                                                       >  1:1830280/1‑150 (MQ=255)
                                                                              cATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAAct                                                                      >  1:1708250/1‑150 (MQ=255)
                                                                              cATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAAct                                                                      >  1:1538338/1‑150 (MQ=255)
                                                                              cATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAAct                                                                      >  1:467980/1‑150 (MQ=255)
                                                                                 cTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAActgct                                                                   >  1:1540067/1‑150 (MQ=255)
                                                                                   tttCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACTAAGGTTAGCAATAACTGCTgg                                                                 <  2:1947625/150‑1 (MQ=255)
                                                                                    ttCTGGTGTTTGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTggg                                                                <  1:546552/150‑1 (MQ=255)
                                                                                     tCTTGTGTTGGTAATAAAGTCCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCTACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGa                                                               <  2:1912450/150‑1 (MQ=255)
                                                                                     tCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGa                                                               <  2:126657/150‑1 (MQ=255)
                                                                                      cTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGaa                                                              >  2:1197816/1‑150 (MQ=255)
                                                                                      cTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTGTTTCTCAACTTCCTGCTTATCCTGCCGATAATTTTTCTTATCTACCTCACAAAGGTTAGCAAAAACTGCTGGGaa                                                              >  2:428746/1‑150 (MQ=255)
                                                                                        ggTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGaaaa                                                            >  1:1235003/1‑150 (MQ=255)
                                                                                        ggTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGaaaa                                                            >  2:1773753/1‑150 (MQ=255)
                                                                                        ggTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGaaaa                                                            >  1:1403915/1‑150 (MQ=255)
                                                                                        ggTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATATTCTTTTTCTGAACTTCGTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTTGCAATAACTGCTGGGaaaa                                                            >  1:1024631/1‑150 (MQ=255)
                                                                                          tgtTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAAtt                                                          <  2:1433060/150‑1 (MQ=255)
                                                                                            ttGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTcc                                                        <  1:352204/150‑1 (MQ=255)
                                                                                             ttgTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCTGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCg                                                       <  2:1680373/148‑1 (MQ=255)
                                                                                               gTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAg                                                     >  2:1851194/1‑150 (MQ=255)
                                                                                               gTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAg                                                     >  1:479693/1‑150 (MQ=255)
                                                                                                tAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGTAAAATTACGAGt                                                    >  2:110686/1‑150 (MQ=255)
                                                                                                tAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGt                                                    >  1:1529727/1‑150 (MQ=255)
                                                                                                 aaaaaaGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGtt                                                   >  1:1663267/1‑150 (MQ=255)
                                                                                                  aaaaaGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCGCAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTa                                                  >  2:1671186/1‑150 (MQ=255)
                                                                                                      aGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAAAGGCTGGGAAAATTACGAGTTATTCg                                              >  1:1296168/1‑150 (MQ=255)
                                                                                                          cAGCCCGTACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTtat                                          <  1:1632526/150‑1 (MQ=255)
                                                                                                          cAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTtat                                          <  1:1051897/150‑1 (MQ=255)
                                                                                                            tcccTGACAAAATCATTGATTCCCTGAATGCCAGCTTAGTGGGGCATTTTCTTTTTCTCTACTTCCTGCTTTTCATGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTatat                                        <  1:1009850/149‑1 (MQ=255)
                                                                                                            gCCCGGACAAAATCATTGATTCCCGGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTatat                                        <  1:271516/150‑1 (MQ=255)
                                                                                                            gCCAGGACAAAATCATTGATTCCCTGAATGCCCGTTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTatat                                        <  1:339621/150‑1 (MQ=255)
                                                                                                               cGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCTTTTTCTTTTTCTCAACTTCCTGCTTTTCTTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGt                                     <  1:708522/150‑1 (MQ=255)
                                                                                                                 taCAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCg                                   <  1:1824212/149‑1 (MQ=255)
                                                                                                                   cAAAATCATTGATTCGCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGcc                                 <  1:179308/150‑1 (MQ=255)
                                                                                                                   cAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGcc                                 >  2:1686248/1‑150 (MQ=255)
                                                                                                                   cAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGcc                                 <  2:1930326/150‑1 (MQ=255)
                                                                                                                   cAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGcc                                 >  1:1881728/1‑150 (MQ=255)
                                                                                                                      aaTCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTAc                              <  2:1865411/150‑1 (MQ=255)
                                                                                                                         cATTGATTCCCTGAATGCCCGCTTATGCGGGAATTTTCTTTTTCTCAACTTCCTGCTTTTACTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATa                           <  2:1189999/150‑1 (MQ=255)
                                                                                                                             gATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGt                       >  1:800324/1‑150 (MQ=255)
                                                                                                                              aTTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACGTCACAATGATTAGAAATAAATGCTGGGAAAATTCAGATTTTGTCGTTAAATTTTCGCATACATAACttt                      >  1:263371/1‑150 (MQ=255)
                                                                                                                               ttCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTa                     >  2:1559221/1‑150 (MQ=255)
                                                                                                                                tcacTGAATGCCAGCTTAGTCGTGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCTATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTACGAGTTAGTCGTTATATTGTCGCCTACATAACGTTAc                    <  1:1001234/147‑1 (MQ=255)
                                                                                                                                tCCCTGAATGCCCGCTTAGTCGGGCATTTTATTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTAc                    <  1:1007960/150‑1 (MQ=255)
                                                                                                                                aaagaaaaTGCCCGACTAAGCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTAc                    >  2:1677357/7‑150 (MQ=255)
                                                                                                                                 cccTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACa                   >  2:1144463/1‑150 (MQ=255)
                                                                                                                                 cccTGAATGCACGATTAGGCGGGCATTTTCGTTTTCTCAACTTCATGCTTTTCCTTCCGATATTTATTCATATCTACCGCACAAAGGATAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACa                   <  2:372322/150‑1 (MQ=255)
                                                                                                                                   gaaaaTGCCCGACTAAGCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACAtt                 >  1:1129807/4‑150 (MQ=255)
                                                                                                                                   cTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACAtt                 >  2:934627/1‑150 (MQ=255)
                                                                                                                                     gAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTaa               >  1:107165/1‑150 (MQ=255)
                                                                                                                                       aTGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAgg             <  1:2017454/150‑1 (MQ=255)
                                                                                                                                              cTTAGTCGGGTATTTTCTTTTTCTCAACTTCCTGCTTTTGCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTATCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGTTAcc      <  1:80185/150‑1 (MQ=255)
                                                                                                                                                tAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTCAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGTTACCaa    >  2:1980610/1‑150 (MQ=255)
                                                                                                                                                  ttCGGTCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGGTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGTTACCAATg  <  2:1849182/149‑1 (MQ=255)
                                                                                                                                                  gTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGTTACCAATg  <  1:825432/150‑1 (MQ=255)
                                                                                                                                                  gTCGGGCATTTTCTTTTTCTCAACTTCCTGATTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAATATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGTTCCCAATg  >  2:1375601/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGTATTCGCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGTTACCAATG  >  NC_000913/794796‑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: