Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,841,978	T→A	10.2%	intergenic (‑114/+138)	yqhD ← / ← PP_2493	NAD(P)‑dependent alcohol dehydrogenase/RadC family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,841,978	0	T	A	10.2%	141.9 / 7.3	49	intergenic (‑114/+138)	yqhD/PP_2493	NAD(P)‑dependent alcohol dehydrogenase/RadC family protein
Reads supporting (aligned to +/- strand):  ref base T (24/20);  new base A (2/3);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

TCGCGTGGGGTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAGGGTTGTTGCTAGGCACTGATGATAGGGCATCGGGAGTCATAGGTT  >  NC_002947/2841832‑2842124                                                                                                                                                   |                                                                                                                                                   tCGCGTCGGGCTGTGGAAGTGAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTAGCTCATTATTTAGTGAAATTTTTATTGCAGAACTTCGTCGGCGCGCATTGCGCACCAAACAGCCGAAGATTGGTTAGAAgcttcc                                                                                                                                                 <  1:799452/150‑3 (MQ=255)   gcgTGGGGTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGGTAGAAgcttgctt                                                                                                                                               <  2:729132/150‑1 (MQ=255)       ggggTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAGGCAAGCTTCTaa                                                                                                                                           <  1:905605/150‑1 (MQ=255)       ggggTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCAAGCTTCTaa                                                                                                                                           <  2:483351/150‑1 (MQ=255)        gggTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCAAGCTTCTaac                                                                                                                                          <  1:519259/150‑2 (MQ=255)                gAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTATATTTTTTTTTTTTTTTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATATCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCa                                                                                                                                  >  2:1056423/1‑150 (MQ=255)                gAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGATGAAAATTGGTTAGAAGCTTTCTTCTAAAGGCAACCa                                                                                                                                  >  2:319965/1‑150 (MQ=255)                      aaaaTTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTgg                                                                                                                            >  2:862876/1‑150 (MQ=255)                         aTTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGACTCTCGGCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGgtgt                                                                                                                         <  1:1241481/150‑1 (MQ=255)                         aTTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGgtgt                                                                                                                         <  1:368745/150‑1 (MQ=255)                                           tACTCCTTAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTcc                                                                                                       >  1:1313733/1‑150 (MQ=255)                                              tCCTTAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCACACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTgg                                                                                                    >  1:766836/1‑150 (MQ=255)                                                cTTAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTaaa                                                                                                                                          >  2:16410/1‑110 (MQ=255)                                                cTTAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTaaa                                                                                                                                          <  1:16410/110‑1 (MQ=255)                                                 ttAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTcc                                                                                                 >  2:137615/1‑150 (MQ=255)                                                 ttAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTcc                                                                                                 >  1:444633/1‑150 (MQ=255)                                                  tAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTccc                                                                                                >  1:1371365/1‑150 (MQ=255)                                                        tAGACCAGTCGTATATTTTGTGGTTAAATTTTTATTGCAGAAGCTCGTCAGCGCGCATTGCGCACAAAATAGCTGAAGATTGGTTAGATGCTTGCTTCTAAAGCGAACCAGCGTGGTGTTGCACTCCAGTGCATTCCTGGTCCCGAGCCa                                                                                          >  2:1156116/1‑150 (MQ=255)                                                          gACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCACACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCaaa                                                                                        <  2:766836/150‑1 (MQ=255)                                                             cAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTgg                                                                                                                       >  1:958135/1‑116 (MQ=255)                                                             cAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTgg                                                                                                                       <  2:958135/116‑1 (MQ=255)                                                              aGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTgg                                                                                                    <  2:165789/134‑1 (MQ=255)                                                              aGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTgg                                                                                                    >  1:165789/1‑134 (MQ=255)                                                               gTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGc                                                                                     <  1:815203/148‑1 (MQ=255)                                                               gTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGc                                                                                     >  2:815203/1‑148 (MQ=255)                                                                   tataTTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTgcc                                                                               >  2:1397012/1‑150 (MQ=255)                                                                     tatTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTaaa                                                                                                                                          >  1:761768/1‑89 (MQ=255)                                                                     tatTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTaaa                                                                                                                                          <  2:761768/89‑1 (MQ=255)                                                                        tttttGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCaa                                                                          <  2:1313733/150‑1 (MQ=255)                                                                               tGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAg                                                                                                                                >  1:1126829/1‑89 (MQ=255)                                                                               tGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAg                                                                                                                                <  2:1126829/89‑1 (MQ=255)                                                                                gAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTa                                                                  >  1:1399187/1‑150 (MQ=255)                                                                                         aTTGCAGAAGCTCGTCGGCGCGCATTGCGCACCACATAGCTGAAGACCGCGCGCATGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGCGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGGTCcgcg                                                         <  1:1156116/150‑1 (MQ=255)                                                                                            gCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGa                                                      >  1:1264259/1‑150 (MQ=255)                                                                                             cAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGAt                                                     >  1:33879/1‑150 (MQ=255)                                                                                               cgatcTCGTCGGCGCGCATTGCGAACCAAATAGATAAAGATTGGTTAGAAGCTTCATTCTAAAGGTAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCTGCCAAGCAGCGTAAGCTCCGCGCGATGt                                                   <  2:299008/146‑1 (MQ=255)                                                                                                  gCTCGTCGGCGCGCATTGCGCACCAAATCGCTGAAGATTGGTTAGACGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCTTGGCTGCCAAGCAGCGTAAGCTCCGCGCGAGGTcag                                                >  1:299008/1‑148 (MQ=255)                                                                                                  gCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGa                                                >  2:1127292/1‑150 (MQ=255)                                                                                                    tcgtcgGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAgg                                              <  2:1399187/150‑1 (MQ=255)                                                                                                       tcgGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAGGgtt                                           >  2:1311079/1‑150 (MQ=255)                                                                                                              gcATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAGGGTTGTTGCTa                                    <  1:1397012/150‑1 (MQ=255)                                                                                                                      cgccAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAGGGTTGTTGCTAGGCACtga                            <  1:163349/148‑1 (MQ=255)                                                                                                                         cAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAGGGTTGTTGCTAGGCACtgatga                         >  2:1005853/1‑150 (MQ=255)                                                                                                                                        cctttAGAAGCAAGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAGGGTTGTTGCTAGGCACTGATGATAGGGCATCGGGAGt          >  1:1116972/4‑150 (MQ=255)                                                                                                                                         ctttAGAAGCAAGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAGGGTTGTTGCTAGGCACTGATGATAGGGCATCGGGAGTc         >  2:766751/3‑150 (MQ=255)                                                                                                                                           ttAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAGGGTTGTTGCTAGGCACTGATGATAGGGCATCGGGAGTCAt       <  1:862876/150‑1 (MQ=255)                                                                                                                                           ttAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAGGGTTGTTGCTAGGCACTGATGATAGGGCATCGGGAGTCAt       <  1:137615/150‑1 (MQ=255)                                                                                                                                            tAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAGGGTTGTTGCTAGGCACTGATGATAGGGCATCGGGAGTCATa      <  2:1264259/150‑1 (MQ=255)                                                                                                                                              gAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAagcagc                                                                     <  2:590605/85‑1 (MQ=255)                                                                                                                                              gAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAagcagc                                                                     >  1:590605/1‑85 (MQ=255)                                                                                                                                                aGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAGGGTTGTTGCTAGGCACTGATGATAGGGCATCGGGAGTCATAGGtt  <  2:444633/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   TCGCGTGGGGTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAGGGTTGTTGCTAGGCACTGATGATAGGGCATCGGGAGTCATAGGTT  >  NC_002947/2841832‑2842124

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