Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,130,163 A→G 9.6% intergenic (+140/‑42) PP_4511 → / → PP_4512 AraC family transcriptional regulator/DMT superfamily permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,130,1630AG 9.6% 618.7 / 39.4 209intergenic (+140/‑42)PP_4511/PP_4512AraC family transcriptional regulator/DMT superfamily permease
Reads supporting (aligned to +/- strand):  ref base A (110/79);  new base G (17/3);  total (127/82)
Fisher's exact test for biased strand distribution p-value = 2.80e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.02e-01

Only 100 of 211 total aligned reads displayed.

CTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTC  >  NC_002947/5130020‑5130312
                                                                                                                                                |                                                                                                                                                     
cTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAAACCGCTCCAACAGGATGCCCCCGCCGGAAAGTTGTGCAAAACCGGTAGGAGCGGGTTTACCCGCGCAAGGGCCg                                                                                                                                                 >  2:5853943/1‑150 (MQ=39)
 tGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGGAAAACCGCGGCAAAAGGATGCCCCCCCCAGAAATTTTTTGAAAAACGGTAGGAGCGGGTTTTACCCCGAAACGGCCgg                                                                                                                                                >  1:1084008/1‑150 (MQ=255)
      gcagTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGGAAAACCGCGCCTAAAGGCTGCCCCCCCCGGAAGGTTGTGCAAAAACAGGAGGTGCGGGGTTTCCCCCCCAAAGGCCCgcccgg                                                                                                                                            >  2:2109263/1‑145 (MQ=255)
       cagACTTGCCACATAATTTTTTTGCATGAGCCTTAACAATTTGCATTTCCCGGCCCCTACGCGGGTAAACCCGCTCAAACAGGATCCCCCCCCCTTAATTTTGTCAAAAACCTCTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGc                                                                                                                                          <  2:7175650/150‑1 (MQ=19)
         gTCTTGCGACTGATTGTGTGTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTGCGCGGGTAAACCCGCTCCTCCAGGATGCCCCCGCCTGTCTGTTGTGCAAACCCTGTCGGCGGGGGTTTACCCGCGAAAGGGCCGGCACTGCCt                                                                                                                                        <  2:6652561/150‑1 (MQ=37)
         gTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCt                                                                                                                                        <  2:6036796/150‑1 (MQ=255)
              gCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGATTGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTagt                                                                                                                                   <  2:5160009/150‑2 (MQ=255)
                 aCAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAAACCGCACCTACACGATTCCCCCCCCGGAATGTTGTGCAAAACCAGTAGGAGCGGGTTTAACCGCCAAAGGGCCGGCACAGCCGGTAGAGAt                                                                                                                                >  2:2737319/1‑150 (MQ=255)
                          gtTTGCCTGCGGCTTTACACTTGGCCTTGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACCGGATGCCCCCGCCTGATTGTTGTGAAAAACCTGTAGGAGCGGGTTAACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCaa                                                                                                                       <  2:2711939/150‑1 (MQ=38)
                                     ggtttACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCTCCTTATTGTTGTGCAAAACCTGTAGGTGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCAt                                                                                                            <  2:6119209/148‑1 (MQ=255)
                                       tttACACTGTGCCTTTGCCGGCCCCTTCGCGGTTAAACCCGATCCTACAGGATGCCCCCGCCTTATTGTTGTGCAAAACCTGTAGGTGGGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGa                                                                                                          <  2:3383034/150‑1 (MQ=38)
                                        ttACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAAACCGCGCCTACAGGATTCCCCCCCCAGAAAGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCCCGAAAGGGCCGGCACTGCCGGTAGAGATCTCCAACAAAGACCCACCATGaa                                                                                                         >  2:1120386/1‑150 (MQ=37)
                                         tACACTTTTCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGACTGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAac                                                                                                        <  2:5430378/150‑1 (MQ=255)
                                              tttgcctttgccGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATTCCCCCGCCTGAATGTTGTGCAAAAACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGCACTGCCGGTAGAGATCTCCCTCAAGGACCCACCCAGAACcaccg                                                                                                   >  2:6773218/1‑150 (MQ=37)
                                              tttgcctttgccGGCCCCTTCGCGGGTAAACCCCCACCAAAAGGGTGCCCCCGCCGGAATGTTTTGCAAAACCGGGAGGAGCGGGTTTACCCGCCAAAGGGCCGGCACAGCCAGTAGTGGTCTCCAACAAAGAACCACCAAGAACcaccg                                                                                                   >  2:4252877/1‑150 (MQ=255)
                                                         cGGCCCCTTCGCGGGCAACCCCGCTCCAACTGTCGCCCCCCGCCTTATTGTTGTCCAAACCTTGTTGGCGCGGGTTCACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCg                                                                                        <  1:6458402/150‑1 (MQ=255)
                                                          ggCCCCTTCGCGGGTAAACCCCCTCCAACAGGATGCCCCCCCCGGAATGTTGTGCAAAACCGGTAGGAGCGGGTTTTCCCGCGAAAGGGCCGGCCCGGCCGGGAGAGATCTCCCTCAAGGCCCCACCATGAACCACCGCACCGCCCTCgg                                                                                       >  2:4906337/1‑150 (MQ=255)
                                                             cccTTCGCGGGTAAAACCGCGCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGTAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGcccaccgccc                                                                                    >  2:7409474/1‑143 (MQ=39)
                                                              ccTTCGCGGGTAAAACCGCTCCTACACGGTGCCCCCGCCGGAAAGTTTTGGAAAACCGGTAGGAGCGGGTTTACCCGCGAAAGGGGCGGCACAGCCTGTAGAGATCTCCATCAAGGACCCCCCAAGTACCACCCCACCCCCCGCGgcgcc                                                                                   >  2:7151724/1‑149 (MQ=25)
                                                               cTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAAACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTATAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCAc                                                                                  >  1:7382287/1‑150 (MQ=255)
                                                                 TCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTTTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCACGGCGCACTg                                                                                <  2:6261984‑M1/150‑1 (MQ=255)
                                                                      GGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGTGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACAt                                                                           >  2:3510489‑M1/1‑150 (MQ=255)
                                                                        gtaaacccgctcctacaggttttgcacaacattcaggcgggggcatcctgtaggagcgggtttacccgcgaaGGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCAc                                                                                  <  2:6695218‑M2/70‑1 (MQ=255)
                                                                        gtaaacccgctcctacaggttttgcacaacattcaggcgggggcatcctgtaggagcgggtttacccgcgaaGGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCAc                                                                                  >  1:6695218‑M2/73‑142 (MQ=255)
                                                                        TAAACCCGCGCCTACCGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTCGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTg                                                                         >  2:3189156‑M1/1‑150 (MQ=255)
                                                                           ACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACCTTGgcg                                                                      <  1:6788105‑M1/150‑1 (MQ=255)
                                                                           ACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGgcg                                                                      <  1:1113944‑M1/150‑1 (MQ=255)
                                                                            CCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAATGGACGGCACTGCCTGTAGAGATCTCCATCAAGGAGCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGgcgc                                                                     <  2:5016331‑M1/150‑1 (MQ=255)
                                                                            CCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGgcgc                                                                     <  2:911198‑M1/150‑1 (MQ=255)
                                                                             CCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCCCTGCCTGTAGAGATCTCCCGCACGGACCCACCAGGAACCACCGCACCGCCCGCGGCGCACTGCACAGTGgcgcg                                                                    <  2:5760232‑M1/150‑1 (MQ=255)
                                                                             CCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGgcgcg                                                                    <  2:4266707‑M1/150‑1 (MQ=255)
                                                                              ccgctcctacaggttttgcacaacattcaggcgggggcatcctgtaggagcgggtttacccgcgaaGGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGgcgcg                                                                    >  2:6054848‑M2/67‑150 (MQ=255)
                                                                               GCTCCTCCCGGTTGCCCCCGCCTTAATGTGGTCCAAAACCTGTCGGCGGGGTTTTCCCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCt                                                                  >  1:4992108‑M1/1‑150 (MQ=255)
                                                                                gctcctacaggttttgcacaacattcaggcgggggcatcctgtaggagcgggtttacccgcgaaAGGGCCGGCACTGCCTGTAGAGATCTCCa                                                                                                                           <  1:6121218‑M2/29‑1 (MQ=255)
                                                                                gctcctacaggttttgcacaacattcaggcgggggcatcctgtaggagcgggtttacccgcgaaAGGGCCGGCACTGCCTGTAGAGATCTCCa                                                                                                                           >  2:6121218‑M2/65‑93 (MQ=255)
                                                                                gctcctacaggttttgcacaacattcaggcgggggcatcctgtaggagcgggtttacccgcgaaAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCt                                                                  >  2:1919992‑M2/65‑150 (MQ=255)
                                                                                  tcctacaggttttgcacaacattcaggcgggggcatcctgtaggagcgggtttacccgcgaaAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGt                                                                >  1:4378454‑M2/63‑150 (MQ=255)
                                                                                   tcctacaggttttgcacaacattcaggcgggggcatcctgtaggagcgggtttacccgcgaAGGGGCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGtt                                                               >  2:5513130‑M2/62‑150 (MQ=255)
                                                                                   CTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAGCCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGttc                                                              <  2:4242559‑M1/150‑1 (MQ=255)
                                                                                    TACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGttct                                                             <  1:6942209‑M1/150‑1 (MQ=255)
                                                                                    TACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGttct                                                             <  1:7063057‑M1/150‑1 (MQ=255)
                                                                                    TACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGttct                                                             <  1:4071365‑M1/150‑1 (MQ=255)
                                                                                     tacaggttttgcacaacattcaggcgggggcatcctgtaggagcgggtttacccgcgaaGGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGttct                                                             >  1:6163614‑M2/60‑150 (MQ=255)
                                                                                     tacaggttttgcacaacattcaggcgggggcatcctgtaggagcgggtttacccgcgaaGGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGttct                                                             >  1:6593798‑M2/60‑150 (MQ=255)
                                                                                     tacaggttttgcacaacattcaggcgggggaatcctgtaggagcgggtttacccgcgaaGGGGCCGGCACTGCATGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGttct                                                             >  1:6593729‑M2/60‑150 (MQ=255)
                                                                                     ACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGttctt                                                            <  2:4567154‑M1/150‑1 (MQ=255)
                                                                                      CAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGttcttc                                                           <  1:1650439‑M1/150‑1 (MQ=255)
                                                                                      CAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGttcttc                                                           <  1:4963320‑M1/150‑1 (MQ=255)
                                                                                       AGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCAGCAAGGACCCACCGTGAACCTCCGCACCGCCCTCGGCGCACTGCACTTTGGCGCGCTGTGCTTCg                                                          <  1:1371681‑M1/150‑1 (MQ=255)
                                                                                         GATGCCCCCGCCTTTATTTTGGGCCACTCTTGAAGGAGCGGCTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGc                                                        >  2:6817784‑M1/1‑150 (MQ=255)
                                                                                         GATGCCCCCGCCTGAATGTTGTGCAATACCTGTAGGAGCGGTTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGc                                                        >  2:1084008‑M1/1‑150 (MQ=255)
                                                                                         GATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGc                                                        <  1:1539547‑M1/150‑1 (MQ=255)
                                                                                         GATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCGGTTCTTCGGc                                                        <  2:4815334‑M1/150‑1 (MQ=255)
                                                                                            tttgcacaacattcaggcgggggcatcctgtaggagcgggtttacccgcgaaGGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCt                                                      >  1:6610982‑M2/53‑150 (MQ=255)
                                                                                            tttgcacaacattcaggcgggggcatcctgtaggagcgggtttacccgcgaaGGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCt                                                      >  1:5073704‑M2/53‑150 (MQ=255)
                                                                                            tttgcacaacattcaggcgggggcatcctgtaggagcgggtttacccgcgaaAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGc                                                        >  2:2561729‑M2/53‑148 (MQ=255)
                                                                                            tttgcacaacattcaggcgggggcatcctgtaggagcgggtttacccgcgaaAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGc                                                        <  1:2561626‑M2/96‑1 (MQ=255)
                                                                                             cccccGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTa                                                    >  2:7246373/1‑150 (MQ=255)
                                                                                              ccccGCCTGCATTTTGCGCCAACCCTGACGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTAc                                                   <  2:1774203/150‑1 (MQ=255)
                                                                                               gcacaacattcaggcgggggcatcctgtaggagcgggtttacccgcgaaAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTAc                                                   <  2:4378454‑M2/101‑1 (MQ=255)
                                                                                               cccGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTAcc                                                  >  1:3787067/1‑150 (MQ=255)
                                                                                               cccGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTAcc                                                  >  2:2884411/1‑150 (MQ=255)
                                                                                                ccGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCg                                                 >  1:7265703/1‑150 (MQ=255)
                                                                                                ccGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCg                                                 >  1:7381356/1‑150 (MQ=255)
                                                                                                  gCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGc                                               >  1:5477875/1‑150 (MQ=255)
                                                                                                  gCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGc                                               >  1:6885955/1‑150 (MQ=255)
                                                                                                   ccTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCg                                              >  2:3848323/1‑150 (MQ=255)
                                                                                                    cTGAATGGTGTGCAAAACCTGTAGGTGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGt                                             <  2:3883226/150‑1 (MQ=255)
                                                                                                     tGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGTACGGAACTTCCTGTATAGATCTCAATCAAGGACACAACATGAACCACCGCACCGCCATCGGCGCACTGCACATTGTCGCGATGTTCTTCGGCCTTACCGGCGTc                                            >  2:1114620/1‑150 (MQ=255)
                                                                                                      gAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCt                                           >  1:4457167/1‑150 (MQ=255)
                                                                                                        aTGTTGTGAAAAACCTGTAGGTGGGGTTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTc                                         <  1:4509215/150‑1 (MQ=255)
                                                                                                             gtgCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCaa                                    >  2:3927331/1‑150 (MQ=255)
                                                                                                              tgCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCaaa                                   >  1:5652286/1‑150 (MQ=255)
                                                                                                               gCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAAc                                  >  2:6457915/1‑150 (MQ=255)
                                                                                                                 aaaaCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCgg                                                                                       >  2:3452023/1‑95 (MQ=255)
                                                                                                                 aaaaCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCgg                                                                                       <  1:3452023/95‑1 (MQ=255)
                                                                                                                 aaaaCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTg                                >  2:2504512/1‑150 (MQ=255)
                                                                                                                  aaaCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTgg                               >  1:5934705/1‑150 (MQ=255)
                                                                                                                   aaCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGc                              >  1:4661315/1‑150 (MQ=255)
                                                                                                                    aCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGcc                             >  2:3320109/1‑150 (MQ=255)
                                                                                                                     ccctgtcggcgggggtttacccgcgaaAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGcc                             <  1:850308‑M2/123‑1 (MQ=255)
                                                                                                                       ctgtaggagcgggtttacccgcgaaAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGccac                           >  2:1385071‑M2/26‑150 (MQ=255)
                                                                                                                        tgtaggagcgggtttacccgcgaaAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGccacc                          >  1:198919‑M2/25‑150 (MQ=255)
                                                                                                                         gtaggagcgggtttacccgcgaaAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGcca                            >  1:2286171‑M2/24‑147 (MQ=255)
                                                                                                                             gagcgggtttacccgcgaaAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAg                                                                                                                      <  1:275582‑M2/34‑1 (MQ=255)
                                                                                                                             gagcgggtttacccgcgaaAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAgcgc                     >  1:5342745‑M2/20‑150 (MQ=255)
                                                                                                                               gcgggtttacccgcgaaAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCa                   >  2:4259334‑M2/18‑150 (MQ=255)
                                                                                                                               gcgggtttacccgcgaaAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCa                   >  1:4931197‑M2/18‑150 (MQ=255)
                                                                                                                               gcgggtttacccgcgaaAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCa                   >  2:1276578‑M2/18‑150 (MQ=255)
                                                                                                                                   gtttacccgcgaaAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGccc               >  2:4181349‑M2/14‑150 (MQ=255)
                                                                                                                                   gtttacccgcgaaAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACAGTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGccc               >  2:6192993‑M2/14‑150 (MQ=255)
                                                                                                                                     ttacccgcgaaAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCgg             <  2:1539547‑M2/139‑1 (MQ=255)
                                                                                                                                       acccgcgaaAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCg           >  1:5321557‑M2/10‑150 (MQ=255)
                                                                                                                                        cccgcgaaAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGa          <  1:4228955‑M2/142‑1 (MQ=255)
                                                                                                                                         ccgcgaaAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGat         >  1:3238206‑M2/8‑150 (MQ=255)
                                                                                                                                          cgcgaaAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGatc        <  2:4728372‑M2/144‑1 (MQ=255)
                                                                                                                                          cgcgaaAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGatc        >  2:3672948‑M2/7‑150 (MQ=255)
                                                                                                                                           gcgaaAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGatca       >  2:5606887‑M2/6‑150 (MQ=255)
                                                                                                                                             gaaGGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGatcatc     <  1:3915734‑M2/147‑1 (MQ=255)
                                                                                                                                                aGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTc  >  2:3631015/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
CTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGCGCTGTTCTTCGGCCTTACCGGCGTCTTCGGCAAACTGGCCACCAGCGCCAGCCCGGCGATCATCGTC  >  NC_002947/5130020‑5130312

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: