Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,431,130 C→T 26.6% D29D (GAC→GAT PP_3042 → terminase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,431,1300CT26.6% 22.0 / 5.0 15D29D (GAC→GATPP_3042terminase large subunit
Reads supporting (aligned to +/- strand):  ref base C (5/6);  new base T (2/2);  total (7/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

ATTGGCCAGAAACAATCCACTAAACAGAGCTGAGCCTATGCCCACCGGATACGCAAACGGTGCAGAGGTGTACCGCGAAGCGTTTTGCCGAGGGCTGACGCCCGACCCTGATCTATGGGTCGACGAGTGGGCAGATGAGTACATGCGAATCCCGCGTGGTACCGGTGCCGCTGAGCCTGGCAAGTACCGCACCGCGCGTACGCCGTATGCCCGAGAACCTATGCGTTGCCTGT  >  NC_002947/3431007‑3431239
                                                                                                                           |                                                                                                             
aTTGGCCAGAAACAATCCACTAAACAGAGCTGAGCCTATGCCCACCGGATACGCAAACGGTGCAGAGGTGTACCGCGAAGCGTTTTGCCGAGGGCTGACGCCCGACCCTGATCTATGGGTCGACGAGTGGGCAGAt                                                                                                   >  2:193811/1‑136 (MQ=255)
               tCCACTAAACAGAGCTGAGCCTATGCCCACCGGATACGCAAACGGTGCAGAGGTGTACCGCGAAGCGTTTTGCCGAGGGCTGACGCCAGACCCTGATCTATGGGTCGACGAGTGGGCAGATGAGTACATGCGAATc                                                                                    <  2:122203/136‑1 (MQ=255)
                   atAAAGAGAGCTGAGCCTATGCCCACCGGATACGCAAACGGTGCAGAGGTGTACCGCGATGCGTTTTGCCGAGGGCTGACGCCCGACCCTGATCTTTGGGTCGACGAGTGGGCAGATGAGTACATGCGAATCCcgc                                                                                <  1:317460/135‑1 (MQ=255)
                                        cccACCGGATACGCAAACGGTGCAGAGGTGTACCGCGAAGCGTTTTGCCGAGGGCTGACGCCCGACCCTGATCTATGGGTCGATGAGTGGGCAGATGAGTACATGCGAATCCCGCGTGGTACCGGTGCCGCTGAGc                                                           >  2:435988/1‑136 (MQ=255)
                                        cccACCGGATACGCAAACGGTGCAGAGGTGTACCGCGAAGCGTTTTGCCGAGGGCTGACGCCCGACCCTGATCTATGGGTCGATGAGTGGGCAGATGAGTACATGCGAATCCCGCGTGGTACCGGTGCCGCTGAGc                                                           >  2:438534/1‑136 (MQ=255)
                                                          ggTGCAGAGGTGTACCGCGAAGCGTTTTGCCGAGGGCTGACGCCCGACCCTGATCTATGGGTCGACGAGTGGGCAGATGAGTACATGCGAATCCCGCGTGGTACCGGTGCCGCTGAGCCTGGCAAGTaccgcaccg                                         >  2:333788/1‑136 (MQ=255)
                                                                       aCCGCGAAGCGTTTTGCCGAGGGCTGACGCCCGACCCTGATCTATGGGTCGACGAGTGGGCa                                                                                                      >  1:366691/1‑62 (MQ=255)
                                                                       aCCGCGAAGCGTTTTGCCGAGGGCTGACGCCCGACCCTGATCTATGGGTCGACGAGTGGGCa                                                                                                      <  2:366691/62‑1 (MQ=255)
                                                                         cgcgAAGCGTTTTGCCGAGGGCTGACGCCCGACCCTGATCTATGGGTCGACGAGTGGGCAGATGAGTACATGCGAATCCCGCGTGGTACCGGTGCCGCTGAGCCt                                                         >  1:377590/1‑105 (MQ=255)
                                                                         cgcgAAGCGTTTTGCCGAGGGCTGACGCCCGACCCTGATCTATGGGTCGACGAGTGGGCAGATGAGTACATGCGAATCCCGCGTGGTACCGGTGCCGCTGAGCCt                                                         <  2:377590/105‑1 (MQ=255)
                                                                            gAAGCGTTTTGCCGAGGGCTGACGCCCGACCCTGATCTATGGGTCGATGAGTGGGCAGATGAGTACATGCGAATCCCGCGTGGTACCGGTGCCGCTGAGCCTGGCAAGTACCGCACCGCGCGTACGCCGTATGccc                       <  1:435988/136‑1 (MQ=255)
                                                                            gAAGCGTTTTGCCGAGGGCTGACGCCCGACCCTGATCTATGGGTCGATGAGTGGGCAGATGAGTACATGCGAATCCCGCGTGGTACCGGTGCCGCTGAGCCTGGCAAGTACCGCACCGCGCGTACGCCGTATGccc                       <  1:438534/136‑1 (MQ=255)
                                                                                        cGAGGGCTGACGCCCGACCCTGATCTATGGGTCGACGAGTGGGCAGATGAGTACATGCGAATCCCGCGTGGTACCGGTGCCGCTGAGCCt                                                         <  1:96651/90‑1 (MQ=255)
                                                                                        cGAGGGCTGACGCCCGACCCTGATCTATGGGTCGACGAGTGGGCAGATGAGTACATGCGAATCCCGCGTGGTACCGGTGCCGCTGAGCCt                                                         >  2:96651/1‑90 (MQ=255)
                                                                                                 aCGCCCGACCCTGATCTATGGGTCGACGAGTGGGCAGATGAGTACATGCGAATCCCGCGTGGTACCGGTGCCGCTGAGCCTGGCAAGTACCGCACCGCGCGTACGCCGTATGCCCGAGAACCTATGCGTTGCCTGt  <  1:193811/136‑1 (MQ=255)
                                                                                                                           |                                                                                                             
ATTGGCCAGAAACAATCCACTAAACAGAGCTGAGCCTATGCCCACCGGATACGCAAACGGTGCAGAGGTGTACCGCGAAGCGTTTTGCCGAGGGCTGACGCCCGACCCTGATCTATGGGTCGACGAGTGGGCAGATGAGTACATGCGAATCCCGCGTGGTACCGGTGCCGCTGAGCCTGGCAAGTACCGCACCGCGCGTACGCCGTATGCCCGAGAACCTATGCGTTGCCTGT  >  NC_002947/3431007‑3431239

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