Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,261,698	T→A	11.8%	D143V (GAC→GTC)	PP_2859 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,261,698	0	T	A	11.8%	100.0 / 4.5	41	D143V (GAC→GTC)	PP_2859	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (18/18);  new base A (3/2);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTTAACACTCTGTCGGGCATGATGGCGCATCGCTATTAAGCAGCTTAGCAACCGATTGCAGGGCCAGCAATTGGCGTTGCGGCCAATTGCCATCAAGTATTTGGTGCGGTTGATTGTGCCGCCTCAGCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGAACAG  >  NC_002947/3261554‑3261846                                                                                                                                                 |                                                                                                                                                     gTTTAACACTCTGTCGGGCATGATGGCGCATCGCTATTAAGCAGCTTAGCAACCGATTGCAGGGCCAGCAATTGGCGTTGCGGCCAATTGCCATCAAGTATTTGGTGCGGTTGATTGTGCCGCCTCAGCCACGCCTGGCTGTCGTCGaa                                                                                                                                                  <  6:50653/149‑1 (MQ=255)  tttAACACTCTGTCGGGCATGATGGCGCATCGCTATTAAGCAGCTTAGCAACCGATTGCAGGGCCAGCAATTGGCGTTGCGGCCAATTGCCATCAAGTATTTGGTGCGGTTGATTGTGCCGCCTCAGCCACGCCTGGCTGTCGTCGaa                                                                                                                                                  >  5:351575/1‑148 (MQ=255)   ttAACACTCTGTCGGGCATGATGGCGCATCGCTATTAAGCAGCTTAGCAACCGATTGCAGGGCCAGCAATTGGCGTTGCGGCCAATTGCCATCAAGTATTTGGTGCGGTTGATTGTGCCGCCTCAGCCACGCCTGGCTGTCGTCGaaa                                                                                                                                                 >  2:12640/1‑148 (MQ=255)        acTCTGTCGGGCATGATGGCGCATCGCTATTAAGCAGCTTAGCAACCGATTGCAGGGCCAGCAATTGGCGTTGCGGCCAATTGCCATCAAGTATTTGGTGCGGTTGATTGTGCCGCCTCAGCCACGCCTGGCTGTCGTCGAAAAAgcgc                                                                                                                                           <  6:351575/149‑1 (MQ=255)            tGTCGGGCATGATGGCGCATCGGTATTAACCTGCTTAGCAACCGATTGCAGGGCCAGCAATTGGCGTTGCGGCCAATTGCCATCAAGTATTTGGTGCGGTTGATTGTGCCGCCTCAGCCACGCGTGGCTGTCGTCGAAAAAGCGCTGAc                                                                                                                                       <  1:35430/149‑1 (MQ=255)                      gatgGCGCATCGCTATTAAGCAGCTTAGCAACCGATTGCAGGGCCAGCAATTGGCGTTGCGGCCAATTGCCATCAAGTATTTGGTGCGGTTGATTGTGCCGCCTCAGCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAg                                                                                                                             >  6:273328/1‑149 (MQ=255)                      gatgGCGCATCGCTATTAAGCAGCTTAGCAACCGATTGCAGGGCCAGCAATTGGCGATGCGGCCAATTGCCATCAAGTATTTGGTGCGGTTGATTGTGCCGCCTCAGCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAg                                                                                                                             >  1:268049/1‑149 (MQ=255)                       atgGCGCATCGCTATTAAGCAGCTTAGCAACCGATTGCAGGGCCAGCAATTGGCGTTGCGGCCAATTGCCATCAAGTATTTGGTGCGGCTGATTGTGCCGCCTCAGCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAg                                                                                                                             <  4:220396/148‑1 (MQ=255)                        tgGCGCATCGCTATTAAGCAGCTTAGCAACCGATTGCAGGGCCAGCAATTGGCGTTGCGGCCAATTGCCATCAAGTATTTGGTGCGGTTGATTGTGCCGCCTCAGCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGt                                                                                                                           >  1:258115/1‑149 (MQ=255)                        tgGCGCATCGCTATTAAGCAGCTTAGCAACCGATTGCAGGGCCAGCAATTGGCGTTGCGGCCAATTGCCATCAAGTATTTGGTGCGGTTGATTGTGCCGCCTCAGCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGt                                                                                                                           >  1:308113/1‑149 (MQ=255)                          gcgcATCGCTATTAAGCAGCTTAGCAACCGATTGCAGGGCCAGCAATTGGCGTTGCGGCCAATTGCCATCAAGTATTTGGTGCGGTTGATTGTGCCGCCTCAGCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTc                                                                                                                          <  1:12640/148‑1 (MQ=255)                            gcATCGCTATTAAGCAGCTTAGCAACCGATTGCAGGGCCAGCAATTGGCGTTGCGGCCAATTGCCATCAAGTATTTGGTGCGGTTGATTGTGCCGCCTCAGCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCAc                                                                                                                                 >  3:202328/1‑139 (MQ=255)                            gcATCGCTATTAAGCAGCTTAGCAACCGATTGCAGGGCCAGCAATTGGCGTTGCGGCCAATTGCCATCAAGTATTTGGTGCGGTTGATTGTGCCGCCTCAGCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCAc                                                                                                                                 <  4:202328/139‑1 (MQ=255)                              aTCGCTATTAAGCAGCTTAGCAACCGATTGCAGGGCCAGCAATTGGCGTTGCGGCCAATTGCCATCAAGTATTTGGTGCGGTTGATTGTGCCGCCTCAGCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGtt                                                                                                                     <  2:258115/149‑1 (MQ=255)                                       aaGCAGCTTAGCAACCGATTGCAGGGCCAGCAATTGGCGTTGCGGCCAATTGCCATCAAGTATTTGGTGCGGTTGATTGTGCCGCCTCAGCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCt                                                                                                            >  4:41685/1‑149 (MQ=255)                                        agcagcTTAGCAACCGATTGCAGGGCCAGCAATTCGCGATGCGTCCAAATCCCATCAAGTATTTGGTGCGGTTTATTGTGCCGCCTCAGCCACGCCTGGCTGTCGACGAAAAAGCGCTGACGATCACTGAGGAACGGTTGGCAACGCTg                                                                                                           >  7:38241/1‑149 (MQ=255)                                             cTTAGCAACCGATTGCAGGGCCAGCAATTGGCGTTGCGGCCAATTGCCATCAAGTATTTGGTGCGGTTGATTGTGCCGCCTCAGCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCAt                                                                                                      >  2:83806/1‑149 (MQ=255)                                             cTTAGCAACCGATTGCAGGGCCAGCAATTGGCGTTGCGGCCAATTGCCATCAAGTATTTGGTGCGGTTGATTGTGCCGCCTCAGCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCAt                                                                                                      >  7:284728/1‑149 (MQ=255)                                                           gCAGGTCCATCAATTGGCGTTGCGTCGAATTGCCATCAAGTATTTGGTGCGGTTTTTTGTGCCGCCTCAGCCTCGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCATCGCTGCCCATCGGCAATCCAGTc                                                                                         <  8:38241/148‑1 (MQ=255)                                                             aGGGCCAGCAATTGGCGTTGCGGCCAATTGCCATCAAGTATTTGGTGCGGTTGATTGTGCCGCCTCAGCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCAc                                                                                      >  8:4304/1‑149 (MQ=255)                                                                 ccAGCAATTGGCGTTCCTGCCAATTGACATCAAGTATTTCGTGCTGTTCATTGTGCCGCCTCAGCCACGCCTGGCTGTAGTCGAAAAAGCGCTCAAGATCACAGAGGTCGGGTTTGCAACGCTGCCAATCGGCAAACCAGTCCACCCCt                                                                                  >  8:178814/1‑149 (MQ=255)                                                                   aGCAATTGGCGTTGCGGCCAATTGCCATCAAGTATTTGGTGCGGTTGATTGTGCCGCCTCAGCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTg                                                                                <  5:253489/149‑1 (MQ=255)                                                                    gCAATTGGCGTTGCGGCCAATTGCCATCAAGTATTTGGTGCGGTTGATTGTGCCGCCTCAGCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGc                                                                               <  2:268049/149‑1 (MQ=255)                                                                    gCAATTGGCGTTGCGGCCAATTGCCATCAAGTATTTGGTGCGGTTGATTGTGCCGCCTCAGCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGc                                                                               <  1:83806/149‑1 (MQ=255)                                                                     cAATTGGCGTTGCGGCCAATTGCCATCAAGTATTTGGTGCGGTTGATTGTGCCGCCTCAGCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCg                                                                              <  7:4304/149‑1 (MQ=255)                                                                      aaTTGGCGTTGCGGCCAATTGCCATCAAGTATTTGGTGCGGTTGATTGTGCCGCCTCAGCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCgg                                                                             <  5:273328/149‑1 (MQ=255)                                                                      aaTTGGCGTTGCGGCCAATTGCCATCAAGTATTTGGTGCGGTTGATTGTGCCGCCTCAGCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCgg                                                                             <  3:41685/149‑1 (MQ=255)                                                                       aTTGGCGTTGCGGCCAATTGCCATCAAGTATTTGGTGCGGTTGATTGTGCCGCCTCAGCCACGCCTGGCTGTCGACGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCggg                                                                            >  6:199754/1‑149 (MQ=255)                                                                       aTTGGCGTTGCGGCCAATTGCCATCAAGTATTTGGTGCGGTTGATTGTGCCGCCTCAGCCACGCCTGGCTGTCGACGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCggg                                                                            >  6:77591/1‑149 (MQ=255)                                                                           gCGTTGCGGCCAATTGCCATCAAGTATTTGGTGCGGTTGATTGTGCCGCCTCAGCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCt                                                                          >  3:297112/1‑147 (MQ=255)                                                                                                     tttGGTGCGGTTGATTGTGCCGCCTCAGCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCg                                                                              >  1:125590/1‑117 (MQ=255)                                                                                                     tttGGTGCGGTTGATTGTGCCGCCTCAGCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCg                                                                              <  2:125590/117‑1 (MQ=255)                                                                                                                      tgCCGCCTCAGCCACGCCTGGCTGTCGACGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTc                             <  5:77591/149‑1 (MQ=255)                                                                                                                      tgCCGCCTCAGCCACGCCTGGCTGTCGACGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTc                             <  5:199754/149‑1 (MQ=255)                                                                                                                        ccgccTCAGCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTcca                           <  8:284728/149‑1 (MQ=255)                                                                                                                        ccgccTCAGCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTc                             >  5:336034/1‑147 (MQ=255)                                                                                                                           ccTCAGCCACGCCTGGCTGTCGTCGTAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTccagcc                        <  7:178814/149‑1 (MQ=255)                                                                                                                             tCAGCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTccagccag                      >  7:96892/1‑149 (MQ=255)                                                                                                                                gCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCaa                                                                    <  1:223338/100‑1 (MQ=255)                                                                                                                                gCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCaa                                                                    >  2:223338/1‑100 (MQ=255)                                                                                                                                 ccACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTg                  >  5:55399/1‑149 (MQ=255)                                                                                                                                              tcgtcgAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCGATCGGCAAGCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGGCTCTTGCAGGCAGGTTGGGCAGTCGGCgga     >  8:155661/1‑147 (MQ=255)                                                                                                                                              tcgtcgAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGAGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGaa     >  5:330196/1‑149 (MQ=255)                                                                                                                                                gtcgAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGaa     >  7:264184/1‑147 (MQ=255)                                                                                                                                                 tcgAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGcc                                  <  2:81504/117‑1 (MQ=255)                                                                                                                                                 tcgAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGAACat  <  8:96892/149‑2 (MQ=255)                                                                                                                                                 tcgAAAAAGCGCTGACGATCACTGAGGTCGCGTTCGCAACGCTGCCCAACGGCAATCCAGTACACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCACCAGTGCCACTTCCAGCCAGGTTGGGCAGTCGCCGAACAg  >  8:15804/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                     GTTTAACACTCTGTCGGGCATGATGGCGCATCGCTATTAAGCAGCTTAGCAACCGATTGCAGGGCCAGCAATTGGCGTTGCGGCCAATTGCCATCAAGTATTTGGTGCGGTTGATTGTGCCGCCTCAGCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGAACAG  >  NC_002947/3261554‑3261846

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