Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 709,861 C→G 17.9% A201A (GCC→GCG ribF → bifunctional riboflavin kinase/FMN adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947709,8610CG17.9% 63.8 / 10.2 28A201A (GCC→GCGribFbifunctional riboflavin kinase/FMN adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (13/10);  new base G (3/2);  total (16/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-01

GCCAACACGGTCATCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGCCGGT  >  NC_002947/709715‑710005
                                                                                                                                                  |                                                                                                                                                
gccAACACGGTCATCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGc                                                                                                                                               >  2:171133/1‑150 (MQ=255)
                       ccTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCaa                                                                                                                        >  2:116805/1‑150 (MQ=255)
                                   cagcaCCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCaa                                                                                                            <  2:750316/150‑1 (MQ=255)
                                   cagcaCCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCaa                                                                                                            <  1:116805/150‑1 (MQ=255)
                                     gcaCCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAgc                                                                                                          <  1:670979/150‑1 (MQ=255)
                                                  cAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGcc                                                                                             <  2:299342/150‑1 (MQ=255)
                                                   aaGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGGAGCTGGCGGGCCAGCttc                                                                                                                                        <  1:485465/106‑3 (MQ=255)
                                                   aaGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGGAGCTGGCGGGCCAGCttc                                                                                                                                        >  2:485465/1‑104 (MQ=255)
                                                                     gACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTg                                                                          >  2:751325/1‑150 (MQ=255)
                                                                      aCTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTgg                                                                         <  1:259479/150‑1 (MQ=255)
                                                                              cTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATc                                                                 <  2:707191/150‑1 (MQ=255)
                                                                                          ctgctgGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCCACCGCCAACATCCAGCTCAAGAGCCGCCGCGTGCCGCTCAACGGTGTATACCTGGCCAGCATCGTGATCGACGGc                                                     >  2:384001/1‑150 (MQ=255)
                                                                                          ctgctgGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGc                                                     >  1:30570/1‑150 (MQ=255)
                                                                                              tgGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAg                                                 >  1:69748/1‑150 (MQ=255)
                                                                                                    gccCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATc                                                           <  2:156607/134‑1 (MQ=255)
                                                                                                    gccCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATc                                                           >  1:156607/1‑134 (MQ=255)
                                                                                                             gCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCATCGAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCg                                  <  1:384001/150‑1 (MQ=255)
                                                                                                                               tgcACGGCCAGAAGCTGGCCCGCCAGCTTGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCggg                                                                                   <  1:55804/83‑1 (MQ=255)
                                                                                                                               tgcACGGCCAGAAGCTGGCCCGCCAGCTTGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCggg                                                                                   >  2:55804/1‑83 (MQ=255)
                                                                                                                                gcACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGccc               >  2:560832/1‑150 (MQ=255)
                                                                                                                                  aCGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACAGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGAGCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCa             >  1:609989/1‑150 (MQ=255)
                                                                                                                                         cgaGCTGGCGGGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGc      >  2:53643/3‑150 (MQ=255)
                                                                                                                                          gaGCTGGCGGGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCg                                                                                     >  2:509075/2‑70 (MQ=255)
                                                                                                                                          gaGCTGGCGGGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCg                                                                                     <  1:509075/69‑1 (MQ=255)
                                                                                                                                          aaGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTg                                                              >  1:577550/1‑93 (MQ=255)
                                                                                                                                          aaGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTg                                                              <  2:577550/93‑1 (MQ=255)
                                                                                                                                             cTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCTCACCGTTGCCGGt  >  2:128569/1‑150 (MQ=255)
                                                                                                                                             cTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGCCGGt  >  2:183213/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GCCAACACGGTCATCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGCCGGT  >  NC_002947/709715‑710005

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