Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,161,974 T→A 7.3% A68A (GCT→GCA)  yijD → DUF1422 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,161,9740TA 7.3% 185.0 / 5.2 55A68A (GCT→GCA) yijDDUF1422 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (25/26);  new base A (2/2);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01

TTAATGGTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAA  >  NC_000913/4161832‑4162105
                                                                                                                                              |                                                                                                                                   
ttAATGGTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGt                                                                                                                               >  2:479563/1‑149 (MQ=255)
      gTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTAtt                                                                                                                         <  1:406783/149‑1 (MQ=255)
      gTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTAtt                                                                                                                         <  3:44202/149‑1 (MQ=255)
          tttCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCg                                                                                                                     <  1:255581/149‑1 (MQ=255)
               cAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTa                                                                                                                <  1:394446/149‑1 (MQ=255)
                    cTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTActgct                                                                                                           >  3:175285/1‑149 (MQ=255)
                    cTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTActgct                                                                                                           <  1:35023/149‑1 (MQ=255)
                         tAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACa                                                                                                      <  4:496482/149‑1 (MQ=255)
                             tCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCAc                                                                                                  <  1:410175/149‑1 (MQ=255)
                                aTTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGt                                                                                               >  3:303058/1‑149 (MQ=255)
                                       cATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGt                                                                                        >  3:167570/1‑149 (MQ=255)
                                                  ttCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATcc                                                                             <  1:173960/149‑1 (MQ=255)
                                                     ccGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGt                                                                                               >  4:126797/1‑128 (MQ=255)
                                                     ccGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGt                                                                                               <  3:126797/128‑1 (MQ=255)
                                                     ccGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGa                                                                          >  1:460162/1‑149 (MQ=255)
                                                     ccGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGa                                                                          <  3:39018/149‑1 (MQ=255)
                                                       gATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCAGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGata                                                                        >  3:245483/1‑149 (MQ=255)
                                                       gATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCAGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGata                                                                        >  3:242998/1‑149 (MQ=255)
                                                              tCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCAGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGctc                                                                 <  4:242998/149‑1 (MQ=255)
                                                              tCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCAGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGctc                                                                 <  4:245483/149‑1 (MQ=255)
                                                               cccTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTTTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGctct                                                                >  4:104042/1‑149 (MQ=255)
                                                                  tGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAAc                                                             <  2:9511/149‑1 (MQ=255)
                                                                  tGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAAc                                                             <  4:78331/149‑1 (MQ=255)
                                                                       cTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAActtctt                                                        >  4:244024/1‑149 (MQ=255)
                                                                         gACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTcc                                                      >  2:518587/1‑149 (MQ=255)
                                                                            ggTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCgg                                                   >  1:506544/1‑149 (MQ=255)
                                                                            ggTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCgg                                                   >  3:128327/1‑149 (MQ=255)
                                                                                  ctgcctgcATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTAc                                             <  2:516561/149‑1 (MQ=255)
                                                                                     cctgcATCAACGTTATCTTAATCGTACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTActct                                          >  1:217975/1‑149 (MQ=255)
                                                                                      ctgcATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACtctc                                         >  1:427558/1‑149 (MQ=255)
                                                                                          aTCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTc                                     >  1:261246/1‑149 (MQ=255)
                                                                                          aTCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTc                                     >  3:335980/1‑149 (MQ=255)
                                                                                          aTCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTc                                     >  1:482328/1‑149 (MQ=255)
                                                                                          aTCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTc                                     >  1:311657/1‑149 (MQ=255)
                                                                                            cAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCAt                                   <  2:390330/149‑1 (MQ=255)
                                                                                             aaCGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCAtt                                  >  1:259412/1‑149 (MQ=255)
                                                                                             aaCGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCAtt                                  >  3:202228/1‑149 (MQ=255)
                                                                                                gTTATCTTAATCGCACTATGCCGGTAGGCTTGCTGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATg                               <  4:375472/149‑1 (MQ=255)
                                                                                                gTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATg                               <  4:46560/149‑1 (MQ=255)
                                                                                                gTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATg                               <  4:167570/149‑1 (MQ=255)
                                                                                                gTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATg                               >  4:130129/1‑149 (MQ=255)
                                                                                                gTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATg                               <  2:427558/149‑1 (MQ=255)
                                                                                                   aTCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGat                                                                         >  3:398323/1‑104 (MQ=255)
                                                                                                   aTCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGat                                                                         <  4:398323/104‑1 (MQ=255)
                                                                                                      ttAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTc                         >  2:463563/1‑149 (MQ=255)
                                                                                                          tCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGa                     >  3:172053/1‑149 (MQ=255)
                                                                                                             cACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTg                  <  1:463563/149‑1 (MQ=255)
                                                                                                                    ccGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGaa           <  4:124270/149‑1 (MQ=255)
                                                                                                                      ggTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGAAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGa         <  2:218837/149‑1 (MQ=255)
                                                                                                                       gTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGAt        <  1:518587/149‑1 (MQ=255)
                                                                                                                       gTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGAt        <  2:217975/149‑1 (MQ=255)
                                                                                                                         aGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGc      <  3:244024/149‑1 (MQ=255)
                                                                                                                           gCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGt    >  1:145057/1‑149 (MQ=255)
                                                                                                                            cTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTa   >  1:575951/1‑149 (MQ=255)
                                                                                                                             ttGACGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGTGTaa  <  1:522539/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                   
TTAATGGTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTAA  >  NC_000913/4161832‑4162105

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