Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,664,539	C→T	7.9%	F134F (TTC→TTT)	acnA‑II →	aconitate hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,664,539	0	C	T	7.9%	137.1 / 2.5	51	F134F (TTC→TTT)	acnA‑II	aconitate hydratase
Reads supporting (aligned to +/- strand):  ref base C (25/21);  new base T (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.06e-01

AGACTGCGCTGGTCGACCTTGCCGGCCTGCGTGACGCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGA‑CCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAG‑G‑CGTTCAAGAACGTCGACGTGATCCAGCCGGGCAACGGCATCATGCACCAGATCAACCTGGAGAAGATGTCACCGG  >  NC_002947/2664391‑2664684                                                                                                                                                      |                                                                                                                                                    aGACTGCGCTGGTCGACCTTGCCGGCCTGCGTGACGCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGA‑CCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTc                                                                                                                                                     <  7:109/149‑1 (MQ=255)   aCTGCGCTGGTCGACCTTGCCGGCCTGCGTGACGCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGA‑CCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCg                                                                                                                                                    >  7:318871/1‑148 (MQ=255)                aCCTTGCCGGCCTGCGTGACGCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGA‑CCAGTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCgcgc                                                                                                                                       >  2:283717/1‑148 (MQ=255)                 ccTTGCCGGCCTGCGTGACGCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGA‑CCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCa                                                                                                                                     >  3:151991/1‑149 (MQ=255)                  cTTGCCGACCTGCGTGACGCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGA‑CCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCa                                                                                                                                     >  8:119766/1‑148 (MQ=255)                   ttGCCGGCCTGCGTGACGCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGA‑CCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATc                                                                                                                                   >  1:324231/1‑149 (MQ=255)                        ggCCTGCGTGACGCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGA‑CCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCg                                                                                                                                  >  6:345097/1‑145 (MQ=255)                         gCCTGCGTGACGCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGA‑CCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCg                                                                                                                                  >  7:164268/1‑144 (MQ=255)                                  aCGCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGA‑CCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCaa                                                                                                                     <  8:318871/148‑1 (MQ=255)                                      cATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGA‑CCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACg                                                                                                                   >  2:328632/1‑146 (MQ=255)                                         tGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGA‑CCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAg                                                                                                                               <  6:57438/131‑1 (MQ=255)                                         tGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGA‑CCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAg                                                                                                                               >  5:57438/1‑131 (MQ=255)                                           ccGGCAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGA‑CCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGc                                                                                                           <  5:353751/149‑1 (MQ=255)                                           ccGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGA‑CCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGCCCGc                                                                                                           >  3:47566/1‑149 (MQ=255)                                              aCAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGA‑CCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGCCCGCTTc                                                                                                        >  5:62052/1‑149 (MQ=255)                                                  aGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGT‑CCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACt                                                                                                    <  6:138790/149‑1 (MQ=255)                                                    gcggcgACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGA‑CCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTtc                                                                                                  <  4:151991/149‑1 (MQ=255)                                                         gACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGA‑CCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTTGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCaa                                                                                             >  8:294433/1‑149 (MQ=255)                                                         gACCCGGCCCAGGTCAACCCGGTGGTGCCGCTGCAACTGAACGTCGA‑CCACTCGCTGGCCGACGAGTGCGGAGGTATCGACCCGCAGGCCGTTGAAAAGAACCGCCACATCGAAGATCGCCGCACCGCAGACCGCTTCCACtgactcaa                                                                                             >  8:294437/1‑142 (MQ=255)                                                               gCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGA‑CCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGAc                                                                                       <  5:345097/149‑1 (MQ=255)                                                                 ccAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGA‑CCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGAc                                                                                       >  4:133933/1‑147 (MQ=255)                                                                 ccAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGA‑CCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGAc                                                                                       <  3:133933/147‑1 (MQ=255)                                                                     gTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGA‑CCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCaag                                                                                   >  6:353881/1‑147 (MQ=255)                                                                      tCAACCCGGTGGTGCCGGTGCAACTGATCGTCGA‑CCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCaagaag                                                                                <  8:164268/149‑1 (MQ=255)                                                                       cAACCCGGTGGTGCCGGTGCAACTGATCGTCGA‑CCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGTCGCAACGAAGACCGCTTCCACTTCATCAACTGGAcc                                                                                      <  5:243949/142‑1 (MQ=255)                                                                       cAACCCGGTGGTGCCGGTGCAACTGATCGTCGA‑CCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGTCGCAACGAAGACCGCTTCCACTTCATCAACTGGAcc                                                                                      >  6:243949/1‑142 (MQ=255)                                                                        aaCCCGGTGGTGCCGGTGCAACTGATCGTCGA‑CCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAG‑G‑c                                                                            <  2:207931/149‑1 (MQ=255)                                                                        aaCCCGGTGGTGCCGGTGCAACTGATCGTCGA‑CCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAG‑G‑c                                                                            >  5:87366/1‑149 (MQ=255)                                                                         aCCCGGTGGTGCCGGTGCAACTGATCGTCGA‑CCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAG‑G‑Cg                                                                           >  8:313328/1‑149 (MQ=255)                                                                          cccGGTGGTGCCGGTGCAACTGATCGTCGA‑CCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAG‑G‑CGt                                                                          >  3:223477/1‑149 (MQ=255)                                                                           ccGGTGGTGCCGGTGCAACTGATCGTCGATCCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCaagaag                                                                                >  1:27034/1‑145 (MQ=255)                                                                           ccGGTGGTGCCGGTGCAACTGATCGTCGA‑CCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCaagaag                                                                                <  2:27034/144‑1 (MQ=255)                                                                              gtggtgCCGGTGCAACTGATCGTCGA‑CCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAG‑G‑CGTTc                                                                        >  2:272784/1‑147 (MQ=255)                                                                                 gtgCCGGTGCAACTGATCGTCGA‑CCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAG‑G‑CGTTCAAGaa                                                                   <  1:272784/149‑1 (MQ=255)                                                                                       gTGCAACTGATCGTGGA‑CCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCTCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCTCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAG‑G‑CGTTCAAGAACGTCGa                                                             <  7:119766/149‑1 (MQ=255)                                                                                              ttttcgtggA‑CCAATCGCTGGCCGTCGAGTTCGGTGTTTTCGACCCGCAGGCCTTTGAAAAGTACCGCGCCATCGAAGATCGCCGCAACGAAGACGGCTTCCACTTCATCAACTGGAAGAATATG‑G‑CGTTCAAGAACGTCGACGTGATc                                                      <  7:294437/146‑1 (MQ=255)                                                                                                aTCGTCGA‑CCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTTGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAG‑G‑CGTTCAAGAACGTCGACGTGATc                                                      <  7:294433/147‑1 (MQ=255)                                                                                                                   gCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAG‑G‑CGTTCAAGAACGTCGACGTGATCCAGCCGGGCAACGGcatcat                                  >  1:355734/1‑149 (MQ=255)                                                                                                                   gCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAG‑G‑CGTTCAAGAACGTCGACGTGATCCAGCCGGGCAACGGcatca                                   <  7:313328/148‑1 (MQ=255)                                                                                                                              cGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAG‑G‑CGTTCAAGAAc                                                                  <  3:157204/106‑1 (MQ=255)                                                                                                                              cGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAG‑G‑CGTTCAAGAAc                                                                  >  4:157204/1‑106 (MQ=255)                                                                                                                                gtggtTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAG‑G‑CGTTCAAGAACGTCGACGTGATCCAGCCGGGCaa                                           <  2:356221/127‑1 (MQ=255)                                                                                                                                gtggtTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAG‑G‑CGTTCAAGAACGTCGACGTGATCCAGCCGGGCaa                                           >  1:356221/1‑127 (MQ=255)                                                                                                                                gtggtTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAG‑G‑CGTTCAAGAACGTCGACGTGATCCAGCCGGGCAACGGCATCATGCACCAGATCAAc                     <  2:355734/149‑1 (MQ=255)                                                                                                                                   gtTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAG‑G‑CGTTCAAGAACGTCGACGTGATCCAGCCGGGCAACGGCATCATGCACCAGATCAACCTg                  <  1:328632/149‑1 (MQ=255)                                                                                                                                      ttgACCCGCACGGCTTAGAAAAACACCCCGCCAACTAATATCGCCCCAACCCAAACCCGTTACACATCAAAAAAATGCCCAAGAAGTGTCGTTCAAGAACGTCGACGTGATCCAGCCGGGCAACGGCATCATGCACCAGATCAACCTg                  <  1:283717/146‑1 (MQ=255)                                                                                                                                          cccGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAG‑G‑CGTTCAAGAACGTCGACGTGATCCAGCCGGGCAACGGCATCATGCACCAGATCAACCTGGagaaga           <  6:62052/149‑1 (MQ=255)                                                                                                                                                   cTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAG‑G‑CGTTCAAGAACGTCGACGTGATCCAGCCGGGCAACGGCATCATGc                                <  2:316906/119‑1 (MQ=255)                                                                                                                                                   cTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAG‑G‑CGTTCAAGAACGTCGACGTGATCCAGCCGGGCAACGGCATCATGc                                >  1:316906/1‑119 (MQ=255)                                                                                                                                                   cTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAG‑G‑CGTTCAAGAACGTCGACGTGATCCAGCCGGGCAACGGCATCATGCACCAGATCAACCTGGAGAAGATGTCACCgc  >  2:249709/1‑148 (MQ=255)                                                                                                                                                    ttCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAG‑G‑CGTTCAAGAACGTCGACGTGATCCAGCCGGGCAACGGCATCATGCACCAGa                          <  7:83/124‑1 (MQ=255)                                                                                                                                                    ttCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAG‑G‑CGTTCAAGAACGTCGACGTGATCCAGCCGGGCAACGGCATCATGCACCAGa                          >  8:83/1‑124 (MQ=255)                                                                                                                                                      |                                                                                                                                                    AGACTGCGCTGGTCGACCTTGCCGGCCTGCGTGACGCCATTGCCGACAAAGGCGGCGACCCGGCCCAGGTCAACCCGGTGGTGCCGGTGCAACTGATCGTCGA‑CCACTCGCTGGCCGTCGAGTGCGGTGGTTTCGACCCGCAGGCCTTCGAAAAGAACCGCGCCATCGAAGATCGCCGCAACGAAGACCGCTTCCACTTCATCAACTGGACCAAGAAG‑G‑CGTTCAAGAACGTCGACGTGATCCAGCCGGGCAACGGCATCATGCACCAGATCAACCTGGAGAAGATGTCACCGG  >  NC_002947/2664391‑2664684

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