Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,847,184 C→A 12.9% A351A (GCG→GCT)  ydjE ← putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,847,1840CA12.9% 91.3 / 2.5 31A351A (GCG→GCT) ydjEputative transporter
Reads supporting (aligned to +/- strand):  ref base C (12/15);  new base A (2/2);  total (14/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.52e-01

AACCAGCGCACAGAGCAATAACATAACACCAAGTACCATAAACACCGTGATCGACCCATAATGTGTTAATAATGCCGCAACGCCATAGGGCGTGAAGACTGCGACGATCCGTCCGACGGCATTAACGAAACCCGAACCGCGCAGGCGTAAATGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGCAAACATACATGTATAAAAAGAAGATCATCACCAGTCCATAGATTAAAATCGCCCACTCTGTAGTCTGAATTGAATAGATATAGCCTAACACGGCAATAATAATAAGTAAGGTGGAGCCAAATAACCGACGAGGAAAATGATCAATAATTAATGCCGCAATAAATATTCCTACCGGAGCGCCAATCATAATAACAGCGGTCATTAATATTGATTTATCGACATCAATGCCGGAGTTAACAAATATGGTCGGTATCCATAC  >  CP009273/1846941‑1847393
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                 
acccAGCGAACAGAGAAATAACATAACACCAAGTACCATAACCACCGTGATCGACCCATAATGTGTTAATAATGCCGCAACGCCATAGGGCGTGAAGACTGCGTCGATCCGTCCGACGGCATTAACGAAACCCGAACCGCGCAGGCGTAAATGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGCAAACATACATGTATAAAAAGAAGATCATCACCAGTCCATAGATTAAAATCGCCCACTCTGTAGTCTGAATTGAATAGata                                                                                                                                                                                     <  1:310715/272‑1 (MQ=255)
                                             cGTGATCGACCCATAATGTGTTAATAATGCCGCAACGCCATAGGGCGTGAAGACTGCGACGATCCGTCCGACGGCATTAACGAAACCCGAACCGCGCAGGCGTAAATGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGCAAACATACATGTATAAAAAGACGATCATCACCAGTCCATAGATTAAAATCg                                                                                                                                                                                                                  >  1:306517/1‑200 (MQ=255)
                                              gTGATCGACCCATAATGTGTTAATAATGCCGCAACGCCATAGGGCGTGAAGACTGCGACGATCCGTCCGACGGCATTAACGAAACCCGAACCGCGCAGGCGTAAATGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGCAAACATACATGTATAAAAAGAAGATCATCACCAGTCCATAGATTAAAATCGCCCACTCTGTAGTCTGAATTGAATAGATATAGCCTAACACGGCAATAATNAAAAGTAAAATATAACCAAATAAcct                                                                                                                                      >  2:332255/1‑274 (MQ=255)
                                                 aTCGACCCATAATGTGTTAATAATGCCGCAACGCCATAGGGCGTGAAGACTGCGACGATCCGTCCGACGGCATTAACGAAACCCGAACCGCGCAGGCGTAAATGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGCAAACATACATGTATAAAAAGAAGATCATCACCAGTCCATAGATTAAAATCGCCCa                                                                                                                                                                                                              <  2:306508/200‑1 (MQ=255)
                                                                gtTAATAATGCCGCAACGCCATAGGGCGTGAAGACTGCGACGATCCGTCCGACGGCATTAACGAAACCCGAACCGCGCAGGCGTAAATGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGCAAACATACATGTATAAAAAGAAGATCATCACCAGTCCATAGATTAAAATCGCCCACTCTGTAGTCTGAATTGAATAGATATAGCCTAACACGGCAATAATAATAAGTAAGGTGGAGCCAAATAACCGACGAGGAAAATGATCaat                                                                                                                    >  1:81336/1‑275 (MQ=255)
                                                                             cAACGCCATAGGGCGTGAAGACTGCGACGATCCGTCCGACGGCATTAACGAAACCCGAACCGCGCAGGCGTAAATGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGCAAACATACATGTATAAAAAGAAGATCATCACCAGTCCATAGATTAAAATCGcc                                                                                                                                                                                                                <  2:340161/170‑1 (MQ=255)
                                                                              aaCCCCATAGGGCGTGAAGACTGCGACGACCCGTCCGACGGCATTAACGAAACCCGAACCGGGCAGGCGTAAATGCGTTGGNCAAAGCGCNGGGATATAAACCGCCGACGCGAAGCAAACATACATGTATAAAAAGAAGATCATCACCAGTCCATAGATTAAAATCGCCCACTCTGTAGTCTGAATTGAATAGATATAGCCTAACACGGCAATAATAATAAGTAAGGTGGAGCCAAATACCCGACGAAGAAAATGATCAATAATTAATGCCGCa                                                                                                       <  1:332266/274‑1 (MQ=255)
                                                                                          cGTGAAGACTGCGACGATCCGTCCGACGGCATTAACGAAACCCGAACCGCGCAGGCGTAAATGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGCAAACATACATGTATAAAAAGAAGATCATCACCAGTCCATAGATTAAAATCGCCCa                                                                                                                                                                                                              >  1:252833/1‑159 (MQ=255)
                                                                                          cGTGAAGACTGCGACGATCCGTCCGACGGCATTAACGAAACCCGAACCGCGCAGGCGTAAATGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGCAAACATACATGTATAAAAAGAAGATCATCACCAGTCCATAGATTAAAATCGCCCACTCTGTAGTCtgaat                                                                                                                                                                                               >  2:490806/1‑174 (MQ=255)
                                                                                            tGAAGACTGCGACGATCCGTCCGACGGCATTAACGAAACCCGAACCGCGCAGGCGTAAATGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGCAAACATACATGTATAAAAAGAAGATCATCACCAGTCCATAGATTAAAATCGCCCACTCTGTAGTCtgaat                                                                                                                                                                                               >  2:284044/1‑172 (MQ=255)
                                                                                              aaGACTGCGACGATCCGTCCGACGGCATTAACGAAACCCGAACCGCGCAGGCGTAAATGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGCAAACATACATGTATAAAAAGAAGATCATCACCAGTCCATAGATTAAAATCGCCCActct                                                                                                                                                                                                          <  2:252826/159‑1 (MQ=255)
                                                                                              aaGACTGCGACGATCCGTCCGACGGCATTAACGAAACCCGAACCGCGCAGGCGTAAATGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGCAAACATACATGTATAAAAAGAAGATCATCACCAGTCCATAGATTAAAATCGCCCACTCTGTAGTCtgaattgaa                                                                                                                                                                                           <  1:490819/174‑1 (MQ=255)
                                                                                                gACTGCGACGATCCGTCCGACGGCATTAACGAAACCCGAACCGCGCAGGCGTAAATGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGCAAACATACATGTATAAAAAGAAGATCATCACCAGTCCATAGATTAAAATCGCCCACTCTGTAGTCtgaattgaa                                                                                                                                                                                           <  1:284052/172‑1 (MQ=255)
                                                                                                        ataTCCGTCAGACGGAATTAACGAACCACGAAACGCGTAGGCGTAAATGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGCAAACATACATGTATAAAAAGTAGATCATACCCATTCCATAGATTAACATCGCCCACTCTGTAGTCTGAATTGAATAGATATAGCCTAACACGGCAATAATAATAAGTAAGGTGGAGCCAAATAACCGACGAGGAAAATGATCAATAATTAATGCCGCAATAAATATTCCTACCGGAGCGCCAAt                                                                             <  1:107199/272‑1 (MQ=255)
                                                                                                                                           cgcAGGCGTAAATGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGCAAACATACATGTATAAAAAGAAGATCATCACCAGTCCATAGATTAAAATCGCCCACTCTGTAGTCTGAATTGAATAGatat                                                                                                                                                                                    >  2:95689/1‑136 (MQ=255)
                                                                                                                                               ggCGTAAATGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGCAAACATACATGTATAAAAAGAAGATCATCACCAGTCCATAGATTAAAATCGCCCACTCTGTAGTCTGAATTGAATAGATATAGcc                                                                                                                                                                                <  1:95691/136‑1 (MQ=255)
                                                                                                                                                     aaTGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGCAAACATACATGTATAAAAAGAAGATCATCACCAGTCCATAGATTAAAATCGCCCACTCTGTAGTCTGAATTGAATAGATATAGCCTAACACgg                                                                                                                                                                        >  2:308817/1‑138 (MQ=255)
                                                                                                                                                         cGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGCAAACATACATGTATAAAAAGAAGATCATCACCAGTCCATAGATTAAAATCGCCCACTCTGTAGTCTGAATTGAATAGATATAGCCTAACACGGCaat                                                                                                                                                                    <  1:308826/138‑1 (MQ=255)
                                                                                                                                                                cAAAGCTCCGGGATATAAACCGCCGACGCGAAGCAAACATACATGTATAAAAAGAAGATCATCACCAGTCCATAGATTAAAATCGc                                                                                                                                                                                                                 <  1:66143/86‑1 (MQ=255)
                                                                                                                                                                               taAACCTCCGACGCGAAGCAACCATACATGTATAAAACGAATATCATCACCAGTCCATAGATTAAAATCGCCCACTCTGTAGTCTGAATTGAATAGATATACCCTAACACg                                                                                                                                                                         >  1:114276/1‑111 (MQ=255)
                                                                                                                                                                                   ccgccgAAGAGAAGCAAACATACATGTATAAAAAGAAGATCATCACCAGTCCATAGATTAAAATCGCCCACTCTGTAGGCGGAATTGAATAGATATAGCCTAACAAGGCaa                                                                                                                                                                     <  2:114270/111‑1 (MQ=255)
                                                                                                                                                                                                             tataAAAAGAAGATCATCACCAGTCCATAGATTAAAATAGCCCACTCTg                                                                                                                                                                                                         >  1:106197/1‑49 (MQ=255)
                                                                                                                                                                                                               taAAAAGAAGATCATCACCAGTCCATAGATTAAAATCGCCCACTCTGTAGTCt                                                                                                                                                                                                   >  1:201419/1‑53 (MQ=255)
                                                                                                                                                                                                                 aaaaGAAGATCATCACCAGTCCATAGATTAAAATAGCCCACTCTGTAGt                                                                                                                                                                                                     <  2:106193/49‑1 (MQ=255)
                                                                                                                                                                                                                  aaaGAAGATCATCACCCGTCCATAGATTAAAATCGCCCACTCTGTAGTCTNAATTGAATNGATATAGCCTAACACGGCAATAAGAATAAGTAAGGTGGAGCCAAATAACCGACGAGGAAAATGATCAATAATTAATGCCGCAATAAATATTCCTACCGGAGCGCCAATCATAATAACAGCGGTCATTAATATTGATTTATCGACATCACTGCCGGAGTTAACAAATATGGTCGGTATCCATAc  <  1:331978/243‑1 (MQ=255)
                                                                                                                                                                                                                   aagaagATCATCACCAGTCCATAGATTAAAATCGCCCACTCTGTAGTCtgaat                                                                                                                                                                                               <  2:201413/53‑1 (MQ=255)
                                                                                                                                                                                                                         atcatcACCAGTCCATAGATTAAAATCGCCCACTCTGTAGTCtgaat                                                                                                                                                                                               >  2:368121/1‑47 (MQ=255)
                                                                                                                                                                                                                         atcatcACCAGTCCATAGATTAAAATAGCCCACTCTGTAGTCtgaat                                                                                                                                                                                               >  2:224813/1‑47 (MQ=255)
                                                                                                                                                                                                                             tcaCCAGTCCATAGATTAAAATAGCCCACTCTGTAGTCtgaattgaa                                                                                                                                                                                           <  1:224820/47‑1 (MQ=255)
                                                                                                                                                                                                                                              aaaaTCGCCCACTCTGTAGTCTGAATTGAATAGATATAGCCTAACACGGCAATAATAATAAGTAAGGTGGAGCCAAATAACCGACGAGGAAAATGATCAATAATTAATGCCGCAATAAATATTCCTACCGGAGCGCCAATCATAATAAc                                                                    >  1:294808/1‑149 (MQ=255)
                                                                                                                                                                                                                                                  tCGCCCACTCTGTAGTCTGAATTGAATAGATATAGCCTAACACGGCAATAATAATAAGTAAGGTGGAGCCAAATAACCGACGAGGAAAATGATCAATAATTAATGCCGCAATAAATATTCCTACCGGAGCGCCAATCATAATAACAGCg                                                                <  2:294800/149‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                 
AACCAGCGCACAGAGCAATAACATAACACCAAGTACCATAAACACCGTGATCGACCCATAATGTGTTAATAATGCCGCAACGCCATAGGGCGTGAAGACTGCGACGATCCGTCCGACGGCATTAACGAAACCCGAACCGCGCAGGCGTAAATGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGCAAACATACATGTATAAAAAGAAGATCATCACCAGTCCATAGATTAAAATCGCCCACTCTGTAGTCTGAATTGAATAGATATAGCCTAACACGGCAATAATAATAAGTAAGGTGGAGCCAAATAACCGACGAGGAAAATGATCAATAATTAATGCCGCAATAAATATTCCTACCGGAGCGCCAATCATAATAACAGCGGTCATTAATATTGATTTATCGACATCAATGCCGGAGTTAACAAATATGGTCGGTATCCATAC  >  CP009273/1846941‑1847393

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑