Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,150,691 G→A 22.2% V222I (GTC→ATC)  msuE → FMN reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,150,6910GA22.2% 30.7 / 4.9 18V222I (GTC→ATC) msuEFMN reductase
Reads supporting (aligned to +/- strand):  ref base G (11/3);  new base A (2/2);  total (13/5)
Fisher's exact test for biased strand distribution p-value = 5.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01

GCCCGCTATTCGCTTTCTTCCAGGCCCACACCCTGCCGTACGGCCTGTACGCCAGTGTCGAGTCTTTCGACGATCAGCGCCTGGCTGACCCTGCGCAGTTCGAACGGATCGAGCGGGTACTCGATACCGTCGGCGCCTTCTTTCATATCCCGGTCGCCCGCGCGGCCTGACTTTCCGAGTTACAGGAGACATCGCCATGAGCCAACCGATCAAATTCGCCTACTGGGTTCCTAACGTCAGCGGTGGCCTCGTTGTC  >  NC_002947/3150563‑3150818
                                                                                                                                |                                                                                                                               
gcCCGCTATTCGCTTTCTTCCAGGCCCACACCCTGCCGTACGGCCTGTACGCCAGTGTCGAGTCTTTCGACGATCAGCGCCTGGCTGACCCTGCGCAGTTCGAACGGATCGAGCGGGTACTCGATACCGTCGgcgc                                                                                                                          <  2:263194/136‑1 (MQ=255)
                              cccTGCCGTACGGCCTGTACGCCAGTGTCGAGTCTTTCGACGATCAGCGCCTGGCTGACCCTGCGCAGTTCGAACGGATCGAGCGGGTACTCGATACCATCGGCGCCTTCTTTCATATCCCGGTCGCCCGCGCGGc                                                                                            >  1:38322/1‑136 (MQ=255)
                              cccTGCCGTACGGCCTGTACGCCAGTGTCGAGTCTTTCGACGATCAGCGCCTGGCTGACCCTGCGCAGTTCGAACGGATCGAGCGGGTACTCGATACCATCGGCGCCTTCTTTCATATCCCGGTCGCCCGCGCGGc                                                                                            >  1:40785/1‑136 (MQ=255)
                                    cGTACGGCCTGTACGCCAGTGTCGAGTCTTTCGACGATCAGCGCCTGGCTGACCCTGCGCAGTTCGAACGGATCGAGCGGGTACTCGATACCGTCGGCGCCTTCTTTCATATCCCGGTCGCCCGCGCGGCCTGACt                                                                                      >  1:450051/1‑136 (MQ=255)
                                    cGTACGGCCTGTACGCCAGTGTCGAGTCTTTCGACGATCAGCGCCTGGCTGACCCTGCGCAGTTCGAACGGATCGAGCGGGTACTCGATACCGTCGGCGCCTTCTTTCATATCCCGGTCGCCCGCGCGGCCTGACt                                                                                      >  2:292852/1‑136 (MQ=255)
                                      tACGGCCTGTACGCCAGTGTCGAGTCTTTCGACGATCAGCGCCTGGCTGACCCTGCGCAGTTCGAACGGATCGAGCGGGTACTCGATACCATCGGCGCCTTCTTTCATATCCCGGTCGCCCGCGCGGCCTGACttt                                                                                    <  2:40785/136‑1 (MQ=255)
                                      tACGGCCTGTACGCCAGTGTCGAGTCTTTCGACGATCAGCGCCTGGCTGACCCTGCGCAGTTCGAACGGATCGAGCGGGTACTCGATACCATCGGCGCCTTCTTTCATATCCCGGTCGCCCGCGCGGCCTGACttt                                                                                    <  2:38322/136‑1 (MQ=255)
                                                aCGCCAGTGTCGAGTCTTTCGACGATCAGCGCCTGGCTGACCCTGCGCAGTTCGAACGGATCGAGCGGGTACTCGATACCGTCGGCGCCTTCTTTCATATCCCGGTCGCCCGCGCGGCCTGACTTTCCGAGTTACa                                                                          >  2:414990/1‑136 (MQ=255)
                                                    cAGTGTCGAGTCTTTCGACGATCAGCGCCTGGCTGACCCTGCGCAGTTCGAACGGATCGAGCGGGTACTCGATACCGTCGGCGCCTTCTTTCATATCCCGGTCGCCCGCGCGGCCTGACTTTCCGAGTTACAGgag                                                                      <  1:126065/136‑1 (MQ=255)
                                                        gtCGAGTCTTTCGACGATCAGCGCCTGGCTGACCCTGCGCAGTTCGAACGGATCGAGCGGGTACTCGATACCGTCGGCGCCTTCTTTCATATCCCGGTCGCCCGCGCGGCCTGACTTTCCGAGTTACAGGAGACAt                                                                  >  2:460165/1‑136 (MQ=255)
                                                               cTTTCGACGATCAGCGCCTGGCTGACCCTGCGCAGTTCGAACGGATCGAGCGGGTACTCGATACCGTCGGCGCCTTCTTTCATATCCCGGTCGCCCGCGCGGCCAGACTTTCCGACTTACAGGAGACATAGCCATg                                                           >  1:326494/1‑136 (MQ=255)
                                                                                             cgcAGTTCGAACGGATCGAGCGGGTACTCGATAGCGTCGGCGCCTTCTTTCATATCCCGGTCGCCCGCGCGGCCTGACTTTCCGAGTTACAGGAGACATCGCCATGAGCCAACCGATCAAATTCGCCTACTGGGtt                             <  1:258273/136‑1 (MQ=255)
                                                                                                aGTTCGAACGGATCGAGCGGGTACTCGATACCGTCGGCGCCTTCTTTCATATCCCGGTCGCCCGCGCGGCCTGACTTTCCGAGTTTCAGGAGACATCGCCATGAGCCAACCGATCAAATTAGCCTACTGGGTTGCt                          >  2:305273/1‑136 (MQ=255)
                                                                                                aGTTCGAACGGATCGAGCGGGTACTCGATACCGTCGGCGCCTTCTTTCATATCCCGGTCGCCCGCGCGGCCTGACTTTCCGAGTTACAGGAGACATCGCCATGAGCCAACCGATCAAATTCGCCTACTGGGTTCCt                          >  1:443805/1‑136 (MQ=255)
                                                                                                               aGCGGGTACTCGATACCGTCGGCGCCTTCTTTCATATCCCGGTCGCCCGCGCGGCCTGACTTTCCGAGTTACAGGAGACATCGCCATGAGCCAACCGATCAAATTCGCCTACTGGGTTCCTAACGTCAGCGGTggg           >  2:378843/1‑135 (MQ=255)
                                                                                                               aGCGGGTACTCGATACCGTCGGCGCCTTCTTTCATATCCCGGTCGCCCGCGCGGCCTGACGTTCCGAGTTACAGGAGACCTCGCCATGAGCCAACCGCTCAAATTCGCCTTCTGGGTTCCTAACGTCAGCGGTGGc           >  1:174880/1‑136 (MQ=255)
                                                                                                                        tCGATACCGTCGGCGCCTTCTTTCATATCCCGGTCGCCCGCGCGGCCTGACTTTCCGAGTTACAGGAGACATCGCCATGAGCCAACCGATCAAATTCGCCTACTGGGTTCCTAACGTCAGCGGTGGCCTCGTTGTc  >  2:316354/1‑136 (MQ=255)
                                                                                                                        tCGATACCGTCGGCGCCTTCTTTCATATCCCGGTCGCCCGCGCGGCCTGACATTACGAGTTAAAGGAGACATAGCCATGAGCCAACCGAACAAATTAGCCTAAAGGGTTACTAACGTAAGCGGTTGCCTCGttttt  >  2:405324/1‑133 (MQ=255)
                                                                                                                                |                                                                                                                               
GCCCGCTATTCGCTTTCTTCCAGGCCCACACCCTGCCGTACGGCCTGTACGCCAGTGTCGAGTCTTTCGACGATCAGCGCCTGGCTGACCCTGCGCAGTTCGAACGGATCGAGCGGGTACTCGATACCGTCGGCGCCTTCTTTCATATCCCGGTCGCCCGCGCGGCCTGACTTTCCGAGTTACAGGAGACATCGCCATGAGCCAACCGATCAAATTCGCCTACTGGGTTCCTAACGTCAGCGGTGGCCTCGTTGTC  >  NC_002947/3150563‑3150818

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: