Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,793,181 T→A 28.6% S222S (TCA→TCT polA ← DNA‑directed DNA polymerase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,793,1810TA28.6% 21.2 / 9.2 14S222S (TCA→TCTpolADNA‑directed DNA polymerase I
Reads supporting (aligned to +/- strand):  ref base T (4/6);  new base A (2/2);  total (6/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 = 1.00e+00

CATCTTGATGAGTCATTAACGTATCTTCAAGTTTTACTTCAATCGGACTATCTACATTAATCGTTGCTAATTCTTTACTCATTAATGCATCTTCTTTGCTATTTTGAAGTTTTTCTTTTAACTTTTTACCTGAAATTTCATCTAAATGTTCATAGACACCTTCTACTGTGTCAAATTGGTTTAGCAATTTTATTGCTGTTTTCTCTCCAACACCTGCAACACCTGGTATATTATCAG  >  CP000730/1793051‑1793287
                                                                                                                                  |                                                                                                          
cATCTTGATGAGTCATTAACGTATCTTCAAGTTTTACTTCAATCGGACTATCTACATTAATCGTTGCTAATTCTTTACTCATTAATGCATCTTCTTTGCTATTTTGAAGTTTTTCTTTTAACTTTTTACCAGAAATTTCAt                                                                                                  >  2:12767/1‑141 (MQ=255)
cATCTTGATGAGTCATTAACGTATCTTCAAGCTTTACTTCAATCGGACTATCTACATTAATCGTTGCTAATTCTTTACTCATTAATGCATCTTCTTTGCTATTTTGAAGTTTTTCTTTTAACTTTTTACCAGAAATTTCAt                                                                                                  >  2:84197/1‑141 (MQ=255)
    ttGATGAGTCATTAACGTATCTTCAAGTTTTACTTCAATCGGACTATCTACATTAATCGTTGCTAATTCTTTACTCATTAATGCATCTTCTTTGCTATTTTGAAGTTTTTCTTTTAACTTTTTACCTGAAATTTCATCTaa                                                                                              <  2:108191/141‑1 (MQ=255)
                tAACGTATCTTCAAGTTTTACTTCAATCGGACTATCTACATTAATCGTTGCTAATTCTTTACTCATTAATGCATCTTCTTTGCTATTTTGAAGTTTTTCTTTTAACTTTTTACCTGAAATTTCATCTAAATGTTCATAGac                                                                                  <  1:62081/141‑1 (MQ=255)
                            aaGTTTTACTTCAATCGGACTATCTACATTAATCGTTGCTAATTCTTTACTCATTAATGCATCTTCTTTGCTATTTTGAAGTTTTTCTTTTAACTTTTTACCTGAAATTTCATCTAAATGTTCATAGACACCTTCTACtgt                                                                      >  2:87535/1‑141 (MQ=255)
                            aaGTTTTACTTCAATCGGACTATCTACATTAATCGTTGCTAATTCTTTACTCATTAATGCATCTTCTTTGCTATTTTGAAGTTTTTCTTTTAACTTTTTACCTGAAATTTCATCTAAATGTTCATAGACACCTTCTACtgt                                                                      >  2:28636/1‑141 (MQ=255)
                                  tACTTCAATCGGACTATCTACATTAATCGTTGCTAATTCTTTACTCATTAATGCATCTTCTTTGCTATTTTGAAGTTTTTCTTTTAACTTTTTACCTGAAATTTCATCTAAATGTTCATAGACACCTTCTACTGTGTCaaa                                                                <  2:4495/141‑1 (MQ=255)
                                                         tAATCTTTGCTAATTCTTTACTCATTAATGCATCTTCTTTGCTATTTTGAAGTTTTTCTTTTAACTTTTTACCTGAAATTTCATCTAAATGTTCATATACACCTTCTACTGTGTCAAATTGGTTTAGCAATTTTATTGCtt                                         >  2:29736/1‑140 (MQ=255)
                                                                   tAATTCTTTACTCATTAATGCATCTTCTTTGCTATTTTGAAGTTTTTCTTTTAACTTTTTACCAGAAATTTCATCTAAATGTTCATAGACACCTTCTACTGTGTCAAATTGGTTTAGCAATTTTATTGCTGTTTTCTCTcc                               <  1:12767/141‑1 (MQ=255)
                                                                   tAATTCTTTACTCATTAATGCATCTTCTTTGCTATTTTGAAGTTTTTCTTTTAACTTTTTACCAGAAATTTCATCTAAATGTTCATAGACACCTTCTACTGTGTCAAATTGGTTTAGCAATTTTATTGCTGTTTTCTCTcc                               <  1:84197/141‑1 (MQ=255)
                                                                             cTCATTAATGCATCTTCTTTGCTATTTTGAAGTTTTTCTTTTAACTTTTTACCTGAAATTTCATCTAAATGTTCATAGACACCTTCTACTGTGTCAAATTGGTTTAGCAATTTTATTGCTGTTTTCTCTCcaacacctgca                     <  1:26218/141‑1 (MQ=255)
                                                                                  tAATGCATCTTCTTTGCTATTTTGAAGTTTTTCTTTTAACTTTTTACCTGAAATTTCATCTAAATGTTCATAGACACCTTCTACTGTGTCAAATTGGTTTAGCAATTTTATTGCTGTTTTCTCTCcaacacctgcaacacc                <  1:28636/141‑1 (MQ=255)
                                                                                  tAATGCATCTTCTTTGCTATTTTGAAGTTTTTCTTTTAACTTTTTACCTGAAATTTCATCTAAATGTTCATAGACACCTTCTACTGTGTCAAATTGGTTTAGCAATTTTATTGCTGTTTTCTCTCcaacacctgcaacacc                <  1:87535/141‑1 (MQ=255)
                                                                                                tGCTATTTTGAAGTTTTTCTTTTAACTTTTTACCTGAAATTTCATCTAAATGTTCATAGACACCTTCTACTGTGTCAAATTGGTTTAGCAATTTTATTGCTGTTTTCTCTCCAACACCTGCAACACCTGGTATATTATCAg  >  1:95482/1‑141 (MQ=255)
                                                                                                                                  |                                                                                                          
CATCTTGATGAGTCATTAACGTATCTTCAAGTTTTACTTCAATCGGACTATCTACATTAATCGTTGCTAATTCTTTACTCATTAATGCATCTTCTTTGCTATTTTGAAGTTTTTCTTTTAACTTTTTACCTGAAATTTCATCTAAATGTTCATAGACACCTTCTACTGTGTCAAATTGGTTTAGCAATTTTATTGCTGTTTTCTCTCCAACACCTGCAACACCTGGTATATTATCAG  >  CP000730/1793051‑1793287

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: