Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 777,208 A→C 6.7% K289N (AAA→AAC)  tolA → membrane anchored protein in TolA‑TolQ‑TolR complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913777,2080AC 6.7% 1159.9 / 11.1 412K289N (AAA→AAC) tolAmembrane anchored protein in TolA‑TolQ‑TolR complex
Reads supporting (aligned to +/- strand):  ref base A (281/91);  new base C (4/35);  total (286/126)
Fisher's exact test for biased strand distribution p-value = 7.15e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 424 total aligned reads displayed.

GCAGAGAAAGCTGCAGCCGACAAAAAAGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAGGGGCCGATATCAAT  >  NC_000913/777059‑777355
                                                                                                                                                     |                                                                                                                                                   
gcagAGAAAGCTGAAGCCGACAAAAAAGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCaaaa                                                                                                                                                     <  2:3883636/150‑1 (MQ=255)
      aaaGCTGCAGCCGACAAAAAAGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGcc                                                                                                                                               >  2:6051406/1‑150 (MQ=255)
          cTGCAGCCGACAAAAAAGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCag                                                                                                                                           >  1:4601596/1‑150 (MQ=255)
           tGCAGCCGACAAAAAAGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCACCAGCGGCCGCaga                                                                                                                                          >  2:4888269/1‑150 (MQ=255)
           tGCAGCCGACAAAAAAGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCaga                                                                                                                                          >  1:9065684/1‑150 (MQ=255)
           tGCAGCCGACAAAAAAGCAGCAGAAAAAGCGGCTGCTGAAAACGCAGCAGCTGATAAGAAAGCGGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCACAAGCGGCCGCaga                                                                                                                                          >  2:8757432/1‑150 (MQ=255)
              aGCCGACACAAAAGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTTATACGAACGCAGCGGCAGAAAAAGCCGCCGCTGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTTAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGc                                                                                                                                       <  2:8080239/150‑1 (MQ=38)
               gCCGACAAAAAAGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCa                                                                                                                                      >  1:5666981/1‑150 (MQ=255)
               gCCGACAAAAAAGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCa                                                                                                                                      >  1:1005577/1‑150 (MQ=255)
                  gACAAAAAAGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAgat                                                                                                                                   >  1:7185605/1‑150 (MQ=255)
                     aaaaaaGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAgatgat                                                                                                                                >  1:7425681/1‑150 (MQ=255)
                     aaaaaaGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAgatgat                                                                                                                                >  1:4030511/1‑150 (MQ=255)
                       aaaaGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGatat                                                                                                                              >  2:8144057/1‑150 (MQ=255)
                       aaaaGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGatat                                                                                                                              >  2:6671779/1‑150 (MQ=255)
                           gcagcaGAAAACGCGGCTGCTGAAAAGGCAGCTGCTTATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTTAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTc                                                                                                                          <  2:2111726/150‑1 (MQ=255)
                             agcagAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAATGCGGCCGCAGAGGCAGATGAGATTTTCgg                                                                                                                        >  2:5839508/1‑150 (MQ=255)
                              gcagAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGt                                                                                                                       >  1:6589400/1‑150 (MQ=255)
                              gcagAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGt                                                                                                                       >  1:8942359/1‑150 (MQ=255)
                                agAAAACGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGa                                                                                                                     <  2:6142570/150‑1 (MQ=255)
                                agAAAAAGCGGCTGCTGAAAAGGCCGCAGCTGAGAAGAAAGCAGCGGCCGAAAAAGCCGCCGCAGACCAAAAAGCGGCAGCGGCAAAAGCTGCAGCGGAAAAAGCCGCTGCAGCACAAGCGGCCGCAGAGGCAGATGATCTTTTCGGTGa                                                                                                                     >  2:2857856/1‑150 (MQ=38)
                                  aaaaaGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGc                                                                                                                   >  1:1542038/1‑150 (MQ=255)
                                           ctgctgAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTg                                                                                                          >  1:3250146/1‑150 (MQ=255)
                                            tgctgAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTgg                                                                                                         >  2:3972782/1‑150 (MQ=255)
                                            tgctgAAAAGGCAGATGCTTATAATAACGATGGGGATTAAAAAGCCGCCGAATACAAAAACGCGGCTGCGGCAAACGCTGCAGGTTAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTgg                                                                                                         <  2:3416362/150‑1 (MQ=255)
                                              ctgAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTa                                                                                                       >  2:626033/1‑150 (MQ=255)
                                                gAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAg                                                                                                     >  2:6579/1‑150 (MQ=255)
                                                gAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATATTCGGTGAGCTAAGCTCTGGTAAg                                                                                                     >  1:7589480/1‑150 (MQ=255)
                                                 aaaaGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGa                                                                                                    >  2:2730836/1‑150 (MQ=255)
                                                   aaGGCAGTAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCCGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGGAAGAAt                                                                                                  >  2:3092530/1‑150 (MQ=255)
                                                   aaGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAAt                                                                                                  >  2:1887840/1‑150 (MQ=255)
                                                   aaGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAAt                                                                                                  <  2:2291362/150‑1 (MQ=255)
                                                      gcagcaGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCa                                                                                               >  2:3791996/1‑150 (MQ=255)
                                                         gcagcTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCg                                                                                            >  1:2982524/1‑150 (MQ=255)
                                                           agcTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGaa                                                                                          >  1:8063556/1‑150 (MQ=255)
                                                                 tAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACggg                                                                                    <  2:748101/150‑1 (MQ=255)
                                                                 tAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACggg                                                                                    >  1:7668571/1‑150 (MQ=255)
                                                                  aacaacccagcGGAAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGAAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACgggg                                                                                   <  2:4290538/143‑1 (MQ=255)
                                                                    gAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGACAACgggggg                                                                                 >  2:8399548/1‑150 (MQ=255)
                                                                          agcGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGc                                                                           >  1:1016453/1‑150 (MQ=255)
                                                                            cGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGa                                                                         >  2:5447265/1‑150 (MQ=255)
                                                                             ggCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGaa                                                                        >  1:8252784/1‑150 (MQ=255)
                                                                              gCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCACAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGaaa                                                                       >  2:1221901/1‑150 (MQ=255)
                                                                              gCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGaaa                                                                       >  1:5193500/1‑150 (MQ=255)
                                                                               cAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAg                                                                      >  1:6845165/1‑150 (MQ=255)
                                                                               cAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAg                                                                      >  1:9208421/1‑150 (MQ=255)
                                                                                aGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAgg                                                                     >  1:4843246/1‑150 (MQ=255)
                                                                                aGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAgg                                                                     >  1:3819975/1‑150 (MQ=255)
                                                                                aGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAgg                                                                     >  1:2609008/1‑150 (MQ=255)
                                                                                  aaaaaGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGa                                                                   >  2:1596138/1‑150 (MQ=255)
                                                                                         cgccgcAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGc                                                            >  1:2033464/1‑150 (MQ=255)
                                                                                             gcAGACAAAAAAGCGGCAGCGGCAAACGCTGCTGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCg                                                        <  2:6424797/150‑1 (MQ=255)
                                                                                             gcAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCg                                                        >  1:3626554/1‑150 (MQ=255)
                                                                                                gACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCt                                                     >  1:9226291/1‑150 (MQ=255)
                                                                                                   aaaaaaGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGcc                                                  >  2:919656/1‑150 (MQ=255)
                                                                                                   aaaaaaGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGcc                                                  >  1:225623/1‑150 (MQ=255)
                                                                                                    aaaatGCGGCTGCGGAAAACGCTGATGCTTAAAACGCGGCTGAAGAAAACGCGGCCGCAGAGGCTGTTGATATTTTCGGTGCGCTAAGCTCTGTTACGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCg                                                 <  1:729788/150‑1 (MQ=255)
                                                                                                     aaaaGCGGCAGCGGCAAAAGCTGCAGCTTAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCgg                                                <  2:5698854/150‑1 (MQ=255)
                                                                                                     aaaaGCGGCAGCGGAAAACGCTGCAGCTCAAAACCCCGCTGCAGCAAACGCGGCCGCAGAGGAAGATGATATTTTCTGTGAGCTAAGCTCTGGTAAGAATGCACCGAACACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCgg                                                <  1:6035448/150‑1 (MQ=255)
                                                                                                       aaGCGGCTGCGGCAAACGCTGCAGCTGAAACCGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAATCGGGGGGTGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGa                                              <  2:3006056/150‑1 (MQ=255)
                                                                                                        aggggcagcggcaAAAGCTGCCGCTGAAAAAGCCGCTGCAGAAAAAGCGGCCGCAGAGGCTGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAg                                             <  2:7251751/147‑1 (MQ=255)
                                                                                                        agcggcagcggcAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAg                                             <  2:5294313/150‑1 (MQ=255)
                                                                                                          cggcaggggcaACCGCGGCCGCTGCAAACGCCGCTGCAGCAAACGCGGCCGCAGCGGCAGATGAGATTTTCGGTGAGCTACGCTCTGGTAAGAATGCCCCGACAACGGGGGGAGGGGCGAACGGGAACACTGCTTCGCCTGCCGGGAGTg                                           <  2:5536462/150‑1 (MQ=255)
                                                                                                             cagcggcaAAAGCTGCAGCTCAAAAAGCCGCTGCAGAAAACGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTa                                        <  1:5929148/150‑1 (MQ=255)
                                                                                                                cggcaAAAGCTGTAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATa                                     >  2:4456463/1‑150 (MQ=255)
                                                                                                                    aaaaGCTGCTGCTTAAAACGCCGCGGCACCAAACGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAATAAGGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTaa                                 <  1:1160666/150‑1 (MQ=255)
                                                                                                                     aaaGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTaaa                                >  2:7177080/1‑150 (MQ=255)
                                                                                                                     aaaGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTaaa                                <  1:5345681/150‑1 (MQ=255)
                                                                                                                       cgCTGCCGCTGAAAAAGCCGCGGCAGCAAACGCGGCCGCAGAGGCAGATTATTTTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTaaaaa                              <  1:8751362/149‑1 (MQ=255)
                                                                                                                          tGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACaa                           >  2:8635788/1‑150 (MQ=255)
                                                                                                                          tGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACaa                           >  2:3624258/1‑150 (MQ=255)
                                                                                                                          ggctgCTCAAAACGCCGCGGCAGCAACCGCGGCCGCAGAGGCTGATGATATTTTGGGTGAGATAAGCTCTGGTAAGAATGCACCGAACACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACaa                           <  1:1924271/146‑1 (MQ=255)
                                                                                                                            cAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGTTGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGCAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATg                         >  2:7365499/1‑150 (MQ=255)
                                                                                                                             aGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATgg                        >  2:462815/1‑150 (MQ=255)
                                                                                                                              gCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGgc                       <  2:5005489/150‑1 (MQ=255)
                                                                                                                              gCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGgc                       >  2:6933409/1‑150 (MQ=255)
                                                                                                                               cTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGgcg                      >  2:7871267/1‑150 (MQ=255)
                                                                                                                                tGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGgcgc                     >  1:5115580/1‑150 (MQ=255)
                                                                                                                                tGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGgcgc                     >  1:6300323/1‑150 (MQ=255)
                                                                                                                                 gAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTACTACTAAAAACAATGGCGCa                    >  1:8542648/1‑150 (MQ=255)
                                                                                                                                 gAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCa                    >  2:2570715/1‑150 (MQ=255)
                                                                                                                                  aaaaaGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCAt                   >  1:4118191/1‑150 (MQ=255)
                                                                                                                                  aaaaaGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCAt                   >  1:666889/1‑150 (MQ=255)
                                                                                                                                   aaaaGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATc                  >  1:2691672/1‑150 (MQ=255)
                                                                                                                                    aaaGCCGCTGCAGCAACAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCa                 <  1:1682478/150‑1 (MQ=255)
                                                                                                                                    aaaGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCa                 >  1:7862517/1‑150 (MQ=255)
                                                                                                                                    aaaGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCa                 >  2:3282075/1‑150 (MQ=255)
                                                                                                                                       gCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAggg              >  2:1952589/1‑150 (MQ=255)
                                                                                                                                       gCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAggg              >  2:8248487/1‑150 (MQ=255)
                                                                                                                                         cGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAGGGGc            >  1:5570624/1‑150 (MQ=255)
                                                                                                                                           cTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAGGGGCCg          <  2:7544158/150‑1 (MQ=255)
                                                                                                                                           cTGCAGAAAACGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAGGGGCCg          <  2:239665/150‑1 (MQ=255)
                                                                                                                                            tGCAGCAAAAGCGGCCGCAGCGGCAGAGGATATTTGCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAGGGGCCGa         <  1:1080541/150‑1 (MQ=255)
                                                                                                                                             gcagcaAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAGGGGCCGat        >  2:4124346/1‑150 (MQ=255)
                                                                                                                                              cagcaAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCATCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAGGGGCCGata       >  1:3434170/1‑150 (MQ=255)
                                                                                                                                              cagcaAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAGGGGCCGata       <  2:1027280/150‑1 (MQ=255)
                                                                                                                                               agcaAACGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAGGGGCCGatat      <  1:1240659/150‑1 (MQ=255)
                                                                                                                                               agcaAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAGGGGCCGatat      >  1:4337243/1‑150 (MQ=255)
                                                                                                                                                gcaAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAGGGGCCGATATc     >  2:9071435/1‑150 (MQ=255)
                                                                                                                                                 caAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAGGGGCCGATATCa    >  1:5023703/1‑150 (MQ=255)
                                                                                                                                                   aaaGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAGGGGCCGATATCAAt  <  2:9009210/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCAGAGAAAGCTGCAGCCGACAAAAAAGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAGGGGCCGATATCAAT  >  NC_000913/777059‑777355

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑