Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,276,169	A→T	11.8%	E61D (GAA→GAT)	PP_2874 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,276,169	0	A	T	11.8%	76.7 / 4.2	34	E61D (GAA→GAT)	PP_2874	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (19/11);  new base T (2/2);  total (21/13)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.81e-01

TGGGCCGACGCACAGCTGCACAGCCGCTGCCACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTGGGGTGTCGACCGCGCCCGGGCGATGGCAACAGGCGCCTTGCCTCTGCCGGAAACCCGTCTGCAGCGGCGCCAGGCATGAAAAAGCCCGGCCTATG  >  NC_002947/3276023‑3276314                                                                                                                                                   |                                                                                                                                                  tGGGCCGACGCACAGCTGCACAGCCGCTGCCACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAaga                                                                                                                                                 <  1:67253/149‑1 (MQ=255)   ggCCGACGCACAGCTGCACAGCCGCTGCCACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGt                                                                                                                                               <  8:48925/149‑1 (MQ=255)                cTGCACAGCCGCTGCCACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCt                                                                                                                                  <  2:198870/149‑1 (MQ=255)                       gCCGCTGCAACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGc                                                                                                                           <  6:35324/149‑1 (MQ=255)                               cACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCAt                                                                                                                   >  1:37973/1‑149 (MQ=255)                                acgacgAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCAt                                                                                                                   >  8:162525/1‑148 (MQ=255)                                              aGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGgt                                                                                                    >  3:35804/1‑149 (MQ=255)                                              aGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGgt                                                                                                    <  4:35804/149‑1 (MQ=255)                                                cGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGgtgt                                                                                                  <  7:162525/149‑1 (MQ=255)                                                    ggCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGATGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTgggg                                                                                              >  4:244624/1‑149 (MQ=255)                                                    ggCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGATGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTgggg                                                                                              >  4:57102/1‑149 (MQ=255)                                                       accacACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTGGGgtgt                                                                                           >  7:222751/1‑149 (MQ=255)                                                         cacaCTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTGGGGTGTc                                                                                          >  3:141878/1‑148 (MQ=255)                                                         cacaCTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTGGGGTGTCg                                                                                         >  7:82556/1‑149 (MQ=255)                                                           cacTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTGGGGTGTCGAc                                                                                       <  5:159225/149‑1 (MQ=255)                                                             cTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTGGGGTGTCGACcg                                                                                     >  2:245121/1‑149 (MQ=255)                                                                gACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTGGGGTGTCGACCgcgc                                                                                  >  7:132006/1‑149 (MQ=255)                                                                    tGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCTCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTGGGGTGTCGACCGCGCCCgg                                                                              <  4:141878/149‑1 (MQ=255)                                                                     gCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTGGGGTGTCGCCCGGGCGCggt                                                                             >  6:142812/1‑148 (MQ=255)                                                                      cAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTGGGGTGTCGACCGCGCCCggg                                                                             <  1:245121/148‑1 (MQ=255)                                                                          tGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTGGGGTGTCGCCCGGGCGCGGTCGaca                                                                        <  5:142812/149‑3 (MQ=255)                                                                                    gcgcCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTGGGGTGTCGCCCGCGCCCGGGCGATGGCAACAGgcg                                                              >  6:215125/1‑149 (MQ=255)                                                                                    gcgcCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTGGGGTGTCGACCGCGCCCGGGCGATGGCAACAGgcg                                                              >  8:242444/1‑149 (MQ=255)                                                                                                     aGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTGGGGTGTCGCCCGCGCCCGGGCGATGGCAACAGGCGCCTTGCCTCTGCCGGa                                              >  2:47429/1‑148 (MQ=255)                                                                                                      gCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGATGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTGGGGTGTCGACCGCGCCCGGGCGATGGCAACAGGCGCCTTGCCTCTGCCGGaaa                                            <  3:57102/149‑1 (MQ=255)                                                                                                      gCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGATGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTGGGGTGTCGACCGCGCCCGGGCGATGGCAACAGGCGCCTTGCCTCTGCCGGaaa                                            <  3:244624/149‑1 (MQ=255)                                                                                                                gggTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTGGGGTGTCGACCGCGCCCGGGCGATGGCAACAGGCGCCTTGCCTCTGCCGGAAACCCGTCTg                                    >  6:127308/1‑147 (MQ=255)                                                                                                                          ccGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTGGGGTGTCGACCGCGCCCGGGCGATGGCAACAGGCGCCTTGCCTCTGCCGGAAACCCGTCTGCAGCGGCGCCAg                        >  4:162132/1‑149 (MQ=255)                                                                                                                            ggggATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTGGGGTGTCGACCGCGCCCGGGCGATGGCAACAGGCGCCTTGCCTCTGCCGGAAACCCGTCTGCAGCGGCGCCAgg                       >  5:51898/1‑148 (MQ=255)                                                                                                                                aggggcgggcGTGGCGGGCAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTGGGGTGTCGACCGCGCCCGGGCGATGGCAACAGGCGCCTTGCCTCTGCCGGAAACCCGCCTGCAGCGGCGCCAGGCATGa                  <  8:132006/147‑1 (MQ=255)                                                                                                                                     cgggcgTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCACGACCCGTGCGCTGGTGTGGGGTGTCGACCGCGCCCGGGCGATGGCAACAGGCGCCTTGCCTCTGCCGGAAACCCGTCTGCAGCGGCGCCAGGCATGAAAAAg             >  8:81829/1‑149 (MQ=255)                                                                                                                                           ttGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTGGGGTGTCGACCGCGCCCGGGCGATGGCAACAGGCGCCTTGCCTCTGCCGGAAACCCGTCTGCAGCGGCGCCAGGCATGAAAAAGCCCGGc       >  8:80730/1‑149 (MQ=255)                                                                                                                                               tGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTGGGGTGTCGACCGCGCCCGGGCGATGGCAACAGGCGCCTTGCCTCTGCCGGAAACCCGTCTGCAGCGGCGCCAGGCATGAAAAAGCCCGGCCTAt   <  7:80730/149‑1 (MQ=255)                                                                                                                                               tGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTGGGGTGTCGACCGCGCCCGGGCGATGGCAACAGGCGCCTTGCCTCTGCCGGAAACCCGTCTGCAGCGGCGCCAGGCATGAAAAAGCCCGGCCTAt   >  2:125947/1‑149 (MQ=255)                                                                                                                                                ggAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTGGGGTGTCGACCGCGCCCGGGCGATGGCAACAGGCGCCTTGCCTCTGCCGGAAACCCGTCTGCAGCGGCGCCAGGCATGAAAAAGCCCGGCCtgca  >  2:153861/1‑146 (MQ=255)                                                                                                                                                   |                                                                                                                                                  TGGGCCGACGCACAGCTGCACAGCCGCTGCCACGACGAAGCACTCAGCGACGGCACCACACTGGACGTGCAGGTGCGCTTGTCGCGCCTGGGGGCGACGCAGCTGTTTCTGGGGTTGTACGCCGGGGATGGGCGGGCGTTGCTGGAAGAGTATTACCCCGCGCTGCCGGGCGAGACCATGACCCGTGCGCTGGTGTGGGGTGTCGACCGCGCCCGGGCGATGGCAACAGGCGCCTTGCCTCTGCCGGAAACCCGTCTGCAGCGGCGCCAGGCATGAAAAAGCCCGGCCTATG  >  NC_002947/3276023‑3276314

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