Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,389,175 Δ1 bp coding (60/849 nt) yfbP → uncharacterized protein YfbP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,389,1720G.97.9% 194.2 / ‑3.0 48coding (57/849 nt)yfbPuncharacterized protein YfbP
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base . (23/24);  total (24/24)
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

TGGTGGGATTAAAACATTTATATAAACAGATAAGCGGATGTTTACATAAATTTTTATGTGAAGGAACATGATGAAACTTATTCCTCGCAGTAGTGATATTTCACCTGGCATTGATGGAATTTGCCCGGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATGCAGCATATGGTAACGGAGACTGTTTCGGTTTGTACTGGCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACG  >  NC_000913/2389046‑2389292
                                                                                                                              |                                                                                                                        
tggtggGATTAAAACATTTATATAAACAGATAAGCGGATGTTTACATAAATTTTTATGTGAAGGAACATGATGAAACTTATTCCTCGCAGTAGTGATATTTCACCTGGCATTGATGGAATTTGCCCgggcctt                                                                                                                    >  2:458411/1‑129 (MQ=255)
 ggtggGATTAAAACATTTATATAAACAGATAAGCGGATGTTTACATAAATTTTTATGTGAAGGAACATGATGAAACTTATTCCTCGCAGTAGTGATATTTCACCTGGCATTGATGGAATTTGCCCGGGCCttt                                                                                                                  >  2:549961/1‑133 (MQ=255)
        ttAAAACATTTATATAAACAGATAAGCGGATGTTTACATAAATTTTTATGTGAAGGAACATGATGAAACTTATTCCTCGCAGTAGTGATATTTCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCg                                                                                                           <  1:577204/133‑1 (MQ=255)
        ttAAAACATTTATATAAACAGATAAGCGGATGTTTACATAAATTTTTATGTGAAGGAACATGATGAAACTTATTCCTCGCAGTAGTGATATTTCACCTGGCAGTGATGGACTTTACCTGGGCCTGTTCCTCCg                                                                                                           >  1:285712/1‑133 (MQ=255)
         tAAAACATTTATATAAACAGATAAGCGGATGTTTACATAAATTTTTATGTGAAGGAACATGATGAAACTTATTCCTCGCAGTAGTGATATTCCACCTGGCATTGATTGAATTTGCCCGGGCCGATTCCTCCGa                                                                                                          >  2:285135/1‑133 (MQ=255)
                 ttATATAAACAGATAAGCGGATGTTTACATAAATTTTTATGTGAAGGAACATGATGAAACTTATTCCTCGCAGTAGTGATATTTCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGAttt                                                                                                  >  1:95363/1‑133 (MQ=255)
                   atataAACAGATAAGCGGATGTTTACATAAATTTTTATGTGAAGGAACATGATGAAACTTATTCCTCGCAGTAGTGATATTTCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTAc                                                                                                <  1:452940/133‑1 (MQ=255)
                         aCAGATAAGCGGATGTTTACATAAATTTTTATGTGAAGGAACATGATGAAACTTATTCCTCGCAGTAGTGATATTTCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACt                                                                                          >  2:151601/1‑133 (MQ=255)
                             aTAAGCGGATGTTTACATAAATTTTTATGTGAAGGAACATGATGAAACTTATTCCTCGCAGTAGTGATATTTCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACa                                                                                      >  2:590578/1‑133 (MQ=255)
                             aTAAGCGGATGTTTACATAAATTTTTATGTGAAGGAACATGATGAAACTTATTCCTCGCAGTAGTGATATTTCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACa                                                                                      >  1:138790/1‑133 (MQ=255)
                             aTAAGCGGATGTTTACATAAATTTTTATGTGAAGGAACATGATGAAACTTATTCCTCGCAGTAGTGATATTTCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACa                                                                                      >  1:621756/1‑133 (MQ=255)
                             aTAAGCGGATGTTTACATAAATTTTTATGTGAAGGAACATGATGAAACTTATTCCTCGCAGTAGTGATATTTCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACa                                                                                      <  2:42370/133‑1 (MQ=255)
                              tAAGCGGATGTTTACATAAATTTTTATGTGAAGGAACATGATGAAACTTATTCCTCGCAGTAGTGATATTTCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAg                                                                                     >  2:345733/1‑133 (MQ=255)
                              tAAGCGGATGTTTACATAAATTTTTATGTGAAGGAACATGATGAAACTTATTCCTCGCAGTAGTGATATTTCACCTGGCATTGATGGAATTGGCCCGGGCATTTACCTCCGAACGGAGTTACTGTACTTACAg                                                                                     >  2:204637/1‑133 (MQ=255)
                                 gCGGATGTTTACATAAATTTTTATGTGAAGGAACATGATGAAACTTATTCCTCGCAGTAGTGATATTTCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATg                                                                                  <  2:100385/133‑1 (MQ=255)
                                 gCGGATGTTTACATAAATTTTTATGTGAAGGAACATGATGAAACTTATTCCTCGCAGTAGTGATATTTCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATg                                                                                  <  2:539044/133‑1 (MQ=255)
                                  cGGATGTTTACATAAATTTTTATGTGAAGGAACATGATGAAACTTATTCCTCGCAGTAGTGATATTTCACCTGGCATTGATGGAATTTGCCGGGGCCTTTTCCTACAAAGGGATTTAACGTACGTACAGAtac                                                                                 >  2:301133/1‑131 (MQ=255)
                                     aTGTTTACATAAATTTTTATGTGAAGGAACATGATGAAACTTATTCCTCGCAGTAGTGATATTTCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATgcagc                                                                              >  2:44244/1‑133 (MQ=255)
                                      tGTTTACATAGATTTGTATGTGACTGAACATGTTGAAACTTCTTCCTCGCAGTAGTGATATTTCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATgcagca                                                                             <  2:285712/133‑1 (MQ=255)
                                      tGTTTACATAAATTTTTATGTGAAGGAACATGATGAAACTTATTCCTCGCAGTAGTGATATTTCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATgcagca                                                                             <  1:151601/133‑1 (MQ=255)
                                       gTTTACATAAATTTTTATGTGAAGGAACATGATGAAACTTATTCCTCGCAGTAGTGATATTTCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATGCAGCat                                                                            >  1:121137/1‑133 (MQ=255)
                                          tACATAAATTTTTATGTGAAGGAACATGATGAAACTTATTCCTCGCAGTAGTGATATTTCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATGCAGCATATg                                                                         >  2:345264/1‑133 (MQ=255)
                                          tACATAAATTTTTATGTGAAGGAACATGATGAAACTTATTCCTCGCAGTAGTGATATTTCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATGCAGCATATg                                                                         >  2:36431/1‑133 (MQ=255)
                                                       aTGTGAAGGAACATGATGAAACTTATTCCTCGCAGTAGTGATATTTCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATGCAGCATATGGTAATGGAGACTg                                                            >  1:262840/1‑133 (MQ=255)
                                                         gtgAAGGAACATGATGAAACTTATTCCTCGCAGTAGTGATATTTCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATGCAGCATATGGTAACGGAGACTGtt                                                          <  1:495499/133‑1 (MQ=255)
                                                              ggAACATGATGAAACTTATTCCTCGCAGTAGTGATATTTCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATGCAGCATATGGTAACGGAGACTGTTTCGGt                                                     <  2:95363/133‑1 (MQ=255)
                                                              ggAACATGATGAAACTTATTCCTCGCAGTAGTGATATTTCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATGCAGCATATGGTAACGGAGACTGTTTCGGt                                                     <  1:345733/133‑1 (MQ=255)
                                                              ggAACATGATGAAACTTATTCCTCGCAGTAGTGATATTTCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATGCAGCATATGGTAACGGAGACTGTTTCGGt                                                     <  1:36431/133‑1 (MQ=255)
                                                              ggAACATGATGAAACATACTCCACGCAGTAGTGATATTTCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATGCAGCATATGGTAACGGAGACTGTTTCGGt                                                     <  1:301133/133‑1 (MQ=255)
                                                                   atgatgAAACTTATTCCTCGCAGTAGTGATATTTCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATGCAGCATATGGTAACGGAGACTGTTTCGGTTTGTa                                                <  1:590578/133‑1 (MQ=255)
                                                                   atgatgAAACTTATTCCTCGCAGTAGTGATATTTCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATGCAGCATATGGTAACGGAGACTGTTTCGGTTTGTa                                                <  1:549961/133‑1 (MQ=255)
                                                                      atgaAACTTATTCCTCGCAGTAGTGATATTTCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATGCAGCATATGGTAACGGAGACTGTTTCGGTTTGTACTg                                             <  1:488785/133‑1 (MQ=255)
                                                                      atgaAACTTATTCCTCGCAGTAGTGATATTTCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATGCAGCATATGGTAACGGAGACTGTTTCGGTTTGTACTg                                             >  1:443673/1‑133 (MQ=255)
                                                                            cTTATTCCTCGCAGTAGTGATATTTCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATGCAGCATATGGTAACGGAGACTGTTTCGGTTTGTACTGGCCTAt                                       <  1:113780/133‑1 (MQ=255)
                                                                                     cGCAGTAGTGATATTTCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATGCAGCATATGGTAACGGAGACTGTTTCGGTTTGTACTGGCCTATTGGTCAGGa                              >  2:376572/1‑133 (MQ=255)
                                                                                      gCAGTAGTGATATTTCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATGCAGCATATGGTAACGGAGACTGTTTCGGTTTGTACTGGCCTATTGGTCAGGAg                             >  1:581833/1‑133 (MQ=255)
                                                                                         gtggtGATATTTCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATGCAGCATATGGTAACGGAGACTGTTTCGGTTTGTACTGGCCTATTGGTCAGGAGCAt                          <  2:581833/130‑1 (MQ=255)
                                                                                           agtGATATTTCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATGCAGCATATGGTAACGGAGACTGTTTCGGTTTGTACTGGCCTATTGGTCAGGAGCATaa                        <  2:443673/133‑1 (MQ=255)
                                                                                           agtGATATTTCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATGCAGCATATGGTAACGGAGACTGTTTCGGTTTGTACTGGCCTATTGGTCAGGAGCATaa                        <  1:1741/133‑1 (MQ=255)
                                                                                           agtGATATTTCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATGCAGCATATGGTAACGGAGACTGTTTCGGTTTGTACTGGCCTATTGGTCAGGAGCATaa                        <  1:126102/133‑1 (MQ=255)
                                                                                           agtGATATCTCACCTGGCATTGTTGGAATTGGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATGCAGCATATGGTAACGGAGACTGTTTCGGTTTGTACTGGCCTATTGGTCAGGAGCATaa                        <  1:285135/133‑1 (MQ=255)
                                                                                            gtGATATTTCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATGCAGCATATGGTAACGGAGACTGTTTCGGTTTGTACTGGCCTATTGGTCAGGAGCATaaa                       <  1:44244/133‑1 (MQ=255)
                                                                                                  tttCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATGCAGCATATGGTAACGGAGACTGTTTCGGTTTGTACTGGCCTATTGGTCAGGAGCATAAACTACCt                 >  1:570581/1‑133 (MQ=255)
                                                                                                   ttCACCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATGCAGCATATGGTAACGGAGACTGTTTCGGTTTGTACTGGCCTATTGGTCAGGAGCATAAACTACCta                <  1:376572/133‑1 (MQ=255)
                                                                                                      aCCTGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATGCAGCATATGGTAACGGAGACTGTTTCGGTTTGTACTGGGCTATTGGGCAGAAGGATAAACAACCTATAg             >  1:183223/1‑133 (MQ=255)
                                                                                                         tGGCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATGCAGCATATGGTAACGGAGACTGTTTCGGTTTGTACTGGCCTATTGGTCAGGAGCATAAACTACCTATAgtgt          >  1:325719/1‑133 (MQ=255)
                                                                                                          ggCATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATGCAGCATATGGTAACGGAGACTGTTTCGGTTTGTACTGGCCTATTGGTCAGGAGCATAAACTACCTATAGtgtg         >  1:337402/1‑133 (MQ=255)
                                                                                                            cATTGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATGCAGCATATGGTAACGGAGACTGTTTCGGTTTGTACTGGCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCg       <  1:458411/133‑1 (MQ=255)
                                                                                                               tGATGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATGCAGCATATGGTAACGGAGACTGTTTCGGTTTGTACTGGCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGaaa    <  1:345264/133‑1 (MQ=255)
                                                                                                                 aTGGAATTTGCCCGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATGCAGCATATGGTAACGGAGACTGTTTCGGTTTGTACTGGCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACg  >  2:11178/1‑133 (MQ=255)
                                                                                                                 aTGGAATTTGCCCGGGCCTTTTCCTCCGAACGCATTTACCGTACTTACAGATGCAGCATATGGTAACGGAGACTGTTTCGGTTTGTACTGGCCTATTGGTCAAGAGCATAAACTGCCGATAGTGAGCGAAACg  >  1:287901/1‑133 (MQ=255)
                                                                                                                              |                                                                                                                        
TGGTGGGATTAAAACATTTATATAAACAGATAAGCGGATGTTTACATAAATTTTTATGTGAAGGAACATGATGAAACTTATTCCTCGCAGTAGTGATATTTCACCTGGCATTGATGGAATTTGCCCGGGGCCTTTTCCTCCGAACGGATTTACCGTACTTACAGATGCAGCATATGGTAACGGAGACTGTTTCGGTTTGTACTGGCCTATTGGTCAGGAGCATAAACTACCTATAGTGTGCGAAACG  >  NC_000913/2389046‑2389292

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