Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,456,659	C→A	10.2%	L115M (CTG→ATG)	leuS →	leucine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,456,659	0	C	A	10.2%	114.3 / 2.5	39	L115M (CTG→ATG)	leuS	leucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base C (16/19);  new base A (2/2);  total (18/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.10e-01

CAGCGCATGCTGGGCAAGAACGTCCTGCAGCCGATGGGCTGGGACGCTTTCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCACCCGCCTGTTCGAAAAAGGCATCATCTACCGCAAGAACGGTACCGTGAACTGGGACCCGGCGGACCAGACC  >  NC_002947/5456512‑5456808                                                                                                                                                    |                                                                                                                                                      cAGCGCATGCTGGGCAAGAACGTCCTGCAGCCGATGGGCTGGGACGCTTTCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAgcctg                                                                                                                                                     <  2:395353/150‑1 (MQ=255)  aGCGCATGCTGGGCAAGAACGTCCTGCAGCCGATGGGCTGGGACGCTTTCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAgcctgg                                                                                                                                                    <  1:487125/150‑1 (MQ=255)   gagcctGCTGGGCAAGAACGGCCTGCAGCCGATGGGCTGGGACGCTTTCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAgcctggg                                                                                                                                                   >  2:491872/6‑150 (MQ=255)    cgcATGCTGGGCAAGAACGTCCTGCAGCCAATGGGCTGGGACGCTTTCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCCCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAgcctgggc                                                                                                                                                  <  2:536520/150‑1 (MQ=255)             ggCAAGAACGTCCTGCAGCCGATGGGCTGGGACGCTTTCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAgcctg                                                                                                                                                     >  1:77670/1‑138 (MQ=255)             ggCAAGAACGTCCTGCAGCCGATGGGCTGGGACGCTTTCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAgcctg                                                                                                                                                     <  2:77670/138‑1 (MQ=255)              gCAAGAACGTCCTGCAGCCGATGGGCTGGGACGCTTTCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCg                                                                                                                                        <  2:497680/150‑1 (MQ=255)               cAAGAACGTCCTGCAGCCGATGGGCTGGGACGCTTTCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCATGGGCCTGGCCATCGa                                                                                                                                       >  1:207029/1‑150 (MQ=255)               cAAGAACGTCCTGCAGCCGATGGGCTGGGACGCTTTCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCATGGGCCTGGCCATCGa                                                                                                                                       >  1:206881/1‑150 (MQ=255)                  gAACGTCCTGCAGCCGATGGGCTGGGACGCTTTCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTg                                                                                                                                    >  2:426711/1‑150 (MQ=255)                          tGCAGCCGATGGGCTGGGACGCTTTCGGCATGCCCGCGGAAAACGCGGAGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTg                                                                                                                            <  2:107585/150‑1 (MQ=255)                                      gCTGGGACGCTTTCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCt                                                                                                                <  1:426711/150‑1 (MQ=255)                                        tGGGACGCTTTCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGc                                                                                                              >  1:487829/1‑150 (MQ=255)                                           gACGCTTTCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCaaa                                                                                                           <  2:487829/150‑1 (MQ=255)                                                tttCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGCACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGa                                                                                                      <  2:688488/150‑1 (MQ=255)                                                 ttCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGAc                                                                                                     <  2:579427/150‑1 (MQ=255)                                                 ttCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCATGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGAc                                                                                                     <  2:207029/150‑1 (MQ=255)                                                 ttCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCATGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGAc                                                                                                     <  2:206881/150‑1 (MQ=255)                                                     gCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTAtt                                                                                                 <  1:492710/150‑1 (MQ=255)                                                       aTGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTAc                                                                                               <  1:674653/150‑1 (MQ=255)                                                             gcgGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGAc                                                                                                                                      <  2:222168/105‑1 (MQ=255)                                                             gcgGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGAc                                                                                                                                      >  1:222168/1‑105 (MQ=255)                                                                          cgATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGt                                                                            <  1:491872/150‑1 (MQ=255)                                                                          cgATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGt                                                                            <  1:341276/150‑1 (MQ=255)                                                                             tGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCa                                                                                                                                <  1:505010/95‑1 (MQ=255)                                                                             tGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCa                                                                                                                                >  2:505010/1‑95 (MQ=255)                                                                             tGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCa                                                                         <  1:763520/150‑1 (MQ=255)                                                                                   acaacGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACGACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGt                                                                            <  1:761307/141‑1 (MQ=255)                                                                                   acaacGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGt                                                                            >  2:761307/1‑141 (MQ=255)                                                                                           gCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTAc                                                                                               >  1:602913/1‑114 (MQ=255)                                                                                           gCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTAc                                                                                               <  2:602913/114‑1 (MQ=255)                                                                                                    aaGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCACCCGCCTGTTCGAAAAAGGcatc                                                  >  2:426313/1‑150 (MQ=255)                                                                                                             tACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGCAGCAGTGGCTGTACACCCGCCTGTTCGAAAAAGCCATCATCTACCGc                                         >  2:81580/1‑150 (MQ=255)                                                                                                                gAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCGTCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCAcc                                                                       >  2:92293/1‑117 (MQ=255)                                                                                                                gAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCa                                                                         <  1:336435/115‑1 (MQ=255)                                                                                                                gAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCa                                                                         >  2:336435/1‑115 (MQ=255)                                                                                                                gAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCAcc                                                                       <  1:92293/117‑1 (MQ=255)                                                                                                                               aTGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCACCCGCCTGTTCGAAAAAGGCATCATCTACCGCAAGAACGGTACCGTGAAc                       >  1:313351/1‑150 (MQ=255)                                                                                                                               aTGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCACCCGCCTGTTCGAAAAAGGCATCATCTACCGCAAGAACGGTACCGTGAAc                       >  1:463746/1‑150 (MQ=255)                                                                                                                               aTGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGACCAGAGGCTGTTCACCCGCCTGTTCGAAAAAGGCATCATCTACCGCACGAACGGTACCGTGacc                       >  2:747418/1‑148 (MQ=255)                                                                                                                                       ccAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTAtt                                                                                                 >  1:218230/1‑68 (MQ=255)                                                                                                                                       ccAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTAtt                                                                                                 <  2:218230/68‑1 (MQ=255)                                                                                                                                       ccAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCACCCGCCTGTTCGAAAAAGGCATCATCTACCGCAAGAACGGTACCGTGAACTGGGAccc               >  1:274960/1‑150 (MQ=255)                                                                                                                                                gagCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCACCCGCCTGTTCGAAAAAg                                                       <  1:298543/101‑1 (MQ=255)                                                                                                                                                gagCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCACCCGCCTGTTCGAAAAAg                                                       >  2:298543/1‑101 (MQ=255)                                                                                                                                                    ctgggcctggCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCACCCGCCTGTTCGAAAAAGGCATCATCTACCGCAAGAACGGTACCGTGAACTGGAACGCGGCGGACCAGAcc  <  2:274960/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                      CAGCGCATGCTGGGCAAGAACGTCCTGCAGCCGATGGGCTGGGACGCTTTCGGCATGCCCGCGGAAAACGCGGCGATGAAAAACAACGTCGCCCCGGCCAAGTGGACGTACGAAAACATCGACTACATGAAGACCCAGCTCAAGAGCCTGGGCCTGGCCATCGACTGGGCACGTGAAGTCACCACCTGCAAACCCGACTATTACCGTTGGGAGCAGTGGCTGTTCACCCGCCTGTTCGAAAAAGGCATCATCTACCGCAAGAACGGTACCGTGAACTGGGACCCGGCGGACCAGACC  >  NC_002947/5456512‑5456808

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