Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	807,424	C→A	11.1%	D205Y (GAC→TAC)	PP_0694 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	807,424	0	C	A	11.1%	116.0 / 2.7	36	D205Y (GAC→TAC)	PP_0694	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (14/18);  new base A (2/2);  total (16/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 = 9.39e-01

GCAATGCGAAGCGGCGGATGCCGTAGCGTTCCTGGCTCATGCGTCAATCTCCTGTTGCTCTTCGCGGATGGCCCGGGCCAGGGCGTCCTCTTCGCTTTCCTCGCCCTGGAAAGGCGTGAGGTCGAGCGGGTCGTCGGTACGGTTGAGTTCTTCGGCGCTTTCTTCCTCGACCAGTACCGGGGTCGGCAGCGGCAGGGCGCTGTGCAGGGTGGCGGCCAGGTCGCGGTC  >  NC_002947/807293‑807520                                                                                                                                    |                                                                                                 gCAATGCGAAGCGGCGGATGCCGTAGCGTTCCTGGCTCATGCGTCAATCTCCTGTTGCTCTTCGCGGATGGCCCGGGCCAGGGCGTCCTCTTCGCTTTCCTCGCCCTGGAAAGGCGTGAGGTCGAGCGGGtcgtcg                                                                                              <  1:423455/136‑1 (MQ=255)          aGCGGCGGATGCCGTAGCGTTCCTGGCTCATGCGTCAATCTCCTGTTGCTCTTCGCGGATGGCCCGGGCCAGGGCGTCCTCTTCGCTTTCCTCGCCCTGGAAAGGCGTGAGGTCGAGCGGGTCGTCGGTACGGTTg                                                                                     >  1:180625/1‑136 (MQ=255)                 gATGCCGTAGCGTTCCTGGCTCATGCGTCAATCTCCTGTTGCTCTTCGCGGATGGCCCGGGCCAGGGCGTCCTCTTCGCTTTCCTCGCCCTGGAAAGGCGTGAGGTCGAGCGGGTCGTCGGTACGGTTGAGttctt                                                                              <  1:264547/136‑1 (MQ=255)                  aTGCCGTAGCGTTCCTGGCTCATGCGTCAATCTCCTGTTGCTCTTCGCGGATGGCCCGGGCCAGGGCGTCCTCTTCGCTTTCCTCGCCCTGGAAAGGCGTGAGGTCGAGCGGGTCGTCGGTACGGTTGAGttcttc                                                                             >  1:66249/1‑136 (MQ=255)                  aTGCCGTAGCGTTCCTGGCTCATGCGTCAATCTCCTGTTGCTCTTCGCGGATGGCCCGGGCCAGGGCGTCCTCTTCGCTTTCCTCGCCCTGGAAAGGCGTGAGGTCGAGCGGGTCGTCGGTACGGTTGAGttcttc                                                                             >  1:33938/1‑136 (MQ=255)                            gTTCCTGGCTCATGCGTCAATCTCCTGTTGCTCTTCGCGGATGGCCCGGGCCAGGGCGTCCTCTTCGCTTTCCTCGCCCTGGAAAGGCGTGAGGTCGAGCGGGTCGTCGGTACGGTTGAGTTCTTCGGCGCTttct                                                                   <  2:83453/136‑1 (MQ=255)                              tCCTGGCTCATGCGTCAATCTCCTGTTGCTCTTCGCGGATGGCCCGGGCCAGGGCGTCCTCTTCGCTTTCCTCGCCCTGGAAAGGCGTGAGGTCGAGCGGGTCGTCGGTACGGTTGAGTTCTTCGGCGCTttcttc                                                                 <  1:377962/136‑1 (MQ=255)                              tCCTGGCTCATGCGTCAATCTCCTGTTGCTCTTCGCGGATGGCCCGGGCCAGGGCGTCCTCTTCGCTTTCCTCGCCCTGGAAAGGCGTGAGGTCGAGCGGGTCGTCGGTACGGTTGAGTTCTTCGGCGCTttcttc                                                                 <  1:339926/136‑1 (MQ=255)                                 tGGCTCATGCGTCAATCTCCTGTTGCTCTTCGCGGATGGCCCGGGCCAGGGCGTCCTCTTCGCTTTCCTCGCCCTGGAAAGGCGTGAGGTCGAGCGGGTCGTCGGTACGGTTGAGTTCTTCGGCGCTTTCTTCCTc                                                              <  2:241031/136‑1 (MQ=255)                                   gCTCATGCGTCAATCTCCTGTTGCTCTTCGCGGATGGCCCGGGCCAGGGCGTCCTCTTCGCTTTCCTCGCCCTGGAAAGGCGTGAGGTCGAGCGGGTCGTCGGTACGGTTGAGTTCTTCGGCGCTTTCTTCCTCGa                                                            >  1:104659/1‑136 (MQ=255)                                   gCTCATGCGTCAATCTCCTGTTGCTCTTCGCGGATGGCCCGGGCCAGGGCGTCCTCTTCGCTTTCCTCGCCCTGGAAAGGCGTGAGGTCGAGCGGGTCGTCGGTACGGATGAGTTCTTCGGCGCTTTCTTCCTCGa                                                            >  1:360899/1‑136 (MQ=255)                                          cGTCAATCTCCTGTTGCTCTTCGCGGATGGCCCGGGCCAGGGCGTCCTCTTCGCTTTCCTCGCCCTGGAAAGGCGTGAGGTCGAGCGGGTCGTCGGTACGGTTGAGTTCTTCGGCGCTTTCTTCCTCGACCAGTAc                                                     >  2:8268/1‑136 (MQ=255)                                           gTCAATCTCCTGTTGCTCTTCGCGGATGGCCCGGGCCAGGGCGTCCTCTTCGCTTTCCTCGCCCTGCAAAGGCGTGAGGTCGAGCGGGTCCTCGGTACGGTTGAGATCTTCGGCGCTATCTTCCTCGACCAGAAcc                                                    >  1:261420/1‑136 (MQ=255)                                           gTCAATCTCCTGTTGCTCTTAGAGGATGGCCCGGGCCAGGGCGTCCTCTTCGCTTTCCTCGCACTGGAAAGGCGTGATGTCGAGCGGGTCGACGGTACTGTTGAGTTCTTCGTCGCTTTCTTCCTCGACCAGTAcc                                                    >  2:33047/1‑136 (MQ=255)                                                 ctcCTGTTGCTCTTCGCGGATGGCCCGGGCCAGGGCGTCCTCTTCGCTTTCCTCGCCCTGGAAAGGCGTGAGGTCGAGCGGGTCGTCGGTACGGTTGAGTTCTTCGGCGCTTTCTTCCTCGACCAGTACCGGGGTc                                              <  2:347649/136‑1 (MQ=255)                                                   ccTGTTGCTCTTCGCGGATGGCCCGGGCCAGGGCGTCCTCTTCGCTTTCCTCGCCCTGGAAAGGCGTGAGGTCGAGCGGGTCGTCGGTACGGTTGAGTTCTTCGGCGCTTTCTTCCTCGACCAGTACCGGGGTcgg                                            <  2:66249/136‑1 (MQ=255)                                                   ccTGTTGCTCTTCGCGGATGGCCCGGGCCAGGGCGTCCTCTTCGCTTTCCTCGCCCTGGAAAGGCGTGAGGTCGAGCGGGTCGTCGGTACGGTTGAGTTCTTCGGCGCTTTCTTCCTCGACCAGTACCGGGGTcgg                                            <  2:33938/136‑1 (MQ=255)                                                        tGCTCTTCGCGGATGGCCCGGGCCAGGGCGTCCTCTTCGCTTTCCTCGCCCTGGAAAGGCGTGAGGTCGAGCGGGTCGTCGGTACGGTTGAGTTCTTCGGCGCTTTCTTCCTCGACCAGTACCGGGGTcggcagcg                                       <  2:360899/136‑1 (MQ=255)                                                        tGCTCTTCGCGGATGGCCCGGGCCAGGGCGTCCTCTTCGCTTTCCTCGCCCTGGAAAGGCGTGAGGTCGAGCGGGTCGTCGGTACGGTTGAGTTCTTCGGCGCTTTCTTCCTCGACCAGTACCGGGGTcggcagcg                                       >  1:189229/1‑136 (MQ=255)                                                          ctctTCGCGGATGGCCCGGGCCAGGGCGTCCTCTTCGCTTTCCTCGCCCTGGAAAGGCGTGAGGTCGAGCGGGTCGTCGGTACGGTTGAGTTCTTCGGCGCTTTCTTCCTCGACCAGTACCGGGGTcggcagcggc                                     >  2:401134/1‑136 (MQ=255)                                                             ttCGCGGATGGCCCGGTCCAGGGCGTCCTCTTCGCTTTCCTCGCCCTGGAAAGGCGTGAGGTCGAGCGGGTCGTCGGTACGGTTGAGTTCTTCGGCGCTTTCTTCCTCGACCAGTACCGGGGTCGGCAGCGGCAgg                                  <  1:35034/136‑1 (MQ=255)                                                             ttCGCGGATGGCCCGGGCCAGGGCGTCCTCTTCGCTTTCCTCGCCCTGGAAAGGCGTGAGGTCGAGCGGGTCGTCGGTACGGTTGAGTTCTTCGGCGCTTTCTTCCTCGACCAGTACCGGGGTCGGCAGCGGCAgg                                  <  1:238254/136‑1 (MQ=255)                                                             ttCGCGGATGGCCCGGGCCAGGGCGTCCTCTTCGCTTTCCTCGCCCTGGAAAGGCGTGAGGTCGAGCGGGTCGTCGGTACGGTTGAGTTCTTCGGCGCTTTCTTCCTCGACCAGTACCGGGGTCGGCAGCGGCAgg                                  <  1:8268/136‑1 (MQ=255)                                                             ttCGCGGATGGCCCGGGCCAGGGCGTCCTCTTCGATTTCCTCGCCCTGGAAATGCGTGAGGTCGAGCGGGTCGTCGGTACAGCTGAGTTCTTCGGCTCTATCTTCCTCGACCAGTACCGGGGTCGGCAGCGGCAgg                                  <  1:33047/136‑1 (MQ=255)                                                                 cgGATGGCCCGGGCCAGGGCGTCCTCTTCGCTTTCCTCGCCCTGGAAAGGCGTGAGGTCGAGCGGGTAGTCGGTACGGTTGAGTTCTTCGGCGCTTTCTTCCTCGACCAGTACCGGGGTCGGCAGCGGCAGGgcgc                              >  2:222378/1‑136 (MQ=255)                                                                 cgGATGGCCCGGGCCAGGGCGTCCTCTTCGCTTTCCTCGCCCTGGAAAGGCGTGAGGTCGAGCGGGTAGTCGGTACGGTTGAGTTCTTCGGCGCTTTCTTCCTCGACCAGTACCGGGGTCGGCAGCGGCAGGgcgc                              >  2:233663/1‑136 (MQ=255)                                                                  ggATGGCCCGGGCCAGGGCGTCCTCTTCGCTTTCCTCGCCCTGGAAAGGCGTGAGGTCGAGCGGGTCGTCGGTACGGTTGAGTTCTTCGGCGCTTTCTTCCTCGACCAGTACCGGGGTCGGCAGCGGCAGGGCGCt                             <  1:150431/136‑1 (MQ=255)                                                                    aTGGCCCGGGCCAGGGCGTCCTCTTCGCTTTCCTCGCCCTGGAAAGGCGTGAGGTCGAGCGGGTCGTCGGTACGGTTGAGTTCTTCGGCGCTTTCTTCCTCGACCAGTACCGGGGTcgg                                            >  2:95632/1‑119 (MQ=255)                                                                    aTGGCCCGGGCCAGGGCGTCCTCTTCGCTTTCCTCGCCCTGGAAAGGCGTGAGGTCGAGCGGGTCGTCGGTACGGTTGAGTTCTTCGGCGCTTTCTTCCTCGACCAGTACCGGGGTcgg                                            <  1:95632/119‑1 (MQ=255)                                                                    aTGGCCCGGGCCAGGGCGTCCTCTTCGCTTTCCTCGCCCTGGAAAGGCGTGAGGTCGAGCGGGTCGTCGGTACGGTTGAGTTCTTCGGCGCTTTCTTCCTCGACCAGTACCGGGGTCGGCAGCGGCAGGGCGCtgt                           <  1:401134/136‑1 (MQ=255)                                                                                aGGGCGTCCTCTTCGCTTTCCTCGCCCTGGAAAGGCGTGAGGTCGAGCGGGTCGTCGGTACGGTTGAGTTCTTCGGCGCTTTCTTCCTCGACCAGTACCGGGGTCGGCAGCGGCAGGGCGCTGTGCAGGGTggcgg               <  1:122583/136‑1 (MQ=255)                                                                                     gTCCTCTTCGCTTTCCTCGCCCTGGAAAGGCGTGAGGTCGAGCGGGTCGTCGGTACGGTTGAGTTCTTCGGCGCTTTCTTCCTCGACCAGTACCGGGGTCGGCAGCGGCAGGGCGCTGTGCAGGGTGGCGGCCAgg          >  2:402517/1‑136 (MQ=255)                                                                                     gTCCTCGTCGCTTTCCTCGCCCAGGAAACGCGTGAGGTCGAGCGGGTCGTCGGTACGGTTGAGTTCTTCGGCGCTTTCTTCCACGACCAGTACCGGGGGCGGCAGAGGCAGGGCGCTGTGCAGGGGGGGGGCCAgg          >  1:395091/1‑136 (MQ=255)                                                                                          ctTCGCTTTCCTCGCCCTGGAAAGGCGTGAGGTCGAGCGGGTCGTCGGTACGGTTGAGTTCTTCGGCGCTTTCTTCCTCGACCAGTACCGGGGTCGGCAGCGGCAGGGCGCTGTGCAGGGTGGCGGCCAGGTcgcg     >  1:7936/1‑136 (MQ=255)                                                                                          ctTCGCTTTCCTCGCCCTGGAAAGGCGTGAGGTCGAGCGGGTAGTCGGTACGGTTGAGTTCTTCGGCGCTTTCTTCCTCGACCAGTACCGGGGTCGGCAGCGGCAGGGCGCTGTGCAGGGTGGCGGCCAGGTcgcg     <  1:233663/136‑1 (MQ=255)                                                                                          ctTCGCTTTCCTCGCCCTGGAAAGGCGTGAGGTCGAGCGGGTAGTCGGTACGGTTGAGTTCTTCGGCGCTTTCTTCCTCGACCAGTACCGGGGTCGGCAGCGGCAGGGCGCTGTGCAGGGTGGCGGCCAGGTcgcg     <  1:222378/136‑1 (MQ=255)                                                                                             cGCTTTCCTCGCCCTGGAAAGGCGTGAGGTCGAGCGGGTCGTCGGTACGGTTGAGTTCTTCGGCGCTTTCTTCCTCGACCAGTACCGGGGTCGGCAGCGGCAGGGCGCTGTGCAGGGTGGCGGCCAGGTCGCGGTc  <  1:52408/136‑1 (MQ=255)                                                                                                                                    |                                                                                                 GCAATGCGAAGCGGCGGATGCCGTAGCGTTCCTGGCTCATGCGTCAATCTCCTGTTGCTCTTCGCGGATGGCCCGGGCCAGGGCGTCCTCTTCGCTTTCCTCGCCCTGGAAAGGCGTGAGGTCGAGCGGGTCGTCGGTACGGTTGAGTTCTTCGGCGCTTTCTTCCTCGACCAGTACCGGGGTCGGCAGCGGCAGGGCGCTGTGCAGGGTGGCGGCCAGGTCGCGGTC  >  NC_002947/807293‑807520

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