Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	624,426	Δ1 bp	9.1%	coding (481/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,426	0	T	.	9.1%	131.8 / 7.1	44	coding (481/702 nt)	PP_0537	transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (21/18);  new base . (2/2);  total (23/21)
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.95e-01

GCTGAGCAGCGGTGTAGCGTTCGACCGTGTAGTTGCGCCGGTCATGGCTG‑TGCAGGTACAGGGTGTCGTGCGGGCCAAGACTGAGGTTGTTTACCCTGAGCATGCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCATG‑CCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATC‑‑AGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCC  >  NC_002947/624279‑624570                                                                                                                                                     |                                                                                                                                                    gCTGAGCAGCGGTGTAGCGTTCGACCGTGTAGTTGCGCCGGTCATGGCTG‑TGCAGGTACAGGGTGTCGTGCGGGCCAAGACTGAGGTTGTTTACCCTGAGCATGCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGt                                                                                                                                                     >  3:16722/1‑148 (MQ=255)   tGAGCAGCGGTGTAGCGTTCGACCGTGTAGTTGCGCCGGTCATGGCTG‑TGCAGGTACAGGGTGTCGTGCGGGCCAAGACTGAGGTTGTTTACCCTGAGCATGCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCAt                                                                                                                                                  >  3:158212/1‑149 (MQ=255)   tGAGCAGCGGTGTAGCGTTCGACCGTGTAGTTGCGCCGGTCATGGCTG‑TGCAGGTAAAGGGTGTCGTGCGGGCCAAGACTGAGGTTGTTTACCCTGAGCATGCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCAt                                                                                                                                                  >  1:73459/1‑149 (MQ=255)     agcagcGGTGTAGCGTTCGACCGTGTAGTTGCGCCGGTCATGGCTGTTGCAGGTACAGGGTTTCGTGCGGGCCAAGACTGAGGTTGTTTACCCTGAGCATGCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCATg                                                                                                                                                 >  1:60121/1‑149 (MQ=255)     agcagcGGTGTAGCGTTCGACCGTGTAGTTGCGCCGGTCATGGCTG‑TGCAGGTACAGGGTGTCGTGCGGGCCAAGACTGAGGTTGTTTACCCTGAGCATGCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCAt                                                                                                                                                  <  4:211622/147‑1 (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                                                              <  6:122792/138‑1 (MQ=255)                                                                                                ccTGAGCATGCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCATG‑CCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATC‑‑AGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGtt                                                              >  5:122792/1‑138 (MQ=255)                                                                                                         gCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCATG‑CCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATC‑‑AGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCtt                                                      >  2:264523/1‑137 (MQ=255)                                                                                                         gCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCATG‑CCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATC‑‑AGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCtt                                                      <  1:264523/137‑1 (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)                                                                                                                                                     |                                                                                                                                                    GCTGAGCAGCGGTGTAGCGTTCGACCGTGTAGTTGCGCCGGTCATGGCTG‑TGCAGGTACAGGGTGTCGTGCGGGCCAAGACTGAGGTTGTTTACCCTGAGCATGCCGTTTTGCAGCAGGGCATAGGTTTCGCCCTCGACCACCTGGGTCATG‑CCCAGGTCGATGGCCAGGGTAGCTTGCAGCGGAATC‑‑AGCGGGAGCATGTTGCTGGCCGGCATTGCCAGGCAGATTGCGTTGCCGGCTTGCACACCCAAGCCCTTCAGGGCTTTTGCCGGCAAGCGCAGGTGCTTGCTGCC  >  NC_002947/624279‑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: ‑