Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 603,658 T→A 9.8% intergenic (+54/‑38) 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,6580TA 9.8% 97.4 / 3.8 41intergenic (+54/‑38)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 (11/25);  new base A (2/2);  total (14/27)
Fisher's exact test for biased strand distribution p-value = 5.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.41e-01

GTGTGCGCAAGATGCGCCTGATGAGTTCGCCGATGA‑A‑GTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTCGTTGTCGAAAGCCTGGTAAGCGGTGCCG  >  NC_002947/603510‑603804
                                                                                                                                                      |                                                                                                                                                  
gtgtGCGCAAGATGCGCCTGATGAGTTCGCCGATGA‑A‑GTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCAAGTCCAGACCCCTGTCGAGGCCGCt                                                                                                                                                    >  7:147015/1‑149 (MQ=255)
   tgCGCAAGATGCGCCTGATGAGTTCGCCGATGA‑A‑GTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGt                                                                                                                                                 <  7:142692/149‑1 (MQ=255)
     cgcAAGATGCGCCTGATGAGTTCGCCGATGA‑AGGTTCAACGCGATATCCGGATTCGATCTGGAAGTTTTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCACTGCTCGAGTCCAGACCCCTGTCGAGGCCCCTTTTa                                                                                                                                                >  1:46766/1‑149 (MQ=255)
        aaGATGCGCCTGATGAGTTCGCCGATGATA‑GTTAAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCGTTGACGCCATATGCTCGAGTCCAGCCCCATGTCGAGGCCGCtt                                                                                                                                                   >  1:222969/1‑143 (MQ=255)
        aaGATGCGCCTGATGAGTTCGCCGATGA‑A‑GTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCGTTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCtt                                                                                                                                                   <  2:222969/142‑1 (MQ=255)
              cgcCTGATGAGTTCGCCGATGA‑A‑GTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCg                                                                                                                                      >  6:89570/1‑149 (MQ=255)
              cgcCTGATGAGTTCGCCGATGA‑A‑GTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCATACAAGCGGCCTc                                                                                                                                       <  4:35480/148‑1 (MQ=255)
              cgcCTGATGAGTTCGCCGATGA‑A‑GTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCATACAAGCGGCCTc                                                                                                                                       <  6:258772/148‑1 (MQ=255)
                       aGTTCGCCGATGA‑A‑GTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGaa                                                                                                                              <  2:46766/148‑1 (MQ=255)
                       aGTTCGCCGATGA‑A‑GTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCAAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGaa                                                                                                                              <  8:147015/148‑1 (MQ=255)
                            gCCGATGA‑A‑GTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGAt                                                                                                                           >  8:20027/1‑146 (MQ=255)
                            gCCGATGA‑A‑GTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATga                                                                                                                         <  2:286593/148‑1 (MQ=255)
                            gCCGATGA‑A‑GTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATga                                                                                                                         <  1:144731/148‑1 (MQ=255)
                                     a‑GTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCgg                                                                                                                <  7:189928/149‑1 (MQ=255)
                                     a‑GTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCgg                                                                                                                <  1:291824/149‑1 (MQ=255)
                                          cAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAAt                                                                                                             <  3:307658/148‑1 (MQ=255)
                                           aaCGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATg                                                                                                            <  6:379752/148‑1 (MQ=255)
                                                    tCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATgag                                                                                                                        >  2:283284/1‑127 (MQ=255)
                                                    tCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATgag                                                                                                                        <  1:283284/127‑1 (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   <  4:108517/149‑1 (MQ=255)
                                                                                                                                                   cGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTCGTTGTCGAAAGCCTGGTAAGCGGTGcc   >  5:240200/1‑149 (MQ=255)
                                                                                                                                                    gCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTCGTTGTCGAAAGCCTGGTAAGCGGTGCCg  >  4:239143/1‑149 (MQ=255)
                                                                                                                                                      ttGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTc                                                                 <  6:272921/84‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
GTGTGCGCAAGATGCGCCTGATGAGTTCGCCGATGA‑A‑GTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTCGTTGTCGAAAGCCTGGTAAGCGGTGCCG  >  NC_002947/603510‑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: ‑