Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,860,168 A→G 13.3% I309I (ATT→ATC ligA ← DNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,860,1680AG13.3% 66.5 / 3.4 30I309I (ATT→ATCligADNA ligase
Reads supporting (aligned to +/- strand):  ref base A (12/14);  new base G (2/2);  total (14/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 = 8.49e-01

GTTGTGCAGGGTGGCGTTGGACACCGTCACACCGGCCACCTTGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCAGGCTGTTACGCCGCTCGCCAATGTCGCGATAGTAGGCCAGGCACTCTTCGATGC  >  NC_002947/4860021‑4860316
                                                                                                                                                    |                                                                                                                                                    
gttgtGCAGGGTGGCGTTGGACACCGTCACACCGGCCACCTTGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAAt                                                                                                                                                    >  1:250751/1‑150 (MQ=255)
 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)
                                                                                                                                                    |                                                                                                                                                    
GTTGTGCAGGGTGGCGTTGGACACCGTCACACCGGCCACCTTGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATCCAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCAGGCTGTTACGCCGCTCGCCAATGTCGCGATAGTAGGCCAGGCACTCTTCGATGC  >  NC_002947/4860021‑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: