Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,419,720 T→C 64.7% D183G (GAT→GGT)  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,7200TC64.7% 8.8 / 9.2 17D183G (GAT→GGT) hisJhistidine/lysine/arginine/ornithine transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (3/3);  new base C (4/7);  total (7/10)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

TTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATA‑CGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGAACGTCTCCTGGGTGGTGCCCTGCAATACGCCTACCCG  >  CP009273/2419578‑2419865
                                                                                                                                              |                                                                                                                                                  
ttCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTCGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGaa                                                                                                                                             <  2:238461/150‑1 (MQ=255)
             ttATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCGTCAATa                                                                                                                                >  2:49603/1‑150 (MQ=255)
                              gCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAACTTGTTATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCACCCTGGAACGCGGCATCAATA‑CGTCCGGCAgtcaggtc                                                                                                              <  2:84020/150‑1 (MQ=255)
                                ccATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCACCCTGGAACGCGGCATCAATA‑CGTCCGGCAgtcaggtcag                                                                                                            <  2:256456/150‑1 (MQ=255)
                                                                  tttAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCCCCGCGATGAGCGACCTCACCCTGGAACGCGGCATCAATA‑CGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGAt                                                                          <  2:261583/150‑1 (MQ=255)
                                                                               ccACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCACCCTGGAACGCGGCATCAATA‑CGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCtttt                                                             >  1:195282/1‑150 (MQ=255)
                                                                                                   aCCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTTACCCTGGAACGCGGCATCAATA‑CGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTAcc                                         <  2:206414/150‑1 (MQ=255)
                                                                                                       aCAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATA‑CGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGAAc                                     <  2:128145/150‑1 (MQ=255)
                                                                                                        cAGGTTGTTTGAGGAATCCTTCGCTGGCAGCGACCTCACCCTGGAACGCGGCATCAATA‑CGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGTCTTTTGGTGCCCAATGTTCATTACCGAACg                                    <  2:214679/150‑1 (MQ=255)
                                                                                                          ggTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCCCCCTGGAACGCGGCATCAATA‑CGTCCGCCAGTCAGGTCAGAACAAATGTTGTCCTGCCCCTGATACGCGACGATTTCAATGCCTGTTGAGGCCCAATGTCCATTACCGAACGtc                                  >  1:23919/1‑150 (MQ=255)
                                                                                                          ggTTGTTTGAGGAAACCTCCGATGCCAGCGACCTCACCCTGGAACGCGGCATCAATA‑CGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCACAGCTTTTTGGTGCCCAATGTCCACTACCGaagcat                                  >  1:23934/1‑146 (MQ=255)
                                                                                                           gttgttTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATA‑CGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGAACGtct                                 >  1:207420/1‑150 (MQ=255)
                                                                                                                           ttCGCTGGCAGCGACCTCACCCTGGAACGCGGCATCAATA‑CGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGAACGTCTCCTGGGTGGTGCCCTg                 <  2:283109/150‑1 (MQ=255)
                                                                                                                            tCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATA‑CGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCTTTACCGAACGTCTCCTGGGTGGTGCCCTGc                <  2:252724/150‑1 (MQ=255)
                                                                                                                              gCTGGCAGCGACCTCACCCTGGAACGCAG‑ATAATGATTTTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGGTCATTACCGACCGTCTCCTGGGTGGTGCCCGGCaa              >  1:346072/1‑150 (MQ=255)
                                                                                                                                  gCAGCGGCCTCACCCTGGAACGCGGCATCAATA‑CGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCACTGCCTTTTGGTGCCCAATGTTCATTACCGAACGTCTCCTGGGTGGTGCCCTGCAATACg          <  2:105175/150‑1 (MQ=255)
                                                                                                                                          cTCATCCTGGAACGCGGCATCAATA‑CGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGGACGTCTCCTGGCTGGTGCCCTGCAATACGCATCCCCg  >  1:173892/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
TTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATA‑CGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGAACGTCTCCTGGGTGGTGCCCTGCAATACGCCTACCCG  >  CP009273/2419578‑2419865

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