Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,960,142	T→A	9.7%	intergenic (+50/‑78)	bolA → / → PP_1758	DNA‑binding transcriptional dual regulator/rhodanese‑like domain‑containing sulfurtransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,960,142	0	T	A	9.7%	124.0 / 4.5	41	intergenic (+50/‑78)	bolA/PP_1758	DNA‑binding transcriptional dual regulator/rhodanese‑like domain‑containing sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base T (19/18);  new base A (2/2);  total (21/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.79e-01

AGATCCATGCCCTGGCCATCCACACCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGG‑CCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGAAGACTACGTCGAGCTGCGCGA  >  NC_002947/1959995‑1960290                                                                                                                                                    |                                                                                                                                                      agaTCCATGCCCTGGCCATCCACACCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAg                                                                                                                                                     >  1:404131/1‑150 (MQ=255) agaTCCATGCCCTGGCCATCCACACCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTAAg                                                                                                                                                     <  1:135116/150‑1 (MQ=255)     ccATGCCCTGGCCATCCACACCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGG‑c                                                                                                                                                >  1:60416/1‑150 (MQ=255)                  aTCCACACCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGG‑CCCCTTGCTTTTTa                                                                                                                                   <  2:404131/150‑1 (MQ=255)                  aTCCACACCTACACCGCCGAAGAGTGGGCCAAGGGCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGAGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGG‑CCCCTTGCTTTTTa                                                                                                                                   <  2:658926/150‑1 (MQ=255)                   tCCACACCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGG‑CCCCTTGCttt                                                                                                                                      <  2:297133/146‑1 (MQ=255)                   tCCACACCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGG‑CCCCTTGCttt                                                                                                                                      >  1:297133/1‑146 (MQ=255)                               ccgccgAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGG‑CCCCTTGCTTTTTATGCAAACGTCaaa                                                                                                                      >  1:579112/1‑150 (MQ=255)                                      agagTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGG‑CCCCTTGCTTTTTATGCAAACGTCAAACCCGGTc                                                                                                               <  1:573303/150‑1 (MQ=255)                                              ccAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGG‑CCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCtt                                                                                                       >  1:422094/1‑150 (MQ=255)                                              ccAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGG‑CCCCTTGCTTTTTATGCAAAAGTCAAACCCGGTCCGCCCCtt                                                                                                       >  2:477167/1‑150 (MQ=255)                                              ccAAGGTCGGCGTTGCACCGGCCTCACCCGTGTGTGCGGGCGGTGGGCTCTGAAACCGTTCTTTTGTTCTGGTAGGATCCGCATCACCCAAAGGCGGCCGCTTAGCGG‑CCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCtt                                                                                                       <  1:452328/150‑1 (MQ=255)                                                  ggTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGG‑CCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGa                                                                                                   <  2:117851/150‑1 (MQ=255)                                                       gCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGG‑CCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCa                                                                                              >  2:573262/1‑150 (MQ=255)                                                        cGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGG‑CCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCaa                                                                                             <  1:476507/150‑1 (MQ=255)                                                          ttGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGG‑CCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAacc                                                                                           >  1:103045/1‑150 (MQ=255)                                                            gCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTGCTGTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGG‑CCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAaccac                                                                                         >  1:623050/1‑150 (MQ=255)                                                              aCCGGCCTCGCCGGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGG‑CCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCt                                                                                       <  1:390490/150‑1 (MQ=255)                                                              aCCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGG‑CCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCt                                                                                       >  2:479178/1‑150 (MQ=255)                                                                       gCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGG‑CCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAg                                                                                                  >  1:566049/1‑130 (MQ=255)                                                                       gCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGG‑CCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAg                                                                                                  <  2:566049/130‑1 (MQ=255)                                                                                                tGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGG‑CCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCaa                                                                                             >  2:255781/1‑110 (MQ=255)                                                                                                tGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGG‑CCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCaa                                                                                             <  1:255781/110‑1 (MQ=255)                                                                                                      cGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGG‑CCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCgcggc                                               >  1:289419/1‑150 (MQ=255)                                                                                                           ttttGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGG‑CCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGt                                          <  1:479178/150‑1 (MQ=255)                                                                                                            tttGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGG‑CCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTa                                         >  2:343017/1‑150 (MQ=255)                                                                                                               gTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGG‑CCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACaa                                      <  2:422094/150‑1 (MQ=255)                                                                                                                ttCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGG‑CCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAg                                     <  2:103045/150‑1 (MQ=255)                                                                                                                ttCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGG‑CCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAg                                     >  1:73733/1‑150 (MQ=255)                                                                                                                              gCATCAGCCAAAGGGGGCCGCTTAGCGG‑CCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGa                       >  2:462218/1‑150 (MQ=255)                                                                                                                                 tCACCCAAAGGGGGCCGCTTAGCGG‑CCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGAAGa                    <  2:60416/150‑1 (MQ=255)                                                                                                                                   aCCCAAAGGGGGCCGCTTAGCGG‑CCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAg                                     <  2:423744/131‑1 (MQ=255)                                                                                                                                   aCCCAAAGGGGGCCGCTTAGCGG‑CCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAg                                     >  1:423744/1‑131 (MQ=255)                                                                                                                                     ccAAAGGGGGCCGCTTAGCGG‑CCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTAc                                        <  2:46935/126‑1 (MQ=255)                                                                                                                                     ccAAAGGGGGCCGCTTAGCGG‑CCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTAc                                        >  1:46935/1‑126 (MQ=255)                                                                                                                                       aaaGGGGGCCGCTTAGCGG‑CCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGAAGACTACGt              <  2:73733/150‑1 (MQ=255)                                                                                                                                         aaggggCCGCTAAGCGGCCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGgcgc                                             >  1:54933/3‑118 (MQ=255)                                                                                                                                         aaggggCCGCTAAGCGGCCCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGgcgc                                             <  2:54933/116‑1 (MQ=255)                                                                                                                                              gCCGCTTAGCGG‑CCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGt                                    >  2:302905/1‑121 (MQ=255)                                                                                                                                              gCCGCTTAGCGG‑CCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGt                                    <  1:302905/121‑1 (MQ=255)                                                                                                                                              gCCGCTAAGCGG‑CCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAATTCGTCACCCTGGAAGACTACGTCGAGCTg       >  2:216336/1‑150 (MQ=255)                                                                                                                                                   ttAGCGG‑CCCCTTGCTTTTTATGCAAACGGCAAACCCGGTCCGCCCCTTACGAGGGCACCCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGCAGACTACGTCGAGCTGcgcgc  >  2:355998/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                      AGATCCATGCCCTGGCCATCCACACCTACACCGCCGAAGAGTGGGCCAAGGTCGGCGTTGCACCGGCCTCGCCCGTGTGTGCGGGCGGCGGGCACTGAAACCGTTCTTTTGTTCTGGTAGAATCCGCATCACCCAAAGGGGGCCGCTTAGCGG‑CCCCTTGCTTTTTATGCAAACGTCAAACCCGGTCCGCCCCTTACGAGGGCAACCACCTGGAGATTCCGTTACATGTCTCAAGCCATCGTCGTCGCGGCGCTGTACAAGTTCGTCACCCTGGAAGACTACGTCGAGCTGCGCGA  >  NC_002947/1959995‑1960290

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