Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,208,661	A→T	7.1%	Q205L (CAG→CTG)	PP_3692 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,208,661	0	A	T	7.1%	168.1 / 3.5	56	Q205L (CAG→CTG)	PP_3692	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (26/26);  new base T (2/2);  total (28/28)
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.89e-01

GGGTCTCGCTCCAAGGACGAGATACGGAGCAACGTATTCGACCCTAACCATACGGACTGCCGGGCTCAATTCAAGGTGTCGATTGACGGCAAATCCGTTGATGTCCACTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAAG‑CGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGAAG  >  NC_002947/4208516‑4208800                                                                                                                                                  |                                                                                                                                             gggTCTCGCTCCAAGGACGAGATACGGAGCAACGTATTCGACCCTAACCATACGGACTGCCGGGCTCAATTCAAGGTGTCGATTGACGGCAAATCCGTTGATGTCCACTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGgc                                                                                                                                           >  5:163522/1‑149 (MQ=255)      tcGCTCCAAGGACGAGATACGGAGCAACGTATTCGACCCTAACCATACGGACTGCCGGGCTCAATTCAAGGTGTCGATTGACGGCAAATCCGTTGATGTCCACTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAg                                                                                                                                      >  2:236469/1‑149 (MQ=255)        gCTCCAAGGACGAGATACGGAGCAACGTATTCGACCCTAACCATACGGACTGCCGGGCTCAATTCAAGGTGTCGATTGACGGCAAATCCGTTGATGTCCACTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAgtg                                                                                                                                    >  1:374562/1‑149 (MQ=255)         cTCCAAGGACGAGATACGGAGCAACGTATTCGACCCTAACCATACGGACTGCCGGGCTCAATTCAAGGTGTCGATTGACGGCAAATCCGTTGATGTCCACTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAgtgt                                                                                                                                   >  5:279758/1‑149 (MQ=255)            cAAGGACGAGATACGGAGCAACGTATTCGACCCTAACCATACGGACTGCCGGGCTCAATTCAAGGTGTCGATTGACGGCAAATCCGTTGATGTCCACTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTc                                                                                                                                 <  2:374562/148‑1 (MQ=255)                gACGAGATACGGAGCAACGTATTCGACCCTAACCATACGGACTGCCGGGCTCAATTCAAGGTGTCGATTGACGGCAAATCCGTTGATGTCCACTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGtt                                                                                                                                  >  8:83154/1‑143 (MQ=255)                gACGAGATACGGAGCAACGTATTCGACCCTAACCATACGGACTGCCGGGCTCAATTCAAGGTGTCGATTGACGGCAAATCCGTTGATGTCCACTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGtt                                                                                                                                  <  7:83154/143‑1 (MQ=255)                gACGAGATACGGAGCAACGTATTCGACCCTAACCATACGGACTGCCGGGCTCAATTCAAGGTGTCGATTGACGGCAAATCCGTTGATGTCCACTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAg                                                                                                                            <  6:252101/149‑1 (MQ=255)                  cGAGATACGGAGCAACGTATTCGACCCTAACCATACGGACTGCCGGGCTCAATTCAAGGTGTCGATTGACGGCAAATCCGTTGATGTCCACTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCTGGCGCaacgcttt                                                                                                                                 <  3:24421/142‑7 (MQ=255)                  cGAGATACGGAGCAACGTATTCGACCCTAACCATACGGACTGCCGGGCTCAATTCAAGGTGTCGATTGACGGCAAATCCGTTGATGTCCACTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCTGGCGCaacgcttt                                                                                                                                 >  4:24421/1‑136 (MQ=255)                          ggAGCAACGTATTCGACCCTAACCATACGGACTGCCGGGCTCAATTCAAGGTGTCGATTGACGGCAAATCCGTTGATGTCCACTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCt                                                                                                                  >  8:209033/1‑149 (MQ=255)                             gCAACGTATTCGACCCTAACCATACGGACTGCCGGGCTCAATTCAAGGTGTCGATTGACGGCAAATCCGTTGATGTCCACTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAg                                                                                                               <  6:163522/149‑1 (MQ=255)                               aaCGTATTCGACCCTAACCATACGGACTGCCGGGCTCAATTCAAGGTGTCGATTGACGGCAAATCCGTTGATGTCCACTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGc                                                                                                              >  8:125459/1‑148 (MQ=255)                                 cGTATTCGACCCTAACCATACGGACTGCCGGGCTCAATTCAAGGTGTCGATTGACGGCAAATCCGTTGATGTCCACTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCaa                                                                                                            >  6:107740/1‑148 (MQ=255)                                       cGACCCTAACCATACGGACTGCCGGGCTCAATTCAAGGTGTCGATTGACGGCAAATCCGTTGATGTCTACTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCTGGCGCaacgcttt                                                                                                                                 >  4:107352/1‑115 (MQ=255)                                       cGACCCTAACCATACGGACTGCCGGGCTCAATTCAAGGTGTCGATTGACGGCAAATCCGTTGATGTCTACTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCTGGCGCaacgcttt                                                                                                                                 <  3:107352/121‑7 (MQ=255)                                        gACCCTAACCATACGGACTGCCGGGCTCAATTCAAGGTGTCGATTGACGGCAAATCCGTTGATGTCCACTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATggg                                                                                                    <  6:279758/149‑1 (MQ=255)                                         aCCCTAACCATACGGACTGCCGGGCTCAATTCAAGGTGTCGATTGACGGCAAATCCGTTGATGTCCACTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATgggg                                                                                                   <  7:125459/149‑1 (MQ=255)                                                ccATACGGACTGCCGGGCTCAATTCAAGGTGTCGATTGACGGCAAATCCGTTGATGTCCACTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTg                                                                                            <  2:352553/149‑1 (MQ=255)                                                  aTACGGACTGCCGGGCTCAATTCAAGGTGTCGATTGACGGCAAATCCGTTGATGTCCACTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTg                                                                                            >  3:126657/1‑147 (MQ=255)                                                       gACTGCCGGGCTCAATTCAAGGTGTCGATTGACGGCAAATCCGTTGATGTCCACTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAAG‑CGTg                                                                                    >  8:94151/1‑149 (MQ=255)                                                             cGGGCTCAATTCAAGGTGTCGATTGACGGCAAATCCGTTGATGTCAACTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAAG‑CGTGCTAt                                                                                >  1:233787/1‑147 (MQ=255)                                                                 ctgaagtCAAGGTGTCGATTGGCGGCAAATCCGTTGATGTCCACTGGGACATCCGAAGCGGCCCGCAGAAAGCGTTGCGCGAGGCGCAAGTGTTCGTGAGTTAAGCCGCTCAGCAAGCGATGGGGCGAACTGAAG‑CGTGGTAGGCGcc                                                                           >  4:244911/7‑148 (MQ=255)                                                                  tCAATTCAAGGTGTCGATTGACGGCAAATCCGTTGATGTCCACTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAg                                                                                                                                             <  2:376270/82‑1 (MQ=255)                                                                  tCAATTCAAGGTGTCGATTGACGGCAAATCCGTTGATGTCCACTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAg                                                                                                                                             >  1:376270/1‑82 (MQ=255)                                                                    aaTTCAAGGTGTCGATTGACGGCAAATCCGTTGATGTCCACTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAAG‑CGTGCTATGAGCCAAg                                                                        <  6:11657/148‑1 (MQ=255)                                                                            gtgtCGATTGACGGCAAATCCGTTGATGTCCACTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAAG‑CGTGCTATTAGCCAAGTCGCATGa                                                                >  8:109205/1‑148 (MQ=255)                                                                                 gATTGACGGCAAATCCGTTGATGTCCACTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATgg                                                                                                     <  3:234110/107‑1 (MQ=255)                                                                                 gATTGACGGCAAATCCGTTGATGTCCACTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATgg                                                                                                     >  4:234110/1‑107 (MQ=255)                                                                                 gATTGACGGCAAATCCGTTGATGTCCACTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAAG‑CGTGCTATGAGCCAAGTCGCATGAAAttt                                                           <  7:209033/148‑1 (MQ=255)                                                                                       cGGCAAATCCGTTGATGTCCACTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAAG‑CGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCgg                                                    <  7:94151/149‑1 (MQ=255)                                                                                                 gTTGATGTCCACTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAAG‑CGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGc                                          >  3:74646/1‑149 (MQ=255)                                                                                                  ttGATGTCCACTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAAGACGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGc                                          <  1:297317/149‑1 (MQ=255)                                                                                                           aCTGGGGCATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAAG‑CGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGAt                                >  6:11635/1‑149 (MQ=255)                                                                                                           aCTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCTATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAAG‑CGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGAt                                <  3:94105/149‑1 (MQ=255)                                                                                                           aCTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAAG‑c                                                                                       <  4:119928/94‑1 (MQ=255)                                                                                                           aCTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAAG‑c                                                                                       >  3:119928/1‑94 (MQ=255)                                                                                                           aCTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAAG‑CGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGAt                                >  2:117906/1‑149 (MQ=255)                                                                                                           aCTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTAGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAAG‑CGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGAt                                <  3:244911/149‑1 (MQ=255)                                                                                                               ggACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAAG‑CGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGcc                             <  1:389486/148‑1 (MQ=255)                                                                                                                 aCATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAAG‑CGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTc                          <  5:107740/149‑1 (MQ=255)                                                                                                                    tCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAAG‑CGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAac                       <  8:50092/149‑1 (MQ=255)                                                                                                                     cccAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAAG‑CGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAaca                      >  2:149806/1‑149 (MQ=255)                                                                                                                      ccAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAAG‑CGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAaca                      <  4:74646/148‑1 (MQ=255)                                                                                                                       cAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAAG‑CGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCg                                                     >  3:398802/1‑116 (MQ=255)                                                                                                                       cAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAAG‑CGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCg                                                     <  4:398802/116‑1 (MQ=255)                                                                                                                       cAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAAG‑CGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAaca                      >  4:332509/1‑147 (MQ=255)                                                                                                                            ggCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAAG‑CGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGtt               <  4:126657/149‑1 (MQ=255)                                                                                                                            ggCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAAG‑CGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGtt               <  1:236469/149‑1 (MQ=255)                                                                                                                                cTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAAG‑CGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACa           <  1:117906/149‑1 (MQ=255)                                                                                                                                   acaAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAAG‑CGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAgg        >  3:304797/1‑149 (MQ=255)                                                                                                                                   acaAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAAG‑CGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAgg        <  1:149806/149‑1 (MQ=255)                                                                                                                                     aaaGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAAG‑CGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGcc                                   <  2:178902/120‑1 (MQ=255)                                                                                                                                     aaaGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAAG‑CGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAgg        >  4:350037/1‑147 (MQ=255)                                                                                                                                       aGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAAG‑CGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACTGAGTTGACAAGGgag     >  3:266297/1‑148 (MQ=255)                                                                                                                                          gTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAAG‑CGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGAAg  >  2:101355/1‑148 (MQ=255)                                                                                                                                                  |                                                                                                                                             GGGTCTCGCTCCAAGGACGAGATACGGAGCAACGTATTCGACCCTAACCATACGGACTGCCGGGCTCAATTCAAGGTGTCGATTGACGGCAAATCCGTTGATGTCCACTGGGACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAAG‑CGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGAAG  >  NC_002947/4208516‑4208800

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