Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,362,215	C→G	10.9%	V49L (GTG→CTG)	tpiA ←	triose phosphate isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,362,215	0	C	G	10.9%	149.6 / 10.0	55	V49L (GTG→CTG)	tpiA	triose phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base C (28/21);  new base G (3/3);  total (31/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 = 7.69e-02

GCTCCGAATGGCCAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAGGGCGACGCATGC  >  NC_002947/5362068‑5362361                                                                                                                                                     |                                                                                                                                                   gCTCCGAATGGCCAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCAcc                                                                                                                                                   >  5:153465/1‑149 (MQ=255)  cTCCGAATGGCCAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCt                                                                                                                                                  <  8:47033/149‑1 (MQ=255)         tGGCCAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCtgattgat                                                                                                                                           <  6:153465/149‑1 (MQ=255)           gCCAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCAGGTGATtgatgg                                                                                                                                         <  4:81780/149‑2 (MQ=255)             cAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGaa                                                                                                                                        >  2:113478/1‑148 (MQ=255)              aaTCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAAc                                                                                                                                       >  6:25349/1‑148 (MQ=255)               aTCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAAc                                                                                                                                       >  6:198797/1‑147 (MQ=255)                tCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAg                                                                                                                                    <  5:25349/149‑1 (MQ=255)                   acacCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGcc                                                                                                                                 >  3:281495/1‑149 (MQ=255)                     acCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCgg                                                                                                                               <  5:245175/149‑1 (MQ=255)                             gCAACCTGCTTCAAACAGCTTCTAGGGAGCAACTTCCCCGGTCAGCGCACACTGTTCGGGTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACATGAGTGATGAAAAAGGCCGGTGGAAATa                                                                                                                       >  2:45516/1‑149 (MQ=255)                              cAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATAc                                                                                                                      >  7:12086/1‑149 (MQ=255)                               aaCCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTTTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACg                                                                                                                     <  5:198797/149‑1 (MQ=255)                                 ccTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCTGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCAc                                                                                                                                                    <  6:204483/116‑1 (MQ=255)                                 ccTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGc                                                                                                                   <  4:281495/149‑1 (MQ=255)                                 ccTGCTTCAACCAGCTGCTCCGGAGCAACTTCCCCTGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCAGCCAGGCCATCAATCAc                                                                                                                                                    >  5:204483/1‑116 (MQ=255)                                       tCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATtgatga                                                                                                                                         >  1:120559/1‑121 (MQ=255)                                       tCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATtgatga                                                                                                                                         <  2:120559/121‑1 (MQ=255)                                         aaCCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCAGGTGATtgatg                                                                                                                                          >  5:145985/1‑118 (MQ=255)                                         aaCCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCAGGTGATtgatg                                                                                                                                          <  6:145985/118‑1 (MQ=255)                                         aaCCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCt                                                                                                                                                  >  4:286847/1‑110 (MQ=255)                                         aaCCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCt                                                                                                                                                  <  3:286847/110‑1 (MQ=255)                                                  ctCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCg                                                                                                  >  4:151548/1‑149 (MQ=255)                                                     gggAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGc                                                                                                <  5:27930/148‑1 (MQ=255)                                                      ggAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCaa                                                                                              <  6:186669/149‑1 (MQ=255)                                                         gCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGc                                                                                           >  1:298030/1‑149 (MQ=255)                                                          cAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGc                                                                                           >  6:289285/1‑148 (MQ=255)                                                          cAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGc                                                                                           >  6:175806/1‑148 (MQ=255)                                                             cTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGAAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTTCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATAAGGCGACTTCCAATCCGCTCGGCTAGGc                                                                                           >  8:286020/1‑145 (MQ=255)                                                                ccccGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTa                                                                                    >  7:247328/1‑149 (MQ=255)                                                                           gcACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCtt                                                                          <  3:151548/148‑1 (MQ=255)                                                                            cACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTg                                                                        >  4:313037/1‑149 (MQ=255)                                                                            cACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATATCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTg                                                                        >  1:233085/1‑149 (MQ=255)                                                                              cccTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGGACAGTAATTAATTTACATGCCAGGCCATCAATCACGTGATTGATGAACAAGGCCGTTGGAAATACGGCGAGTACCAGTCCGATCGGCAAGGCGAGATTAATCAAGCGTTTgg                                                                       >  1:214774/1‑149 (MQ=255)                                                                              cccTGTTCGGGTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGt                                                                      <  1:113478/149‑1 (MQ=255)                                                                               ccTGTTCGGGTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGt                                                                      >  8:189155/1‑148 (MQ=255)                                                                                   ttCGGGTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGtcagct                                                                 >  3:85998/1‑149 (MQ=255)                                                                                       ggTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTCCCAGGCCAACAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGtcagctcagc                                                             >  6:230995/1‑149 (MQ=255)                                                                                         ttGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGa                                                           >  6:73441/1‑149 (MQ=255)                                                                                            tACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACg                                                         <  5:289285/148‑1 (MQ=255)                                                                                            tACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACg                                                         <  5:175806/148‑1 (MQ=255)                                                                                                   gAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGcgcg                                                 <  2:214774/149‑1 (MQ=255)                                                                                                         tgtgCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTAcc                                           >  3:294120/1‑149 (MQ=255)                                                                                                           tgCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTAc                                            <  8:12086/146‑1 (MQ=255)                                                                                                           tgCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGt                                         >  7:97372/1‑149 (MQ=255)                                                                                                              gcC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCa                                      <  5:73441/149‑1 (MQ=255)                                                                                                                  gACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCAttt                                   >  6:80241/1‑149 (MQ=255)                                                                                                                      gTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCa                               <  8:247328/149‑1 (MQ=255)                                                                                                                         aTTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGt                             >  4:161466/1‑148 (MQ=255)                                                                                                                                tACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAg                       <  7:189155/147‑1 (MQ=255)                                                                                                                                  ccTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTa                    >  5:85765/1‑148 (MQ=255)                                                                                                                                          ttcATCAATCAGGTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTc                                                                                 >  3:185021/3‑79 (MQ=39)                                                                                                                                          ttcATCAATCAGGTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTc                                                                                 <  4:185021/77‑1 (MQ=39)                                                                                                                                          gCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAGGGc           >  8:263461/1‑149 (MQ=255)                                                                                                                                             atcaatcaGGTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAGGGCGAc        >  7:250428/1‑149 (MQ=255)                                                                                                                                             atcaatcaCCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAGGGCGAc        <  5:80241/149‑1 (MQ=255)                                                                                                                                                 atcaCCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAGGGCGACGCAt    >  7:255805/1‑149 (MQ=255)                                                                                                                                                    aCCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAGGGCGACGCATGc  <  2:298030/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GCTCCGAATGGCCAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCC‑GACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAGGGCGACGCATGC  >  NC_002947/5362068‑5362361

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