Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 603,661 T→C 9.7% intergenic (+57/‑35) ribAB‑I → / → ribH bifunctional 3,4‑dihydroxy‑2‑butanone 4‑phosphate synthase/GTP cyclohydrolase II/6,7‑dimethyl‑8‑ribityllumazine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947603,6610TC 9.7% 106.2 / 3.2 41intergenic (+57/‑35)ribAB‑I/ribHbifunctional 3,4‑dihydroxy‑2‑butanone 4‑phosphate synthase/GTP cyclohydrolase II/6,7‑dimethyl‑8‑ribityllumazine synthase
Reads supporting (aligned to +/- strand):  ref base T (12/25);  new base C (2/2);  total (14/27)
Fisher's exact test for biased strand distribution p-value = 5.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

TGCGCAAGATGCGCCTGATGAGTTCGCCGATGAA‑GTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTCGTTGTCGAAAGCCTGGTAAGCGGTGCCG  >  NC_002947/603513‑603804
                                                                                                                                                     |                                                                                                                                               
tgCGCAAGATGCGCCTGATGAGTTCGCCGATGAA‑GTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGt                                                                                                                                                 <  7:142692/149‑1 (MQ=255)
  cgcAAGATGCGCCTGATGAGTTCGCCGATGAAGGTTCAACGCGATATCCGGATTCGATCTGGAAGTTTTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCACTGCTCGAGTCCAGACCCCTGTCGAGGCCCCTTTTa                                                                                                                                                >  1:46766/1‑149 (MQ=255)
           cgcCTGATGAGTTCGCCGATGAA‑GTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCg                                                                                                                                      >  6:89570/1‑149 (MQ=255)
           cgcCTGATGAGTTCGCCGATGAA‑GTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCATACAAGCGGCCTc                                                                                                                                       <  4:35480/148‑1 (MQ=255)
           cgcCTGATGAGTTCGCCGATGAA‑GTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCATACAAGCGGCCTc                                                                                                                                       <  6:258772/148‑1 (MQ=255)
                    aGTTCGCCGATGAA‑GTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGaa                                                                                                                              <  2:46766/148‑1 (MQ=255)
                    aGTTCGCCGATGAA‑GTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCAAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGaa                                                                                                                              <  8:147015/148‑1 (MQ=255)
                         gCCGATGAA‑GTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGAt                                                                                                                           >  8:20027/1‑146 (MQ=255)
                         gCCGATGAA‑GTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATga                                                                                                                         <  1:144731/148‑1 (MQ=255)
                         gCCGATGAA‑GTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATga                                                                                                                         <  2:286593/148‑1 (MQ=255)
                                 a‑GTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCgg                                                                                                                <  1:291824/149‑1 (MQ=255)
                                 a‑GTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCgg                                                                                                                <  7:189928/149‑1 (MQ=255)
                                      cAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAAt                                                                                                             <  3:307658/148‑1 (MQ=255)
                                       aaCGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATg                                                                                                            <  6:379752/148‑1 (MQ=255)
                                                tCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATgag                                                                                                                        <  1:283284/127‑1 (MQ=255)
                                                tCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATgag                                                                                                                        >  2:283284/1‑127 (MQ=255)
                                                tCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGa                                                                                                  <  4:332951/149‑1 (MQ=255)
                                                   ggATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGAcc                                                                                                >  8:364975/1‑148 (MQ=255)
                                                     aTTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATc                                                                                             <  6:273442/149‑1 (MQ=255)
                                                         gATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGGCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGaa                                                                                          <  2:243246/148‑1 (MQ=255)
                                                              ggAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTAcc                                                                                    >  3:80439/1‑149 (MQ=255)
                                                               gAAGTTGTAGAATACGTGCCCTCCTAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGTGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTAcc                                                                                    <  3:6246/148‑1 (MQ=255)
                                                               gAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAg                                                                                         <  2:309076/143‑1 (MQ=255)
                                                                        gAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGCGTCCAGCCCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCGCTTGAAGATGAGAATATAGGAATGGCACTGCAGCCCCTCGAACGTACCGTCCTCcccc                                                                          >  5:306044/1‑149 (MQ=255)
                                                                           tACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGCTGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCa                                                                       <  6:285123/149‑1 (MQ=255)
                                                                                  ccTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATTCGGCCACGCTCATGAAGATTAGAATATAGGAATGACCCTGCACACCATAGAAGGTACCATCATCGCCCCCAAAGGTCg                                                                >  3:330787/1‑149 (MQ=255)
                                                                                            aCTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTtg                                                      <  7:20027/149‑1 (MQ=255)
                                                                                              ttGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCGTGTTTGCGGGCTGGCTATTGAAGATGAGAATATCGGAATGAACCGGAAGACCATCGAAGGTTCCATAATCGCCCCCAAAGGTCGCGATGCTTtggt                                                    >  2:57958/1‑149 (MQ=255)
                                                                                                             cGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCaa                                     <  4:80439/149‑1 (MQ=255)
                                                                                                                         gAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTCGtt                           <  7:364975/147‑1 (MQ=255)
                                                                                                                                   cccTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCtt                                                        <  5:266814/108‑1 (MQ=255)
                                                                                                                                   cccTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTCGTTGTCGAAAGCCTg               <  5:5400/149‑1 (MQ=255)
                                                                                                                                   cccTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATAGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTCGTTGTCGAAAGCCTg               >  3:108517/1‑149 (MQ=255)
                                                                                                                                    ccTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTCGTTGTCGAAAGCCTgg              <  8:253172/149‑1 (MQ=255)
                                                                                                                                      gagcGAGGCCGCATACAAGCGGCCTCGCTCTTGAAGATGAGAAAATCGGAATGACCCTGAAGACCATCGAAGGTACATTCATCGCACCCAAAGGTCGCTATGCTTTGGTCGTTCGCCGCTTCAACAGCTTCGTTGTCGAAAGCCAGGt             >  2:283505/4‑148 (MQ=255)
                                                                                                                                        tCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTCGTTGTCGAAAGCCTGGTAAg          <  5:89570/149‑1 (MQ=255)
                                                                                                                                        gcGAGGCCGCATACAAGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGACGCTTCAACAGCTACGTTGTCGAAAGCCTGGTAAg          >  7:215672/2‑149 (MQ=255)
                                                                                                                                          gAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCACTTCAACAGCTTCGTTGTCGAAAGCCTGGTAAGCg        <  1:57958/149‑1 (MQ=255)
                                                                                                                                               cGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTCGTTGTCGAAAGCCTGGTAAGCGGTGcc   >  5:240200/1‑149 (MQ=255)
                                                                                                                                               cGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTCGTTGTCGAAAGCCTGGTAAGCGGTGcc   <  4:108517/149‑1 (MQ=255)
                                                                                                                                                gCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTCGTTGTCGAAAGCCTGGTAAGCGGTGCCg  >  4:239143/1‑149 (MQ=255)
                                                                                                                                                  ttGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTc                                                                 <  6:272921/84‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
TGCGCAAGATGCGCCTGATGAGTTCGCCGATGAA‑GTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTCGTTGTCGAAAGCCTGGTAAGCGGTGCCG  >  NC_002947/603513‑603804

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑