Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 508,601 (C)6→7 coding (570/1293 nt) ybaT → APC family permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273508,5951.C87.5% 42.6 / 0.5 16coding (564/1293 nt)ybaTAPC family permease
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (8/6);  total (9/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

CATTAAAATGATGATCCTGTTATTGCTGATTATTGCCGGTGTCTGGTCGCTGCAACCGGCGCATATTTCCGTCTCTGCGCCCCCCAGCTCCGGTGCGTTCTTCTCCTGTATTGGGATAACTTTCCTTGCCTATGCGGGCTTTGGCATGATGGCGAAC  >  NZ_CP009273/508517‑508673
                                                                               |                                                                              
cATTAAAATGATGATCCTGTTTTTGCTGATTATTGCCGGTGTCTGGTCGCTGCAACCGGCGCATATTTCCGTCTCTGCGCCCCCCAGCTc                                                                     <  2:262997/90‑1 (MQ=255)
            gatCCTGTTATTGCTGATTATTGCCGGTGTCTGGTCGCTGCAACCGGCGCATATTTCCGTCTCTGCGCCCCCCCAGCTCCGGTGCGttct                                                          <  1:274942/90‑1 (MQ=255)
                 tGTTATTGCTGATTATTGCCGGTGTCTGGTCGCTGCAACCGGCGCATATTTCCGTCTCTGCGCCCCCCCAGCTCCGGTGCGTTCTTCTcc                                                     >  1:250731/1‑90 (MQ=255)
                       tGCTGATTATTGCCGGTGTCTGGTCGCTGCAACCGGCGCATATTTCCGTCTCTGCGCCCCCCCAGCTCCGGTGCGTTCTTCTCCTGTAtt                                               >  2:191275/1‑90 (MQ=255)
                           gATTATTGCCGGTGTCTGGTCGCTGCAACCGGCGCATATTTCCGTCTCTGCGCCCCCCCAGCTCCGGTGCGTTCTTCTCCTGTATTGGGa                                           >  1:284216/1‑90 (MQ=255)
                           gATTATTGCCGGTGTCTGGTCGCTGCAACCGGCGCATATTTCCGTCTCTGCGCCCCCCCAGCTCCGGTGCGTTCTTCTCCTGTATTGGGa                                           >  2:73191/1‑90 (MQ=255)
                                     ggTGTCTGGTCGCTGCAACCGGCGCATATTTCCGTCTCTGCGCCCCCCCAGCTCCGGTGCGTTCTTCTCCTGTATTGGGATAACTTTCCt                                 <  1:274360/90‑1 (MQ=255)
                                     ggTGTCTGGTCGCTGCAACCGGCGCATATTTCCGTCTCTGCGCCCCCCAGCTCCGGTGCGTTCTTCTCCTGTATTGGGATAACTTTCCtt                                >  1:27734/1‑90 (MQ=255)
                                       tgtCTGGTCGCTGCAACCGGCGCATATTTCCGTCTCTGCGCCCCCCCAGCTCCGGTGCGTTCTTCTCCTGTATTGGGATAACTTTCCTTg                               >  1:375392/1‑90 (MQ=255)
                                         tCTGGTCGCTGCAACCGGCGCATATTTCCGTCTCTGCGCCCCCCCAGCTCCGGTGCGTTCTTCTCCTGTATTGGGATAACTTTCCTTGcc                             <  1:115895/90‑1 (MQ=255)
                                         tCTGGTCGCTGCAACCGGCGCATATTTCCGTCTCTGCGCCCCCCCAGCTCCGGTGCGTTCTTCTCCTGTATTGGGATAACTTTCCTTGcc                             <  1:254630/90‑1 (MQ=255)
                                               cGCTGCAACCGGCGCATATTTCCGTCTCTGCGCCCCCCCAGCTCCGGTGCGTTCTTCTCCTGTATTGGGATAACTTTCCTTGCCTATGCg                       <  1:191275/90‑1 (MQ=255)
                                                gCTGCAACCGGCGCATATTTCCGTCTCTGCGCCCCCCCAGCTCCGGTGCGTTCTTCTCCt                                                    <  1:94716/60‑1 (MQ=39)
                                                gCTGCAACCGGCGCATATTTCCGTCTCTGCGCCCCCCCAGCTCCGGTGCGTTCTTCTCCt                                                    >  2:94716/1‑60 (MQ=39)
                                                 cTGCAACCGGCGCATATTTCCGTCTCTGCGCCCCCCCAGCTCCGGTGCGTTCTTCTCCTGTATTGGGATAACTTTCCTTGCCTATGCggg                     >  2:259414/1‑90 (MQ=255)
                                                                    ccGTCTCTGCGCCCCCCCAGCTCCGGTGCGTTCTTCTCCTGTATTGGGATAACTTTCCTTGCCTATGCGGGCTTTGGCATGATGGCGAAc  >  2:132437/1‑90 (MQ=255)
                                                                               |                                                                              
CATTAAAATGATGATCCTGTTATTGCTGATTATTGCCGGTGTCTGGTCGCTGCAACCGGCGCATATTTCCGTCTCTGCGCCCCCCAGCTCCGGTGCGTTCTTCTCCTGTATTGGGATAACTTTCCTTGCCTATGCGGGCTTTGGCATGATGGCGAAC  >  NZ_CP009273/508517‑508673

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