Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 603,661 T→C 10.4% 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,6610TC10.4% 128.0 / 5.3 48intergenic (+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 (21/22);  new base C (3/2);  total (24/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

GCGCAAGATGCGCCTGATGAGTTCGCCGATGAAGTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTCGTTGTCGAAAGCCTGGTAAGCGGTGCCGTTGA  >  NC_002947/603514‑603808
                                                                                                                                                   |                                                                                                                                                   
gcgcAAGATGCGCCTGATGAGTTCGCCGATGAAGTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGACCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCAAGACCCCTGTCGAGGCAGCTTGTa                                                                                                                                                    >  6:153995/1‑149 (MQ=255)
           gcCTGATGAGTTCGCCGATGAAGTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGc                                                                                                                                         >  7:276618/1‑149 (MQ=255)
             cTGATGAGTTCGCCGATGAAGTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGctc                                                                                                                                       <  4:177699/149‑1 (MQ=255)
             cTGATGAGTTCGCCGATGAAGTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCATACAAGCGGCCTCgaca                                                                                                                                       <  5:66019/149‑4 (MQ=255)
                atgaGTTCGCCGATGAAGTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGctct                                                                                                                                      <  1:170817/147‑1 (MQ=255)
                atgaGTTCGCCGATGAAGTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCtt                                                                                                                                     <  7:123118/148‑1 (MQ=255)
                    gTTCGCCGATGAAGTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGa                                                                                                                                <  8:276618/149‑1 (MQ=255)
                          cGATGAAGTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCATACAAGCGGCCTc                                                                                                                                           >  8:277157/1‑132 (MQ=255)
                          cGATGAAGTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCATACAAGCGGCCTc                                                                                                                                           <  7:277157/132‑1 (MQ=255)
                           gATGAAGTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGaa                                                                                                                          >  2:56256/1‑148 (MQ=255)
                              gAAGTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAata                                                                                                                        >  6:262020/1‑147 (MQ=255)
                               aaGTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATc                                                                                                                      <  4:302531/148‑1 (MQ=255)
                                     aaCGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGATGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGa                                                                                                               <  3:293133/149‑1 (MQ=255)
                                           atatCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTg                                                                                                          >  2:27798/1‑148 (MQ=255)
                                              tCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGa                                                                                                      >  1:36459/1‑149 (MQ=255)
                                              tCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGa                                                                                                      <  1:27798/149‑1 (MQ=255)
                                                 ggATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCa                                                                                                   <  3:162490/149‑1 (MQ=255)
                                                   aTTCGATCTTGAAGTTTTATAATACGTGCCCTCAGAGTGACTTGAGGCTGAATTTACTCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAATATGAGAATATCGGAATGACCCTGAATACCAt                                                                                                  <  5:153995/148‑1 (MQ=255)
                                                       gATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTa                                                                                                                                                    <  6:4179/94‑1 (MQ=255)
                                                       gATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTa                                                                                                                                                    >  5:4179/1‑94 (MQ=255)
                                                        aTCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGAAACCTCGCTCTTGAAGATGAGAATTTCGGAATGACCCTGAAGACCAACGAAgg                                                                                            >  1:259759/1‑149 (MQ=255)
                                                           tGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGAcc                                                                                                    <  7:33488/138‑1 (MQ=255)
                                                           tGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCCCGTGTTTGCGGCCACGCGGTTGAAGATGAGAATATAGGAATTACCCAGAAGACCaaaaaagg                                                                                            >  8:33488/1‑139 (MQ=255)
                                                             gAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAg                                                                                             <  5:262020/143‑1 (MQ=255)
                                                               aGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCa                                                                                                   <  8:282099/135‑1 (MQ=255)
                                                               aGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCa                                                                                                   >  7:282099/1‑135 (MQ=255)
                                                                         tACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGccccc                                                                            >  5:227305/1‑148 (MQ=255)
                                                                            gTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGctct                                                                                                                                      <  8:92038/87‑1 (MQ=255)
                                                                            gTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGctct                                                                                                                                      >  7:92038/1‑87 (MQ=255)
                                                                                            ttGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAg                                                                                                       <  1:157997/102‑1 (MQ=255)
                                                                                             aatggcggcATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTtg                                                          <  1:271680/143‑1 (MQ=255)
                                                                                               aGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCCAGGCCCCTTGTATGCAGCCAAGCTCTACAAAATGATAATATAGGAATGAAACTAAAGACCAACCAAGGTACCTTCATCGCCCCCAAACGTCCCAATGCATTggtggt                                                     >  6:321463/1‑149 (MQ=255)
                                                                                                ggcggcATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGtt                                                    <  2:36459/149‑1 (MQ=255)
                                                                                                     cATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGGCCCGTGTATTGGGGCGCGCGCATTAAAATGAGAATAAAGGAATGGCCCAGGAGACCCGAGAAAGTAAAATGATAGCCCCCAAACGTTGCGAAGCTTTTGTGGTTGGACg                                               >  6:272656/1‑149 (MQ=255)
                                                                                                      aTTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGGCTCGCTCTTGCAGATGAGAATATAGGGATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTAGGTATGATTTGGTGGTTGGCCGc                                              >  5:145776/1‑149 (MQ=255)
                                                                                                      aTTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCGTGTATGCGGCCTAGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGc                                              >  6:138593/1‑149 (MQ=255)
                                                                                                            gCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTAGCTCTTGAAGATGAGAATATAGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCaa                                         >  5:72697/1‑148 (MQ=255)
                                                                                                              ccTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGt                                                                                           <  1:47524/96‑1 (MQ=255)
                                                                                                              ccTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGt                                                                                           >  2:47524/1‑96 (MQ=255)
                                                                                                                       gAGTCCATACGCCTTTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCTATGGTACCTTCATCGCGCGCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTCGTTg                              <  5:321463/148‑1 (MQ=255)
                                                                                                                           ccAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAATATGAGAATATCGGAATTACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTCGTTGTCGa                          <  2:259759/148‑1 (MQ=255)
                                                                                                                            cAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATAGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTaaa                                         >  7:207553/1‑132 (MQ=255)
                                                                                                                                 cccTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTCGTTGTCGAAAGCCTg                   <  6:227305/149‑1 (MQ=255)
                                                                                                                                    tGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACa                                       >  1:74661/1‑126 (MQ=255)
                                                                                                                                    tGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACa                                       <  2:74661/126‑1 (MQ=255)
                                                                                                                                     agcGAGGCCGCATACAAGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTCGTTGTCGAAAGCCTGGt                 >  3:155245/3‑147 (MQ=255)
                                                                                                                                      gcGAGGCCGCATACAAGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTCGTTGTCGAAAGCCTGGTAAg              >  6:67206/2‑149 (MQ=255)
                                                                                                                                              gCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTCGTTGTCGAAAGCCTGGTAAGCGGTGCCg      >  4:59839/1‑149 (MQ=255)
                                                                                                                                                   tATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTCGTTGTCGAAAGCCTGGTAAGCGGTGCCGTTGa  <  1:56256/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCGCAAGATGCGCCTGATGAGTTCGCCGATGAAGTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTCGTTGTCGAAAGCCTGGTAAGCGGTGCCGTTGA  >  NC_002947/603514‑603808

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