Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,744,610	A→C	7.8%	A298A (GCA→GCC)	paoB →	aromatic aldehyde dehydrogenase FAD‑binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,744,610	0	A	C	7.8%	157.4 / 4.5	51	A298A (GCA→GCC)	paoB	aromatic aldehyde dehydrogenase FAD‑binding subunit
Reads supporting (aligned to +/- strand):  ref base A (25/22);  new base C (2/2);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.68e-01

TGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCGAGACCCGGTTGTGGAGCGCATGCTGGCCGGCCAGGCGCCTTGCGAAGCGCTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGCA‑ATTGTGCGGGCATTGGATGACGCGGCGCGCGGCGCCGCTGACGGGGGGCCACGACCATGATCAAACCCACACCGGTAGCCGGTACGTCCCTGGACCGC  >  NC_002947/3744463‑3744756                                                                                                                                                    |                                                                                                                                                     tGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCACACAAGCCGTGGCGAGACCCGGTTGTGGAGCGCATGCTGGCCGGCCAGGCGCCTTGCGAAGCGCTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCa                                                                                                                                                      >  1:76330/1‑148 (MQ=255)   gCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCGAGACCCGGTTGTGGAGCGCATGCTGGCCGGCCAGGCGCCTTGCGAAGCGCTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCaaaa                                                                                                                                                   >  5:237440/1‑149 (MQ=255)    ccAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCGAGACCCGGTTGTGGAGCGCATGCTGGCCGGCCAGGCGCCTTGCGAAGCGCTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAAc                                                                                                                                                  <  4:167833/149‑1 (MQ=255)    ccAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCGAGACCCGGTTGTGGAGCGCATGCTGGCCGGCCAGGCGCCTTGCGAAGCGCTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAAc                                                                                                                                                  >  3:126317/1‑149 (MQ=255)    ccAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCGAGACACGGTTGTGGAGCGCATGCTGGCCGGCCAGGCGCCTTGCGAAGCGCTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCaaaa                                                                                                                                                   >  6:127062/1‑148 (MQ=255)       ggTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCGAGACCCGGTTGTGGTGCGCATGCTGGCCGGCCAGGCGCCTTGCGAAGCGCTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAACctc                                                                                                                                               >  8:306633/1‑149 (MQ=255)        gTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCGAGACCCGGTTGTGGAGCGCATGCTGGCCGGCCAGGCGCCTTGCGAAGCGCTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAACctct                                                                                                                                              <  5:208811/149‑1 (MQ=255)           aTCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCGAGACCCGGTTGTGGAGCGCATGCTGGCCGGCCAGGCGCCTTGCGAAGCGCTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAACCTCTGa                                                                                                                                            <  8:127539/148‑1 (MQ=255)              cGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCGAGACCCGGTTGTGGAGCGCATGCTGGCCGGCCAGGCGCCTTGCGAAGCGCTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAACCTCTGAGCca                                                                                                                                        <  2:97467/149‑1 (MQ=255)                    gCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCGAGACCCGGTTGTGGAGCGCATGCTGGCCGGCCAGGCGCCTTGCGAAGCGCTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAACCTCTGAGCCA‑CAACgg                                                                                                                                 >  7:194418/1‑149 (MQ=255)                         ccTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCGAGACCCGGTTGTGGAGCGCATGCTGGCCGGCCAGGCGCCTTGCGAAGCGCTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTCa                                                                                                                            <  6:237440/149‑1 (MQ=255)                               gcTGGGCGGCGTTGCTCACAAGCCGTGGCGAGACCCGGTTGTGGAGCGCATGCTGGCCGGCCAGGCGCCTTGCGAAGCGCTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCa                                                                                                                      <  3:16897/149‑1 (MQ=255)                                              tCACAAGCCGTGGCGAGACCCGGTTGTGGAGCGCATGCTGGCCGGCCAGGCGCCTTGCGAAGCGCTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGc                                                                                                               >  7:304574/1‑141 (MQ=255)                                              tCACAAGCCGTGGCGAGACCCGGTTGTGGAGCGCATGCTGGCCGGCCAGGCGCCTTGCGAAGCGCTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGc                                                                                                               <  8:304574/141‑1 (MQ=255)                                                acaAGCCGTGGCGAGACCCGGTTGTGGAGCGCATGCTGGCCGGCCAGGCGCCTTGCGAAGCGCTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAAc                                                                                                          >  5:2971/1‑144 (MQ=255)                                                acaAGCCGTGGCGAGACCCGGTTGTGGAGCGCATGCTGGCCGGCCAGGCGCCTTGCGAAGCGCTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAAc                                                                                                          <  6:2971/144‑1 (MQ=255)                                                acaAGCCGTGGCGAGACCCGGTTGTGGAGCGCATGCTGGCCGGCCAGGCGCCTTGCGAAGCGCTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGc                                                                                                      >  8:91653/1‑148 (MQ=255)                                                                  cGGTTGTGGAGCGCATGCTGGCCGGCCAGGCGCCTTGCGAAGCGCTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGCA‑ATTGTGCGGGCATTGGAt                                                                                  >  1:91657/1‑149 (MQ=255)                                                                    gTTGTGGAGCGCATGCTGGCCGGCCAGGCGCCTTGCGAAGCGCTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGCA‑ATTGTGCGGGCATTGGAt                                                                                  >  1:213816/1‑147 (MQ=255)                                                                    gTTGTGGAGCGCATGCTGGCCGGCCAGGCGCCTTGCGAAGCGCTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGCA‑ATTGTGCGGGCATTGGAt                                                                                  <  2:213816/147‑1 (MQ=255)                                                                    gTTGTGGAGCGCATGCTGGCCGGCCAGGCGCCTTGCGAAGCGCTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGCA‑ATTGTGCGGGCATTGGATGa                                                                                >  3:122140/1‑149 (MQ=255)                                                                          gAGCGCATGCTGGCCGGCCAGGCGCCTTGCGAAGCGCTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGCA‑ATTGTGCGGGCATTGGATGAcgcggc                                                                          >  8:26979/1‑149 (MQ=255)                                                                                 tGCTGGCCGGCCAGGCGCCTTGCGAAGCGCTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGCA‑ATTGTGCGGGCATTGGATGAcgcggcgcgcggc                                                                   >  1:79916/1‑149 (MQ=255)                                                                                 tGCTGGCCGGCCAGGCGCCTTGCGAAGCGCTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGCA‑ATTGTGCGGGCATTGGATGAcgcggcgcgcggc                                                                   >  5:67823/1‑149 (MQ=255)                                                                                   cTGGCCGGCCAGGCGCCTTGCGAAGCGCTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGCA‑ATTGTGCGGGCATTGGATGAcgcggcgcgcggcg                                                                  <  5:127062/148‑1 (MQ=255)                                                                                   cTGGCCGGCCAGGCGCCTTGCGAAGCGCTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGCA‑ATTGTGCGGGCATTGGATGACgcggcgcgcggcgc                                                                 <  4:63140/149‑1 (MQ=255)                                                                                                   cTTGCGAAGCGCTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGCA‑ATTGTGCGGGCATTGGATGACGCGGCGCGCGGCGCCGCTGACggggg                                                     <  4:132885/145‑1 (MQ=255)                                                                                                   cTTGCGAAGCGCTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGCA‑ATTGTGCGGGCATTGGATGACGCGGCGCGCGGCGCCGCTGACggggg                                                     <  4:137974/145‑1 (MQ=255)                                                                                                   cTTGCGAAGCGCTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGCA‑ATTGTGCGGGCATTGGATGACGCGGCGCGCGGCGCCGCTGACggggg                                                     >  3:137974/1‑145 (MQ=255)                                                                                                   cTTGCGAAGCGCTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGCA‑ATTGTGCGGGCATTGGATGACGCGGCGCGCGGCGCCGCTGACggggg                                                     >  3:132885/1‑145 (MQ=255)                                                                                                      gCGAAGCGCTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGCA‑ATTGTGCGGGCATTGGATGACGCGGCGCGCGGCGCCGCTGACGGGGGGCCAc                                                <  7:91653/147‑1 (MQ=255)                                                                                                           gcgcTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGCA‑ATTGTGCGGGCATTGGATGACGCGGCGCGCGGCGCCGCTGACGGGGGGCCACGACCATg                                         <  4:126317/149‑1 (MQ=255)                                                                                                            cgcTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGTTGCAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGCA‑ATTGTGCGGGCATTGGATGACGCGGCGCGCGGCGCCGCTGACGGGGGGCCa                                                 <  6:104003/140‑1 (MQ=255)                                                                                                            cgcTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGCA‑ATTGTGCGGGCATTGGATGACGCGGCGCGCGGCGCCGCTGACGGGGGGCCa                                                 >  5:104003/1‑140 (MQ=255)                                                                                                             gcTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGCA‑ATTGTGCGGGCATTGGATGACGCGGCGCGCGGCGCCGCTGACGGGGGGCCACGACCATGAt                                       <  2:5094/149‑1 (MQ=255)                                                                                                               tGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGCA‑ATTGTGCGGGCATTGGATGACGCGGCGCGCGGCGCCGCTGACGGGGGGCCACGACCATGATc                                      >  8:197268/1‑148 (MQ=255)                                                                                                                   cACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAACCTCTGCGCCA‑CAACGGCTTc                                                                                                                             <  7:291328/58‑1 (MQ=255)                                                                                                                   cACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTc                                                                                                                             >  8:291328/1‑58 (MQ=255)                                                                                                                    aCCCGGGCTTCAAATAGTTTGCTGCAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGCA‑ATTGTGGGGGCATTGGATGACGCGGCGCGCGGCGCCGCTGACGGGGGGCCACGACCATGATCAAAcc                                 <  8:194418/148‑1 (MQ=255)                                                                                                                      ccGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAACCTCTGAGCCAGCAACGGCTTCAAGGTCACGCTCGCCCAACGTGCA‑ATTGTGCGGGCATTGGATGACGCGGCGCGCGGCGCCGCTGACGGGGGGCCACGACCATGATCAAACCca                               >  1:244997/1‑149 (MQ=255)                                                                                                                            tGCAAAAACGTTGCTGCAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGCA‑ATTGTGCGGGCATTGGATGACGCGGCGCGCGGCGCCGCTGACGGGGGGCCACGACCATGATCAAACCCACACCgg                         >  5:325374/1‑148 (MQ=255)                                                                                                                            tGCAAAAACGTTGCTGCAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGCA‑ATTGTGCGGGCATTGGATGACGCGGCGCGCGGCGCCGCTGACGGGGGGCCACGACCATGATCAAACCCACACCgg                         >  6:156499/1‑148 (MQ=255)                                                                                                                              cAAAAACGTTGCTGCAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGCA‑ATTGTGCGGGCATTGGATGACGCGGCGCGCGGCGCCGCTGACGGGGGGCCACGACCATGATCAAACCCACACCGGt                        <  2:91657/147‑1 (MQ=255)                                                                                                                                 aaaCGTTGCTGCAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGCA‑ATTGTGCGGGCATTGGATGACGCGGCGCGCGGCGCCGCTGACGGGGGGCCACGACCATGATCAAACCCACACCGGTAGCCg                   >  6:121879/1‑149 (MQ=255)                                                                                                                                         ctgcAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGCA‑ATTGTGCGGGCATTGGATGACGCGGCGCGCGGCGCCGCTGACGGGGGGCCACGACCATGATCAAACCCACACCGGTAGCCGGTACGTcc           <  7:306633/149‑1 (MQ=255)                                                                                                                                         ctgcAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGCA‑ATTGTGCGGGCATTGGATGACGCGGCGCGCGGCGCCGCTGACGGGGGGCCACGACCATGATCAAACCCACACCGGTAGCCGGTACGTcc           >  5:37785/1‑149 (MQ=255)                                                                                                                                         ctgcAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGCA‑ATTGTGCGGGCATTGGATGACGCGGCGCGCGGCGCCGCTGACGGGGGGCCACGACCATGATCAAACCCACACCGGTAGCCGGTACGTcc           >  5:295135/1‑149 (MQ=255)                                                                                                                                          tgcAAGGTGCCAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGCA‑ATTGTGCGGGCATTGGATGACGCGGCGCGCGGCGCCGCTGACGGGGGGcc                                                  >  8:180147/1‑109 (MQ=255)                                                                                                                                          tgcAAGGTGCCAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGCA‑ATTGTGCGGGCATTGGATGACGCGGCGCGCGGCGCCGCTGACGGGGGGcc                                                  >  8:185158/1‑109 (MQ=255)                                                                                                                                          tgcAAGGTGCCAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGCA‑ATTGTGCGGGCATTGGATGACGCGGCGCGCGGCGCCGCTGACGGGGGGcc                                                  <  7:180147/109‑1 (MQ=255)                                                                                                                                          tgcAAGGTGCCAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGCA‑ATTGTGCGGGCATTGGATGACGCGGCGCGCGGCGCCGCTGACGGGGGGcc                                                  <  7:185158/109‑1 (MQ=255)                                                                                                                                           gcAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGCA‑ATTGTGCGGGCATTGGATGACGCGGCGCGCGGCGCCGCTGACGGGGGGCCACGACCATGATCAAACCCACACCGGTAGCCGGTACGTCCCt         >  3:103085/1‑149 (MQ=255)                                                                                                                                                   cAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGCAGATTGTGCGGGCATTGGATGACGCGGCGCGCGGCGCCGCTGACGGGGGGCCACGACCATGATCAAACCCACACCGGTAGCCGGTACGTCCCTGGACcgc  >  1:258850/1‑149 (MQ=255)                                                                                                                                                   cAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGCA‑ATTGTGCGGGCATTGGATGACGCGGCGCGCGGCGCCGCTGACGGGGGGCCACGACCATGATCAAACCCACACCGGTAGCCGGTACGTCCCTGGACcgc  <  2:79916/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     TGGCCAGGTGATCCGGCAAGCTCGCCTGGCGCTGGGCGGCGTTGCTCACAAGCCGTGGCGAGACCCGGTTGTGGAGCGCATGCTGGCCGGCCAGGCGCCTTGCGAAGCGCTGTTCACCCGGGCTGCAAAAACGTTGCTGCAAGGTGCAAAACCTCTGAGCCA‑CAACGGCTTCAAGGTCACGCTCGCCCAACGTGCA‑ATTGTGCGGGCATTGGATGACGCGGCGCGCGGCGCCGCTGACGGGGGGCCACGACCATGATCAAACCCACACCGGTAGCCGGTACGTCCCTGGACCGC  >  NC_002947/3744463‑3744756

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