Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	624,430:1	+A	10.3%	coding (477/702 nt)	PP_0537 ←	transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	624,430	1	.	A	10.3%	139.2 / 10.4	39	coding (477/702 nt)	PP_0537	transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (18/17);  new base A (2/2);  total (20/19)
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

AGCAGCGGTGTAGCGTTCGACCGTGTAGTTGCGCCGGTCATGGCTG‑TGCAGGTACAGGGTGTCGTGCGGGCCAAGACTGAGGTTGTTTACCCTGAGCATGCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCATG‑CCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATC‑‑AGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCC  >  NC_002947/624283‑624570                                                                                                                                                      |                                                                                                                                               agcagcGGTGTAGCGTTCGACCGTGTAGTTGCGCCGGTCATGGCTGTTGCAGGTACAGGGTTTCGTGCGGGCCAAGACTGAGGTTGTTTACCCTGAGCATGCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCATg                                                                                                                                                 >  1:60121/1‑149 (MQ=255) agcagcGGTGTAGCGTTCGACCGTGTAGTTGCGCCGGTCATGGCTG‑TGCAGGTACAGGGTGTCGTGCGGGCCAAGACTGAGGTTGTTTACCCTGAGCATGCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCATG‑c                                                                                                                                               <  5:142718/149‑1 (MQ=255)       ggTGTAGCGTTCGACCGTGTAGTTGCGCCGGTCATGGCTG‑TGCAGGTACAGGGTGTCGTGCGGGCCAAGACTGAGGTTGTTTACCCTGAGCATGCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCATG‑CCCAGGt                                                                                                                                         <  4:158212/149‑1 (MQ=255)        gtgtAGCGTTCGACCGTGTAGTTGCGCCGGTCATGGCTG‑TGCAGGTACAGGGTGTCGTGCGGGCCAAGACTGAGGTTGTTTACCCTGAGCATGCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGG‑CATGACCCAGGtg                                                                                                                                        <  1:143623/149‑2 (MQ=255)           tAGCGTTCGACCGTGTAGTTGCGCCGGTCATGGCTG‑TGCAGGTACAGGGTGTCGTGCGGGCCAAGACTGAGGTTGTTTACCCTGAGCATGCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCATG‑CCCAGGTCGAt                                                                                                                                     <  2:60121/149‑1 (MQ=255)             gCGTTCGACCGTGTAGTTGCGCCGGTCATGGCTG‑TGCAGGTACAGGGTGTCGTGCGGGCCAAGACTGAGGTTGTTTACCCTGAGCATGCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCATG‑CCCAGGTCGATg                                                                                                                                    <  4:16722/148‑1 (MQ=255)                                 gcCGGTCATGGCTG‑TGCAGGTACAGGGTGTCGTGCGGGCCAAGACTGAGGTTGTTTACCCTGAGCATGCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCATG‑CCCAGGTCGATGGCCAGGGTAGCTTGCAGCgg                                                                                                                >  8:41955/1‑148 (MQ=255)                                      tCATGGCTG‑TGCAGGTACAGGGTGTCGTGCGGGCCAAGACTGAGGTTGTTTACCCTGAGCATGCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCATG‑CCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATC‑‑Ag                                                                                                        >  7:50271/1‑149 (MQ=255)                                      tCATGGCTG‑TGCAGGTACAGGGTGTCGTGCGGGCCAAGACTGAGGTTGTTTAACCTGAGCATGCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCATG‑CCCAGGTCGATGGCCAGGGTAGCTTGCAGCCGAATC‑‑Ag                                                                                                        >  7:50270/1‑149 (MQ=255)                                          ggCTG‑TGCAGGTACAGGGTGTCGTGCGGGCCAAGACTGAGGTTGTTTACCCTGAGCATGCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCATG‑CCCAGGTCGATGGCCAGGGTAGc                                                                                                                         <  6:270875/130‑1 (MQ=255)                                          ggCTG‑TGCAGGTACAGGGTGTCGTGCGGGCCAAGACTGAGGTTGTTTACCCTGAGCATGCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCATG‑CCCAGGTCGATGGCCAGGGTAGc                                                                                                                         >  5:270875/1‑130 (MQ=255)                                                             tgtCGTGCGGGCCAAGACTGAGGTTGTTTACCCTGAGCATGCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCATG‑CCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATC‑‑AGCGGGAGCATg                                                                                              <  3:62473/137‑1 (MQ=255)                                                             tgtCGTGCGGGCCAAGACTGAGGTTGTTTACCCTGAGCATGCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCATG‑CCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATC‑‑AGCGGGAGCATg                                                                                              >  4:62473/1‑137 (MQ=255)                                                             tgtCGTGCGGGCCAAGACTGAGGTTGTTTACCCTGAGCATGCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCATG‑CCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATC‑‑AGCGGGAGCATGTTGCTGGCCgg                                                                                   >  7:133814/1‑148 (MQ=255)                                                               tCGTGCGGGCCAAGACTGAGGTTGTTTACCCTGAGCATGCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCATG‑CCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATC‑‑AGCGGGAGCATGTTGCTGGCCGGCAt                                                                                >  6:209429/1‑149 (MQ=255)                                                                           aGACTGAGGTTGTTTACCCTGAGCATGCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCATG‑CCAAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATC‑‑AGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGAt                                                                    >  7:13315/1‑149 (MQ=255)                                                                           aGACTGAGGTTGTTTACCCTGAGCATGCCGTTTGGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCATG‑CCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATC‑‑AGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGAt                                                                    <  7:41955/149‑1 (MQ=255)                                                                                  ggTTGTTTACCCTGAGCATGCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCATG‑CCGAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATC‑‑AGCGGGAGCATGTTGCTGGCCGGCAt                                                                                >  6:125799/1‑130 (MQ=255)                                                                                  ggTTGTTTACCCTGAGCATGCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCATG‑CCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATC‑‑AGCGGGAGCATGTTGCTGGCCGGCAt                                                                                <  5:125799/130‑1 (MQ=255)                                                                                       tttACCCTGAGCATGCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCATG‑CCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATC‑‑AGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGc                                                        >  7:251041/1‑149 (MQ=255)                                                                                           cccTGAGCATGCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCATG‑CCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATC‑‑AGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTg                                                     <  5:209429/148‑1 (MQ=255)                                                                                            ccTGAGCATGCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCATG‑CCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATC‑‑AGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGtt                                                              >  5:122792/1‑138 (MQ=255)                                                                                            ccTGAGCATGCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCATG‑CCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATC‑‑AGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGtt                                                              <  6:122792/138‑1 (MQ=255)                                                                                                     gCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCATG‑CCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATC‑‑AGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCtt                                                      <  1:264523/137‑1 (MQ=255)                                                                                                     gCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCATG‑CCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATC‑‑AGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCtt                                                      >  2:264523/1‑137 (MQ=255)                                                                                                     gCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCATG‑CCCAGGTCGATAGCCAGGGTAGCTTGCAGCGGAATC‑‑AGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGc                                          >  6:236987/1‑149 (MQ=255)                                                                                                       cGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCATG‑CCCAGGTCGATAGCCAGGGTAGCTTGCAGCGGAATC‑‑AGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGccc                                        <  5:236987/149‑1 (MQ=255)                                                                                                               agcagGGCATAGGTTTCGCCCTCGACCACCTGGGTCATG‑CCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATC‑‑AGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGc                                >  6:72578/1‑149 (MQ=255)                                                                                                                   gggCATAGGTTTCGCCCTCGACCACCTGGGTCATG‑CCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATC‑‑AGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCtttt                            <  8:13315/149‑1 (MQ=255)                                                                                                                   gggCATAGGTTTCGCCCTCGACCACCTGGGTCATG‑CCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATC‑‑AGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTGAAGGGCttgg                            <  2:73459/149‑3 (MQ=255)                                                                                                                       aTAGGTTTCGCCCTCGACCACCTGGGTCATG‑CCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATCTGGGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGc                          >  1:84057/1‑149 (MQ=255)                                                                                                                            tttCGCCCTCGACCACCTGGGTCATG‑CCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATC‑‑AGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAg                   >  2:60529/1‑149 (MQ=255)                                                                                                                              tCGCCCTCGACCACCTGGGTCATG‑CCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATC‑‑AGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAgc                  <  8:133814/148‑1 (MQ=255)                                                                                                                                  ccTCGACCACCTGGGTCATG‑CCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATC‑‑AGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAgg             <  8:50271/149‑1 (MQ=255)                                                                                                                                  ccTCGACCACCTGGGTCATG‑CCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATC‑‑AGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAgg             >  5:262636/1‑149 (MQ=255)                                                                                                                                  ccTCGACCACCTGGGGCATG‑CCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATC‑‑AGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCTGGGCTTTTGCCGGCAAGCGCAgg             <  8:50270/149‑1 (MQ=255)                                                                                                                                     cGACCACCTGGGTCATG‑CCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATC‑‑AGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGtgc          >  6:250133/1‑149 (MQ=255)                                                                                                                                       tcgaccTGGG‑CATGACCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATC‑‑AGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAggg                                 >  7:14389/4‑124 (MQ=255)                                                                                                                                       tcgaccTGGG‑CATGACCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATC‑‑AGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAggg                                 <  8:14389/121‑1 (MQ=255)                                                                                                                                          accTGGG‑CATGACCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATC‑‑AGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGtgcttgct     >  2:19130/1‑149 (MQ=255)                                                                                                                                               ggTCATG‑CCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATC‑‑AGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGcc  <  6:262636/147‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                               AGCAGCGGTGTAGCGTTCGACCGTGTAGTTGCGCCGGTCATGGCTG‑TGCAGGTACAGGGTGTCGTGCGGGCCAAGACTGAGGTTGTTTACCCTGAGCATGCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCATG‑CCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATC‑‑AGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCC  >  NC_002947/624283‑624570

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