Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,376,784 G→A A251T (GCG→ACG)  ycjR → 3‑dehydro‑D‑guloside 4‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,376,7840GA94.4% 41.6 / NA 18A251T (GCG→ACG) ycjR3‑dehydro‑D‑guloside 4‑epimerase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (9/8);  minor base C (1/0);  total (10/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

GGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCGGATAACCATCGCTACCAGCCGGGCAGCGGCACCCTGGATTTCCACGCGCTGTTTGAACAGCTGCGCGCGGATAACTATCAGGGCTATGTAGTGTACGAAGGGCGTATCCGGGCGGAAGATCCTGCCCAGGCGTACCGTGATTCGTTGGCCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAG  >  NC_000913/1376565‑1376959
                                                                                                                                                                                                                           |                                                                                                                                                                               
ggCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCATGCGCTGCATGACAACCGCGTCCTGCTCGGTCATGTGCATATTGCGGATAACCATCTCTACCAGCCGGTCAGCGGCACCCTGGATTTCCACGCGCGGTTTGAACAGCTGCTCACGGGTAACTATCAGGGCTATGTAGTGTACGAAGGGCGTATCCGGGCGGAAGATCCTGCCCAGACGTACCGTGATTCGTTGGCCTGGTTGc                                                                                                                                                      <  1:353323/247‑1 (MQ=255)
    aTTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCGGATAACCATCGCTACCAGCCGGGCAGCGGCACCCTGGATTTCCACGCGCTGTTTGAACAGCTGCGCGCGGATAACTATCAGGGCTATGTAGTGTACGAAGGGCGTATCCGGGCGGAAGATCCTGCCCAGACGTACCGTGATTCGTTGGCCTGGTTGCGTACCTgn                                                                                                                                              <  1:320109/251‑2 (MQ=255)
                     atCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCGGATAACCATCGCTACCAGCCGGGCAGCGGCACCCTGGATTTCCACGCGCTGTTTGAACAGCTGCGCGCGGATAACTATCAGGGCTATGTAGTGTACGAAGGGCGTATCCGGGCGGAAGATCCTGCCCAGACGTACCGTGATTCGTTGGCCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGa                                                                                                                              >  1:35062/1‑250 (MQ=255)
                          agaagaTAACCGGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCGGATAACCATCGCTACCAGCCGGGCAGCGGCACCCTGGATTTCCACGCGCTGTTTGAACAGCTGCGCGCGGATAACTATCAGGGCTATGTAGTGTACGAAGGGCGGATCCGGGTGGAAGATTCTGTCCAGACGTACCGTGATTCGTTGGCCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTg                                                                                                                        <  1:370533/251‑1 (MQ=255)
                                                    aTGACAACCGCGACCTGCTCGGTCATGTGCATATTGCGGATAACCATCGCTACCAGCCGGGCAGCGGCACCCTGGATTTCCACGCGCTGTTTGAACAGCTGCGCGCGGATAACTATCAGGGCTATGTAGTGTACGAAGGGCGTATCCGGGCGGAAGATCCTGCCCAGACGTACCGTGATTCGTTGGCCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTc                                                                                              <  1:136366/251‑1 (MQ=255)
                                                                         gTCATGTGCCTATTGCGGATAACCATCGCTACCAGCCGGGCAGCGGCACCCTGGATTTCCACGCGCTGTTTGAACAGCTGCGCGCGGATAACTATCAGGGCTATGTAGTGTACGAAGGGCGTATCCGGGCGGAAGATCCTGCCCAGACGTACCGTGATTCGTTGGCCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACGAGCTCTCCGCTGCGGGTCGCGATAATAgg                                                                                   >  1:386235/1‑241 (MQ=255)
                                                                                    atTGCGGATAACCATCGCTACCAGCCGGGCAGCGGCACCCTGGATTTCCACGCGCTGTTGGAACAGCTGCGCGCGGATAACTTTCAGGGCTATGTAGTGTACGAAGGGCGTATCCGGGTGGAAGATTCTGTCCAGACGTACCGTGATTCGTTGGCCGGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGc                                                                                                                       <  1:437507/194‑1 (MQ=255)
                                                                                          gATAACCATCGCTACCAGCCGGGCAGCGGCACCCTGGATTTCCACGCGCTGTTTGAACAGCTGCGCGCGGATAACTATCCGGGCTATGTAGTGTACGAAGGGCGTATCCGGGCGGCAGATCCTGCCCAGACGTACCGTGATTCGTTGGCCTGGTTGCGTACCTg                                                                                                                                               >  1:291258/1‑164 (MQ=255)
                                                                                              aCCATCGCTACCAGCCGGGCAGCGGCACCCTGGATTTCCACGCGCTGTTTGAACAGCTGCGCGCGGATAACTATCAGGGCTATGTAGTGTACGAAGGGCGTATCCGGGCGGAAGATCCTGCCCAGACGTACCGTGATTCGTTGGCCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAg                                                             >  1:13776/1‑242 (MQ=255)
                                                                                                         cAGCCGGGCAGCGGCACCCTGGATTTCCACGCGCTGTTTGAACAGCTGCGCGCGGATAACTATCAGGGCTATGTAGTGTACGAAGGGCGTATCCGGGCGGAAGATCCTGCCCAGACGTACCGTGATTCGTTGGCCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTactactn                                         <  1:305260/251‑2 (MQ=255)
                                                                                                                          ccTGGATTTCCACGCGCTGTTTGAACAGCTGCGCGCGGCTAACTATCAGGGCTATGTAGTGTACGAAGGGCGTATCCGGGCGGAAGATCCTGCCCAGACGTACCGTGATTCGTTGGCCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCCGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTgg                        >  1:174538/1‑251 (MQ=255)
                                                                                                                           cTGGATTTCCACGCGCTGTTTGAACAGCTGCGCGCGGATACCTATCAGGGCTATGTAGTGTCCGAAGGGCGTCTCAGGGCGGAAGATCCTGCCCAGACGAACCGTGATTCGTTGGCCTGGTTGCGt                                                                                                                                                    >  1:321984/1‑126 (MQ=255)
                                                                                                                              gATTTCCACGCGCTGTTTGAACAGCTGCGCGCGGATAACTATCAGGGCTATGTAGTGTACGAAGGGCGTATCCGGGCGGAAGATCCTGCCCAGACGTACCGTGATTCGTTGGCCTGGTTGCGTACCTGCTACGCGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGCTAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGtac                                              >  1:488987/1‑225 (MQ=255)
                                                                                                                                   ccACGCGCTGTTTGAACAGCTGCGCGCGGATAACTATCAGGGCTATGTAGTGTACGAAGGGCGTATCCGGGCGGAAGATCCTGCCCAGCCGTACCGTGATTCGTTGGCCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCg                                                                                                       >  1:413084/1‑163 (MQ=255)
                                                                                                                                                gAACAGCTGCGCGCGGATAACTATCAGGGCGATGTAGTGTACGAAGGGCGTATCCGGGCGGAAGATCCGGCCCAGACGTACCGTGATTCGTTCGCCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAATGCAATGGCAGGCTCCCCGCTTCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTGCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAg  <  1:225922/251‑1 (MQ=255)
                                                                                                                                                           cgcgGATAACTATCAGGGCTATGTAGTGTACGAAGGGCGTATCCGGGCGGAAGATCCTGCCCAGACGTACCGTGATTCGTTGGCCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCnctcn                                                                                                        >  2:318758/1‑137 (MQ=255)
                                                                                                                                                           cgcgGATAACTATCAGGGCTATGTAGTGTACGAAGGGCGTATCCGGGCGGAAGATCCTGCCCAGACGTACCGTGATTCGTTGGCCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACtggtggn              <  1:318758/228‑2 (MQ=255)
                                                                                                                                                                           ggCTATGTAGTGTACGAAGGGCGTATCCGGGCGGAAGATCCTGCCCAGACGTACCGTGATTCGTTGGCCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATCGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGc              >  1:490011/1‑212 (MQ=255)
                                                                                                                                                                                                                           |                                                                                                                                                                               
GGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCGGATAACCATCGCTACCAGCCGGGCAGCGGCACCCTGGATTTCCACGCGCTGTTTGAACAGCTGCGCGCGGATAACTATCAGGGCTATGTAGTGTACGAAGGGCGTATCCGGGCGGAAGATCCTGCCCAGGCGTACCGTGATTCGTTGGCCTGGTTGCGTACCTGCTAAGAGGTCTTTGTGAAAAGTGCAATGACAAGCTCTCCGCTGCGGGTCGCGATAATAGGCGCAGGCCAGGTGGCGGATAAGGTTCATGCTTCGTACTACTGCACCCGCAACGATCTGGAACTGGTGGCTGTCTGTGACAG  >  NC_000913/1376565‑1376959

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