Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,911,371 C→A P366T (CCG→ACG)  yebQ → putative transporter YebQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,911,3710CA96.5% 89.5 / ‑3.9 29P366T (CCG→ACG) yebQputative transporter YebQ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (18/10);  minor base G (1/0);  total (19/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

GGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACCACACCATTATTACCTCCGCGCCTCG  >  NC_000913/1911293‑1911447
                                                                              |                                                                            
gggTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGACg                                                                            <  2:326182/81‑1 (MQ=255)
     gttCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGACgatgat                                                                       >  1:326203/1‑81 (MQ=255)
     gttCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGACgatgat                                                                       >  1:326210/1‑81 (MQ=255)
         atcatGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGACGATGATCTTa                                                                   >  2:326232/1‑81 (MQ=255)
          tcatGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGACGATGATCTTAt                                                                  >  2:326221/1‑81 (MQ=255)
          tcatGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGACGATGATCTTAt                                                                  >  2:326208/1‑81 (MQ=255)
                     gggCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGACGATGATCTTATGTGGTGCTGGa                                                       >  1:326216/1‑81 (MQ=255)
                           ttttCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGACGATGATCTTATGTGGTGCTGGATTTgg                                                  <  1:326184/80‑1 (MQ=255)
                           ttttCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGACGATGATCTTATGTGGTGCTGGATTTg                                                   <  2:326183/79‑1 (MQ=255)
                           ttttCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGACGATGATCTTATGTGGTGCTGGATTTGGc                                                 >  1:326233/1‑81 (MQ=255)
                                ccTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGACGATGATCTTATGTGGTGCTGGATTTGGCTTAtt                                            >  2:326228/1‑81 (MQ=255)
                                 cTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGACGATGATCTTATGTGGTGCTGGATTTGGCTTATTc                                           <  2:326185/81‑1 (MQ=255)
                                 cTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGACGATGATCTTATGTGGTGCTGGATTTGGCTTATTc                                           >  2:326212/1‑81 (MQ=255)
                                   ggTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGACGATGATCTTATGTGGTGCTGGATTTGGCTTATTc                                           >  1:326227/1‑79 (MQ=255)
                                    gTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGACGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAg                                        >  1:326217/1‑81 (MQ=255)
                                       ctgctgCCCGCGTCACCTGCGGATATCAATATTATCTGGACGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCa                                     >  2:326213/1‑81 (MQ=255)
                                         gctgcCCGCGTCACCTGCGGATATCAATATTATCTGGGCGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCAcc                                   >  2:326229/1‑81 (MQ=255)
                                          ctgcCCGCGTCACCTGCGGATATCAATATTATCTGGACGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCAccc                                  <  2:326186/81‑1 (MQ=255)
                                             cccGCGTCACCTGCGGATATCAATATTATCTGGACGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCTCCCAAt                               <  2:326187/81‑1 (MQ=255)
                                                gcgTCACCTGCGGATATCAATATTATCTGGACGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAAc                            >  1:598753/1‑81 (MQ=255)
                                                  gTCACCTGCGGATATCAATATTATCTGGACGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACca                          >  1:326225/1‑81 (MQ=255)
                                                     aCCTGCGGATATCAATATTATCTGGACGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACCacac                       >  2:326223/1‑81 (MQ=255)
                                                          cGGATATCAATATTATCTGGACGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACCACACCat                    <  2:326188/79‑1 (MQ=255)
                                                              tatCAATATTATCTGGACGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACCACACCATTATTAc              >  1:326231/1‑81 (MQ=255)
                                                                 cAATATTATCTGGACGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACCACACCATTATTACCTc           >  1:326226/1‑81 (MQ=255)
                                                                  aaTATTATCTGGACGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACCACACCATTATTACCTcc          <  2:326189/81‑1 (MQ=255)
                                                                    tattatCTGGACGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCGGTCACCCAATAACCACACCATTATTACCTCcgc        <  1:326190/81‑1 (MQ=255)
                                                                        atCTGGACGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACCACACCATTATTACCTCCGCGcc     <  2:326191/80‑1 (MQ=255)
                                                                          cTGGACGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACCACACCATTATTACCTCCGCGCCTcg  >  2:326218/1‑81 (MQ=255)
                                                                              |                                                                            
GGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACCACACCATTATTACCTCCGCGCCTCG  >  NC_000913/1911293‑1911447

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