Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,104,698	C→T	7.2%	E62E (GAG→GAA)	PP_1879 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,104,698	0	C	T	7.2%	162.6 / 2.2	55	E62E (GAG→GAA)	PP_1879	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (29/22);  new base T (2/2);  total (31/24)
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.70e-01

AAGCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCG  >  NC_002947/2104554‑2104844                                                                                                                                                 |                                                                                                                                                   aaGCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCa                                                                                                                                                >  8:25011/1‑149 (MQ=255)        tCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAGTTCCACCCgtgg                                                                                                                                         <  3:219607/149‑4 (MQ=255)           aCCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAgggg                                                                                                                                      >  5:184931/1‑149 (MQ=255)                  ccAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAg                                                                                                                               >  2:67309/1‑149 (MQ=255)                   cAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAgg                                                                                                                              >  1:46544/1‑149 (MQ=255)                   cAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAgg                                                                                                                              <  7:217614/149‑1 (MQ=255)                          aaTAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGccc                                                                                                                        <  6:184931/148‑1 (MQ=255)                              gCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGTCCACGGGTGGAGTTCCACCCgt                                                                                                                                           <  8:218756/125‑2 (MQ=255)                              gCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGTCCACGGGTGGAGTTCCACCCgt                                                                                                                                           >  7:218756/1‑124 (MQ=255)                                  ttCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAgg                                                                                                                                        >  8:211557/1‑124 (MQ=255)                                  ttCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAgg                                                                                                                                        <  7:211557/124‑1 (MQ=255)                                                   tGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGTAATAGGTTGCCCGTAACTCGACCGTGTTGGGTTTccc                                                                                               <  4:73712/148‑1 (MQ=255)                                                    gCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTccc                                                                                               <  8:149500/147‑1 (MQ=255)                                                    gCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTccc                                                                                               <  6:227397/147‑1 (MQ=255)                                                      tGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACACGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCt                                                                                           >  7:276060/1‑149 (MQ=255)                                                        aaGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTc                                                                                         <  1:67309/149‑1 (MQ=255)                                                              ttcttGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCACGGAACTCGACCGTGTTGGGTTTCACCCCTTCAACATc                                                                                   >  7:64563/1‑149 (MQ=255)                                                               tcttGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCg                                                                                  >  3:97443/1‑149 (MQ=255)                                                                cttGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGt                                                                                 >  2:263912/1‑149 (MQ=255)                                                                          aGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCt                                                                       <  1:314349/149‑1 (MQ=255)                                                                           gTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCACGGAACTCGACCGTGTTGGGTTTCACCCATTCAACATCGTCCACCACGc                                                                        >  6:91247/1‑147 (MQ=255)                                                                              cGTCATTCCAGTAGATGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTTCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGt                                                                   <  5:91247/149‑1 (MQ=255)                                                                              cGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTgg                                                                    <  4:97443/148‑1 (MQ=255)                                                                                  aTTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTg                                                                     >  6:190158/1‑143 (MQ=255)                                                                                  aTTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTg                                                                     <  5:190158/143‑1 (MQ=255)                                                                                    tCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGtt                                                             >  1:295080/1‑149 (MQ=255)                                                                                         tAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTggggg                                                        >  2:28210/1‑149 (MQ=255)                                                                                         tAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTggggg                                                        <  2:46544/149‑1 (MQ=255)                                                                                           gagGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGa                                                      <  3:259712/149‑1 (MQ=255)                                                                                                  aGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAaga                                                 >  8:129960/1‑147 (MQ=255)                                                                                                   ggTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAagagag                                              <  8:64563/149‑1 (MQ=255)                                                                                                          cTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGAGCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGt                                       >  6:246564/1‑149 (MQ=255)                                                                                                          cTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAAATCGACAGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGt                                       >  5:266346/1‑149 (MQ=255)                                                                                                           tCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGAc                                                     >  1:142833/1‑134 (MQ=255)                                                                                                           tCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGAc                                                     <  2:142833/134‑1 (MQ=255)                                                                                                           tCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTc                                      >  7:190029/1‑149 (MQ=255)                                                                                                               cTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGt                                  >  4:209375/1‑149 (MQ=255)                                                                                                                aaTTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCACAGGGGGGAATAGGTTGCCCTGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGtt                                 <  3:266409/148‑1 (MQ=255)                                                                                                                 aTTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGtt                                 >  3:233960/1‑148 (MQ=255)                                                                                                                  tttCGGGTTGGGAATTTGCCCACGGGTGGAACTCCAGACGGAGGGGGGAATAGGTGGCCCGGAACTCGACCGTGTGGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTTGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGtt                                 <  8:276060/147‑1 (MQ=255)                                                                                                                         gTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTg                        <  2:295080/149‑1 (MQ=255)                                                                                                                          tGCGAATTTGCACACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTg                        >  8:96224/1‑148 (MQ=255)                                                                                                                             gAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGTAATAGGGTGCCCGGAACTCGACCGTGTTGTGTTTCCCCCCTTCAACATCGTCCACCACGCTCGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAg                              <  5:57955/139‑1 (MQ=255)                                                                                                                             gAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGt                    >  7:162591/1‑149 (MQ=255)                                                                                                                             gAATTTGCCCACGGGTGGAACTCAACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTCGGTACGGTTGGTGGGACAATAGAGGACCAGTCCACTTCAg                              >  6:57955/1‑139 (MQ=255)                                                                                                                                   gCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGc               <  7:25011/148‑1 (MQ=255)                                                                                                                                   gCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGc               <  7:129960/148‑1 (MQ=255)                                                                                                                                   gCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCg              >  6:104643/1‑149 (MQ=255)                                                                                                                                     ccACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGt            >  3:164806/1‑149 (MQ=255)                                                                                                                                      tgcGGGTGGAGTTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGt            >  7:44525/3‑148 (MQ=255)                                                                                                                                       aCGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCgg             <  6:80128/146‑1 (MQ=255)                                                                                                                                       aCGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGtag          >  5:80128/1‑147 (MQ=255)                                                                                                                                       aCGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTa          >  6:268156/1‑149 (MQ=255)                                                                                                                                       aCGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTAGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGATCAgggct                          >  4:222815/1‑131 (MQ=255)                                                                                                                                               aaCTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAgcg  >  1:310097/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                   AAGCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGCG  >  NC_002947/2104554‑2104844

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: ‑