Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,632,086	A→C	6.9%	E327A (GAG→GCG)	PP_5590 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,632,086	0	A	C	6.9%	156.4 / 2.2	55	E327A (GAG→GCG)	PP_5590	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (22/29);  new base C (2/2);  total (24/31)
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.96e-01

CGACCCCGGCCACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTCGACGAAGACACCTGGGCTGCC  >  NC_002947/3631939‑3632234                                                                                                                                                    |                                                                                                                                                      cGACCCCGGCCACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAg                                                                                                                                                      <  3:212797/149‑1 (MQ=255)          ccACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAGCGCCTCGGCGttgctg                                                                                                                                             <  1:296647/149‑1 (MQ=255)                        aTGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCAc                                                                                                                                 <  1:35128/147‑1 (MQ=255)                        aTGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCAc                                                                                                                                 >  2:35128/1‑147 (MQ=255)                          gAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTg                                                                                                                              >  4:25059/1‑148 (MQ=255)                                cGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAAccc                                                                                                                         <  2:186811/147‑1 (MQ=255)                                   tGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACg                                                                                                                    >  7:186555/1‑149 (MQ=255)                                         gTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATcc                                                                                                              >  3:73267/1‑149 (MQ=255)                                             gtgtTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTg                                                                                                          >  5:241523/1‑149 (MQ=255)                                                  cAACGGTCGCCTCGACCACTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAg                                                                                                                                                      <  2:346775/100‑1 (MQ=255)                                                  cAACGGTCGCCTCGACCACTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAg                                                                                                                                                      >  1:346775/1‑100 (MQ=255)                                                    aCGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTg                                                                                                  <  6:241523/149‑1 (MQ=255)                                                          gCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCCTCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCgt                                                                                            <  6:4058/149‑1 (MQ=255)                                                          gCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCgt                                                                                            <  3:221657/149‑1 (MQ=255)                                                          gCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCgt                                                                                            <  1:299581/149‑1 (MQ=255)                                                          gCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCCACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCgt                                                                                            <  5:352921/149‑1 (MQ=255)                                                            cTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCgtgt                                                                                          <  3:25059/149‑1 (MQ=255)                                                             tCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCgtgtg                                                                                         >  2:12600/1‑149 (MQ=255)                                                             tCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCgtgtg                                                                                         >  4:268490/1‑149 (MQ=255)                                                               gACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTCG‑GAGCTGCACCGtgtgtg                                                                                       >  5:196535/1‑149 (MQ=255)                                                                aCCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGT‑GAGCTGCACCGTGTgtgg                                                                                      <  8:186555/149‑1 (MQ=255)                                                                  cGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAGCGCCTCGGCGTTGCTgg                                                                                                                                            >  4:9515/1‑94 (MQ=39)                                                                       cGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACATCGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCGTGTGTGGGtggtgg                                                                               >  4:53691/1‑149 (MQ=255)                                                                       cGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCGTGTGTGggtggt                                                                                 <  4:73267/147‑1 (MQ=255)                                                                           tGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCGTGTGTGGGTGGTGGAGGc                                                                           <  1:12600/149‑1 (MQ=255)                                                                                gCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACTAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGc                                                                      <  2:85278/149‑1 (MQ=255)                                                                                  tcgtcggCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCAc                                                                                               <  8:56813/117‑1 (MQ=255)                                                                                    gtggGAAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCCCCGCAAGCAACGCCGAGGCGCTGCTGGCCCCCCATAACCGGAACCCCGACGCCAGCCGCGGG‑CAGCTGCACCGAGTGTGGGTGGAGGAGGCGACGCTga                                                                   >  8:130577/1‑148 (MQ=255)                                                                                       ggCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTgaagaa                                                               <  3:268490/149‑1 (MQ=255)                                                                                          aagaagGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATgg                                                            >  6:39035/1‑149 (MQ=255)                                                                                                 agaTGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCGTGTgtgggt                                                                                    >  8:212324/1‑118 (MQ=255)                                                                                                 agaTGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCGTGTgtgggt                                                                                    <  7:212324/118‑1 (MQ=255)                                                                                                      tACCTCCCCTACAACGCCAACAAGGTGTTCACCGAGAGCAACGCCGATGCGCTGCTGGCCCCCCATCACCTGAACCCCGACCCCATCCGCTGG‑TAGCTGCAACGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCAtg                                                 <  7:130577/148‑1 (MQ=255)                                                                                                      tACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCAtgt                                                <  6:196535/149‑1 (MQ=255)                                                                                                         aTCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGAt                                             <  7:176234/149‑1 (MQ=255)                                                                                                          tCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATg                                            >  7:119075/1‑149 (MQ=255)                                                                                                           ccccTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCgtgtgt                                                                                        <  5:95970/104‑1 (MQ=255)                                                                                                           ccccTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCgtgtgt                                                                                        >  6:95970/1‑104 (MQ=255)                                                                                                                aCAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCcgcg                                       <  7:297759/148‑1 (MQ=255)                                                                                                                aCAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCcgcg                                       >  1:121189/1‑148 (MQ=255)                                                                                                                aCAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCcgcg                                       >  8:297759/1‑148 (MQ=255)                                                                                                                  atcGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAg                                    <  3:53691/147‑1 (MQ=255)                                                                                                                  aaCGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAg                                    >  5:229003/1‑149 (MQ=255)                                                                                                                        aacaaGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCTATCCGCTGG‑GAGCTGCACCGTGTgtgg                                                                                      <  1:9994/93‑1 (MQ=255)                                                                                                                        aacaaGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCGTGTgtgg                                                                                      >  2:9994/1‑93 (MQ=255)                                                                                                                         acaaGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCAc                                 >  7:162022/1‑145 (MQ=255)                                                                                                                         acaaGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCAc                                 <  8:162022/145‑1 (MQ=255)                                                                                                                                     ccGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGACCCCCGACGCCATCCGCTGGATAGCTGCACCGTGTGTGGGTGGTGGAGGCGCCGCTGAAGAATGGCCAGCGCCAGGTGATGCCGCGTTGCACCTACTACCTcgac                    >  1:286427/1‑147 (MQ=255)                                                                                                                                     ccGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTCGACGaa                 <  2:121189/149‑1 (MQ=255)                                                                                                                                     ccGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTCGACGaa                 <  5:39035/149‑1 (MQ=255)                                                                                                                                     ccGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTCGACGaa                 <  2:315794/149‑1 (MQ=255)                                                                                                                                       gCCCGCACCGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTCGACGAAGa               <  2:286427/149‑1 (MQ=255)                                                                                                                                         cagcagCGCCTCGGCGTTGCTGGCGCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTga                                                                   <  3:148816/95‑1 (MQ=38)                                                                                                                                         cagcagCGCCTCGGCGTTGCTGGCGCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑CAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTga                                                                   >  4:148816/1‑95 (MQ=38)                                                                                                                                              aCGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTCGACGAAGACACCTgg        >  8:8894/1‑149 (MQ=255)                                                                                                                                                    aGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTCGACGAAGACACCTGGGCTGcc  <  7:146204/149‑1 (MQ=255)                                                                                                                                                    aGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTCGACGAAGACACCTGGGCTGcc  >  8:146204/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                      CGACCCCGGCCACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCCATCACCTGAACCCCGACGCCATCCGCTGG‑GAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTCGACGAAGACACCTGGGCTGCC  >  NC_002947/3631939‑3632234

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