Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,360,846	(C)5→6	coding (569/1206 nt)	rhmD ←	L‑rhamnonate dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,360,841	1	.	C	97.9%	166.6 / ‑2.7	48	coding (574/1206 nt)	rhmD	L‑rhamnonate dehydratase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (22/25);  total (22/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 = 9.15e-01

AGGGCGCGCAAGCGTGGGCCAGTTTGGTCGCATAGTTCACGTCCTGACTCATCCAGCAGTCGAGCATTAACCAGAAATCCTCACCGCATTTTTCACGCATATCCGCGACCATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGG‑CCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCGGACGCGCACCTGTGGCGTAGAACTGAATCTCATCACGAACAGCGCCGCCTAAAAGTTTATAAACCGGAAGCCCGACCACTTTGCCG  >  NC_000913/2360691‑2360983                                                                                                                                                        |                                                                                                                                                aGGGCGCGCAAGCGTGGGCCAGTTTGGTCGCATAGTTCACGTCCTGACTCATCCAGCAGTCGAGCATTAACCAGAAATCCTCACCGCATTTTTCACGCATATCCGCGACCATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTgg                                                                                                                                                  >  2:234648/1‑151 (MQ=255)  gggCGCGCAAGCGTGGGCCAGTTTGGTCGCATAGTTCACGTCCTGACTCATCCAGCAGTCGAGCATTAACCAGAAATCCTCACCGCATTTTTCACGCATATCCGCGACCATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGG‑c                                                                                                                                                >  2:582610/1‑151 (MQ=255)        gcAAGCGTGGGCCAGTTTGGTCGCATAGTTCACGTCCTGACTCATCCAGCAGTCGAGCATTAACCAGAAATCCTCACCGCATTTTTCACGCATATCCGCGACCATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGG‑ccc‑ccca                                                                                                                                         <  1:883710/151‑3 (MQ=255)             cGTGGGCCAGTTTGGTCGCATAGTTCACGTCCTGACTCATCCAGCAGTCGAGCATTAACCAGAAATCCTCACCGCATTTTTCACGCATATCCGCGACCATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCg                                                                                                                                     >  1:582713/1‑151 (MQ=255)              gTGGGCCAGTTTGGTCGCATAGTTCACGTCCTGACTCATCCAGCAGTCGAGCATTAACCAGAAATCCTCACCGCATTTTTCACGCATATCCGCGACCATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCGt                                                                                                                                    >  1:388607/1‑151 (MQ=255)                     aGTTTGGTCGCATAGTTCACGTCCTGACTCATCCAGCAGTCGAGCATTAACCAGAAATCCTCACCGCATTTTTCACGCATATCCGCGACCATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCAtt                                                                                                                             <  1:84339/151‑1 (MQ=255)                      gTTTGGTCGCATAGTTCACGTCCTGACTCATCCAGCAGTCGAGCATTAACCAGAAATCCTCACCGCATTTTTCACGCATATCCGCGACCATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCAttt                                                                                                                            <  1:101929/151‑1 (MQ=255)                          ggTCGCATAGTTCACGTCCTGACTCATCCAGCAGTCGAGCATTAACCAGAAATCCTCACCGCATTTTTCACGCATATCCGCGACCATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGcc                                                                                                                        >  1:354788/1‑151 (MQ=255)                             cGCATAGTTCACGTCCTGACTCATCCAGCAGTCGAGCATTAACCAGAAATCCTCACCGCATTTTTCACGCATATCCGCGACCATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCAcc                                                                                                                     <  2:482378/151‑1 (MQ=255)                                   gTTCACGTCCTGACTCATCCAGCAGTCGAGCATTAACCAGAAATCCTCACCGCATTTTTCACGCATATCCGCGACCATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATaaa                                                                                                               <  2:364146/151‑1 (MQ=255)                                     tCACGTCCTGACTCATCCAGCAGTCGAGCATTAACCAGAAATCCTCACCGCATTTTTCACGCATATCCGCGACCATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGc                                                                                                             >  1:319645/1‑151 (MQ=255)                                              gACTCATCCAGCAGTCGAGCATTAACCAGAAATCCTCACCGCATTTTTCACGCATATCCGCGACCATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCAtctc                                                                                                      <  1:232088/149‑1 (MQ=255)                                                  cATCCAGCAGTCGAGCATTAACCAGAAATCCTCACCGCATTTTTCACGCATATCCGCGACCATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCtttt                                                                                                  <  1:600593/149‑1 (MQ=255)                                                  cATCCAGCAGTCGAGCATTAACCAGAAATCCTCACCGCATTTTTCACGCATATCCGCGACCATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTg                                                                                                 <  2:804382/150‑1 (MQ=255)                                                    tCCAGCAGTCGAGCATTAACCAGAAATCCTCACCGCATTTTTCACGCATATCCGCGACCATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGcc                                                                                               >  2:242213/1‑150 (MQ=255)                                                        gcagTCGAGCATTAACCAGAAATCCTCACCGCATTTTTCACGCATATCCGCGACCATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATc                                                                                          <  2:248351/151‑1 (MQ=255)                                                            tCGAGCATTAACCAGAAATCCTCACCGCATTTTTCACGCATATCCGCGACCATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCgg                                                                                       >  1:46525/1‑150 (MQ=255)                                                              gAGCATTAACCAGAAATCCTCACCGCATTTTTCACGCATATCCGCGACCATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCGGAcg                                                                                    <  2:149871/151‑1 (MQ=255)                                                                 cATTAACCAGAAATCCTCACCGCATTTTTCACGCATATCCGCGACCATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCGGACgcgc                                                                                 >  2:526670/1‑151 (MQ=255)                                                                  aTTAACCAGAAATCCTCACCGCATTTTTCACGCATATCCGCGACCATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCGGACgcgc                                                                                 >  2:884870/1‑150 (MQ=255)                                                                    tAACCAGAAATCCTCACCGCATTTTTCACGCATATCCGCGACCATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCGGACGCGCAcc                                                                              >  2:287564/1‑151 (MQ=255)                                                                           aaaTCCTCACCGCATTTTTCACGCATATCCGCGACCATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCGGAcg                                                                                    >  2:803965/1‑138 (MQ=255)                                                                           aaaTCCTCACCGCATTTTTCACGCATATCCGCGACCATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCGGAcg                                                                                    <  1:804030/138‑1 (MQ=255)                                                                                        aTTTTTCACGCATATCCGCGACCATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGG‑CCCACCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCGGACGCGCACCTGTGGCGTAGAACTGAAtct                                                          <  2:310904/151‑1 (MQ=255)                                                                                             tCACGCATATCCGCGACCATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCGGACGCGCACCTGTGGCGTAGAACTGAATCtcatc                                                      >  2:823483/1‑150 (MQ=255)                                                                                              cACGCATATCCGCGACCATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCGGACGCGCACCTGTGGCGTAGAACTGAATCtcatca                                                     <  1:810916/150‑1 (MQ=255)                                                                                               aCGCATATCCGCGACCATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCGGACGCGCACCTGTGGCGTAGAACTGAATCTCATCACg                                                   >  2:518237/1‑151 (MQ=255)                                                                                                    tatCCGCGACCATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCGGACGCGCACCTGTGGCGTAGAACTGAATCTCATCACGAAc                                                >  2:804245/1‑149 (MQ=255)                                                                                                     atCCGCGACCATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCGGACGCGCACCTGTGGCGTAGAACTGAATCTCATCACGAACAgc                                             <  1:570703/151‑1 (MQ=255)                                                                                                       ccGCGACCATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCGGACGCGCACCTGTGGCGTAGAACTGAATCTCATCACGAACAgcg                                            <  1:462800/150‑1 (MQ=255)                                                                                                          cgACCATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCGGACGCGCACCTGTGGCGTAGAACTGAATCTCATCACGAACAGcgccgc                                        >  1:877856/1‑151 (MQ=255)                                                                                                           gACCATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCGGACGCGCACCTGTGGCGTAGAACTGAATCTCATCACGAACAGCgccgcc                                       >  1:7315/1‑151 (MQ=255)                                                                                                           gACCATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCGGACGCGCACCTGTGGCGTAGAACTGAATCTCATCACGAACAGCgccgcc                                       <  1:803406/151‑1 (MQ=255)                                                                                                            aCCATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCGGACGCGCACCTGTGGCGTAGAACTGAATCTCATCACGAACAGCGCCGCCt                                      <  2:323613/151‑1 (MQ=255)                                                                                                             ccATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCGGACGCGCACCTGTGGCGTAGAACTGAATCTCATCACGAACAGCGCCGCCTa                                     <  1:526720/151‑1 (MQ=255)                                                                                                                tAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCGGACGCGCACCTGTGGCGTAGAACTGAATCTCATCACGAACAGCGCCGCCTaaa                                   >  1:683769/1‑150 (MQ=255)                                                                                                                   cGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCGGACGCGCACCTGTGGCGTAGAACTGAATCTCATCACGAACAGCGCCGCCTAAAAGtt                               <  2:540326/151‑1 (MQ=255)                                                                                                                      cATCTTTGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCGGACGCGCACCTGTGGCGTAGAACTGAATCTCATCACGAACAGCGCCGCCTAAAAGTTtat                            <  2:377774/151‑1 (MQ=255)                                                                                                                       aTCTTTGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCGGACGCGCACCTGTGGCGTAGAACTGAATCTCATCACGAACAGCGCCGCCTAAAAGTTtat                            >  2:117625/1‑150 (MQ=255)                                                                                                                           ttGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCGGACGCGCACCTGTGGCGTAGAACTGAATCTCATCACGAACAGCGCCGCCTAAAAGTTTATaa                          <  2:328667/148‑1 (MQ=255)                                                                                                                           ttGCGGATCCCCGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCGGACGCGCACCTGTGGCGTAGAACTGAATCTCATCACGAACAGCGCCGCCTAAAAGTTTATAAAc                        <  1:193506/150‑1 (MQ=255)                                                                                                                                   ccccGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCGGACGCGCACCTGTGGCGTAGAACTGAATCTCATCACGAACAGCGCCGCCTAAAAGTTTATAAACCGGAAGccc               >  1:619093/1‑151 (MQ=255)                                                                                                                                   ccccGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCGGACGCGCACCTGTGGCGTAGAACTGAATCTCATCACGAACAGCGCCGCCTAAAAGTTTATAAACCGGAAGccc               <  1:203936/151‑1 (MQ=255)                                                                                                                                   ccccGCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCGGACGCGCACCTGTGGCGTAGAACTGAATCTCATCACGAACAGCGCCGCCTAAAAGTTTATAAACCGGAAGcc                <  2:396509/150‑1 (MQ=255)                                                                                                                                       gCATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCGGACGCGCACCTGTGGCGTAGAACTGAATCTCATCACGAACAGCGCCGCCTAAAAGTTTATAAACCGGAAGCCCGAcc           <  1:804310/151‑1 (MQ=255)                                                                                                                                        cATCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCGGACGCGCACCTGTGGCGTAGAACTGAATCTCATCACGAACAGCGCCGCCTAAAAGTTTATAAACCGGAAGCCCGAcc           >  2:703387/1‑150 (MQ=255)                                                                                                                                         aTCGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCGGACGCGCACCTGTGGCGTAGAACTGAATCTCATCACGAACAGCGCCGCCTAAAAGTTTATAAACCGGAAGCCCGACCAc         >  2:176653/1‑151 (MQ=255)                                                                                                                                           cGCCATCATGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCGGACGCGCACCTGTGGCGTAGAACTGAATCTCATCACGAACAGCGCCGCCTAAAAGTTTATAAACCGGAAGCCCGACCACt        <  2:534742/150‑1 (MQ=255)                                                                                                                                              catcatGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCGGACGCGCACCTGTGGCGTAGAACTGAATCTCATCACGAACAGCGCCGCCTAAAAGTTTATAAACCGGAAGCCCGACCACTTTGc    >  1:55041/1‑151 (MQ=255)                                                                                                                                                 catGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCGGACGCGCACCTGTGGCGTAGAACTGAATCTCATCACGAACAGCGCCGCCTAAAAGTTTATAAACCGGAAGCCCGACCACTTTGCCg  >  2:100486/1‑150 (MQ=255)                                                                                                                                                 catGTGGCCCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCGGACGCGCACCTGTGGCGTAGAACTGAATCTCATCACGAACAGCGCCGCCTAAAAGTTTATAAACCGGAAGCCCGACCACTTTGCCg  <  1:341473/150‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                AGGGCGCGCAAGCGTGGGCCAGTTTGGTCGCATAGTTCACGTCCTGACTCATCCAGCAGTCGAGCATTAACCAGAAATCCTCACCGCATTTTTCACGCATATCCGCGACCATAGCGGCATCTTTGCGGATCCCCGCATCGCCATCATGTGG‑CCC‑CCAGTGCGTCGGCATTTTGCCACCGATAAAGCCCATCTCTTTTGCCAGATCCGGACGCGCACCTGTGGCGTAGAACTGAATCTCATCACGAACAGCGCCGCCTAAAAGTTTATAAACCGGAAGCCCGACCACTTTGCCG  >  NC_000913/2360691‑2360983

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