Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,882,025 G→C 8.3% intergenic (+42/+63) rpiA → / ← PP_5151 ribose‑5‑phosphate isomerase A/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,882,0250GC 8.3% 622.8 / 7.3 200intergenic (+42/+63)rpiA/PP_5151ribose‑5‑phosphate isomerase A/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (67/113);  new base C (15/5);  total (82/118)
Fisher's exact test for biased strand distribution p-value = 1.53e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 203 total aligned reads displayed.

ATCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTCCCCGCTTCATCCATCGCAATCCCCTC  >  NC_002947/5881876‑5882174
                                                                                                                                                      |                                                                                                                                                     
aTCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCg                                                                                                                                                       <  1:3749645/150‑1 (MQ=255)
  cAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGc                                                                                                                                                     <  2:145637/150‑1 (MQ=255)
    aCGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCt                                                                                                                                                   >  1:4802387/1‑150 (MQ=255)
      gCCATCGTCGGCGTGGTGACCAACGGCCTGGTCGCGGCGCGGCCGGCGGATCTGCTGCTGCGGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAATCCCTTGGCGGGTAAACCCGCGCCGACAAGCTACCCCACCGGCCggg                                                                                                                                                 >  1:2750738/1‑150 (MQ=255)
            gtcgGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCTTCGCGGGGTAAACCCCCACGACAAACGTACCCCCCCGCCCGGGACCCgg                                                                                                                                           >  1:4474595/1‑149 (MQ=255)
              cgGTGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCCGAAGGCCGAGTACAACCCTTCGCGGGTCAACCCGCTCCTACACGGTACCCCACCGGCCTGGAACCGTgg                                                                                                                                         <  2:358410/150‑1 (MQ=255)
              cgGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTgg                                                                                                                                         <  2:2738779/150‑1 (MQ=255)
               ggCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGc                                                                                                                                        <  1:1071755/150‑1 (MQ=255)
                    ggtgACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAATGTACCCCACCGGCCTGGAACCGTGGCGATcc                                                                                                                                   <  1:1128062/150‑1 (MQ=255)
                    ggtgACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATcc                                                                                                                                   <  1:4367705/150‑1 (MQ=255)
                    ggtgACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATcc                                                                                                                                   <  2:1114943/150‑1 (MQ=255)
                     gtgACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATccc                                                                                                                                  >  2:3642327/1‑150 (MQ=255)
                       ggatgaagcggggaccatttacctggtcagtgagccgaacctgttctatgtgttcaagcagccgactgactgatcgcctgtcagggcctgTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTg                                                                                                                                >  2:1360052‑M1/91‑150 (MQ=255)
                         atgaagcggggaccatttacctggtcagtgagccgaacctgttctatgtgttcaagcagccgactgactgatcgcctgtcagggcctgTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTa                                                                                                                              >  1:3928160‑M1/89‑150 (MQ=255)
                         atgaagcggggaccatttacctggtcagtgagccgaacctgttctatgtgttcaagcagccgactgactgatcgcctgtcagggcctgTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTa                                                                                                                              >  2:1618727‑M1/89‑150 (MQ=255)
                          tgaagcggggaccatttacctggtcagtgagccgaacctgttctatgtgttcaagcagccgactgactgatcgcctgtcagggcctgTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAg                                                                                                                             >  1:3670779‑M1/88‑150 (MQ=255)
                          cAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTATAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtag                                                                                                                             >  2:4849712‑M2/1‑145 (MQ=255)
                          cAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGGTAAACCCCACCAAAAAAGTAACCCAACCCCCGGGAACCGGGGCGGTACCAGTAg                                                                                                                             >  1:4290912/1‑150 (MQ=255)
                          cAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGGAAAACCGCTCCAACAACGGACCCCACCCGCCAGGAACCGTGGCGATACCagtag                                                                                                                             <  1:3987654‑M2/150‑6 (MQ=255)
                            aCGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCGGGAACCGTGGCGATCCCgggagga                                                                                                                           <  2:392774‑M2/150‑8 (MQ=255)
                            aCGGCCTGTTCGCGGCGCGGCCGGCGGAT‑‑‑CTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGGAAAACCCCGCCAACAACGGACCCCACCCGCCAGGGACCCGGGCGATACCGGTAGGAGCg                                                                                                                        >  1:2116269/1‑150 (MQ=255)
                              ggCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagc                                                                                                                         >  2:5136045‑M2/1‑141 (MQ=255)
                                ccTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgg                                                                                                                       >  2:4688395‑M2/1‑139 (MQ=255)
                                ccTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgg                                                                                                                       >  1:1354035‑M2/1‑139 (MQ=255)
                                 cTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCagtaggagcggg                                                                                                                      <  2:2792571‑M2/150‑13 (MQ=255)
                                      cgcggcgcggCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgggtttac                                                                                                                 >  2:395678‑M2/1‑133 (MQ=255)
                                       gcggcgcggcCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCGCCTACAAGGGACCCCACCCGCCTGGGACCGTGGCGATCCCaggaggagcggggttacc                                                                                                                <  2:3051749‑M2/150‑19 (MQ=255)
                                       gcggcgcggcCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAAACCCCGCCAACAAAGTAACCCACCCGCCGGGGACCCTGGCGAAACCcgtaggagcggggttaac                                                                                                                <  1:3326161‑M2/150‑19 (MQ=255)
                                        cggcgcggcCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAACACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGCACCGTGGCGATCCCtgtaggagcgggtttaccc                                                                                                               >  2:2376788‑M2/1‑131 (MQ=255)
                                          gcgcggcCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgggtttacccgc                                                                                                             >  2:248420‑M2/1‑129 (MQ=255)
                                             cggccggcGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgggtttacccgcgaa                                                                                                          >  2:672567‑M2/1‑126 (MQ=255)
                                              ggccggcGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAAc                                                                                                         <  1:3650601/150‑1 (MQ=255)
                                               gccggcGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACa                                                                                                        <  1:4613662/150‑1 (MQ=255)
                                                     gtaTCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACACGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCt                                                                                                  <  2:4740815/148‑1 (MQ=255)
                                                      gATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTTAACCCGCTCCTACAAGGGACCCCACCCGCCGGGAACCGTGGCGATCCCGGTAGGGGCGGGGTTTCCCGCGAAAAGGcccag                                                                                                 >  2:1554698/1‑148 (MQ=255)
                                                        tCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTAAAAGGTACCCCACCCGCCAGGAACCCTGGCGATCCCGGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGAc                                                                                               >  2:5316783/1‑150 (MQ=255)
                                                        tCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAAACCCCTCCGACCAGGTACCCCACCCGCCTGGAAACGTGGCGATCCCGGTTGGAGCGGGGTTTCCCCCGAACAGGCCCGGAc                                                                                               >  2:529062/1‑150 (MQ=255)
                                                         ctgctgctgctgGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCATGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACa                                                                                              <  2:1035063/150‑1 (MQ=255)
                                                           gctgctgctgGGCACCTCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAgg                                                                                            <  2:5109365/150‑1 (MQ=255)
                                                            ctgctgctgGGCACCGGTGATGGCGTGAAGAGCCTGATGGCGGTTTAAAACCCTTGGGGGTTAAACCCGCTCCTACCCGGTACCCCACGGGCCTGGAACCGTGGCGATCCCTTTTGGAGGGGTTTTACCCGCGAACAGGCCCTGACAGGc                                                                                           <  1:781433/150‑1 (MQ=255)
                                                             tgctgctgGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTTAAACCCCACCAACAAGGTTCCCCACCCGCCTGGAACCCTGGCGATACCGGTAGGAGCGGGTTTACCCCCGAACAGGCCCAGACAGGCg                                                                                          >  1:1756252/1‑150 (MQ=255)
                                                              gctgctgGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAAACCGCGCCAAAAAGGTAACCCACCCGCCGGGCACCCTGGCGATACCGGTAGGAGCGGGTTTACCCCCCAACAGGCCCGGACACGCGa                                                                                         >  1:1796476/1‑150 (MQ=255)
                                                               ctgctgGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGATTAAAACCCTTCGCGGGTAAACCCGCTCCTACACGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGAt                                                                                        <  1:1275188/150‑1 (MQ=255)
                                                                 ggtgGGCCCCGCTGATGGGGTGAAGCGCCTGAAGGCCGAGTAAAACCCTGCGCGGGTAAACCCGCTCCTACAAGGTCCCCCCCCGGCCTGGAACGGTGGCGATCCGTGTTGGTGGGGTTTTACCCGCGAACAGGCCCTGACAGGCGAtca                                                                                      <  1:3294607/148‑1 (MQ=255)
                                                                 gctgGGCCCCGCTTACGGGTTGAAGAGCCTGACGGCCGTGTAAAACCCTTCGGGGGTAAACCCGCTCCTCCACGGTACCCCCCCGGCCTGGAACCGTGGCGATCCCTGTAGGTGGGGGTTTACCCGCGAACAGGCCCTGACAGGCGAtca                                                                                      <  1:1984464/150‑1 (MQ=255)
                                                                 gctgGGCACCGCTGAAGGCGTGAAGAGCCGGAAGGCCGAGTAAAACCCTGCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGAtca                                                                                      <  2:4026390/150‑1 (MQ=255)
                                                                   tgGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACACGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGAtcagt                                                                                    <  2:5240495/150‑1 (MQ=255)
                                                                       cACCGCTGACGGCGTGAAGCGCCTGACGGCCGATTAAACCCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGAGCCCTGTAGGCGCGGTTTTACCCGCGAACAGGCCCTGACAGGCGATcagtcagt                                                                                <  1:1160809/150‑1 (MQ=255)
                                                                        aCCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTTAAACCGCGCCAACAAGGTACCCCACCCGCCAGGGACCGTGGCGATCCCTGTAGGCGCGGGTTTTCCCGCGGACAGGCCCAGACAGGCGATCcgtcagtc                                                                               >  2:2824449/1‑150 (MQ=255)
                                                                           gCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGgc                                                                            <  2:3681546/150‑1 (MQ=255)
                                                                              gagggCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTCCGCGGGTAAACCCGCTCCACCACGGTACCCCCCCGGCCTGGAACCGGGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGgctgc                                                                         <  1:3335614/147‑1 (MQ=255)
                                                                              gAAGGCGTGAAGAGCCTGAAGGCCGATTAAACCCCTTCGCGGTTAAACCCGCTCCCACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGGGGGTTAACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGgctgc                                                                         <  1:5334590/150‑1 (MQ=255)
                                                                               aaGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTCCACGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGgctgct                                                                        <  1:5262822/150‑1 (MQ=255)
                                                                                  gCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGa                                                                     <  1:4579282/150‑1 (MQ=255)
                                                                                   cGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGaa                                                                    <  2:4591781/150‑1 (MQ=255)
                                                                                    gTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCAGGAAACCTGGCGATCCCAGTAGGAGCGGGTTTACCCGCGAACACGCCCTGACAGGCGAACAGTCACTAGGCTGCTTGAac                                                                   >  2:4622068/1‑150 (MQ=255)
                                                                                      gAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAacac                                                                 <  2:5420483/150‑1 (MQ=255)
                                                                                       aaGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTTAACCCGCTCCAAAAATGTTACCCACCGGGCTGGTAACGTGGCGGTCCCGGTAGGGTCGGGTTTTCCCGCGGAAAGGCCCTGGGGGGCGATGAGTCAGTGGGCGGCGTGTAgaca                                                                >  2:637746/1‑150 (MQ=255)
                                                                                       aaGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAAcaca                                                                <  1:2321478/150‑1 (MQ=255)
                                                                                         gagCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATa                                                              <  2:2423593/150‑1 (MQ=255)
                                                                                            ccTGAAGGCCAAGTAAAACCCTTCGCGGGTAAACCCCCGCCAACAAGGTACCCCACCCGCCAGGAAACCGGGCGCTCCCGGTAGGAGCGGGTTTACCCCCGAACAGGCCCAGACAGGCGATCAGTAAGTCGGCTGCTTTAACACAAAGaa                                                           >  1:5147047/1‑150 (MQ=255)
                                                                                              tGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACa                                                         <  2:1206343/150‑1 (MQ=255)
                                                                                               gAAGGCCGAGTAAAACCCTTCGCGGGTTAAACCCCTCCTACAAGGTCCCCCACCGGCCGGGAACCCTGGCGATCCCGGTAGGAGCGGGTTTACCCGCGAACACGCCCTGACAGGCGATCAGTCAGTCGGCGGCGTGTACACATAGAAccc                                                        >  2:3424618/1‑148 (MQ=255)
                                                                                               gAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACACGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAg                                                        <  2:3916575/150‑1 (MQ=255)
                                                                                               atcgcctgtcagggcctgTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTg                                                                                                                                <  1:1076829‑M1/60‑1 (MQ=255)
                                                                                               atcgcctgtcagggcctgTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTg                                                                                                                                >  2:1076829‑M1/19‑78 (MQ=255)
                                                                                                   gCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGTTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTc                                                    <  1:3571515/150‑1 (MQ=255)
                                                                                                   gCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTc                                                    <  2:1964617/150‑1 (MQ=255)
                                                                                                     cGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCgg                                                  <  2:4369968/150‑1 (MQ=255)
                                                                                                     cGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCgg                                                  <  2:3293811/150‑1 (MQ=255)
                                                                                                      gtttAAAACCGTTCGCGGGTAAACCCGCTCCTAAACGGTACCCCACCGGCCTGGAACCTTGGCGTTCCCTGTAGTAGGGTTTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGc                                                 <  1:287450/147‑1 (MQ=255)
                                                                                                      gattAAAACCCTTCGGGGTTAAACCCGCTCCTACACGTTACCCCACCGGCCTGGAACCGTGGCTATCCCTGTTGGTGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGc                                                 <  1:1588243/147‑1 (MQ=255)
                                                                                                      gattAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGc                                                 <  1:3619647/147‑1 (MQ=255)
                                                                                                      gAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGTAACCGTGGCGATCCCTGTTGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGc                                                 <  1:4348527/150‑1 (MQ=255)
                                                                                                      gAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGATcacg                                                 <  2:5043908/150‑4 (MQ=255)
                                                                                                             aCCCTTCGGGGGTAAACCCGCTCCTACATGGTACCCCACCGGCCTGGAACCGTGGCTATCCCTGTAGGTGGGGTTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGa                                          <  1:2364604/150‑1 (MQ=255)
                                                                                                               ccTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCAGACAGGCGATCAGTCAGTCGGCTGCTTGAACACCTAGAACAGGTTCGGCTCACTGAcc                                        >  2:5622031/1‑150 (MQ=255)
                                                                                                                cTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCa                                       <  1:3712669/150‑1 (MQ=255)
                                                                                                                       GGGTAACCCCGCTCCTACACGGTACCCCACCGGCCTGTAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTaaa                                 >  1:3000300‑M1/1‑150 (MQ=255)
                                                                                                                        GGTAAACCCGCTCATACAATTTACCCCACCGGCCTGGAACCGTGGATATCCGTGTTGTCGCGGTTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCTGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTTACCAGGTAAAt                                >  1:2192599‑M1/1‑150 (MQ=255)
                                                                                                                         GTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATg                               >  1:4986723‑M1/1‑150 (MQ=255)
                                                                                                                              cccGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcgaacag                          <  1:4347241‑M2/150‑105 (MQ=255)
                                                                                                                               ccGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcgaacagg                         <  2:3288573‑M2/150‑106 (MQ=255)
                                                                                                                               CCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTGGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCagtcagtc                                                                               <  1:891537‑M1/96‑1 (MQ=255)
                                                                                                                               CCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAAcaca                                                                >  2:4291141‑M1/1‑111 (MQ=255)
                                                                                                                               CCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCCGGTAAATGGTcccc                         <  2:4126594‑M1/150‑1 (MQ=255)
                                                                                                                                 GCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTCCCCGc                       >  1:3580476‑M1/1‑150 (MQ=255)
                                                                                                                                   tCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcgaacaggccgt                     <  1:5618457‑M2/150‑110 (MQ=255)
                                                                                                                                     cTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcgaacaggccgttg                   <  1:2782365‑M2/150‑112 (MQ=255)
                                                                                                                                         AAGGTCCCCCACCGGCCTGGAACCGTGGCTATCCCTGTAGGAGCGGTTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTCCCCGCTTCATCCa               >  1:2047070‑M1/1‑150 (MQ=255)
                                                                                                                                           ggTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcgaacaggccgttggtcacc             <  1:1527440‑M2/150‑118 (MQ=255)
                                                                                                                                            gTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGATTACCCGCGAACAGGCCCTGACAGGCGAACACTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTCCCCGCTTCATCCATCg            <  2:4173304/150‑1 (MQ=255)
                                                                                                                                             tACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTCCCCGCTTCATCCATCGc           <  2:2194442/150‑1 (MQ=255)
                                                                                                                                                cccACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTCCCCGCTTCATCCATCGCAAt        >  1:154481/1‑150 (MQ=255)
                                                                                                                                                cccACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTCCCCGCTTCATCCATCGCAAt        <  2:4467770/150‑1 (MQ=255)
                                                                                                                                                  cACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGAtcagtca                                                                                  >  1:4068937/1‑74 (MQ=255)
                                                                                                                                                  cACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTCCCCGCTTCATCCATCGCAATcc      >  1:4698389/1‑150 (MQ=255)
                                                                                                                                                  cACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTCCCCGCTTCATCCATCGCAATcc      >  2:3321366/1‑150 (MQ=255)
                                                                                                                                                     cGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTCCCCGCTTCATCCATCGCAATCCCCt   >  2:4366402/1‑150 (MQ=255)
                                                                                                                                                      ggCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTCCCCGCTTCATCCATCGCAATCCCCTc  <  2:2103773/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
ATCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTCCCCGCTTCATCCATCGCAATCCCCTC  >  NC_002947/5881876‑5882174

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