Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,726,241	T→C	7.7%	S139G (AGT→GGT)	PP_2387 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,726,241	0	T	C	7.7%	147.3 / 2.8	52	S139G (AGT→GGT)	PP_2387	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (24/24);  new base C (2/2);  total (26/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.50e-01

ACAGCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAACATCATCAATTCCTGGGTGGCGTTCAGCCCGCAGCCGGCCAGGCTCAGCCAGACTTCGCCCGTGAGCATCGACAGGTGCGAGGCCACTGT  >  NC_002947/2726094‑2726382                                                                                                                                                    |                                                                                                                                               acaGCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACt                                                                                                                                                >  6:287708/1‑148 (MQ=255)    gCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAg                                                                                                                                             >  6:131862/1‑148 (MQ=255)                gaGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATgg                                                                                                                               <  4:114670/149‑1 (MQ=255)                  ggtggtGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGcc                                                                                                                             <  5:287708/149‑1 (MQ=255)                    tggtgCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGcc                                                                                                                             >  8:48261/1‑147 (MQ=255)                       tgCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACCGAGCACCACCA‑CATGGCCGAGtt                                                                                                                        >  1:275968/1‑149 (MQ=255)                                   cATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCgg                                                                                                             <  7:201475/148‑1 (MQ=255)                                   cATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCAACTCGGCCatgtggtg                                                                                                            <  5:131862/149‑8 (MQ=255)                                    aTGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTg                                                                                                           >  4:86879/1‑149 (MQ=255)                                         ttCGCTAAAACTCTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCtggt                                                                                                      >  4:126826/1‑149 (MQ=255)                                            gCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTggg                                                                                                   >  2:67431/1‑149 (MQ=255)                                             cTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGt                                                                                                  >  4:242634/1‑149 (MQ=255)                                              tAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACCGAGc                                                                                                                                            <  8:25552/107‑1 (MQ=255)                                              tAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACCGAGc                                                                                                                                            >  7:25552/1‑107 (MQ=255)                                              tAAAACAGTTGAGAAAGGCCTGGTACGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTAAACACACGGCTTCGTGGGACAGCCTGTTGAGGCCACTGAGCACCACCA‑CATCGCCGAGTTGGCCAAGGCGGTGCTGGTGGGt                                                                                                  >  7:73447/1‑148 (MQ=255)                                                    tGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGaa                                                                                            >  8:257835/1‑148 (MQ=255)                                                      ttGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAAcat                                                                                         >  6:184418/1‑149 (MQ=255)                                                      ttGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAAc                                                                                           <  1:167183/147‑1 (MQ=255)                                                        gagaAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGGTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGTCCACTGAGCACCACCA‑CAAGGCCGAGTTGGAAAAGGCGGTTCTGGTGGGTCAGGAAcatac                                                                                       >  8:280659/1‑147 (MQ=255)                                                        cgccaaGGCCTACCAGGCCTTGGCGAGCTCCGGCTCGTTGGTACCTGCGAGGACTTCACACAGGGCATCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGCTGCTGGTGGGTCAGGAAcatca                                                                                       >  1:127718/5‑149 (MQ=255)                                                          gaAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCACCATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAAcatcat                                                                                      <  1:67431/149‑1 (MQ=255)                                                           caaGGCCTACCAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAACatcatc                                                                                     >  8:10343/2‑148 (MQ=255)                                                              ggCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAACATCATCAAtt                                                                                 >  4:54350/1‑149 (MQ=255)                                                                 cTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAACATCATCAATTCCt                                                                              >  5:150058/1‑149 (MQ=255)                                                                    gTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAACATCATCAATTCCTggg                                                                           <  7:48261/149‑1 (MQ=255)                                                                    gTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACCGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAACATCATCAATTCCTgg                                                                            <  2:275968/148‑1 (MQ=255)                                                                     tAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAACATCATCAATTCCTGGGt                                                                          <  3:242634/149‑1 (MQ=255)                                                                      aGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAACATCATCAATTCCTgg                                                                            <  4:129124/146‑1 (MQ=255)                                                                      aGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAACATCATCAATTCCTgg                                                                            >  3:129124/1‑146 (MQ=255)                                                                         ccTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAACATCATCAATTCCTGGGTGGCg                                                                      >  6:133408/1‑149 (MQ=255)                                                                             ggCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAACATCATCAATTCCTGGGTGGCGtt                                                                    <  5:88370/147‑1 (MQ=255)                                                                                 tgctggGGGTGGTTGGTACCTGCGAGTAGTTCACAAAGGGCTTCGTGGGCCAGTCTGTTGATGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAACATCATCAATTCCTGGGTGGCGTTCAGccc                                                              <  8:73447/149‑1 (MQ=255)                                                                                       gggTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAACATCATCAATTCCTGGGTGGCGTTCAGCCCGCAGcc                                                        <  6:150058/149‑1 (MQ=255)                                                                                         gTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAACATCATCAATTCCTGGGTGGCGTTCAGCCCGCAGCCgg                                                      <  7:10343/149‑1 (MQ=255)                                                                                            gTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAACATCATCAATTCCTGGGTGGCGTTCAGCCCGCAGCCGGc                                                     >  8:234558/1‑147 (MQ=255)                                                                                            gTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAACATCATCAATTCCTGGGTGGCGTTCAGCCCGCAGCCGGCCa                                                   <  5:184418/149‑1 (MQ=255)                                                                                                gTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAACATCATCAATTCCTGGGTGGCGTTCAGCCCGCAGCCGGCCAGGCt                                               >  2:100666/1‑149 (MQ=255)                                                                                                          ggACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAACATCATCAATTc                                                                                >  8:287732/1‑106 (MQ=255)                                                                                                          ggACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAACATCATCAATTc                                                                                <  7:287732/106‑1 (MQ=255)                                                                                                                 acacaGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAACATCATCAATTCCTGGGTGGCGTTCAGCCCGCAGCCGGCCAGGCTCAGCCAGACTTCGccc                               <  7:257835/148‑1 (MQ=255)                                                                                                                  cacaGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAACATCATCAATTCCTGGGTGGCGTTCAGCCCGCAGCCGGCCAGGCTCAGCCAGACTTCGCCCGt                             >  3:28064/1‑149 (MQ=255)                                                                                                                   acaGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAACATCATCAATTCCTGGGTGGCGTTCAGCCCGCAGCCGGCCAGGCTCAGCCAGACTTCGCCCGt                             >  4:117071/1‑148 (MQ=255)                                                                                                                    caGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAACATCATCAATTCCTGGGTGGCGTTCAGCCCGCAGCCGGCCAGGCTCAGCCAGACTTCGCCCGTGa                           <  7:234558/149‑1 (MQ=255)                                                                                                                     aGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAACATCATCAATTCCTGGGTGGCGTTCAGCCCGCAGCCGGCCAGGCTCAGCCAGACTTCGCCCGTGAg                          <  3:86879/149‑1 (MQ=255)                                                                                                                     aGGGCTTCGGGGGCCACGCTGGTGAGGCCACTGAGCACCACCT‑CATGGCCGAGGTCGCCAAGGCGGTGCTGGTGGGTCAGGAACATCATCAATTCCTGGGTGGCGTTCAGCCCGCAGCCGGCCAGGCTCAGCCAGACTTCGCCCGTg                            <  7:287527/147‑1 (MQ=255)                                                                                                                         cTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAACATCATCAATTCCTGGGTGGCGTTCAGCCCGCAGCCGGCCAGGCTCAGCCAGACTTCGCCCGTGAGCATc                      <  1:100666/149‑1 (MQ=255)                                                                                                                            cGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAACATCATCAATTCCTGGGTGGCGTTCAGCCCGCAGCCGGCCAGGCTCAGCCAGACTTCGCCCGTGAGCATCGa                    >  5:159929/1‑148 (MQ=255)                                                                                                                               gggCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAACATCATCAATTCCTGGGTGGCGTTCAGCCCGCAGCCGGCCAGGCTCAGCCAGACTTCGCCCGTGAGCATCGACa                  <  3:126826/147‑1 (MQ=255)                                                                                                                                  ccAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAACATCATCAATTCCTGGGTGGCGTTCAGCCCGCAGCCGGCCAGGCTCAGCCAGACTTCGCCCGTGAGCATCGACAGGTGc             >  5:203388/1‑149 (MQ=255)                                                                                                                                    aGGCTGTTGCGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAACATCATCAATTCCTGGGTGGCGTTCAGCCCGCAGCCGGCCAGGCTCAGCCAGACTTCGCCCGTGAGCATCGACAGGTGCGa           <  2:186126/149‑1 (MQ=255)                                                                                                                                             aGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAACATCATCAATTCCTGGGTGGCGTTCAGCCCGCAGCCGGCCAGGCTCAGCCAGACTTCGCCCGTGAGCATCGACAGGTGCGAGGCCACtgt  <  3:40149/149‑1 (MQ=255)                                                                                                                                               gCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAACATCATCAATTCCTGGGTGGCGTTCAGCCCGCAGCCGGCCAGGCTCAGCCAGACTTCGCCCGTGAGCATCGACAGGTGCGAGGCCa      <  8:222009/143‑1 (MQ=255)                                                                                                                                               gCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAACATCATCAATTCCTGGGTGGCGTTCAGCCCGCAGCCGGCCAGGCTCAGCCAGACTTCGCCCGTGAGCATCGACAGGTGCGAGGCCa      >  7:222009/1‑143 (MQ=255)                                                                                                                                                    |                                                                                                                                               ACAGCGACTGCCGACGAGGTGGTGCGCAGGGCTGCATGCCTTCGCTAAAACTGTTGAGAAAGGCCTGGTAGGCCTTGGCGTGCTGCGGGTCGTTGGTACCTGCGAGGACTTCACACAGGGCTTCGTGGGCCAGGCTGTTGAGGCCACTGAGCACCACCA‑CATGGCCGAGTTGGCCAAGGCGGTGCTGGTGGGTCAGGAACATCATCAATTCCTGGGTGGCGTTCAGCCCGCAGCCGGCCAGGCTCAGCCAGACTTCGCCCGTGAGCATCGACAGGTGCGAGGCCACTGT  >  NC_002947/2726094‑2726382

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