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

AGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGAACGTCTCCTGGGTGGTGCCCTGCAATACGCCTA  >  CP009273/2419573‑2419861
                                                                                                                                                   |                                                                                                                                             
aGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATc                                                                                                                                              >  1:50836/1‑149 (MQ=255)
            aGTTCGTTATCTTCTTTACGCAGGCCCTGGCCGGCCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCACCCTGGaa                                                                                                                                        <  2:197548/143‑1 (MQ=255)
            aGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGc                                                                                                                                  <  1:456043/149‑1 (MQ=255)
                 gTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCACCCTGGAACGCGGCATCaa                                                                                                                             <  1:522392/149‑1 (MQ=255)
                    aTCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCCGCGACCTACCCCTGGAACGCGGcct                                                                                                                                >  2:402208/1‑141 (MQ=255)
                      cttcttTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGt                                                                                                                        >  1:461848/1‑149 (MQ=255)
                       ttcttTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCaa                                                                                                                             <  2:205520/143‑1 (MQ=255)
                       ttcttTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCACCCTGGAACGCGGCATCAATACGTc                                                                                                                       >  1:62080/1‑149 (MQ=255)
                          tttACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCACCCTGGAACGCGGCATCAATACGTCCgg                                                                                                                    <  1:331188/149‑1 (MQ=255)
                                   gCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCACCCTGGAACGCGGCATCAATACGTCCGGCAgtcaggt                                                                                                           <  1:257600/149‑1 (MQ=255)
                                    cccATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAgt                                                                                                                <  2:45048/143‑1 (MQ=255)
                                         gCCGGTCCCTACGCCAAACAGTTCTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGATTCCTTCGCGGGCAGCGACCTCACCCGGGAACGCGGCATCCATACGTCCGGCAGTCAGGTCAGAAt                                                                                                     <  1:125912/149‑1 (MQ=255)
                                                ccTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCACCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAAtgtt                                                                                              <  1:452830/149‑1 (MQ=255)
                                                  tACGCCAAACAGTCTTTCATCTTTAACAGACGGGCCACCGAATTTGTACTCTTTACCGACAGGTTGTTTGAGGAAGCCTTCGCTGGCAGCGACCTAACCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATtaa                                                                                                  <  2:360876/143‑1 (MQ=255)
                                                   aCGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCACCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAAt                                                                                                 <  2:477103/143‑1 (MQ=255)
                                                         aaCAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCACCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAAGGTTGTc                                                                                           >  2:62100/1‑143 (MQ=255)
                                                          aCAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCACCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGcccc                                                                                    <  1:447327/149‑1 (MQ=255)
                                                                    aTCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCACTGATACgaga                                                                          >  1:273775/1‑149 (MQ=255)
                                                                       tttAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGa                                                                       <  1:134720/149‑1 (MQ=255)
                                                                              gACGGCCCCACGAATTTGTAATCTTACCCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAAt                                                                <  1:274782/149‑1 (MQ=255)
                                                                               aCGGTCCACCGAATCTGTAATCTTGCCCGCTACGCTTCTTCAGAAAATCCTCGCTGGCCGCGACCTCCCCCTGGAACGCTGCATCAATACGGCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATg                                                               <  1:332138/149‑1 (MQ=255)
                                                                               aCGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGAGGAGTTCAATg                                                               >  1:477851/1‑149 (MQ=255)
                                                                               aCGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCACCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATg                                                               <  1:232719/149‑1 (MQ=255)
                                                                                   gCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCaa                                                                 >  2:29659/1‑143 (MQ=255)
                                                                                            tttttcaTCTTCCCCGACAGCTTGTTTGAGGAAACCTCCGCTGGCAGCGACCGCATCCTGGAACGCGGCATCAATACGCCCGGCAGTCAGTTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTg                                                        <  2:86396/137‑1 (MQ=255)
                                                                                                tAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCACCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATg                                              >  1:306595/1‑149 (MQ=255)
                                                                                                          cGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCa                                          >  2:278176/1‑143 (MQ=255)
                                                                                                            aCAGGTTTGTTGAGGAAGCCTTCTCGGTCAGCGACCGCATCCTGGAAGGCGTCATCAATACGTCCGGCAGTCAGGGCAGAATAAATGTTGGCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGaa                                  <  1:435183/149‑1 (MQ=255)
                                                                                                            aCAGGTTGTTTGAGGAAACCTTAGCTGTCAGCGACCTAATCTTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCACTGTTCATTACCGaa                                  >  1:106046/1‑149 (MQ=255)
                                                                                                                         ggAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGAACGTCTCCTGggtg                     <  1:328121/149‑1 (MQ=255)
                                                                                                                         ggAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGAACGTCTCCTGggtg                     <  1:328139/149‑1 (MQ=255)
                                                                                                                            aaCCTTCGCTGGCAGCGACCTCACCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGAACGTCTCCTGGgtggtg                  <  1:304433/149‑1 (MQ=255)
                                                                                                                                ttCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGAACGTCTCCTGGGTGGTGCCCt              >  1:444931/1‑149 (MQ=255)
                                                                                                                                  cGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGAACGTCTCCTGGgtggtg                  <  2:336916/143‑1 (MQ=255)
                                                                                                                                      ggCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGAACGTCTCCTGGGTGGTGCCCTGCAATa        <  1:224141/149‑1 (MQ=255)
                                                                                                                                         aGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGAACGTCTCCTGGGTGGTGCCCTGCa           >  2:383009/1‑143 (MQ=255)
                                                                                                                                                  aTCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCGGAGTAAGTGTTGTGCTGCCCCTGGTACGAGAGGATTTCAATGCCTTTTGGTGCCCAATGTTGATTACCGACCGTCTCCTGGGTGGTGGCCTGCCATCCGCCtg  >  2:133231/1‑142 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
AGTGCTTCGCGCAGTTCGTTATCTTCTTTACGCAGGCCCATGCCGGTCCCTACGCCAAACAGTTTTTCATCTTTAACAGACGGGCCACCGAATTTGTAATCTTTACCGACAGGTTGTTTGAGGAAACCTTCGCTGGCAGCGACCTCATCCTGGAACGCGGCATCAATACGTCCGGCAGTCAGGTCAGAATAAATGTTGTCCTGCCCCTGATACGAGACGATTTCAATGCCTTTTGGTGCCCAATGTTCATTACCGAACGTCTCCTGGGTGGTGCCCTGCAATACGCCTA  >  CP009273/2419573‑2419861

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: ‑