Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,191,701 T→C 41.9% L375L (TTA→CTA)  icd → e14 prophage; isocitrate dehydrogenase, specific for NADP+

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,191,7010TC41.9% 4.8 / 8.8 12L375L (TTA→CTA) icde14 prophage; isocitrate dehydrogenase, specific for NADP+
Reads supporting (aligned to +/- strand):  ref base T (4/3);  new base C (3/2);  total (7/5)
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.37e-01

GAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAA  >  CP009273/1191584‑1191846
                                                                                                                     |                                                                                                                                                 
gAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCaa                                                                                                                    >  2:271899/1‑149 (MQ=255)
 aaGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCaat                                                                                                                   >  2:178264/1‑148 (MQ=25)
      acccacGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGaa                                                                                                              >  1:105899/1‑149 (MQ=255)
         cccGGTACCACGGCGAAATATGCCGGTCAGGACAAACTAAATCCTGGCTCTATTATTCTCTCCGCGGAGATGATGCTGCGCCATATGGGTAGTAATGAAGCGGCTGACCTGATTGTTACAGGTATGGAAGGCGCAATCaatgccaagat                                                                                                           <  1:40435/147‑10 (MQ=18)
                             tGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGCTGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTCACTGTTAACGGCATGGAAGGCGCAATCAACGCGAACACCGTAACCTATGACTTCGAGc                                                                                       <  1:106441/149‑1 (MQ=255)
                                        aCAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGCACGTCGgatggaa                                                                            >  2:158847/1‑148 (MQ=11)
                                             gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGgcgc                                                                       >  2:1801/1‑149 (MQ=11)
                                                    cTGGCTCTATTATTCTATCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGAATGGCGCTAAActg                                                                >  2:195569/1‑149 (MQ=255)
                                                    cTGGCTCTATTATGCTCCCCGCTGAGACGATGCTGCGCCATATGGGTTTGACTGAAGCGGCTGACCTTATTGTTAAAGGTATGGAAGGCGCAATCGATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGctg                                                                <  2:177434/149‑1 (MQ=11)
                                                                         cTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGt                                           <  2:105899/149‑1 (MQ=255)
                                                                                                      aCCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCg              >  2:218049/1‑149 (MQ=255)
                                                                                                                  gACTTAATAGTTAAAGGTATGGAAGGCACAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTaaa  <  1:39072/149‑1 (MQ=255)
                                                                                                                     |                                                                                                                                                 
GAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAA  >  CP009273/1191584‑1191846

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