Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,172,868 G→C G190A (GGC→GCC)  yegT → nucleoside permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,172,8680GC96.0% 77.7 / ‑4.5 25G190A (GGC→GCC) yegTnucleoside permease
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base C (17/7);  total (18/7)
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

TCCCGCTGCTGATTACCGCCGGAAGTTCTGCTCTGCTCGGTGTGTTTGCGTTTTTCCTGCCCGACACGCCACCAAAAAGCACCGGCAAAATGGATATTAAAGTCATGCTCGGCCTGGATGCGCTGATCCTGCTGCGCGATAAAAACTTCCTCGTCTTTTTCTTCTGTT  >  NZ_CP009273/2172784‑2172951
                                                                                    |                                                                                   
tCCCGCTGCTGATTACCGCCGGAAGTTCTGCTCTGCTCGGTGTGTTTGCGTTTTTCCTGCCCGACACGCCACCAAAAAGCACCGCCaaaa                                                                                <  2:661702/90‑1 (MQ=255)
  ccGCTGCTGATTACCGCCGGAAGTTCTGCTCTGCTCGGTGTGTTTGCGTTTTTCCTGCCCGACACGCCACCAAAAAGCACCGCCAAAATg                                                                              >  1:695394/1‑90 (MQ=255)
               ccgccgGAAGTTCTGCTCTGCTCGGTGTGTTTGCGTTTTTCCTGCCCGACACGCCACCAAAAAGCACCGCCAAAATGGATATTAAAGTCa                                                                 <  1:115294/90‑1 (MQ=255)
                 gccgGAAGTTCTGCTCTGCTCGGTGTGTTTGCGTTTTTCCTGCCCGACACGCCACCAAAAAGCACCGCCAAAATGGATATTAAAGTCATg                                                               >  2:266764/1‑90 (MQ=255)
                   cgGAAGTTCTGCTCTGCTCGGTGTGTTTGCGTTTTTCCTGCCCGACACGCCACCAAAAAGCACCGCCAAAATGGATATTAAAGTCATGCt                                                             <  2:347962/90‑1 (MQ=255)
                    ggAAGTTCTGCTCTGCTCGGTGTGTTTGCGTTTTTCCTGCCCGACACGCCACCAAAAAGCACCGCCAAAATGGATATTAAAGTCATGCTc                                                            <  1:179801/90‑1 (MQ=255)
                    ggAAGTTCTGCTCTGCTCGGTGTGTTTGCGTTTTTCCTGCCCGACACGCCACCAAAAAGCACCGCCAAAATGGATATTAAAGTCATGCTc                                                            <  2:202155/90‑1 (MQ=255)
                                  gctcGGTGTGTTTGCGTTTTTCCTGCCCGACACGCCACCAAAAAGCACCGGCAAAATGGATATTAAAGTCATGCTCGGCCTGGATGCGCt                                              >  1:232091/1‑90 (MQ=255)
                                  gctcGGTGTGTTTGCGTTTTTCCTGCCCGACACGCCACCAAAAAGCACCGCCAAAATGGATATTAAAGTCATGCTCGGCCTGGATGCGCt                                              >  1:53492/1‑90 (MQ=255)
                                  gctcGGTGTGTTTGCGTTTTTCCTGCCCGACACGCCACCAAAAAGCACCGCCAAAATGGATATTAAAGTCATGCTCGGCCTGGATGCGCt                                              >  1:58160/1‑90 (MQ=255)
                                  gctcGGTGTGTTTGCGTTTTTCCTGCCCGACACGCCACCAAAAAGCACCGCCAAAATGGATATTAAAGTCATGCTCGGCCTGGATGCGCt                                              >  2:19699/1‑90 (MQ=255)
                                  gctcGGTGTGTTTGCGTTTTTCCTGCCCGACACGCCACCAAAAAGCACCGCCAAAATGGATATTAAAGTCATGCTCGGCCTGGATGCGCt                                              >  1:421760/1‑90 (MQ=255)
                                       gtgtgtTTGCGTTTTTCCTGCCCGACACGCCACCAAAAAGCACCGCCAAAATGGATATTAAAGTCATGCTCGGCCTGGATGCGCTGATcc                                         >  1:389775/1‑90 (MQ=255)
                                       gtgtgtTTGCGTTTTTCCTGCCCGACACGCCACCAAAAAGCACCGCCAAAATGGATATTAAAGTCATGCTCGGCCTGGATGCGCTGATcc                                         >  2:350307/1‑90 (MQ=255)
                                            tttGCGTTTTTCCTGCCCGACACGCCACCAAAAAGCACCGCCAAAATGGATATTAAAGTCATGCTCGGCCTGGATGCGCTGATCctgctg                                    >  1:429390/1‑90 (MQ=255)
                                                 gTTTTTCCTGCCCGACACGCCACCAAAAAGCACCGCCAAAATGGATATTAAAGTCATGCTCGGCCTGGATGCGCTGATCCTGCTGCGCGa                               >  1:64240/1‑90 (MQ=255)
                                                 gTTTTTCCTGCCCGACACGCCACCAAAAAGCACCGCCAAAATGGATATTAAAGTCATGCTCGGCCTGGATGCGCTGATCCTGCTGCGCGa                               >  1:448955/1‑90 (MQ=255)
                                                 gTTTTTCCTGCCCGACACGCCACCAAAAAGCACCGCCAAAATGGATATTAAAGTCATGCTCGGCCTGGATGCGCTGATCCTGCTGCGCGa                               <  1:700958/90‑1 (MQ=255)
                                                 gTTTTTCCTGCCCGACACGCCACCAAAAAGCACCGCCAAAATGGATATTAAAGTCATGCTCGGCCTGGATGCGCTGATCCTGCTGCGCGa                               >  2:604625/1‑90 (MQ=255)
                                                          gCCCGACACGCCACCAAAAAGCACCGCCAAAATGGATATTAAAGTCATGCTCGGCCTGGATGCGCTGATCCTGCTGCGCGATAAAAACtt                      >  1:556168/1‑90 (MQ=255)
                                                                         aaaaaGCACCGCCAAAATGGATATTAAAGTCATGCTCGGCCTGGATGCGCTGATCCTGCTGCGCGATAAAAACTTCCTCGTCTTTttctt       >  2:659459/1‑90 (MQ=255)
                                                                            aaGCACCGCCAAAATGGATATTAAAGTCATGCTCGGCCTGGATGCGCTGATCCTGCTGCGCGATAAAAACTTCCTCGTCTTTTTCTTCTg    >  2:485056/1‑90 (MQ=255)
                                                                              gCACCGCCAAAATGGATATTAAAGTCATGCTCGGCCTGGATGCGCTGATCCTGCTGCGCGATAAAAACTTCCTCGTCTTTTTCTTCTGtt  >  1:595081/1‑90 (MQ=255)
                                                                               cACCGCCAAAATGGATATTAAAGTCATGCTCGGCCt                                                       >  1:707561/1‑36 (MQ=38)
                                                                               cACCGCCAAAATGGATATTAAAGTCATGCTCGGCCt                                                       <  2:707561/36‑1 (MQ=38)
                                                                                    |                                                                                   
TCCCGCTGCTGATTACCGCCGGAAGTTCTGCTCTGCTCGGTGTGTTTGCGTTTTTCCTGCCCGACACGCCACCAAAAAGCACCGGCAAAATGGATATTAAAGTCATGCTCGGCCTGGATGCGCTGATCCTGCTGCGCGATAAAAACTTCCTCGTCTTTTTCTTCTGTT  >  NZ_CP009273/2172784‑2172951

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