Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,327,814	C→T	8.9%	A261V (GCG→GTG)	PP_3797 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,327,814	0	C	T	8.9%	124.0 / 2.3	45	A261V (GCG→GTG)	PP_3797	membrane protein
Reads supporting (aligned to +/- strand):  ref base C (21/20);  new base T (2/2);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.77e-01

CTGGTTCCTGATCCAGGCGATCCTGTCTGACAACCACATCACCATTTCCAGCCGACCCATCGTGCCGGTGATCGCCCGTGAGGAGGTGCCGCAAACGCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGACGAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGACCTCCATCACGCTCCCCGCCACTGCCTACAGCCACGTCACGCTGTCTGGCCCGGGCTGGTACCCGTTGATGTTCCAGAGTGCTTCGGTGAATCCGTACACACGCAAT  >  NC_002947/4327677‑4327962                                                                                                                                          |                                                                                                                                                     cTGGTTCCTGATCCAGGCGATCCTGTCTGACAACCACATCACCATTTCCAGCCGACCCATCGTGCCGGTGATCGCCCGTGAGGAGGTGCCGCAAACGCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGACGAAGCGGCGCGCCTTg                                                                                                                                           >  3:149417/1‑149 (MQ=255)     ttCCTGATCCAGGCGATCCTGTCTGACAACCACATCACCATTTCCAGCCGACCCATCGTGCCGGTGATCGCCCGTGAGGAGGTGCCGCAAACGCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGACGAAGCGGCGCGCCTTGCCgg                                                                                                                                       >  1:16408/1‑149 (MQ=255)         tGATCCAGGCGATCCTGTCTGACAACCACATCACCATTTCCAGCCGACCCATCGTGCCGGTGATCGCCCGTGAGGAGGTGCCGCAAACGCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGACGAAGCGGCGCGCCTTGCCGGTCTg                                                                                                                                   >  4:175482/1‑149 (MQ=255)          gATCCAGGCGATCCTGTCTGACAACCACATCACCATTTCCAGCCGACCCATCGTGCCGGTGATCGCCCGTGAGGAGGTGCCGCAAACGCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGACGAAGTGGCGCGCCTTGCCGGTCTgg                                                                                                                                  >  1:348714/1‑149 (MQ=255)            tCCAGGCGATCCTGTCTGACAACCACATCACCATTTCCAGCCGACCCATCGTGCCGGTGATCGCCCGTGAGGAGGTGCCGCAAACGCCAGATGGCAGCCCCGCCTCGCACATCGACCTGGACGAAGCGGCGCGCCTTGCCGGTCTGGc                                                                                                                                 >  5:52066/1‑148 (MQ=255)                            tGACAACCACATCACCATTTCCAGCCGACCCATCGTGCCGGTGACCGCCCGTGAGGAGGTGCCGCAAACGCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGACGAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAg                                                                                                                >  8:212514/1‑149 (MQ=255)                               cAACCACATCACCATTTCCAGCCGACCCATCGTGCCGGTGATCGCCCGTGAGGAGGTGCCGCAAACGCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGACGAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTg                                                                                                              <  8:203800/148‑1 (MQ=255)                                  ccACATCACCATTTCCAGCCGACCCATCGTGCCGGTGATCGCCCGTGAGGAGGTGCCGCAAACGCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGACGAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGACCt                                                                                                          <  3:175482/149‑1 (MQ=255)                                  ccACATCACCATTTCCAGCCGACCCATCGTGCCGGTGATCGCCCGTGAGGAGGTGCCGCAAACGCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGACGAAGCGGCGCGCCTTGCCGGTATGGCGATCCCGGGGTTGGAAGTGACCt                                                                                                          <  1:16053/149‑1 (MQ=255)                                    acaTCACCATTTCCAGCCGACACATCGTGCCGGTGATCGCCCGTGAGGAGGTGCCGCAAACGCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGCCGAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGACCTc                                                                                                         >  6:197048/1‑148 (MQ=255)                                         aCCATTTCCAGCCGACCCATCGTGCCGGTGACCGCCCGTGAGGAGGTGCCGCAAACGCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGACGAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGACCTCCATCa                                                                                                    <  7:212514/148‑1 (MQ=255)                                             tttCCAGCCGACCCATCGTGCCGGTGATCGCCCGTGAGGAGGTGCCGCAAACGCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGACGAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGACCTCCATCACGCTc                                                                                               >  1:302454/1‑149 (MQ=255)                                                ccAGCCGACCCATCGTGCCGGTGATCGCCCGTGAGGAGGTGCCGCAAACGCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGACGAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGAcc                                                                                                           <  5:197566/134‑1 (MQ=255)                                                ccAGCCGACCCATCGTGCCGGTGATCGCCCGTGAGGAGGTGCCGCAAACGCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGACGAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGAcc                                                                                                           >  6:197566/1‑134 (MQ=255)                                                 cggCCGACCCATCGTGCCGGTGATCGCCCGTGAGGAGGTGCCGCAAACGCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGACGAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGACCTCCATCACGCTCCCCg                                                                                           >  3:40288/3‑149 (MQ=255)                                                  aGCCGACCCATCGTGCCGGTGATCGCCCGTGAGGAGGTGCCGCAAACGCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGACGAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGACCTCCATCACGCTCCCCGc                                                                                          >  7:37563/1‑149 (MQ=255)                                                   gCCGACCCATCGTGCCGGTGATCGCCCGTGAGGAGGTGCCGCAAACGCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGACGAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGACCTCCATCACGCTCCCCGcc                                                                                         >  1:291974/1‑149 (MQ=255)                                                        cccATCGTGCCGGTGATCGCCCGTGAGGAGGTGCCGCAAACGCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGACGAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGACCTCCATCACGCTCCCCGCCACTGc                                                                                    <  4:149417/149‑1 (MQ=255)                                                        cccATCGTGCCGGTGATCGCCCGTGAGGAGGTGCCGCAAACGCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGACGAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGACCTCCATCACGCTCCCCGCCACTGc                                                                                    >  4:251716/1‑149 (MQ=255)                                                           aTCGTGCCGGTGATCGCCCGTGAGGAGGTGCCGCAAACGCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGACGAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGACCTCCATCACGCTTCCCGCCACt                                                                                      <  6:41182/144‑1 (MQ=255)                                                           aTCGTGCCGGTGATCGCCCGTGAGGAGGTGCCGCAAACGCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGACGAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGACCTCCATCACGCTTCCCGCCACt                                                                                      >  5:41182/1‑144 (MQ=255)                                                            tCGTGCCGGTGATCGCCCGTGAGGAGGTGCCGCAAACGCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGACGAAGTGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGACCTCCATCACGCTCCCCGCCACTGCCTAc                                                                                <  2:348714/149‑1 (MQ=255)                                                                 ccGGTGATCGCCCGTGAGGAGGTGCCGCAAACGCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGACGAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGACCTCCATCACGCTCCCCGCCACTGCCTACAGcc                                                                            <  2:291974/148‑1 (MQ=255)                                                                  cGGTGATCGCCCGTGAGGAGGTGCCGCAAACGCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGACGAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGACCTCCATCACGCTCCCCGCCACTGCCTACAGCCAc                                                                          >  6:148193/1‑149 (MQ=255)                                                                     tGATCGCCCGTGAGGAGGTGCCGCAAACGCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGACGAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGACCTCCATCACGCTCCCCGCCACTGCCTACAGCCACGTc                                                                       <  8:37563/149‑1 (MQ=255)                                                                      gATCGCCCGTGAGGAGGTGCCGCAAACGCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGACGAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGACCTCCATCACGCTCCCCGCCACTGCCTACAGCCACGTCa                                                                      >  2:95194/1‑149 (MQ=255)                                                                          gCCCGTGAGGAGGTGCCGCAAACGCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGACGAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGACCTCCATCACGCTCCACGCCACTGCCTACAGCCACGTCACgag                                                                  >  6:204941/1‑147 (MQ=255)                                                                           cccGTGAGGAGGTGCCGCAAACGCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGACGAAGTGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGACCTCCATCACGCTCCCCGCCACTGCCTACAGCCACGTCACGCTg                                                                 >  1:246445/1‑149 (MQ=255)                                                                            ccGTGAGGAGGTGCCGCAAACGCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGACGAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGACCTCCATCACGCTCCCCGCCACTGCCTACAGCCACGTCACGCTGt                                                                >  8:230017/1‑149 (MQ=255)                                                                                      gTGCCGCAAACGCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGACGAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTCGGAAGTGACCTCCATCACGCTCCCCGCCACTGCCTACAGCCACGTCACGCTGTCTGGCCCggg                                                      <  7:230017/149‑1 (MQ=255)                                                                                       tGCCGCAAACGCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGACGAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGACCTCCATCACGCTCCCCGCCACTGCCTACAGCCACGTCACGCTGTCTGGCCCGGGc                                                     >  8:198912/1‑149 (MQ=255)                                                                                                cGCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGACGAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCg                                                                                                                          <  8:333285/71‑1 (MQ=255)                                                                                                cGCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGACGAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCg                                                                                                                          >  7:333285/1‑71 (MQ=255)                                                                                                cGCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGACGAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGACCTCCATCACGCTCCCCGCCACTGCCTACAGCCACGTCACGCTGTCTGGCCCGGGCTGGTACCCg                                            <  7:198912/149‑1 (MQ=255)                                                                                                 gCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGACGAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGACCTCCATCACGCTCCCCGCCACTGCCTACAGCCACGTCACGCTGTCTGGCCCGGGCTGGTACCCGt                                           <  2:16408/149‑1 (MQ=255)                                                                                                 gCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGACGAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGACCTCCATCACGCTCCCCGCCACTGCCTACAGCCACGTCACGCTGTCTGGCCCGGGCTGGTACCCGt                                           <  5:134211/149‑1 (MQ=255)                                                                                                        ggCAGCCCCGCCTCGCGCATCGACCTGGACGAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGACCACCATCACGCTCACCGCCACTGCCTACAGCCACCTCACGCTGTCTCGCCCGGGCTGGTACCCGTTCATGt                                     >  4:277732/1‑148 (MQ=255)                                                                                                             ccccGCCTCGCGCATCGACCTGGACGAATCGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGACCTCCATCACGCTCCCCGCCACTGCCTACAGCCACGTCACGCTGTCTGGCCCGGGCTGGTACCCGTTGATGTTCCaga                               <  3:277732/149‑1 (MQ=255)                                                                                                              cccGCCTCGCGCATCGACCTGGACGAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGACCTCCATCACGCTCCCCGCCACTGCCTACAGCCACGTCACGCTGTCTGGCCCGGGCTGGTACCCGTTGATGTTCCagag                              <  5:148193/149‑1 (MQ=255)                                                                                                                       cgcATCGACCTGGACGAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGACCTCCATCACGCTCCCCGCCACTGCCTACAGCCACGTCACGCTGTCTGGCCCGGGCTGGTACCCGTTGATGTTCCAGAGTGCTTCGGt                     <  5:204941/149‑1 (MQ=255)                                                                                                                           tCGACCTGGACGAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGACCTCCATCACGCTCCCCGCCACTGCCTACAGCCACGTCACGCTGTCTGGCCCGGGCTGGTACCCGTTGATGTTCCAGAGTGCTTCGGTGAAt                 <  3:251716/149‑1 (MQ=255)                                                                                                                             gACCTGGACGAAGTGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGACCTCCATCACGCTCCCCGCCACTGCCTACAGCCACGTCACGCTGTCTGGCCCGGGCTGGTACCCGTTGATGTTCCAGAGTGCTTCGGTGAATc                <  2:246445/148‑1 (MQ=255)                                                                                                                                cTGGACGAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGACCTCCATCACGCTCCCCGCCACTGCCTACAGCCACGTCACGCTGTCTGGCCCGGGCTGGTACCCGTTGATGTTCCAGAGTGCTTCGGTGAATCCGTa            <  2:302454/149‑1 (MQ=255)                                                                                                                                      gAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGACCTCCATCACGCTCCCCGCCACTGCCTACAGCCACGTCACGCTGTCTGGCCCGGGCTGGTACCCGTTGATGTTCCAGAGTGCTTCGGTGAATCCGt             >  3:264471/1‑142 (MQ=255)                                                                                                                                      gAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGACCTCCATCACGCTCCCCGCCACTGCCTACAGCCACGTCACGCTGTCTGGCCCGGGCTGGTACCCGTTGATGTTCCAGAGTGCTTCGGTGAATCCGt             <  4:264471/142‑1 (MQ=255)                                                                                                                                          cggcgCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGACCTCCATCACGCTCCCCGCCACTGCCTACAGCCACGTCACGCTGTCTGGCCCGGGCTGGTACCCGTTGATGTTCCAGAGTGCTTCGGTGAATCCGTACACACGCAAt  <  6:52066/149‑1 (MQ=255)                                                                                                                                          |                                                                                                                                                     CTGGTTCCTGATCCAGGCGATCCTGTCTGACAACCACATCACCATTTCCAGCCGACCCATCGTGCCGGTGATCGCCCGTGAGGAGGTGCCGCAAACGCCAGATGGCAGCCCCGCCTCGCGCATCGACCTGGACGAAGCGGCGCGCCTTGCCGGTCTGGCGATCCCGGGGTTGGAAGTGACCTCCATCACGCTCCCCGCCACTGCCTACAGCCACGTCACGCTGTCTGGCCCGGGCTGGTACCCGTTGATGTTCCAGAGTGCTTCGGTGAATCCGTACACACGCAAT  >  NC_002947/4327677‑4327962

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