Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 603,661 T→C 11.1% 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,6610TC11.1% 117.6 / 5.0 45intergenic (+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 (17/23);  new base C (3/2);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGCAAGATGCGCCTGATGAGTTCGCCGATGAAGTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTCGTTGTCGAAAGCCTGGTAAGCGGTGCCGTT  >  NC_002947/603513‑603806
                                                                                                                                                    |                                                                                                                                                 
tgCGCAAGATGCGCCTGATGAGTTCGCCGATGAAGTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCGAGACACCTGTCGACGCCGCTTGt                                                                                                                                                   >  4:186099/1‑149 (MQ=255)
tgCGCAAGATGCGCCTGATGAGTTCGCCGATGAAGTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTACTCGAGTCCAGACCCCTGTCGAGGCCGCTTGt                                                                                                                                                   <  6:181451/149‑1 (MQ=255)
            gcCTGATGAGTTCGCCGATGAAGTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCATGTATGAGGCCTccc                                                                                                                                       >  4:197104/1‑147 (MQ=255)
                gatgaGTTCGCCGATGAAGTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGTTTGTATGCGGCCTCGCTCtt                                                                                                                                   >  6:280077/1‑149 (MQ=255)
                            gATGAAGTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAat                                                                                                                       <  7:137834/149‑1 (MQ=255)
                                  gTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGa                                                                                                                 <  8:288262/149‑1 (MQ=255)
                                    tCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGGCGAGGCCGCATACAAGCGGCCTCgacag                                                                                                                                    <  6:8618/128‑5 (MQ=255)
                                    tCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGGCGAGGCCGCATACAAGCGGCCTCgacag                                                                                                                                    >  5:8618/1‑124 (MQ=255)
                                       aCGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGTTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGAc                                                                                                            <  5:280077/149‑1 (MQ=255)
                                           gATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTTGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCt                                                                                                         >  2:5485/1‑148 (MQ=255)
                                           gATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCATCCGAGTGACTTGTGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTg                                                                                                        <  3:186099/149‑1 (MQ=255)
                                                       cGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCt                                                                                                                                          >  3:179545/1‑103 (MQ=255)
                                                       cGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCt                                                                                                                                          <  4:179545/103‑1 (MQ=255)
                                                        gATCTGGAAGTTGTAGAATACGTGCCCTCCGATTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGaa                                                                                            <  8:271467/148‑1 (MQ=255)
                                                        gATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGCCCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAAAATCGGAATGACCCGGAAGACCCTCGAAg                                                                                           >  6:243540/1‑149 (MQ=255)
                                                        gATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCg                                                                                              >  6:61338/1‑146 (MQ=255)
                                                        gATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGaa                                                                                            <  5:61338/148‑1 (MQ=255)
                                                        gATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGaa                                                                                            <  5:205058/148‑1 (MQ=255)
                                                                    gTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTGGTATGCGGCCTCGCTGTTGAAGATGAGAAAATCGGAATCAACCTGAAGACCCTCGAAGGTACCTTCATCg                                                                               >  5:231133/1‑149 (MQ=255)
                                                                    gTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCGTGTATTCGGCCGAGCTCTTGAAGATGAGAATATAGGAAATACCCTGAAGACCCTCGAAGGTACcatcatag                                                                               >  2:87851/1‑147 (MQ=255)
                                                                     tAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCAt                                                                                                <  5:126951/131‑1 (MQ=255)
                                                                      aGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGc                                                                              <  3:197104/148‑1 (MQ=255)
                                                                                  cTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCt                                                                                                         <  1:206963/109‑1 (MQ=255)
                                                                                  cTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCt                                                                                                         >  2:206963/1‑109 (MQ=255)
                                                                                   tCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGccccc                                                                          <  4:272558/139‑1 (MQ=255)
                                                                                   tCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGccccc                                                                          >  3:272558/1‑139 (MQ=255)
                                                                                         tGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGGCGAGGCCGCATACAAGCGGCCTCg                                                                                                                                        >  1:189284/1‑71 (MQ=37)
                                                                                         tGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGGCGAGGCCGCATACAAGCGGCCTCg                                                                                                                                        <  2:189284/71‑1 (MQ=37)
                                                                                            cTTGAGGCGGCATTGGCGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTtgg                                                       <  5:236668/149‑1 (MQ=255)
                                                                                            cTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTtgg                                                       <  3:145872/149‑1 (MQ=255)
                                                                                                aGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTtggtgg                                                    <  4:42780/148‑1 (MQ=255)
                                                                                                     gcATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGGGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGCCCATCGAACGGACCTTCATCCCCCCCCAACGTCGCGATGCTTTGGTTGTTGGcc                                              >  6:60326/1‑149 (MQ=255)
                                                                                                                 tctGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTgg                                                <  3:203453/135‑1 (MQ=255)
                                                                                                                 tctGCTCGAGTCCAGACCCCTGTCGAGGCCGCGTGTATGCGGCCTCGCTCTTGAAGATGAGAATATAGGAATGACCCTGAAGACCATAGAAGGTACCTTCAAAGCCCCCAAAGGTCGCGATGCTTTGGTGGTTgg                                                >  4:203453/1‑135 (MQ=255)
                                                                                                                     cTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCtt                                 <  5:60326/146‑1 (MQ=255)
                                                                                                                      tCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCGCTTGAAGATGAGAATATAGGAATGAACCAGAAGAACATCGAAGGTACCTTAATAGCCCCCAAAGGTAGCTATGCGTTTGTGGTTGGCCGCTTCAACAGCTTCGtt                             >  4:78717/1‑149 (MQ=255)
                                                                                                                      tCGAGGCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCGCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGGCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTCGtt                             <  5:243540/149‑1 (MQ=255)
                                                                                                                       cGAGTCCAGACCCCTGTCGAGGCCGCTTGTATTCGGCCTAGCTCTTGAAGATTAGAATATAGGAATTACCCTGAAGACCATCGAAGGTACCTTAATAGCCCCCAAAGGTCg                                                                  >  6:303012/1‑111 (MQ=255)
                                                                                                                       cGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCg                                                                  <  5:303012/111‑1 (MQ=255)
                                                                                                                        gAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTCGtt                             <  3:157015/147‑1 (MQ=255)
                                                                                                                                  cccTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGc                                            <  5:201110/122‑1 (MQ=255)
                                                                                                                                   ccTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGCCCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTCGTTGTCGAAAGCCTgg                >  7:234024/1‑149 (MQ=255)
                                                                                                                                    cTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTCGTTGTCGAAAGCCTgg                >  7:1345/1‑148 (MQ=255)
                                                                                                                                    agagcGAGGCCGCATACAAGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTCGTTGTCGAAAGCCTgg                >  2:121543/5‑148 (MQ=255)
                                                                                                                                               gCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTCGTTGTCGAAAGCCTGGTAAGCGGTGcgg    >  5:283467/1‑147 (MQ=255)
                                                                                                                                                cTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTCGTTGTCGAAAGCCTGGTAAGCGGTGCCGt   <  1:87851/149‑1 (MQ=255)
                                                                                                                                                 ttGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTGGTTGTCGAAAGCCTGGTAAGCGGTGCCGtt  <  6:231133/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TGCGCAAGATGCGCCTGATGAGTTCGCCGATGAAGTTCAACGCGATATCCGGATTCGATCTGGAAGTTGTAGAATACGTGCCCTCCGAGTGACTTGAGGCGGCATTGACGCCCTCTGCTCGAGTCCAGACCCCTGTCGAGGCCGCTTGTATGCGGCCTCGCTCTTGAAGATGAGAATATCGGAATGACCCTGAAGACCATCGAAGGTACCTTCATCGCCCCCAAAGGTCGCTATGCTTTGGTGGTTGGCCGCTTCAACAGCTTCGTTGTCGAAAGCCTGGTAAGCGGTGCCGTT  >  NC_002947/603513‑603806

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