Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,798,666 A→T 12.4% L216* (TTA→TAA)  waaU ← putative ADP‑heptose:LPS heptosyltransferase 4

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,798,6660AT12.4% 160.0 / 3.1 63L216* (TTA→TAA) waaUputative ADP‑heptose:LPS heptosyltransferase 4
Reads supporting (aligned to +/- strand):  ref base A (31/23);  new base T (6/2);  total (37/26)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCTTCAGTATTTCTATTGCTCGTGTAC  >  NC_000913/3798529‑3798807
                                                                                                                                         |                                                                                                                                             
tttCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGa                                                                                                                                          >  1:557202/1‑143 (MQ=255)
 ttCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGAt                                                                                                                                         >  1:279166/1‑143 (MQ=255)
    aaGTATTGGTATTGTCAATAGATCTTGCGGTAGCCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCTCCTAATGGATtag                                                                                                                                      >  2:489902/1‑142 (MQ=255)
    aaGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACTTAATGGATtaa                                                                                                                                      >  2:208639/1‑143 (MQ=255)
       tattggtattgTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATATGCTCAAACGTAAGACGGCATATTTTATTTGCCACTATTGGATtactaa                                                                                                                                   >  2:711754/1‑143 (MQ=255)
        attggtattgTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTaataat                                                                                                                                  <  2:539157/143‑1 (MQ=255)
         ttggtattgTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATg                                                                                                                                 >  1:757838/1‑143 (MQ=255)
          tggtattgTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGa                                                                                                                                >  1:184860/1‑143 (MQ=255)
           ggtattgTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCCCCTTATGGATTTATTATggc                                                                                                                               >  1:532782/1‑141 (MQ=255)
                  tCAATTAATCTTGCGGTAACCCGTTAATTATAATCCGATAATTTTCAAATTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCtt                                                                                                                        <  2:779929/143‑1 (MQ=255)
                   cAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATGTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTg                                                                                                                       <  1:463652/143‑1 (MQ=255)
                    aaTAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGt                                                                                                                      <  1:852885/143‑1 (MQ=255)
                    aaTAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGGTTTCTCTTCTTGATGTATAACTTTTATTTGCTCAAACGTAAGACGGCGTATTTTCTTTTCATCTTATGGATGAATAGTGTCTATTcgcgt                                                                                                                      >  2:874529/1‑139 (MQ=255)
                     aTAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTTCTCCTTTTGGATTTATAATGTCTATttttttt                                                                                                                     >  2:903807/1‑140 (MQ=255)
                      tAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGtct                                                                                                                    >  2:725315/1‑141 (MQ=255)
                         aTCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTAtctc                                                                                                                 <  2:441331/143‑1 (MQ=255)
                          tCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTcc                                                                                                                >  2:689241/1‑143 (MQ=255)
                            ttGCGGTATCCCGGTAAATATAATCTGATATTTTTCAAAGTGTGTTTTCACTTCTTGATTTATAACTTTTATTTGCTCAAACGTCAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCaa                                                                                                              <  2:464019/143‑1 (MQ=255)
                            ttGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGTACCTTGTGGATTAATAATGACTACTCTTGTCTCTCCaa                                                                                                              >  2:658703/1‑143 (MQ=255)
                                ggTAAGCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTGGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATaaa                                                                                                          <  1:112728/143‑1 (MQ=255)
                                ggTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATaaa                                                                                                          >  1:720220/1‑143 (MQ=255)
                                ggTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACGTTTATTTGCTCACACGTAAGACGGCATATTTTTTTTGCACCTAATGGGTTAATAATGACTATTCTTGTATATCCactcaa                                                                                                          >  1:86360/1‑138 (MQ=255)
                                   aaCCCGTTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCt                                                                                                      <  1:342390/143‑1 (MQ=255)
                                         gTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCACTAAATTAttttat                                                                                                 >  2:625188/1‑143 (MQ=255)
                                          tAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCttttatt                                                                                                <  2:932517/143‑1 (MQ=255)
                                            aaTATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCttttatttt                                                                                              <  1:824371/143‑1 (MQ=255)
                                             atatAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCttttatttta                                                                                             >  1:346487/1‑143 (MQ=255)
                                               ataATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATc                                                                                           >  1:883655/1‑143 (MQ=255)
                                                      gATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACa                                                                                    >  1:545618/1‑143 (MQ=255)
                                                         aaTTTTCAAAGTGTGTTTTCACTCCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCCATCAATTCTTTTATTTTATCTTCAACctct                                                                                 >  1:871645/1‑143 (MQ=255)
                                                             ttCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTg                                                                             <  1:913312/143‑1 (MQ=255)
                                                                  aGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTaaa                                                                        <  1:144884/143‑1 (MQ=255)
                                                                       gtTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCaa                                                                   <  1:447540/143‑1 (MQ=255)
                                                                                  ttGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATTACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTtgg                                                        <  2:763452/143‑1 (MQ=255)
                                                                                    gATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCACTAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTtggtc                                                      >  1:339285/1‑142 (MQ=255)
                                                                                    gATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTtggtgg                                                     >  1:359232/1‑143 (MQ=255)
                                                                                       atataACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTTTTTTTTTTTTATTAATAATGACTATTCTTGTATCTTTTTTAATTTCTTTTATTTTATCTTCATCATCTACTGGTTATTTCTAGTCTTATTTGGTGGaa                                                   >  2:803570/1‑143 (MQ=255)
                                                                                        tataACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTTATGAATTAATAATGAATATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTACATGCAAGTATTATTTGGTGGaaa                                                  >  2:39243/1‑143 (MQ=255)
                                                                                        tataACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGaaa                                                  <  2:852806/143‑1 (MQ=255)
                                                                                        tataACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGaaa                                                  >  1:943220/1‑143 (MQ=255)
                                                                                         ataACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGaaaa                                                 <  2:545618/143‑1 (MQ=255)
                                                                                            aCTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAAttt                                              <  1:761513/143‑1 (MQ=255)
                                                                                             cTTTTATTTGCTCAAACGTACGACGGCAGGTGTTATTAATGGCTTAGTTTTTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTa                                             <  1:933077/143‑1 (MQ=255)
                                                                                             cTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTGTGCACCTAATGGATGAATCCTGAACATCTTTTTATCTCCAATAAATTCGTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTa                                             <  1:938351/143‑1 (MQ=255)
                                                                                              ttttATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAACTCGTAGTTGGTGGAAAAtttcc                                            >  1:160664/1‑141 (MQ=255)
                                                                                              tcttATTTGCTCAATTGTAAGACTGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTAc                                            <  1:164107/141‑1 (MQ=255)
                                                                                                    ttGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAACTCGTCGTTGGTGGAAAATTTCCCCTccc                                      >  1:16693/1‑141 (MQ=255)
                                                                                                        tCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATa                                  >  1:559737/1‑143 (MQ=255)
                                                                                                             cGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAatat                             <  1:91530/143‑1 (MQ=255)
                                                                                                               tAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGc                           <  2:315354/143‑1 (MQ=255)
                                                                                                                 aGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTTCAATAAATTCTTTTATTTTATCTTATACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCtt                         >  2:109587/1‑143 (MQ=255)
                                                                                                                  gACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTATTGTATCTACAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATTCAACTCGTAGTTGGTGGAAAATTTACCCTGACCATAACTATggtat                        >  2:649285/1‑139 (MQ=255)
                                                                                                                       cATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTATTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCTTCAGTAt                   >  1:500443/1‑143 (MQ=255)
                                                                                                                        catatttttttGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTaaa                                                                        <  2:73732/85‑1 (MQ=255)
                                                                                                                        catatttttttGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTaaa                                                                        >  1:73732/5‑89 (MQ=255)
                                                                                                                           ttttttttGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCTTCAGTATTTCt               >  1:775491/1‑143 (MQ=255)
                                                                                                                             ttttttGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCGTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATTTGCTTCAGTATTTatat             >  1:148011/1‑143 (MQ=255)
                                                                                                                               ttttGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCTTCAGTATTTCTATTg           >  1:211655/1‑143 (MQ=255)
                                                                                                                                  tttaCCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTCCTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATTCTTCAGTATTTCTATTGCTc        >  2:830485/4‑143 (MQ=255)
                                                                                                                                  tGCACCTTATGTATTAATAACGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACACCTACTGGTAAATGCAAGTCGTAGTGGGTGGAAAAGTTACCCTCACCATAAATATGCTTCAGTATTTCTATTGCTc        <  2:292032/143‑1 (MQ=255)
                                                                                                                                    cACCTGATGGATTAATAATCACTATTCTTGTAGCTCCAATAAATTCTTTCATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCTTCAGTATTTCTATTGCTCgt      <  1:60013/143‑1 (MQ=255)
                                                                                                                                       cTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATGTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGCAAAATTTACCCTCACCATAAATATGATTCAGTATTTCTATTGCTccggta   >  2:739908/1‑138 (MQ=255)
                                                                                                                                        tAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATACATATGCTTCAGTATTTCTATTGCTCGTGTAc  <  1:657437/143‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                             
TTTCAAGTATTGGTATTGTCAATAAATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGTAGTTGGTGGAAAATTTACCCTCACCATAAATATGCTTCAGTATTTCTATTGCTCGTGTAC  >  NC_000913/3798529‑3798807

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: