Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 303,093 (G)6→5 coding (103/615 nt) yagU → inner membrane protein that contributes to acid resistance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913303,0880G.98.2% 225.2 / ‑1.9 54coding (98/615 nt)yagUinner membrane protein that contributes to acid resistance
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base . (27/26);  total (27/27)
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

TCCTGAAGAGGTGTATAACATAATGGTACAAGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACC  >  NC_000913/302949‑303211
                                                                                                                                           |                                                                                                                           
tCCTGAAGAGGTGTATAACATAATGGTACAAGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTg                                                                                                                             >  1:376815/1‑140 (MQ=255)
    gAAGAGGTGTATAACATAATGGTACAAGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTggggg                                                                                                                         <  2:376012/140‑1 (MQ=255)
     aaGAGGTGTATAACATAATGGTACAAGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTgggggc                                                                                                                        >  1:606957/1‑139 (MQ=255)
          gtgtATAACATAATGGTATAAGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCACAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAg                                                                                                                  >  1:593096/1‑140 (MQ=255)
          gtgtATAACATAATGGTACAAGCAGGGTGGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAg                                                                                                                  >  2:116240/1‑140 (MQ=255)
          gtgtATAACATAATGGTACAAGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAg                                                                                                                  >  1:1798/1‑140 (MQ=255)
          gtgtATAACATAATGGTACAAGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAg                                                                                                                  >  2:274894/1‑140 (MQ=255)
          gtgtATAACATAATGGTACAAGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAg                                                                                                                  >  1:207038/1‑140 (MQ=255)
          gtgtATAACATAATGGTACAAGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAg                                                                                                                  >  2:722552/1‑140 (MQ=255)
            gtATAACATAATGGTACAAGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGtt                                                                                                                >  1:242690/1‑140 (MQ=255)
                aaCATAATGGTACAAGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCAt                                                                                                            <  2:606957/140‑1 (MQ=255)
                  cATAATGGTACAAGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTg                                                                                                          >  1:576946/1‑140 (MQ=255)
                    tAATGGTACAAGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTgcc                                                                                                        >  2:168969/1‑140 (MQ=255)
                    tAATGGTACAAGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTgcc                                                                                                        >  2:495423/1‑140 (MQ=255)
                     aaTGGTACAAGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTgccg                                                                                                       >  2:561880/1‑140 (MQ=255)
                     aaTGGTACAAGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTgccg                                                                                                       <  1:561880/140‑1 (MQ=255)
                     aaTGGTACAAGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGCTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTgccg                                                                                                       <  1:206712/140‑1 (MQ=255)
                       tGGTACAAGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTgccgcc                                                                                                     >  2:403766/1‑140 (MQ=255)
                          tACAAGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTGCCGCCACg                                                                                                  >  1:933663/1‑140 (MQ=255)
                          tACAAGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTGCCGCCACg                                                                                                  <  1:722552/140‑1 (MQ=255)
                             aaGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTGCCGCCACGTAg                                                                                               >  1:417208/1‑140 (MQ=255)
                               gCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTgccgcc                                                                                                     <  1:814465/132‑1 (MQ=255)
                               gCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTgccgcc                                                                                                     >  2:814465/1‑132 (MQ=255)
                               gCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTGCCGCCACGTAGcc                                                                                             <  1:35007/140‑1 (MQ=255)
                               gCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTGCCGCCACGTAGcc                                                                                             >  1:273921/1‑140 (MQ=255)
                                  gggTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCgg                                                                                          >  1:924819/1‑140 (MQ=255)
                                    gTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTg                                                                                        <  2:273921/140‑1 (MQ=255)
                                       gATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTGGat                                                                                     <  2:225618/140‑1 (MQ=255)
                                       gATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTGGat                                                                                     <  2:933663/140‑1 (MQ=255)
                                       gATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTGGat                                                                                     <  1:116240/140‑1 (MQ=255)
                                       gATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCATGAAGTGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTGGat                                                                                     <  1:219602/140‑1 (MQ=255)
                                        atatGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTGGata                                                                                    >  1:621644/1‑140 (MQ=255)
                                             aaTATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTGGATATGttt                                                                               >  1:321868/1‑140 (MQ=255)
                                               tatatTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTGGat                                                                                     >  2:17239/1‑132 (MQ=255)
                                               tatatTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTGGat                                                                                     <  1:17239/132‑1 (MQ=255)
                                                    ttGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTGGATATGTTTAATGCAg                                                                        <  1:84089/140‑1 (MQ=255)
                                                            aCTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTGGATATGTTTAATGCAGCGTGTGGc                                                                >  2:533113/1‑140 (MQ=255)
                                                             cTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTGGATATGTTTAATGCAGCGTGTGGcc                                                               <  1:513135/140‑1 (MQ=255)
                                                                cACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTGGATATGTTTAATGCAGCGTGTGGCCCgg                                                            <  1:422426/140‑1 (MQ=255)
                                                                cACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTGGATATGTTTAATGCAGCGTGTGGCCCgg                                                            <  2:924819/140‑1 (MQ=255)
                                                                 aCCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGa                                                           >  2:208934/1‑140 (MQ=255)
                                                                     aaCCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCa                                                       <  2:242690/140‑1 (MQ=255)
                                                                       ccGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCAtt                                                     >  1:604124/1‑140 (MQ=255)
                                                                        cGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTa                                                    >  1:218414/1‑140 (MQ=255)
                                                                                  aTGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTg                                          <  2:321868/140‑1 (MQ=255)
                                                                                   tGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGc                                         <  2:1798/140‑1 (MQ=255)
                                                                                     gTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAg                                       <  2:218414/140‑1 (MQ=255)
                                                                                         tgctgcATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCa                                   >  2:278526/1‑140 (MQ=255)
                                                                                                ttcattGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAattgat                            >  1:29745/1‑140 (MQ=255)
                                                                                                   attGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGc                         <  2:417208/140‑1 (MQ=255)
                                                                                                           gATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTa                  <  2:376815/140‑1 (MQ=255)
                                                                                                            aTTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAAt                <  1:533113/140‑1 (MQ=255)
                                                                                                             ttGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGc                         <  1:114266/130‑1 (MQ=255)
                                                                                                             ttGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGc                         >  2:114266/1‑130 (MQ=255)
                                                                                                                     gttgttTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAAt       <  1:208934/140‑1 (MQ=255)
                                                                                                                      ttgttTCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATc      <  2:207038/140‑1 (MQ=255)
                                                                                                                          ttCCGCATTCGTGAAGTGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCAcc  <  1:495423/140‑1 (MQ=255)
                                                                                                                                           |                                                                                                                           
TCCTGAAGAGGTGTATAACATAATGGTACAAGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACC  >  NC_000913/302949‑303211

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