Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,148,465 C→G 9.5% R135G (CGC→GGC)  msuE → FMN reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,148,4650CG 9.5% 121.4 / 5.2 42R135G (CGC→GGC) msuEFMN reductase
Reads supporting (aligned to +/- strand):  ref base C (18/20);  new base G (2/2);  total (20/22)
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.35e-01

GAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCT  >  NC_002947/3148324‑3148604
                                                                                                                                             |                                                                                                                                           
gAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTc                                                                                                                                      >  2:283499/1‑149 (MQ=255)
       aCTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTc                                                                                                                                      >  5:217042/1‑142 (MQ=255)
       aCTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTc                                                                                                                                      <  6:217042/142‑1 (MQ=255)
             tGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCCCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTg                                                                                                                         <  1:283499/149‑1 (MQ=255)
                  aTCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTc                                                                                                                      >  4:173378/1‑147 (MQ=255)
                  aTCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTAGGACGCGCGATACACCgg                                                                                                                             >  3:93798/1‑140 (MQ=255)
                  aTCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTAGGACGCGCGATACACCgg                                                                                                                             <  4:93798/140‑1 (MQ=255)
                   tCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAg                                                                                                                   <  7:331697/149‑1 (MQ=255)
                                cGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGAc                                                                                                       <  5:44419/148‑1 (MQ=255)
                                     gCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCAc                                                                                                                <  8:96328/134‑1 (MQ=255)
                                     gCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCAc                                                                                                                >  7:96328/1‑134 (MQ=255)
                                      ccTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGAc                                                                                                       >  1:305458/1‑142 (MQ=255)
                                      ccTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGAc                                                                                                       <  2:305458/142‑1 (MQ=255)
                                       cTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGa                                                                                               >  7:199726/1‑149 (MQ=255)
                                        tGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGAt                                                                                              >  1:186292/1‑149 (MQ=255)
                                          gCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGAt                                                                                              <  2:186292/147‑1 (MQ=255)
                                           ccGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCAc                                                                                           >  2:263920/1‑149 (MQ=255)
                                             ggTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCa                                                                                         <  3:35367/149‑1 (MQ=255)
                                                  gCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTAGGACGCGCGATACACCgggc                                                                                                                           >  8:234427/1‑108 (MQ=38)
                                                  gCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTAGGACGCGCGATACACCgggc                                                                                                                           <  7:234427/110‑3 (MQ=38)
                                                        gcgcGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAgg                                                                              >  7:113616/1‑149 (MQ=255)
                                                                       ggCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCgtggt                                                                 <  8:113616/147‑1 (MQ=255)
                                                                       ggCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCgtggt                                                                 <  1:263920/147‑1 (MQ=255)
                                                                            aCCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTgcc                                                          >  7:322653/1‑149 (MQ=255)
                                                                              cTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTgccgc                                                        >  2:153520/1‑149 (MQ=255)
                                                                               tGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTgccgc                                                        >  6:346612/1‑148 (MQ=255)
                                                                                   aCCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCg                                                   <  1:153520/149‑1 (MQ=255)
                                                                                                         cAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCtt                             >  1:295238/1‑149 (MQ=255)
                                                                                                           aGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTtga                           <  2:295238/149‑1 (MQ=255)
                                                                                                              cGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTtgatga                        <  3:173378/149‑1 (MQ=255)
                                                                                                                    gCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCggtg                                                >  8:65245/1‑119 (MQ=255)
                                                                                                                    gCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCggtg                                                <  7:65245/119‑1 (MQ=255)
                                                                                                                    gCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATgatcgat                   >  3:249536/1‑148 (MQ=255)
                                                                                                                     cTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATgatcgat                   <  8:199726/147‑1 (MQ=255)
                                                                                                                      tACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCAc                >  5:272225/1‑149 (MQ=255)
                                                                                                                           gTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCGGTGGTagcgagc                                       <  8:132704/121‑1 (MQ=255)
                                                                                                                           gTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGtt           <  4:249536/149‑1 (MQ=255)
                                                                                                                              ggCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCgt                                                                    >  2:24179/1‑89 (MQ=255)
                                                                                                                              ggCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTgc         <  5:346612/148‑1 (MQ=255)
                                                                                                                                    ccGGTGGATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCt  <  3:180872/149‑1 (MQ=255)
                                                                                                                                         gtATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATgatcgat                   <  2:18470/127‑1 (MQ=255)
                                                                                                                                         gtATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATgatcgat                   >  1:18470/1‑127 (MQ=255)
                                                                                                                                             |                                                                                                                                           
GAGGTGCACTGGGTGCGCATCGCCGAACTGTCCGCCAGCCTGGCCGGTTCGCTGGAGCGCGATAGCGCCTCGGCTGACCTGCAACCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCT  >  NC_002947/3148324‑3148604

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