Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,860,171 C→T 14.3% A308A (GCG→GCA ligA ← DNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,860,1710CT14.3% 64.9 / 4.1 28A308A (GCG→GCAligADNA ligase
Reads supporting (aligned to +/- strand):  ref base C (10/14);  new base T (2/2);  total (12/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01

TTGTGCAGGGTGGCGTTGGACACCGTCACACCGGCCACCTTGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCAGGCTGTTACGCCGCTCGCCAATGTCGCGATAGTAGGCCAGGCACTCTTCGATGC  >  NC_002947/4860022‑4860316
                                                                                                                                                      |                                                                                                                                                 
ttgtGCAGGGTGGCGTTGGACACCGTCACACCGGCCACCTTGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATc                                                                                                                                                   >  1:186672/1‑150 (MQ=255)
  gtgCAGGGTGGCGTTGGACACCGTCACACCGGCCACCTTGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGc                                                                                                                                                 >  1:123136/1‑150 (MQ=255)
   tgCAGGGTGGCGTTGGAGACCGTCACACCGGCCACCTTGACCGGCTTCAGCCGCGCCACAGGCGTTACCGCACCCGTACGGCCAACCTGGAACTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGcc                                                                                                                                                <  2:525967/150‑1 (MQ=255)
   tgCAGGGTGGCGTTGGACACCGTCACACCGGCCACCTTGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGcc                                                                                                                                                <  1:353177/150‑1 (MQ=255)
              gTTGGACACCGTCACACCGGCCACCTTGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTc                                                                                                                                     >  2:43825/1‑150 (MQ=255)
                    cacCGTCACACCGGCCACCTTGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCGATTGCCCacaa                                                                                                                                           >  2:452595/1‑135 (MQ=255)
                    cacCGTCACACCGGCCACCTTGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCGATTGCCCacaa                                                                                                                                           <  1:452595/138‑4 (MQ=255)
                          cacaCCGGCCACCTTGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACTGGCGCGGaa                                                                                                                         <  2:586942/150‑1 (MQ=255)
                             caaggCAACCTTGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGTCACTCCCCCAGCGCGGCTCACGGGCGCGGAAccc                                                                                                                      >  2:100295/4‑150 (MQ=255)
                             acCGGCCACCTTGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGcc                                                                                                                      >  1:74195/1‑150 (MQ=255)
                                ggCCACCTTGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATCACAGACCGCGGGCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAg                                                                                                                   >  1:536962/1‑150 (MQ=255)
                                ggCCACCTTGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGACACCAACGAACTCGACATCAACCACCGCGTTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAg                                                                                                                   >  1:537185/1‑150 (MQ=255)
                                     ccTTGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTcgc                                                                                                              >  1:342316/1‑150 (MQ=255)
                                                ttCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCGCGAGCTGggcttcc                                                                                                   <  2:139191/150‑5 (MQ=255)
                                                        cgcCACCGGCGTTACTGCACCCGTACGGCCCACCTGGACCTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCAAGCTCGCGCTGGGCGCCCAGGCTGtt                                                                                           <  2:130244/150‑1 (MQ=255)
                                                                 cGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGaa                                                                                  <  1:441997/150‑1 (MQ=255)
                                                                          aCCCGTACGGCCAACCTGGAACTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTAcccc                                                                         <  2:186672/150‑1 (MQ=255)
                                                                                aCGGCCAACCTGGAACTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGAt                                                                   >  2:494206/1‑150 (MQ=255)
                                                                                   gCCAACCTGGAACTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGAtctc                                                                >  2:524968/1‑150 (MQ=255)
                                                                                          tGGAACTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGc                                                         >  1:37406/1‑150 (MQ=255)
                                                                                            gAACTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCAg                                                       <  2:537185/150‑1 (MQ=255)
                                                                                            gAACTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCAg                                                       <  2:536962/150‑1 (MQ=255)
                                                                                             aaCTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCAgg                                                      <  1:524968/150‑1 (MQ=255)
                                                                                              aCTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCGATTGCCCacaa                                                                                                                                           <  2:385391/64‑4 (MQ=255)
                                                                                              aCTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCGATTGCCCacaa                                                                                                                                           >  1:385391/1‑61 (MQ=255)
                                                                                                                     cAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCAGGCTGTTACGCCGCTCGCCAATGTc                               <  1:115215/150‑1 (MQ=255)
                                                                                                                                     gggAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCAGGCTGTTACGCCGCTCGCCAATGTCGCGATAGTAGGCCAgg               <  2:250751/150‑1 (MQ=255)
                                                                                                                                                ggCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGAGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCAGGCTGTTACGCCGCTCGCCAATGTCGCGATAGTAGGCCAGGCACTCTTCGAt    <  1:62296/150‑1 (MQ=255)
                                                                                                                                                  cAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCAGGCTGTTACGCCGCTCGCCAATGTCGCGATAGTAGGCCAGGCACTCTTCGATGc  <  2:37406/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
TTGTGCAGGGTGGCGTTGGACACCGTCACACCGGCCACCTTGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCAGGCTGTTACGCCGCTCGCCAATGTCGCGATAGTAGGCCAGGCACTCTTCGATGC  >  NC_002947/4860022‑4860316

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