Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,419,933 G→T 66.4% A112D (GCC→GAC)  hisJ ← histidine/lysine/arginine/ornithine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,419,9330GT66.4% 7.8 / 2.6 11A112D (GCC→GAC) hisJhistidine/lysine/arginine/ornithine transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base T (4/3);  total (6/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 = 1.00e+00

CCAATGTTCATTACCGAACGTCTCCTGGGTGGTGCCCTGCAATACGCCTACCCGTTTGCCTTTCAGCGACTCGACTGTCGGCTGAATGTCAGAATTTTTCGCCACCACCAAACGAGAATCGGCAGCGTACAGTTTGTCGGTGAAGGCTATTTCTTGCTGACGTTTTTCCGTAATGGAAAGCGATGACATGATGGCGTCAATCTTCTTCGCTTTTAAGGACGGGATTAACGCATCCAGCGGATTTTCGACAAACGTACAT  >  CP009273/2419812‑2420070
                                                                                                                         |                                                                                                                                         
ccAATGTTCATTACCGAACGTCTCCTGGGTGGTGCCCTGCAATACGCCTACCCGTTTGCCTTTCAGCGACTCGACTGTCGGCTGAATGGCAGAATTTTTCGCCACCACCAACCGAGAATGGTCAGCGCACGGTTGATCGGTGAGGGCttt                                                                                                               >  2:297358/1‑148 (MQ=255)
     gTTCATTACCGAACGTCTCCTGGGTGGTGCCCTGCAATACGCCTACCCGTTTGCCTTTCAGCGACTCGACTGTCGGCTGAATGTCAGAATTTTTCGCCACCACCAAACGAGAATCGTCAGCGTACAGTTTGTCGGTGAAGGCTATTTCtt                                                                                                          <  2:346072/150‑1 (MQ=255)
                           ggtggtGCCCTGCAATACGCCTACCCGTTTGCCTTTCAGCGACTCGACTGTCGGCTGAATGTCAGAATTTTTCGCCACCACCAAACGAGAATCGTCAGCGTACAGTTTGTCGGTGAAGGCTATTTCTTGCTGACGTTTTTCCGTAATGGa                                                                                    >  2:151550/1‑150 (MQ=255)
                                tgCCCTGCAATACGCCTACCCGTTTGCCTTTCAGCGACTCGACTGTCGCCTGAATGTCAGAATTTTTCGCCACCCACAAACGAGAATCGTCAGCGTACAGTTTGTCGGTGCAGGCTATTTCTTGCCGACGTTTTTCGGTTATGGtaaagg                                                                               >  1:235909/1‑147 (MQ=255)
                                   ccTGCAATACGCCTACCCGTTTGCCTTTCAGCGACTCGACTGTCGGCTGAATGTCAGAATTTTTCGCCACCACCAAACGAGAATCGTCAGCGTACAGTTTGTCGGTGAAGGCTATTTCTTGCTGACGTTTTTCCGTAATGGAAAGCGATg                                                                            <  2:195282/150‑1 (MQ=255)
                                        aaTACGCCTACCTGTTTGCCTTTCAGCGACTCGGCTGTCCGCTGAATGTCAGAATTTTTCGCCACCACCAAACGAGAATCGTCAGCGTACAGTTTGTCGGTGAAGGCTATTTCTTGCTGACGTTTTTCCGTAATGGAAAGCGATGACatg                                                                       <  1:26895/150‑1 (MQ=255)
                                                                                tcTGAATGTCGTAAGTTTTCGCCACCACCAAACGACAATCGGCAGCGTACAGTTTGGCGGTGAAGGCTATTTCTTGCTGACGTTTTTCCGTAATGGAAAGCGATGACATGATGGCGTCAATCTTCTTCGCTTTTAAGGACGGGATTAACg                               <  1:49603/149‑1 (MQ=255)
                                                                                     aTGTCAGAATTTTTCGCCACCACCAAACGAGAATCGGCAGCGTACAGTTTGTCGGTGAAGGCTATTTCTTGCTGACGTTTTTCCGTAATGGAAAGCGATGACATGATGGCGTCAATCTTCTTCGCTTTTAAGGACGGGATTAACGCATcc                          >  2:304496/1‑150 (MQ=255)
                                                                                                 ttCGCCACCACCAAACGAGAAAAAGCAGCGTACAGTCAGTCGGTGAAGGCTATTTCATGCTACCGTTTTAACGTACAGGAAAGCGACGACATCGTGGCGTCAAGCATCCTCCCTTTTAAGGACGGGATTAACGAATCCCGCgggcttcgc              >  1:101324/1‑143 (MQ=255)
                                                                                                  tCGCCACCACCAAACGAGAATCGGCAGCGTACAGTTTGTCGGTGAAGGCTATTTCTTGCTGACGTTTTTCCGTAATGGAAAGCGATGACATGATGGCGTCAATCTTCTTCGCTTTTAAGGACGGGATTAACGCATCCAGCGGATTTTCGa             <  2:207420/150‑1 (MQ=255)
                                                                                                             aaaCGAGAATCGTCAGCGTACAGTTTGTCGGTGAAAGCTATTTCTTGCTGACGTTTTTCCGTAATGGAAAGCGATGACATGATGGCGTCAATCTTCTTCGCTTTTAAGAACGGGATTAACGCATCCAGCGGATTTTCGACAAACGTAcct  >  2:180646/1‑148 (MQ=255)
                                                                                                                         |                                                                                                                                         
CCAATGTTCATTACCGAACGTCTCCTGGGTGGTGCCCTGCAATACGCCTACCCGTTTGCCTTTCAGCGACTCGACTGTCGGCTGAATGTCAGAATTTTTCGCCACCACCAAACGAGAATCGGCAGCGTACAGTTTGTCGGTGAAGGCTATTTCTTGCTGACGTTTTTCCGTAATGGAAAGCGATGACATGATGGCGTCAATCTTCTTCGCTTTTAAGGACGGGATTAACGCATCCAGCGGATTTTCGACAAACGTACAT  >  CP009273/2419812‑2420070

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