Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,805,527 T→A 9.5% D186V (GAC→GTC)  ydiY ← acid‑inducible putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,805,5270TA 9.5% 150.3 / 5.4 52D186V (GAC→GTC) ydiYacid‑inducible putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base T (27/20);  new base A (2/3);  total (29/23)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTCTGGAGGTTCAGAGTTCCATGTGACGTTGTAAGCCACCTTCAGTCCAAAATGTTCATTAATAGCAACATTTAATGCGCTTTCGGAGTTCAGCGTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCCGGAGGCATAGCCCAGCGGCTGGGTTTCACTGGCATTGTCGGTATATTTGTCGTAACGCACACCTGGACCGAATTCAAAGCGGAAACTGTGTACCGGGCCATTCAAAAACTGACGACCA  >  NC_000913/1805381‑1805670
                                                                                                                                                  |                                                                                                                                               
aCTCTGGAGGTTCAGAGTTCCATGTGACGTTGTAAGCCACCTTCAGTCCAAAATGTTCATTAATAGCAACATTTAATGCGCTTTCGGAGTTCAGCGTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTgtc                                                                                                                                                <  8:36557/148‑1 (MQ=255)
      gAGGTTCAGCGTTGCATGTGACGTTGTAAGCCACCCTCAGTCAAAAATGTTCATTAATAGCAACATTTAATGCGATTACGGAGTTCAGCGTTGTGTCTTCCGCGCCAACGACCGAAACACCCTGCGTAAATTTACCGGTGACAGTCaa                                                                                                                                          <  6:79302/148‑1 (MQ=255)
       aGGTTCAGAGTTCCATGTGACGTTGTAAGCCACCTTCAGTCCAAAATGTTCATTAATAGCAACATTTAATGCGCTTTCGGAGTTCAGCGTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAAc                                                                                                                                         >  4:230649/1‑148 (MQ=255)
       aGGTTCAGAGTTCCATGTGACGTTGTAAGCCACCTTCAGTCCAAAATGTTCATTAATAGCAACATTTAATGCGCTTTCGGAGTTCAGCGTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAAc                                                                                                                                         >  7:290710/1‑148 (MQ=255)
              gagTTCCATGTGACGTTGTAAGCCACCTTCAGTCCAAAATGTTCATTAATAGCAACATTTAATGCGCTTTCGGAGTTCAGCGTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATg                                                                                                                                  <  3:155788/148‑1 (MQ=255)
               agTTCCATGTGACGTTGTAAGCCACCTTCAGTCCAAAATGTTCATTAATAGCAACATTTAATGCGCTTTCGGAGTTCAGCGTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATg                                                                                                                                  >  3:203500/1‑147 (MQ=255)
                  tCCATGTGACGTTGTAAGCCACCTTCAGTCCAAAATGTTCATTAATAGCAACATTTAATGCGCTTTCGGAGTTCAGCGTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGACAGTCAACTGCCATGCAt                                                                                                                               >  3:256836/1‑147 (MQ=255)
                  tCCATGTGACGTTGTAAGCCACCTTCAGTCCAAAATGTTCATTAATAGCAACATTTAATGCGCTTTCGGAGTTCAGCGTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGACAGTCAACTGCCATGCAt                                                                                                                               >  3:260003/1‑147 (MQ=255)
                             ttGTAAGCCACCTTCAGTCCAAAATGTTCATTAATAGCAACATTTAATGCGCTTTCGGAGTTCAGCGTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCCGGAgg                                                                                                                   <  8:290710/148‑1 (MQ=255)
                              tGTACGCCACCTTCATTCCAAAAAGTACATTAATAGCACCATTTAATTCGCTTTCGGAGTTCAGCGATGTGTATTACGAGCCAAAGACCGAAACACACTGCCTAAATTTAGCGTTGTCAGTCAACTGCCATGCAGAAGCCCCGCAgg                                                                                                                   >  2:244608/1‑147 (MQ=255)
                              tGTAAGCCACCTTCAGTCCAAAATGTTCATTAATAGCAACATTTAATGCGCTTTCGGAGTTCAGCGTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCCGGAgg                                                                                                                   >  7:167281/1‑147 (MQ=255)
                                        cTTCAGTCCAAAATGTTCATTAATAGCAACATTTAATGCGCTTTCGGAGTTCAGCGTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCTTTGTCAGTCAACTGCCATGCATAAGCCCCGGa                                                                                                                     <  5:93463/135‑1 (MQ=255)
                                        cTTCAGTCCAAAATGTTCATTAATAGCAACATTTAATGCGCTTTCGGAGTTCAGCGTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCCGGAGGCATAGCCCAGc                                                                                                        >  3:136481/1‑148 (MQ=255)
                                        cTTCAGTCCAAAATGTTCATTAATAGCAACATTTAATGCGCTTTCGGAGTTCAGCGTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCATTGTCAGTCAACTGCCATGCATAAGCCCCGGa                                                                                                                     >  6:93463/1‑135 (MQ=255)
                                          tCAGTCCAAAATGTTCATTAATAGCAACATTTAATGCGCTTTCGGAGTTCAGCGTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCCGGa                                                                                                                     >  1:123103/1‑133 (MQ=255)
                                          tCAGTCCAAAATGTTCATTAATAGCAACATTTAATGCGCTTTCGGAGTTCAGCGTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCCGGa                                                                                                                     <  2:123103/133‑1 (MQ=255)
                                               ccAAAATGTTCATTAATAGCAACATTTAATGCGCTTTCGGAGTTCAGCGTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCCGGAGGCATAGCCCAGCGGCTggg                                                                                                 >  6:165088/1‑148 (MQ=255)
                                                cAAAATGTTCATTAATAGCAACATTTAATGCGCTTTCGGAGTTCAGCGTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCCGGAGGCATAGCCCAGCGGCTggg                                                                                                 <  4:203500/147‑1 (MQ=255)
                                                  aaaTGTTCATTAATAGCAACATTTAATGCGCTTTCGGAGTTCAGCGTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGcccc                                                                                                                        <  1:138014/122‑1 (MQ=255)
                                                  aaaTGTTCATTAATAGCAACATTTAATGCGCTTTCGGAGTTCAGCGTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGcccc                                                                                                                        >  2:138014/1‑122 (MQ=255)
                                                    aTGTTCATTAATAGCAACATTTAATGCGCTTTCGGAGTTCAGCGTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCCGGAGGCATAGCCCAGCGGCTGGGTTTCa                                                                                            <  3:230649/148‑1 (MQ=255)
                                                        tCATTAATAGCAACATTTAATGCGCTTTCGGAGTTCAGCGTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCCGGAGGCATAGCCCAGCGGCTGGGTTTCACTgg                                                                                        >  7:26048/1‑148 (MQ=255)
                                                            tAATAGCAACATTTAATGCGCTTTCGGAGTTCAGCGTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCCGGAGGCATAGCCCAGCGGCTGGGTTTCa                                                                                            >  7:151731/1‑140 (MQ=255)
                                                            tAATAGCAACATTTAATGCGCTTTCGGAGTTCAGCGTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCCGGAGGCATAGCCCAGCGGCTGGGTTTCa                                                                                            <  8:151731/140‑1 (MQ=255)
                                                              aTAGCAACATTTAATGCGCTTTCGGAGTTCAGCGTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCCgg                                                                                                                      >  5:18009/1‑112 (MQ=255)
                                                              aTAGCAACATTTAATGCGCTTTCGGAGTTCAGCGTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCCGGAGGCATAGCCCAGCGGCTGGGTTTCACTGGCATTg                                                                                   <  8:167281/147‑1 (MQ=255)
                                                              aTAGCAACATTTAATGCGCTTTCGGAGTTCAGCGTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCAgg                                                                                                                      <  6:18009/112‑1 (MQ=255)
                                                                 gCAACATTTAATGCGCTTTCGGAGTTCAGCGTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCCGGAGGCATAGCCCAGCGGCTGGGTTTCACTGGCATTGTCgg                                                                               >  3:280541/1‑148 (MQ=255)
                                                                  cAACATTTAATGCGCTTTCGGAGTTCAGCGTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCCGGAGGCATAGCCCAGCGGCTGGGTTTCACTGGCATTGTCgg                                                                               >  3:150225/1‑147 (MQ=255)
                                                                    aCATTTAATGCGCTTTCGGAGTTCAGCGTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCCGGAGGCATAGCCCAGc                                                                                                        >  7:334439/1‑120 (MQ=255)
                                                                    aCATTTAATGCGCTTTCGGAGTTCAGCGTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCCGGAGGCATAGCCCAGc                                                                                                        <  8:334439/120‑1 (MQ=255)
                                                                           aTGCGCTTTCGGAGTTCAGCGTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCCGGAGGCATAGCCCAGCGGCTGGGTTTCACTGGCATTGTCGGTATATTTGt                                                                      <  5:165088/147‑1 (MQ=255)
                                                                           aTGCGCTTTCGGAGTTCAGCGTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGACAGTCAACTGCCATGCATAAGCCCCGGAGGCATAGCCCAGCGGCTGGGTTTCACTGGCATTGTCGGTATATTTGt                                                                      <  4:260003/147‑1 (MQ=255)
                                                                           aTGCGCTTTCGGAGTTCAGCGTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGACAGTCAACTGCCATGCATAAGCCCCGGAGGCATAGCCCAGCGGCTGGGTTTCACTGGCATTGTCGGTATATTTGt                                                                      <  4:256836/147‑1 (MQ=255)
                                                                                          tCAGCGTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCCGGAGGCATAGCCCAGCGGATGGGTTTCACTGGCATTGTCGGTATATTTGTCGTAACGCACACCTg                                                       >  4:245277/1‑147 (MQ=255)
                                                                                              cGTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCCGGAGGCATAGCCCAGCGGCTGGGTTTCACTGGCATTGTCGGTATATTTGTCGTAACGCACACCTGGAcc                                                   >  6:218416/1‑147 (MQ=255)
                                                                                               gTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCCGGAGGCATAGCCCAGCGGCTGGGTTTCACTGGCATTGTCGGTATATTTGTCGTAACGCACACCTGGACCg                                                  >  6:245535/1‑147 (MQ=255)
                                                                                                   tgtCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCCGGAGGCATAGCCCAGCGGCTGGGTTTCACTGGCATTGTCGGTATATTTGTCGTAACGCACACCt                                                        >  5:168134/1‑137 (MQ=255)
                                                                                                   tgtCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCCGGAGGCATAGCCCAGCGGCTGGGTTTCACTGGCATTGTCGGTATATTTGTCGTAACGCACACCt                                                        <  6:168134/137‑1 (MQ=255)
                                                                                                        tCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCCGGAGGCATAGccc                                                                                                           <  3:66150/81‑1 (MQ=255)
                                                                                                        tCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCCGGAGGCATAGccc                                                                                                           >  4:66150/1‑81 (MQ=255)
                                                                                                         ccGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCCGGAGGCATAGCCCAGCGGCTGGGTTTCACTGGCATTGTCGGTATATTTGTCGTa                                                                  <  4:260348/121‑1 (MQ=255)
                                                                                                         ccGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCCGGAGGCATAGCCCAGCGGCTGGGTTTCACTGGCATTGTCGGTATATTTGTCGTa                                                                  >  3:260348/1‑121 (MQ=255)
                                                                                                          cgcgCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCCGGAGGCATAGCCCAGCGGCTGGGTTTCACTGGCATTGTCGGTATATTTGTCGTAACGCACACCTGGACCGAATTCAAAGCgg                                      <  4:136481/148‑1 (MQ=255)
                                                                                                          cgcgCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCCGGAGGCATAGCCCAGCGGCTGGGTTTCACTGGCATTGTCGGTATATTTGTCGTAACGCACACCTGGACCGAATTCAAAGCgg                                      <  5:245535/148‑1 (MQ=255)
                                                                                                          cgcgCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCCGGAGGCATAGCCCAGCGGCTGGGTTTCACTGGCATTGTCGGTATATTTGTCGTAACGCACACCTGGACCGAATTCAAAGCg                                       <  5:218416/147‑1 (MQ=255)
                                                                                                                   gACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCCGGAGGCATAGCCCAGCGGCTGGGTTTCACTGGCATTGTCGGTATATTTGTCGTAACGCACACCTGGACCGAATTCAAAGCGGAAACTgtgt                             <  4:150225/148‑1 (MQ=255)
                                                                                                                       gAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCCGGAGGCATAGCCCAGCGGCTGGGTTTCACTGGCATTGTCGGTATATTTGTCGTAACGCACACCTGGACCGAATTCAAAGCGGAAACTGTGTACCg                         >  5:138282/1‑148 (MQ=255)
                                                                                                                        aaaCACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCCGGAGGCATAGCCCAGCGGCTGGGTTTCACTGGCATTGTCGGTATATTTGTCGTAACGCACACCTGGACCGAATTCAAAGCGGAAACTGTGTACCgg                        >  1:90666/1‑148 (MQ=255)
                                                                                                                                   gTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCCGGAGGCATAGCCCAGCGGCTGGGTTTCACTGGCATTGTCGGTATATTTGTCGTAACGCACACCTGGACCGAATTCAAAGCGGAAACTGTGTACCGGGCCATTCaaaa             <  8:26048/148‑1 (MQ=255)
                                                                                                                                             cGTTGTCAGTCAACTGCCATGCATAAGCCCCGGAGGCATAGCCCCGCGGCTGGGTTTCACTGGCATTGTCGGTATATTTGTCGTAACGCACACCTGGGCCGAATTCAACGCGGAAACTGTGTACCGGGCCATTCAAAAACTgacgac    >  5:278967/1‑147 (MQ=255)
                                                                                                                                             cGTTGTCAGTCAACTGCCATGCATAAGCCCCGGAGGCATAGCCCAGCGGCTGGGTTTCACTGGCATTGTCGGTATATTTGTCGTAACGCACACCTGGACCGAATTCAAAGCGGAAACTGTGTACCGGGCCATTCAAAAACTgacgac    >  6:161910/1‑147 (MQ=255)
                                                                                                                                              gTTGTCAGTCAACTGCCATGCATAAGCCCCGGAGGCATAGCCCAGCGGCTGGGTTTCACTGGCATTGTCGGTATATTTGTCGTAACGCACACCTGGACCGAATTCAAAGCGGAAACTGTGTACCGGGCCATTCAAAAACTGACGACCa  >  2:137221/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
ACTCTGGAGGTTCAGAGTTCCATGTGACGTTGTAAGCCACCTTCAGTCCAAAATGTTCATTAATAGCAACATTTAATGCGCTTTCGGAGTTCAGCGTTGTGTCTTCCGCGCCAAAGACCGAAACACCCTGCGTAAATTTAGCGTTGTCAGTCAACTGCCATGCATAAGCCCCGGAGGCATAGCCCAGCGGCTGGGTTTCACTGGCATTGTCGGTATATTTGTCGTAACGCACACCTGGACCGAATTCAAAGCGGAAACTGTGTACCGGGCCATTCAAAAACTGACGACCA  >  NC_000913/1805381‑1805670

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