Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 146,238 A→G K386K (AAA→AAG)  yadE → putative polysaccharide deacetylase lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913146,2380AG98.0% 169.5 / NA 50K386K (AAA→AAG) yadEputative polysaccharide deacetylase lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (24/25);  minor base C (0/1);  total (24/26)
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

AATTTAATCCGCATGTCTGGTATCTTTCGTATCCGTTTGGCGGATTTAATGACAACGCCGTGAAGGCAGCAAACGATGCCGGATTTCACCTGGCGGTGACAACCATGAAAGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTAAAACGACTTTATATCTTAAGAACGGATTCGCTGGAGACGATGTCGCGGCTGGTGAGTAACCAGCCGCAGGGATAACAATCAAGCAACCTGTACCGGAATCGCTTTCGCGGTACGTTTCATTTCATTGTCGCCTTCAAAATAGGCGACG  >  NC_000913/146090‑146383
                                                                                                                                                    |                                                                                                                                                 
aaTTTAATCCGCATGTCTGGTATCTTTCGTATCCGTTTGGCGGATTTAATGACAACGCCGTGAAGGCAGCAAACGATGCCGGATTTCACCTGGCGGTGACAACCATGAAAGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTAAAGCg                                                                                                                                                 <  2:687655/151‑1 (MQ=255)
    tAATCCGCATGTCTGGTATCTTTCGTATCCGTTTGGCGGATTTAATGACAACGCCGTGAAGGCAGCAAACGATGCCGGATTTCACCTGGCGGTGACAACCATGAAAGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTAAAGCGACt                                                                                                                                              <  1:191860/150‑1 (MQ=255)
    tAATCCGCATGTCTGGTATCTTTCGTATCCGTTTGGCGGATTTAATGACAACGCCGTGAAGGCAGCAAACGATGCCGGATTTCACCTGGCGGTGACAACCATGAAAGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTAAAGCGACt                                                                                                                                              <  2:320537/150‑1 (MQ=255)
     aaTCCGCATGTCTGGTATCTTTCGTATCCGTTTGGCGGATTTAATGACAACGCCGTGAAGGCAGCAAACGATGCCGGATTTCACCTGGCGGTGACAACCATGAAAGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTAAAGCGACttt                                                                                                                                            <  2:242796/151‑1 (MQ=255)
      aTCCGCATGTCTGGTATCTTTCGTATCCGTTTGGCGGATTTAATGACAACGCCGTGAAGGCAGCAAACGATGCCGGATTTCACCTGGCGGTGACAACCATGAAAGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTAAAGCGACttt                                                                                                                                            >  1:327800/1‑150 (MQ=255)
       tCCGCATGTCTGTTATCTTTCGTATCCGTTTGGCGGATTTAATGACAACGCCGTGAAGGCAGCAAACGATGCCGGATTTCACCTGGCGGTGACAACCATGCAAGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTAAAGCGACGTta                                                                                                                                           <  1:513013/150‑1 (MQ=255)
           cATGTCTGGTATCTTTCGTATCCGTTTGGCGGATTTAATGACAACGCCGTGAAGGCAGCAAACGATGCCGGATTTCACCTGGCGGTGACAACCATGAAAGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTAAAGCGACTTTATATCt                                                                                                                                      <  1:743125/151‑1 (MQ=255)
           cATGTCTGGTATCTTTCGTATCCGTTTGGCGGATTTAATGACAACGCCGTGAAGGCAGCAAACGATGCCGGATTTCACCTGGCGGTGACAACCATGAAAGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTAAAGCGACTTTATATCt                                                                                                                                      <  2:351811/151‑1 (MQ=255)
           cATGTCTGGTATCTTTCGTATCCGTTTGGCGGATTTAATGACAACGCCGTGAAGGCAGCAAACGATGCCGGATTTCACCTGGCGGTGACAACCATGAAAGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTAAAGCGACTTTATATCt                                                                                                                                      <  1:525/151‑1 (MQ=255)
            aTGTCTGGTATCTTTCGTATCCGTTTGGCGGATTTAATGACAACGCCGTGAAGGCAGCAAACGATGCCGGATTTCACCTGGCGGTGACAACCATGAAAGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTAAAGCGACTTTATATCtt                                                                                                                                     <  1:286833/151‑1 (MQ=255)
              gTCTGGTATCTTTCGTATCCGTTTGGCGGATTTAATGACAACGCCGTGAAGGCAGCAAACGATGCCGGATTTCACCTGGCGGTGACAACCATGAAAGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTAAAGCGACTTTATATCTTaa                                                                                                                                   <  2:174564/151‑1 (MQ=255)
                 tGGTATCTTTCGTATCCGTTTGGCGGATTTAATGACAACGCCGTGAAGGCAGCAAACGATGCCGGATTTCACCTGGCGGTGACAACCATGAAAGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTAAAGCGACTTTATATCTTAAGa                                                                                                                                 >  2:150344/1‑150 (MQ=255)
                   gTATCTTTCGTATCCGTTTGGCGGATTTAATGACAACGCCGTGAAGGCAGTAAACGATGCCGGATTTCACCTGGCGGTGACAACCATGAAAGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTAAAGCGACTTTATATCTTAAGAACg                                                                                                                              <  1:816843/151‑1 (MQ=255)
                     aTCTTTCGTATCCGTTTGGCGGATTTAATGACAACGCCGTGAAGGCAGCAAACGATGCCGGATTTCACCTGGCGGTGACAACCATGAAAGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTAAAGCGACTTTATATCTTAAGAACgg                                                                                                                             >  2:358426/1‑150 (MQ=255)
                      tCTTTCGTATCCGTTTGGCGGATTTAATGACAACGCCGTGAAGGCAGCAAACGATGCCGGATTTCACCTGGCGGTGACAACCATGAAAGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTAAAGCGACTTTATATCTTAAGAACGGAt                                                                                                                           >  1:356820/1‑151 (MQ=255)
                         ttCGTATCCGTTTGGCGGATTTAATGACAACGCCGTGAAGGCAGCAAACGATGCCGGATTTCACCTGGCGGTGACAACCATGAAAGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTAAAGCGACTTTATATCTTAAGAACGGATTc                                                                                                                         <  1:460260/150‑1 (MQ=255)
                                  gTTTGGCGGATTTAATGACGACGCCGTGAAGGCAGCAAACGATGCCGGATTTCCCCTGGCGTTGACAACCATGAAAGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTACAGCGACTTTATATCTTAAGAACGGATTCGCTGGAGAc                                                                                                                <  1:325866/150‑1 (MQ=255)
                                      ggCGGATTTAATGACAACGCCGTGAAGGCAGCAAACGATGCCGGATTTCACCTGGCGGTGACAACCATGAAAGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTAAAGCGACTTTATATCTTAAGAACGGATTCGCTGGAGACGATGt                                                                                                           >  1:408221/1‑151 (MQ=255)
                                       gCGGATTTAATGACAACGCCGTGAAGGCAGCAAACGATGCCGGATTTCACCTGGCGGTGACAACCATGAAAGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTAAAGCGACTTTATATCTTAAGAACGGATTCGCTGGAGACGATGTc                                                                                                          >  2:792645/1‑151 (MQ=255)
                                                       cGCCGTGAAGGCAGCAAACGATGCCGGATTTCACCTGGCGGTGACAACCATGAAAGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTAAAGCGACTTTATATCTTAAGAACGGATTCGCTGGAGACGATGTCGCGGCTGGTGAGTaa                                                                                           >  2:208653/1‑150 (MQ=255)
                                                        gCCGTGAAGGCAGCAAACGATGCCGGATTTCACCTGGCGGTGACAACCATGAAAGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTAAAGCGACTTTATATCTTAAGAACGGATTCGCTGGAGACGATGTCGCGGCTGGTGAGTAAcc                                                                                         <  2:1063/151‑1 (MQ=255)
                                                         ccGTGAAGGCAGCAAACGATGACGGATTTCGCCTGGCGGTGACAACCATGAAAGGCAACGTAAAACCGGGGGAAAATCCGTAGTTCCTAAACCGACTGTACATCGTAAGATCGGATGCGATATATCCGATGTCGAGGCTGGTGAGTAACCa                                                                                        <  1:563082/151‑1 (MQ=255)
                                                          cGTGAAGGCAGCAAACGATGCCGGATTTCACCTGGCGGTGACAACCATGAAAGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTAAAGCGACTTTATATCTTAAGAACGGATTCGCTGGAGACGATGTCGCGGCTGGTGAGTAAcc                                                                                         >  1:153331/1‑149 (MQ=255)
                                                               aGGCAGCAAACGATGCCGGATTTCACATGGCGGTGACAACCATGAAAGGCAAAGCAAAACCGGGGGATAATCCGTTGGTACTAAAGAGAGTTTATATCTTAAGAACGGATTCGCTGAAGACGATGTCGAGGCTGGTGAGTaa                                                                                           >  1:494064/1‑142 (MQ=255)
                                                                      aaaCGATGCCGGATTTCACCTGGCGGTGACAACCATGAAAGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTAAAGCGACTTTATATCTTAAGAACGGATTCGCTGGAGACGATGTCGCGGCTGGTGAGTAACCAGCCGCAGGGATa                                                                            >  1:22411/1‑150 (MQ=255)
                                                                      aaaCGATGCCGGATTTCACCTGGCGGTGACAACCATGAAAGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTAAAGCGACTTTATATCTTAAGAACGGATTCGCTGGAGACGATGTCGCGGCTGGTGAGTAACCAGCCGCAGGGATa                                                                            >  2:509266/1‑150 (MQ=255)
                                                                        aCGATGCCGGATTTCACCTGGCGGTGACAACCATGAAAGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTAAAGCGACTTTATATCTTAAGAACGGATTCGCTGGAGACGATGTCGCGGCTGGTGAGTAACCAGCCGCAGGGATAAc                                                                          >  1:505840/1‑150 (MQ=255)
                                                                               cGGATTTCACCTGGCGGTGACAACCATGAAAGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTAAAGCGACTTTATATCTTAAGAACGGATTCGCTGGAGACGATGTCGCGGCTGGTGAGTAACCAGCCGCAGGGATAACAATCaaac                                                                  >  1:725836/1‑149 (MQ=255)
                                                                                    ttCACCTGGCGGTGACAACCATGAAAGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTAAAGCGACTTTATATCTTAAGAACGGATTCGCTGGAGACGATGTCGCGGCTGGTGAGTAACCAGCCGCAGGGATAACAATCAAGCAAcc                                                              <  2:845338/150‑1 (MQ=255)
                                                                                          tGGCGGTGACAACCATGAAAGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTAAAGCGACTTTATATCTTAAGAACGGATTCGCTGGAGACGATGTCGCGGCTGGTGAGTAACCAGCCGCAGGGATAACAATCAAGCAACCTGTACCg                                                       >  2:579022/1‑151 (MQ=255)
                                                                                            gCGGTGACAACCATGAAAGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTAAAGCGACTTTATATCTTAAGAACGGATTCGCTGGAGACGATGTCGCGGCTGGTGAGTAACCAGCCGCAGGGATAACAATCAAGCAACCTGTACCgg                                                      >  1:506047/1‑150 (MQ=255)
                                                                                             cGGTGACAACCATGAAAGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTAAAGCGACTTTATATCTTAAGAACGGATTCGCTGGAGACGATGTCGCGGCTGGTGAGTAACCAGCCGCAGGGATAACAATCAAGCAACCTGTACCGGaa                                                    <  1:393801/151‑1 (MQ=255)
                                                                                                   cAACCATGAAAGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTAAAGCGACTTTATATCTTAAGAACGGATTCGCTGGAGACGATGTCGCGGCTGGTGAGTAACCAGCCGCAGGGATAACAATCAAGCAACCTGTACCGGAATCGCtt                                              >  2:435369/1‑151 (MQ=255)
                                                                                                           aaaGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTAAAGCGACTTTATATCTTAAGAACGGATTCGCTGGAGACGATGTCGCGGCTGGTGAGTAACCAGCCGCAGGGATAACAATCAAGCAACCTGTACCGGAATCGCTTTCGCGGTa                                      <  2:343497/151‑1 (MQ=255)
                                                                                                           aaaGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTAAAGCGACTTTATATCTTAAGAACGGATTCGCTGGAGACGATGTCGCGGCTGGTGAGTAACCAGCCGCAGGGATAACAATCAAGCAACCTGTACCGGAATCGCTTTCGCGGTa                                      <  2:733899/151‑1 (MQ=255)
                                                                                                             aGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTAAAGCGACTTTATATCTTAAGAACGGATTCGCTGGAGACGATGTCGCGGCTGGTGAGTAACCAGCCGCAGGGATAACAATCAAGCAACCTGTACCGGAATCGCTTTCGCGGTACg                                    >  2:821771/1‑151 (MQ=255)
                                                                                                             aGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTAAAGCGACTTTATATCTTAAGAACGGATTCGCTGGAGACGATGTCGCGGCTGGTGAGTAACCAGCCGCAGGGATAACAATCAAGCAACCTGTACCGGAATCGCTTTCGCGGTACg                                    >  1:692450/1‑151 (MQ=255)
                                                                                                                 aaaGTAAAACCGGGGGATAATCCGTTGTTACTAAAGCGACTTTATATCTTAAGAACGGATTCGCTGGAGACGATGTCGCGGCTGGTGAGTAACCAGCCGCAGGGATAACAATCAAGCAACCTGTACCGGAATCGCTTTCGCGGTACGtttc                                <  1:288799/151‑1 (MQ=255)
                                                                                                                  aaGTAAAACCGGGGGATAATCCGTTGTTACTAAAGCGACTTTATATCTTAAGAACGGATTCGCTGGAGACGATGTCGCGGCTGGTGAGTAACCAGCCGCAGGGATAACAATCAAGCAACCTGTACCGGAATCGCTTTCGCGGTACGtttca                               <  1:605194/151‑1 (MQ=255)
                                                                                                                    gTAAAACCGGGGGATAATCCGTTGTTACTAAAGCGACTTTATATCTTAAGAACGGATTCGCTGGAGACGATGTCGCGGCTGGTGAGTAACCAGCCGCAGGGATAACAATCAAGCAACCTGTACCGGAATCGCTTTCGCGGTACGtttcatt                             >  2:761485/1‑151 (MQ=255)
                                                                                                                      aaaaCCGGGGGATAATCCGTTGTTACTAAAGCGACTTTATATCTTAAGAACGGATTCGCTGGAGACGATGTCGCGGCTGGTGAGTAACCAGCCGCAGGGATAACAATCAAGCAACCTGTACCGGAATCGCTTTCGCGGTACGtttcatttc                           <  1:761428/151‑1 (MQ=255)
                                                                                                                          ccGGGGTATAATCAGTTAGTACTAAAGCGACTTTATATCTTAAGAACGGATTAGCTGGAGACGATGTCGCGGCTGGTGAGTAACAAGCCGCAGGGATAACAATCAAGCAACCTGTACCGGAATCGCTTTCGCGGAACGTTtcatttcatt                        >  1:609918/1‑150 (MQ=255)
                                                                                                                          ccGGGGGATAATCCGTTGTTACTAAAGCGACTTTATATCTTAAGAACGGATTCGCTGGAGACGATGTCGCGGCTGGTGAGTAACCAGCCGCAGGGATAACAATCAAGCAACCTGTACCGGAATCGCTTTCGCGGTACGTTTCATTTCATTg                       <  1:299971/151‑1 (MQ=255)
                                                                                                                                 aTAATCCGTTGTTACTAAAGCGCCTTTATATCTTAAGAACGGATTCGCTGGAGACGATGTCGCGGCTGGTGAGTAACCAGCCGCAGGGATAACAATCAAGCAACCTGTACCGGAATCGCTTTCGCGGTACGTTTCATTTCATTGTCGCCtt                >  2:562796/1‑151 (MQ=255)
                                                                                                                                   aaTCCGTTGTTACTAAAGCGACTTTATATCTTAAGAACGGATTCGCTGGAGACGATGTCGCGGCTGGTGAGTAACCAGCCGCAGGGATAACAATCAAGCAACCTGTACCGGAATCGCTTTCGCGGTACGTTTCATTTCATTGTCGCCTTCa              <  2:375026/151‑1 (MQ=255)
                                                                                                                                          tgttACTAAAGCGACTTTATATCTTAAGAACGGATTCGCTGGAGACGATGTCGCGGCTGGTGAGTAACCAGCCGCAGGGATAACAATCAAGCAACCTGTACCGGAATCGCTTTCGCGGTACGTTTCATTTCATTGTCGCCTTCAAAATAg        <  1:56347/150‑1 (MQ=255)
                                                                                                                                             tACTAAAGCGACTTTATATCTTAAGAACGGATTCGCTGGAGACGATGTCGCGGCTGGTGAGTAACCAGCCGCAGGGATAACAATCAAGCAACCTGTACCGGAATCGCTTTCGCGGTACGTTTCATTTCATTGTCGCCTTCAAAATAGGCg     <  2:375051/150‑1 (MQ=255)
                                                                                                                                             tACTAAAGCGACTTTATATCTTAAGAACGGATTCGCTGGAGACGATGTCGCGGCTGGTGAGTAACCAGCCGCAGGGATAACAATCAAGCAACCTGTACCGGAATCGCTTTCGCGGTACGTTTCATTTCATTGTCGCCTTCAAAATAGGCg     >  1:204683/1‑150 (MQ=255)
                                                                                                                                               cTAAAGCGACTTTATATCTTAAGAACGGATTCGCTGGAGACGATGTCGCGGCTGGTGAGTAACCAGCCGCAGGGATAACAATCAAGCAACCTGTACCGGAATCGCTTTCGCGGTACGTTTCATTTCATTGTCGCCTTCAAAATAGGCGACg  >  2:543473/1‑151 (MQ=255)
                                                                                                                                               cTAAAGCGACTTTATATCTTAAGAACGGATTCGCTGGAGACGATGTCGCGGCTGGTGAGTAACCAGCCGCAGGGATAACAATCAAGCAACCTGTACCGGAATCGCTTTCGCGGTACGTTTCATTTCATTGTCGCCTTCAAAATAGGCGACg  >  2:150959/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AATTTAATCCGCATGTCTGGTATCTTTCGTATCCGTTTGGCGGATTTAATGACAACGCCGTGAAGGCAGCAAACGATGCCGGATTTCACCTGGCGGTGACAACCATGAAAGGCAAAGTAAAACCGGGGGATAATCCGTTGTTACTAAAACGACTTTATATCTTAAGAACGGATTCGCTGGAGACGATGTCGCGGCTGGTGAGTAACCAGCCGCAGGGATAACAATCAAGCAACCTGTACCGGAATCGCTTTCGCGGTACGTTTCATTTCATTGTCGCCTTCAAAATAGGCGACG  >  NC_000913/146090‑146383

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