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

TTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGAACGTCTCCTGGGTGGTGCCCTGCAATACGCCTACCCG  >  CP009273/2419578‑2419865
                                                                                                                                              |                                                                                                                                                  
ttCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTCGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGaa                                                                                                                                             <  2:238461/150‑1 (MQ=255)
             ttATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCGTCAATa                                                                                                                                >  2:49603/1‑150 (MQ=255)
                              gCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAACTTGTTATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCACCCTGGAACGCGGCATCAATACGTCCGGCAgtcaggtc                                                                                                              <  2:84020/150‑1 (MQ=255)
                                ccATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCACCCTGGAACGCGGCATCAATACGTCCGGCAgtcaggtcag                                                                                                            <  2:256456/150‑1 (MQ=255)
                                                                  tttAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCCCCGCGATGAGCGACCTCACCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGAt                                                                          <  2:261583/150‑1 (MQ=255)
                                                                               ccACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCACCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCtttt                                                             >  1:195282/1‑150 (MQ=255)
                                                                                                   aCCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTTACCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTAcc                                         <  2:206414/150‑1 (MQ=255)
                                                                                                       aCAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGAAc                                     <  2:128145/150‑1 (MQ=255)
                                                                                                        cAGGTTGTTTGAGGAATCCTTCGCTGGCAGCGACCTCACCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGTCTTTTGGTGCCCAATGTTCATTACCGAACg                                    <  2:214679/150‑1 (MQ=255)
                                                                                                          ggTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCCCCCTGGAACGCGGCATCAATACGTCCGCCAGTCAGGTCAGAACAAATGTTGTCCTGCCCCTGATACGCGACGATTTCAATGCCTGTTGAGGCCCAATGTCCATTACCGAACGtc                                  >  1:23919/1‑150 (MQ=255)
                                                                                                          ggTTGTTTGAGGAAACCTCCGATGCCAGCGACCTCACCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCACAGCTTTTTGGTGCCCAATGTCCACTACCGaagcat                                  >  1:23934/1‑146 (MQ=255)
                                                                                                           gttgttTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGAACGtct                                 >  1:207420/1‑150 (MQ=255)
                                                                                                                           ttCGCTGGCAGCGACCTCACCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGAACGTCTCCTGGGTGGTGCCCTg                 <  2:283109/150‑1 (MQ=255)
                                                                                                                            tCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCTTTACCGAACGTCTCCTGGGTGGTGCCCTGc                <  2:252724/150‑1 (MQ=255)
                                                                                                                              gCTGGCAGCGACCTCACCCTGGAACGCAGATAATGATTTTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGGTCATTACCGACCGTCTCCTGGGTGGTGCCCGGCaa              >  1:346072/1‑150 (MQ=255)
                                                                                                                                  gCAGCGGCCTCACCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCACTGCCTTTTGGTGCCCAATGTTCATTACCGAACGTCTCCTGGGTGGTGCCCTGCAATACg          <  2:105175/150‑1 (MQ=255)
                                                                                                                                          cTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGGACGTCTCCTGGCTGGTGCCCTGCAATACGCATCCCCg  >  1:173892/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
TTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGAACGTCTCCTGGGTGGTGCCCTGCAATACGCCTACCCG  >  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: