Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 603,663 T→A 9.5% intergenic (+59/‑33) 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,6630TA 9.5% 112.0 / 3.6 42intergenic (+59/‑33)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/26);  new base A (2/2);  total (14/28)
Fisher's exact test for biased strand distribution p-value = 5.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

CGCCTGATGAGTTCGCCGATGAAGTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTCGTTGTCGAAAGCCTGGTAAGCGGTGCCGTTGATG  >  NC_002947/603524‑603810
                                                                                                                                           |                                                                                                                                                   
cgcCTGATGAGTTCGCCGATGAAGTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCg                                                                                                                                            >  6:89570/1‑149 (MQ=255)
cgcCTGATGAGTTCGCCGATGAAGTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCATACAAGCGGCCTc                                                                                                                                             <  6:258772/148‑1 (MQ=255)
cgcCTGATGAGTTCGCCGATGAAGTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCATACAAGCGGCCTc                                                                                                                                             <  4:35480/148‑1 (MQ=255)
         aGTTCGCCGATGAAGTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGaa                                                                                                                                    <  2:46766/148‑1 (MQ=255)
         aGTTCGCCGATGAAGTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCAAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGaa                                                                                                                                    <  8:147015/148‑1 (MQ=255)
              gCCGATGAAGTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGAt                                                                                                                                 >  8:20027/1‑146 (MQ=255)
              gCCGATGAAGTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATga                                                                                                                               <  1:144731/148‑1 (MQ=255)
              gCCGATGAAGTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATga                                                                                                                               <  2:286593/148‑1 (MQ=255)
                      aGTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCgg                                                                                                                      <  1:291824/149‑1 (MQ=255)
                      aGTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCgg                                                                                                                      <  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         <  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)
                                                                                                                                          aTGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTCGTTGTCGAAAGCCTGGTAAGCGGTGCCGTTGATg  >  2:328465/1‑149 (MQ=255)
                                                                                                                                          aTGCGGCCTCGCGCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCGCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTCGTTGTCGAAAGCCTGGTAAGCGGTGCCGTTGAt   <  4:330787/148‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
CGCCTGATGAGTTCGCCGATGAAGTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTCGTTGTCGAAAGCCTGGTAAGCGGTGCCGTTGATG  >  NC_002947/603524‑603810

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: ‑