Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,882,009 T→A 6.5% intergenic (+26/+79) 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,0090TA 6.5% 686.2 / 3.0 237intergenic (+26/+79)rpiA/PP_5151ribose‑5‑phosphate isomerase A/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (78/135);  new base A (18/3);  total (98/139)
Fisher's exact test for biased strand distribution p-value = 2.23e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 240 total aligned reads displayed.

AGCTGGAAGCGCAGATCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTCCCCGCTTCAT  >  NC_002947/5881862‑5882158
                                                                                                                                                   |                                                                                                                                                     
tgcaggcgacggtgccgcaggcggaggggattgcgatggatgaagcggggaccatttacctggtcagtgagccgaacctgttctatgtgttcaagcagccgactgactgatcgcctgtcagggcctgTTCGCGGGTAAACCCGCTCCTAc                                                                                                                                                     >  2:2448333‑M1/128‑150 (MQ=255)
aGCTAGAAGCGCAGATCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTAC                                                                                                                                                     <  1:4534251‑M1/150‑1 (MQ=255)
aGCTAGAAGCGCAGATCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTAC                                                                                                                                                     <  1:4534657‑M1/150‑1 (MQ=255)
  cTGGAAGCGCAGATCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACaa                                                                                                                                                   <  2:3455980/150‑1 (MQ=255)
  cTGGAAGCGCAGATCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACaa                                                                                                                                                   <  2:4678139/150‑1 (MQ=255)
  cTGGAAGCGCAGATCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACaa                                                                                                                                                   >  2:2909765/1‑150 (MQ=255)
  cTGGAAGCGCAGATCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACaa                                                                                                                                                   <  1:5452401/150‑1 (MQ=255)
  cTGGAAGCGCAGATCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCGCCTACaa                                                                                                                                                   >  1:3635559/1‑150 (MQ=255)
   tGGAAGCGCAGATCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAg                                                                                                                                                  >  2:2736779/1‑150 (MQ=255)
              aTCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCg                                                                                                                                       <  1:3749645/150‑1 (MQ=255)
                     ccATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGa                                                                                                                                <  2:1377348/150‑1 (MQ=255)
                            cgGTGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCCGAAGGCCGAGTACAACCCTTCGCGGGTCAACCCGCTCCTACACGGTACCCCACCGGCCTGGAACCGTgg                                                                                                                         <  2:358410/150‑1 (MQ=255)
                            cgGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTgg                                                                                                                         <  2:3938860/150‑1 (MQ=255)
                             ggCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGc                                                                                                                        <  1:1071755/150‑1 (MQ=255)
                                  ggtgACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAATGTACCCCACCGGCCTGGAACCGTGGCGATcc                                                                                                                   <  1:1128062/150‑1 (MQ=255)
                                  ggtgACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATcc                                                                                                                   <  1:4893243/150‑1 (MQ=255)
                                  ggtgACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATcc                                                                                                                   <  1:4367705/150‑1 (MQ=255)
                                   gtgACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATccc                                                                                                                  >  2:3642690/1‑150 (MQ=255)
                                       ccAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgta                                                                                                              >  2:5455797‑M2/1‑146 (MQ=255)
                                       atgaagcggggaccatttacctggtcagtgagccgaacctgttctatgtgttcaagcagccgactgactgatcgcctgtcagggcctgTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTa                                                                                                              >  2:1618727‑M1/89‑150 (MQ=255)
                                        tgaagcggggaccatttacctggtcagtgagccgaacctgttctatgtgttcaagcagccgactgactgatcgcctgtcagggcctgTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAg                                                                                                             >  1:3671029‑M1/88‑150 (MQ=255)
                                        cAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTATAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtag                                                                                                             >  2:4849712‑M2/1‑145 (MQ=255)
                                        cAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtag                                                                                                             >  2:4849713‑M2/1‑145 (MQ=255)
                                        cAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGGAAAACCGCTCCAACAACGGACCCCACCCGCCAGGAACCGTGGCGATACCagtag                                                                                                             <  1:3987654‑M2/150‑6 (MQ=255)
                                          aCGGCCTGTTCGCGGCGCGGCCGGCGGAT‑‑‑CTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGGAAAACCCCGCCAACAACGGACCCCACCCGCCAGGGACCCGGGCGATACCGGTAGGAGCg                                                                                                        >  1:2116269/1‑150 (MQ=255)
                                            ggCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAgggatcgccacggttccaggccggtggggtaccttgtaggagc                                                                                                         <  2:3032536‑M1/150‑44 (MQ=255)
                                            ggCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagc                                                                                                         >  2:5135790‑M2/1‑141 (MQ=255)
                                            ggCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagc                                                                                                         >  2:5136045‑M2/1‑141 (MQ=255)
                                            gcggggaccatttacctggtcagtgagccgaacctgttctatgtgttcaagcagccgactgactgatcgcctgtcagggcctgTTCGCGGGTAAACCCGCTCCTACAAGGGACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGc                                                                                                         >  1:77620‑M1/84‑150 (MQ=255)
                                              ccTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAgggatcgccacggttccaggccggtggggtaccttgtaggagcgg                                                                                                       <  2:3288573‑M1/150‑46 (MQ=255)
                                              ccTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgg                                                                                                       >  2:1977193‑M2/1‑139 (MQ=255)
                                               cTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCagtaggagcggg                                                                                                      <  2:2792571‑M2/150‑13 (MQ=255)
                                                    cgcggcgcggCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgggtttac                                                                                                 >  2:395678‑M2/1‑133 (MQ=255)
                                                     gcggcgcggcCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTGAAACCGCGCCAACAAAGGAACCCCCCCGCCAGGGACCCAGGCCATACCGGTAGGGGCGGGGTTagc                                                                                                >  1:3326662/1‑148 (MQ=255)
                                                      cggcgcggcCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAACACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGCACCGTGGCGATCCCtgtaggagcgggtttaccc                                                                                               >  2:2376788‑M2/1‑131 (MQ=255)
                                                        gcgcggcCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgggtttacccgc                                                                                             >  2:248420‑M2/1‑129 (MQ=255)
                                                        gcgcggcCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgggtttacccgc                                                                                             >  1:5113433‑M2/1‑129 (MQ=255)
                                                           cggccggcGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgggtttacccgcgaa                                                                                          >  2:672567‑M2/1‑126 (MQ=255)
                                                                   ggATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAAACCGCGCCTACAAGGTAACCCACCGGCCTGGAACCGTGGCGATACCGGTAGGAGCGGGTTTACCCGCGAACAGGCCCt                                                                                  >  2:2749963/1‑150 (MQ=255)
                                                                      tCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTTAAACCCCAACTACAAGGGACCCCACCGGCCGGGAACCCTGGCGCGCCCAGTAGGAGCCGGTTTTCCCCCGGCCAGGcccggcc                                                                               >  2:3288773/1‑146 (MQ=255)
                                                                         gctgctgctgGGCACCTCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAgg                                                                            <  2:5109365/150‑1 (MQ=255)
                                                                         gctgctgcggGGCACCGCTGAAGGCGTTAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAgg                                                                            <  2:2341858/150‑1 (MQ=255)
                                                                           tgctgctgGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTCCAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCg                                                                          <  2:1798938/150‑1 (MQ=255)
                                                                            gctgctgGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAAACCGCGCCAAAAAGGTAACCCACCCGCCGGGCACCCTGGCGATACCGGTAGGAGCGGGTTTACCCCCCAACAGGCCCGGACACGCGa                                                                         >  1:1796476/1‑150 (MQ=255)
                                                                               gctgGGCACCGCTGAAGGCGTGAAGAGCCGGAAGGCCGAGTAAAACCCTGCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGAtca                                                                      <  2:4026390/150‑1 (MQ=255)
                                                                               gcggggCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGGGGGTTTACCCGCGAACAGGCCCTGACAGGCGAtca                                                                      <  1:3181502/147‑1 (MQ=255)
                                                                                   ggCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGAtcagtca                                                                  <  2:147190/150‑1 (MQ=255)
                                                                                     cACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAAACCGCTCCAACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATcagtcagt                                                                >  2:4605783/1‑150 (MQ=255)
                                                                                      aCCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCAACAAGGTACCCCACCCGCCTGGAACCCTGGCGATCCCAGTAGGAGCGGGTTTACCCGCGAACAGGCCCAGACAGGCGATCagtcagtc                                                               >  2:4156303/1‑150 (MQ=39)
                                                                                        cGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCgg                                                             <  2:3692063/150‑1 (MQ=255)
                                                                                           tGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCCGCCTGGAACCGTGGCGATCCCGGTAGGGGCGGGTTTACCCGCGAACAGGCCCGGACAGGCGCTCAGTCAGTCGgcgg                                                          >  2:5309196/1‑148 (MQ=255)
                                                                                            gAAGGCGTGAAGAGCCTGAAGGCCGATTAAACCCCTTCGCGGTTAAACCCGCTCCCACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGGGGGTTAACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGgctgc                                                         <  1:5334590/150‑1 (MQ=255)
                                                                                            gAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGgctgc                                                         <  1:4133973/150‑1 (MQ=255)
                                                                                             acggCGTGAAGAGCCTGACGGCCGAGTAAAACCCTCCGGGGGTAAACCCGCTCCTACACGGTACCCCACCGGCCTGGACCCGGGGCGATCCCTGTAGGAGGGGTTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGgctgct                                                        <  1:194334/148‑1 (MQ=255)
                                                                                             aaGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTCCACGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGgctgct                                                        <  1:5262822/150‑1 (MQ=255)
                                                                                              aGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACACGGTACCCCACCGGCCTGGAACCGTGGCGAGCCCTGTCGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCtt                                                       <  1:1190245/150‑1 (MQ=255)
                                                                                                gCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTTGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGa                                                     >  1:3303422/1‑150 (MQ=255)
                                                                                                gCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGa                                                     <  1:4579282/150‑1 (MQ=255)
                                                                                                 cGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGaa                                                    <  2:4710091/150‑1 (MQ=255)
                                                                                                   tGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAaca                                                  <  2:5266212/150‑1 (MQ=255)
                                                                                                    gAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAacac                                                 <  2:5420483/150‑1 (MQ=255)
                                                                                                     aaGTGCCTGAAGGCCGAGTAAAACCCTTCGGGGTTAAACCCGCTCCTACCCGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAAcaca                                                <  1:276680/150‑1 (MQ=255)
                                                                                                     aaGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAAcaca                                                <  1:2321647/150‑1 (MQ=255)
                                                                                                     aaGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCCCTCCAACAAGGTACCCCACCGGCCGGGCACCGTGGCGGAGCCCGGAGGTGCGGGGTTAACCCCCAACAGGCCCAGACAGGCCATCACTCAGGAGGCGGCGTGAAcaca                                                >  2:638160/1‑150 (MQ=255)
                                                                                                     aaGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAAACCGCTCCTAAAACGTACCCCACCGGCCAGGAACCGTGGCGATCCCGGTAGGAGCGGGTTTAACCGCCAACAGGCCCAGAAAGGCGATCAGTCAGTAGGCTGCTTGAAcaca                                                >  2:572597/1‑150 (MQ=255)
                                                                                                       gagCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATa                                              <  2:2423593/150‑1 (MQ=255)
                                                                                                            tGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGTTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACa                                         <  2:1206566/150‑1 (MQ=255)
                                                                                                             atcgcctgtcagggcctgTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTg                                                                                                                <  1:1076829‑M1/60‑1 (MQ=255)
                                                                                                                 gCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGTTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTc                                    <  1:3571515/150‑1 (MQ=255)
                                                                                                                 gCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTc                                    <  2:4301749/150‑1 (MQ=255)
                                                                                                                   cGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCgg                                  <  2:4369968/150‑1 (MQ=255)
                                                                                                                    gtttAAAACCGTTCGCGGGTAAACCCGCTCCTAAACGGTACCCCACCGGCCTGGAACCTTGGCGTTCCCTGTAGTAGGGTTTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGc                                 <  1:287450/147‑1 (MQ=255)
                                                                                                                    gattAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGc                                 <  1:3619647/147‑1 (MQ=255)
                                                                                                                    gAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGTAACCGTGGCGATCCCTGTTGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGc                                 <  1:4348527/150‑1 (MQ=255)
                                                                                                                         aaaCCCTTCGCGGGTAAACCCGCTCATACAAGGTACCCCACCGGCCTGGAACGGTGGCGATCCCTGTCGGAGGGGTTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACt                            <  1:1554698/150‑1 (MQ=255)
                                                                                                                             ccTTCGCGGGTAAACCCGCGCCAACAAGGTACCCCACCGGCCAGGAAACGTGGCGATCCCAGTAGGAGCGGGTTTAACCGCGAACAGGCCCAGACAGGCGATCAGTCACTAGGCTGCTTGAAAACAAATAACAGGTTAGGCTCACAGAcc                        >  2:1065152/1‑150 (MQ=255)
                                                                                                                              cTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCa                       <  1:3712669/150‑1 (MQ=255)
                                                                                                                                TCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCGGTAGGAGCGGGGTTGCCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCGTGTACACCTCGTACAGGTTCGGCTCACTGACCAgg                     <  2:3324229‑M1/150‑1 (MQ=255)
                                                                                                                                     gtTAAACCCGCTCCTACATGGTACCCCCCCGGCCGGGACCGGTGGCGATCCGTGTTGGTGGGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAAt                >  1:4696267‑M1/3‑150 (MQ=255)
                                                                                                                                     ggTAACCCCGCTCCTCCTgggatcgccacggtcccaggcgggtggggtcccttgtaggcgcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgc                <  2:3051749‑M1/150‑133 (MQ=255)
                                                                                                                                     ggTAACCCCGCTCCTACCgggatcgcccgggttcccggcgggtggggtaccttgtaggtgcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgc                <  2:2123646‑M1/150‑133 (MQ=255)
                                                                                                                                     GGTAAACCCGCTCATACAATTTACCCCACCGGCCTGGAACCGTGGATATCCGTGTTGTCGCGGTTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCTGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTTACCAGGTAAAt                >  1:2192599‑M1/1‑150 (MQ=255)
                                                                                                                                          ACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTcc           <  1:863954‑M1/150‑1 (MQ=255)
                                                                                                                                           cccGCTCCTACAgggatcgccacggttccaggccggtggggtaccttgtaggagcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcgaacag          >  2:4507424‑M1/1‑12 (MQ=255)
                                                                                                                                           cccGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcgaacag          <  1:4347241‑M2/150‑105 (MQ=255)
                                                                                                                                           CCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTccc          <  1:1410358‑M1/150‑1 (MQ=255)
                                                                                                                                            ccGCTCCTACAgggatcgccacggttccaggccggtggggtaccttgtaggagcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcgaacagg         >  2:4688395‑M1/1‑11 (MQ=255)
                                                                                                                                            ccGCTCCTACAgggatcgccacggttccaggccggtggggtaccttgtaggagcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcgaacagg         >  1:5170886‑M1/1‑11 (MQ=255)
                                                                                                                                            CCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTGGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCagtcagtc                                                               <  1:891537‑M1/96‑1 (MQ=255)
                                                                                                                                            CCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAAcaca                                                <  1:4291141‑M1/111‑1 (MQ=255)
                                                                                                                                              GCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTCCCCGc       <  1:3969351‑M1/150‑1 (MQ=255)
                                                                                                                                                tcctccCgggatcgccacggttcccggccggtggggtaccttgtaggagcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcgaacaggccgt     <  2:392774‑M1/150‑144 (MQ=255)
                                                                                                                                                  cTACCgggtgcgcccgggttcccggcgggtggggtcccttgtaggcgggggttaacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcgaacaggccgttg   <  1:2103773‑M1/150‑146 (MQ=255)
                                                                                                                                                  cTACAgggatcgccacggttccaggccggtggggtaccttgtaggagcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcgaacaggccgttg   >  1:2429246‑M1/1‑5 (MQ=255)
                                                                                                                                                  cTACAgggatcgccacggttccaggccggtggggtacctggtaggagcgggtttacccgcgacgggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcgaacaggccgttg   >  2:4849428‑M1/1‑5 (MQ=255)
                                                                                                                                                  cTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCtgtaggagcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcgaacaggccgttg   <  1:2782365‑M2/150‑112 (MQ=255)
                                                                                                                                                  CTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTCCCCGCTTCa   <  1:1285369‑M1/150‑1 (MQ=255)
                                                                                                                                                   tACAgggatcgccacggttccaggccggtggggtaccttgtaggagcgggtttacccgcgaagggttttactcggccttcaggctcttcacgccttcagcggtgcccagcagcagcagatccgccggccgcgccgcgaacaggccgttgg  >  1:3559967‑M1/1‑4 (MQ=255)
                                                                                                                                                   TACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTCCCCGCTTCAt  <  1:1739104‑M1/150‑1 (MQ=255)
                                                                                                                                                   TACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTCCCCGCTTCAt  <  2:1097667‑M1/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AGCTGGAAGCGCAGATCAACGCCATCGTCGGCGTGGTGACCAACGGCCTGTTCGCGGCGCGGCCGGCGGATCTGCTGCTGCTGGGCACCGCTGAAGGCGTGAAGAGCCTGAAGGCCGAGTAAAACCCTTCGCGGGTAAACCCGCTCCTACAAGGTACCCCACCGGCCTGGAACCGTGGCGATCCCTGTAGGAGCGGGTTTACCCGCGAACAGGCCCTGACAGGCGATCAGTCAGTCGGCTGCTTGAACACATAGAACAGGTTCGGCTCACTGACCAGGTAAATGGTCCCCGCTTCAT  >  NC_002947/5881862‑5882158

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: