Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,789,047 C→T A316V (GCT→GTT) 
*345* (TAG→TAA
yibQ →
waaH ←
putative polysaccharide deacetylase
LPS(HepIII)‑glucuronic acid glycosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,789,0470CT98.2% 169.7 / NA 54A316V (GCT→GTT) 
*345* (TAG→TAA
yibQ
waaH
putative polysaccharide deacetylase
LPS(HepIII)‑glucuronic acid glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (26/27);  total (26/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 = 1.00e+00

TGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACATTGGTAATCAGG  >  NC_000913/3788925‑3789176
                                                                                                                          |                                                                                                                                 
tGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAgg                                                                                                                   >  2:517079/1‑139 (MQ=255)
      tGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTAc                                                                                                             >  1:320016/1‑139 (MQ=255)
       gCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTAcg                                                                                                            <  2:320016/139‑1 (MQ=255)
           aTCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTa                                                                                                        >  2:142795/1‑139 (MQ=255)
            tCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTaa                                                                                                       >  1:118759/1‑139 (MQ=255)
            tCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTaa                                                                                                       >  2:238168/1‑139 (MQ=255)
            tCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTaa                                                                                                       >  1:570380/1‑139 (MQ=255)
                 tcttTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAATAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGa                                                                                                  >  2:300746/1‑139 (MQ=255)
                            ttAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACcg                                                                                       >  2:471241/1‑139 (MQ=255)
                            ttAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACcg                                                                                       >  1:307619/1‑139 (MQ=255)
                                    aaGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGt                                                                               >  2:286322/1‑139 (MQ=255)
                                    aaGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGt                                                                               >  1:677573/1‑139 (MQ=255)
                                    aaGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGt                                                                               >  1:533017/1‑139 (MQ=255)
                                        aTCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTcgcg                                                                           <  2:626221/139‑1 (MQ=255)
                                                 agCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCa                                                                  >  1:238281/1‑139 (MQ=255)
                                                 agCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCa                                                                  >  2:177686/1‑139 (MQ=255)
                                                       cTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCg                                                            >  2:391594/1‑139 (MQ=255)
                                                       cTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCg                                                            >  1:155812/1‑139 (MQ=255)
                                                         gATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGAAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGa                                                          <  1:286322/139‑1 (MQ=255)
                                                           tCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGaaa                                                        <  1:158895/139‑1 (MQ=255)
                                                           tCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGaaa                                                        <  1:238168/139‑1 (MQ=255)
                                                           tCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGaaa                                                        <  2:118759/139‑1 (MQ=255)
                                                              tttACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGAc                                                     <  1:538822/139‑1 (MQ=255)
                                                                    tCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTAcgcg                                                                                                          <  1:436922/80‑1 (MQ=255)
                                                                    tCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTAcgcg                                                                                                          >  2:436922/1‑80 (MQ=255)
                                                                    tCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGa                                                                                                  >  1:474775/1‑88 (MQ=255)
                                                                    tCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGa                                                                                                  <  2:474775/88‑1 (MQ=255)
                                                                       aGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAAc                                            >  2:361633/1‑139 (MQ=255)
                                                                         cATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGcc                                                                   <  1:148005/114‑1 (MQ=255)
                                                                         cATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGcc                                                                   >  2:148005/1‑114 (MQ=255)
                                                                           tCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCt                                        <  2:155812/139‑1 (MQ=255)
                                                                             aGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATa                                                                >  2:601179/1‑113 (MQ=255)
                                                                             aGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATa                                                                <  1:601179/113‑1 (MQ=255)
                                                                             aGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTgg                                      >  2:551075/1‑139 (MQ=255)
                                                                                  cAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAAcc                                 <  2:570380/139‑1 (MQ=255)
                                                                                          gggCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCa                         >  2:193026/1‑139 (MQ=255)
                                                                                              aaaCAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCt                     <  1:391594/139‑1 (MQ=255)
                                                                                                aCAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCAt                                                                       >  2:798706/1‑87 (MQ=255)
                                                                                                aCAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCAt                                                                       <  1:798706/87‑1 (MQ=255)
                                                                                                 cAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTgcg                  <  1:177686/139‑1 (MQ=255)
                                                                                                 cAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTgcg                  <  1:142795/139‑1 (MQ=255)
                                                                                                   gCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGca                <  2:238281/139‑1 (MQ=255)
                                                                                                   gCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGca                <  1:685932/139‑1 (MQ=255)
                                                                                                   gCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGca                <  1:471241/139‑1 (MQ=255)
                                                                                                    cccGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGcac               >  1:338171/1‑139 (MQ=255)
                                                                                                    cccGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGcac               <  2:677573/139‑1 (MQ=255)
                                                                                                      cGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTcgcg                                                                           >  2:226252/1‑77 (MQ=255)
                                                                                                      cGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTcgcg                                                                           <  1:226252/77‑1 (MQ=255)
                                                                                                               gAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACATTGGTAATCa    <  1:647133/139‑1 (MQ=255)
                                                                                                               gAAGTTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACATGGGTAATCa    <  2:522906/139‑1 (MQ=255)
                                                                                                                 aGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACATTGGTAATCAgg  <  1:551075/139‑1 (MQ=255)
                                                                                                                 aGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACATTGGTAATCAgg  <  2:533017/139‑1 (MQ=255)
                                                                                                                 aGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACATTGGTAATCAgg  <  1:300746/139‑1 (MQ=255)
                                                                                                                 aGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACACTGGTAATCAgg  <  2:135916/139‑1 (MQ=255)
                                                                                                                          |                                                                                                                                 
TGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACATTGGTAATCAGG  >  NC_000913/3788925‑3789176

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