Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	729,698	A→G	H298R (CAC→CGC)	rhsC →	RhsC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	729,698	0	A	G	90.4%	155.6 / 11.7	73	H298R (CAC→CGC)	rhsC	RhsC
Reads supporting (aligned to +/- strand):  ref base A (6/1);  new base G (40/26);  total (46/27)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GGCATCCGTCTGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTGGCGGCGGTGTATGACC  >  NC_000913/729667‑729787                                |                                                                                           ggCATCCGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑ccc                              >  2:769956/1‑93 (MQ=14) ggCATCCGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑ccc                              <  2:3870388/93‑1 (MQ=14) ggCATCCGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑ccc                              >  2:1838634/1‑93 (MQ=14) ggCATCCGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑ccc                              >  1:5340499/1‑93 (MQ=14)  gCATCCGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑cccc                             >  2:5189826/1‑93 (MQ=17)  gCATCCGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑cccc                             >  2:4256475/1‑93 (MQ=17)  gCATCCGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑cccc                             >  2:1552476/1‑93 (MQ=17)  gCATCCGTCTGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑cccc                             >  2:535133/1‑93 (MQ=32)   cATCCGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGAGG‑GCCcg                            >  1:3546213/1‑93 (MQ=255)    aTCCGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCcgc                           >  2:1045460/1‑93 (MQ=18)    aTCCGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCcgc                           <  2:169047/93‑1 (MQ=18)    aTCCGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCcgc                           <  1:6093501/93‑1 (MQ=18)     tCCGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCcgcg                          >  2:5826458/1‑93 (MQ=18)     tCCGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCcgcg                          >  1:3061001/1‑93 (MQ=18)     tCCGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCcgcg                          >  2:1975887/1‑93 (MQ=18)     tCCGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCcgcg                          >  1:1916296/1‑93 (MQ=18)     tCCGTCTGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCcgcg                          >  1:4600432/1‑93 (MQ=33)       cGTCTGTCTGCCGTGTGGCTGCCGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCgcggc                        <  2:3313707/93‑1 (MQ=18)       cGTCTGTCTGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCgcggc                        >  2:2726025/1‑93 (MQ=18)        gtctgtctGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGGCG‑CCCCgcggcg                       >  2:5701262/1‑93 (MQ=255)        gtctgtctGCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGGCG‑CCCCgcggca                       >  2:1798169/1‑92 (MQ=255)        gtctgtctGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCgcggcg                       >  1:1348442/1‑93 (MQ=33)        gtctgtctGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCgcggcg                       >  1:3149377/1‑93 (MQ=33)        gtctgtctGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGAAG‑GCCCgcggcg                       >  1:2349130/1‑93 (MQ=16)         tctgtctgCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGa                      <  2:3135409/93‑1 (MQ=21)          ctgtctgCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGaa                     <  1:4267020/93‑1 (MQ=21)           tgtctgCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAAc                    <  1:3916737/93‑1 (MQ=18)           tgtctgCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGGGCGCTATGGCTGGGCG‑CCCCcccgcggaa                    >  1:1383298/1‑90 (MQ=11)            gtctgCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACt                   <  1:5364074/93‑1 (MQ=21)            gtctgCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACt                   >  1:3206867/1‑93 (MQ=18)             tctgCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTg                  >  1:1065850/1‑93 (MQ=18)              ctgCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGGCG‑CCCCCCGGCGAACttg                 >  1:249884/1‑91 (MQ=11)              ctgCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTgg                 >  2:6103009/1‑93 (MQ=18)              ctgCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTgg                 >  1:5059267/1‑93 (MQ=18)               tgCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTggc                <  1:1538657/93‑1 (MQ=21)               tgCCGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTggc                <  1:2372438/93‑1 (MQ=21)               tgCCGTGTGGCTGACGCACGACCCGGAATACCCGGCGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑GCCCCCCGCGAACGggc                >  1:3251028/1‑93 (MQ=11)                 ccGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTggcgg              <  1:2984458/93‑1 (MQ=21)                  cGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGGCG‑CCCCCCGGCGAACcggcggg             >  1:1472344/1‑92 (MQ=11)                  cGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACT‑CCCGGGGGCGAACTggcggg             >  2:5323369/1‑92 (MQ=11)                  cGTGTGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTggcggc             <  2:1290368/93‑1 (MQ=25)                   gtgtGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑GCCCGCCGCGGACCTgcggcg            >  2:3567353/1‑93 (MQ=11)                   gtgtGGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTGgcggcg            >  1:899545/1‑93 (MQ=25)                    tgtgGCTGCCGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTGGcggcgg           <  1:2986515/93‑1 (MQ=255)                    tgtgGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTGGcggcgg           >  2:4101648/1‑93 (MQ=255)                    tgtgGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTGGcggcgg           <  1:5828216/93‑1 (MQ=255)                    tgtgGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTGGcggcgg           <  1:5846031/93‑1 (MQ=255)                     gtgGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTGGcggcggg          >  1:2678695/1‑92 (MQ=255)                     gtgGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTGGCGGCGgt          <  2:575230/93‑1 (MQ=255)                     gtgGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTGGCGGCGgt          <  2:6177442/93‑1 (MQ=255)                      tgGCTGACGCGCGACCCGGAATCCCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTGGCGGCGgtg         <  1:3368112/93‑1 (MQ=255)                      tgGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTGGCGGCGgtg         <  2:4570919/93‑1 (MQ=37)                      tgGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTGGCGGCGgtg         <  2:2743402/93‑1 (MQ=255)                      tgGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTGGCGGCGgtg         <  2:339531/93‑1 (MQ=37)                      tgGCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCCCCGCGAACTGGcggcgggg         >  1:3660907/1‑91 (MQ=11)                       ggCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTGGcggcggggg        >  2:1451817/1‑90 (MQ=255)                       ggCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTGGCGGCGgtgt        <  2:3483621/93‑1 (MQ=255)                       ggCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTGGCGGCGgtgt        >  1:6173668/1‑93 (MQ=255)                       ggCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCCCGGCCAACCGGCGGCGgcgt        >  2:2818190/1‑93 (MQ=11)                        gCTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTGGCGGCGGTGTa       <  1:78154/93‑1 (MQ=37)                         cTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTGGCGGCGGTGTAt      <  1:5740258/93‑1 (MQ=255)                         cTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTGGCGGCGGTGTAt      <  2:2383867/93‑1 (MQ=255)                         cTGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTGGCGGCGGTGTAt      <  1:5425153/93‑1 (MQ=255)                          tGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGGGCTATGGCTGGGCG‑CCCCGCGGCGCACTGGCGGCGGGGTATg     >  2:6153286/1‑93 (MQ=11)                          tGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTGGCGGCGGTGTATg     <  2:5897440/93‑1 (MQ=37)                          tGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTGGCGGCGGTGTATg     >  1:2154606/1‑93 (MQ=255)                          tGACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTGGCGGCGGTGTATg     >  2:6153287/1‑93 (MQ=255)                           gACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGTTG‑CCCCGCGGCGAACGGGCGGGGGGGTagtg    >  2:491598/1‑90 (MQ=11)                           gACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTGGCGGCGGTGTATGa    >  2:5886707/1‑93 (MQ=255)                           gACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTGGCGGCGGTGTATGa    >  2:4060360/1‑93 (MQ=255)                           gACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTGGCGGCGGTGTATGa    >  1:1883798/1‑93 (MQ=255)                           gACGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTGGCGGCGGTGTATGa    >  2:2130427/1‑93 (MQ=37)                           gACGCGCGACCCGAAATACCCGGAGAATTTACCTGCCGCGCCGCTCGAGCGCTATGGCTGGACGACCCCCCGGGGAA‑TGGCGGCGGGGTATGa    >  1:300911/1‑93 (MQ=11)                            aCGCGCGACCCGGAATACCCGGAGAATTTGCCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTGGCGGCGGTGTATGAc   <  2:1845094/93‑1 (MQ=255)                            aCGCGCGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTGGCGGCGGTGTATGAc   <  1:4898286/93‑1 (MQ=255)                             cgcgcgACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTGGCGGCGGTGTATGAcc  <  1:4442691/93‑1 (MQ=255)                             cGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTGGCGGCGGTGTATGAcc  <  2:2457496/93‑1 (MQ=38)                                |                                                                                           GGCATCCGTCTGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACG‑CCCCGCGGCGAACTGGCGGCGGTGTATGACC  >  NC_000913/729667‑729787

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