Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	132,875	A→G	T421A (ACC→GCC)	acnB →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	132,875	0	A	G	96.4%	175.5 / ‑2.7	55	T421A (ACC→GCC)	acnB	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (25/28);  total (26/29)
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.99e-01

TCGCTGAGAGCGATCGCGGCTTCTCGCTGGCGCAAAAAATGGTAGGCCGTGCCTGTGGCGTGAAAGGCATTCGTCCGGGCGCGTACTGTGAACCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGACCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGCAGTCTTTCTGCCACACCGCGGCGTATCCGAAGCCAGTTGACGTGAACACGCACCACACGCTGCCGGACTT  >  NC_000913/132738‑133002                                                                                                                                          |                                                                                                                                tCGCTGAGAGCGATCGCGGCTTCTCGCTGGCGCAAAAAATGGTAGGCCGTGCCTGTGGCGTGAAAGGCATTCGTCCGGGCGCGTACTGTGAACCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGcc                                                                                                                               >  2:192915/1‑140 (MQ=255)     tGAGAGCGATCGCGGCTTCTCGCTGGCGCAAAAAATGGTAGGCCGTGCCTGTGGCGTGAAAGGCATTCGTCCGGGCGCGTACTGTGAACCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGtg                                                                                                                           <  2:800912/140‑1 (MQ=255)       agagCGATCGCGGCTTCTCGCTGGCGCAAAAAATGGTAGGCCGTGCCTGTGGCGTGAAAGGCATTCGTCCGGGCGCGTACTGTGAACCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGtgat                                                                                                                         <  2:467699/140‑1 (MQ=255)        gagCGATCGCGGCTTCTCGCTGGCGCAAAAAATGGTAGGCCGTGCCTGTGGCGTGAAAGGCATTCGTCCGGGCGCGTACTGTGAACCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGtgatg                                                                                                                        <  2:956855/140‑1 (MQ=255)         agCGATCGCGGCTTCTCGCTGGCGCAAAAAATGGTAGGCCGTGCCTGTGGCGTGAAAGGCATTCGTCCGGGCGCGTACTGTGAACCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGtgatga                                                                                                                       <  2:1162051/140‑1 (MQ=255)             aTCGCGGCTTCTCGCTGGCGCAAAAAATGGTAGGCCGTGCCTGTGGCGTGAAAGGCATTCGTCCGGGCGCGTACTGTGAACCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGTGATGAACTg                                                                                                                   >  1:564004/1‑140 (MQ=255)              tCGCGGCTTCTCGCTGGCGCAAAAAATGGTAGGCCGTGCCTGTGGCGTGAAAGGCATTCGTCCGGGCGCGTACTGTGAACCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGtg                                                                                                                           <  2:944986/131‑1 (MQ=255)              tCGCGGCTTCTCGCTGGCGCAAAAAATGGTAGGCCGTGCCTGTGGCGTGAAAGGCATTCGTCCGGGCGCGTACTGTGAACCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGtg                                                                                                                           >  1:944986/1‑131 (MQ=255)              tCGCGGCTTCTCGCTGGCGCAAAAAATGGTAGGCCGTGCCTGTGGCGTGAAAGGCATTCGTCCGGGCGCGTACTGTGAACCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGtg                                                                                                                           >  1:568060/1‑131 (MQ=255)              tCGCGGCTTCTCGCTGGCGCAAAAAATGGTAGGCCGTGCCTGTGGCGTGAAAGGCATTCGTCCGGGCGCGTACTGTGAACCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGtg                                                                                                                           <  2:568060/131‑1 (MQ=255)              tCGCGGCTTCTCGCTGGCGCAAAAAATGGTAGGCCGTGCCTGTGGCGTGAAAGGCATTCGTCCGGGCGCGTACTGTGAACCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGTGATGAACTGa                                                                                                                  <  1:747214/140‑1 (MQ=255)              tCGCGGCTTCTCGCTGGCGCAAAAAATGGTAGGCCGTGCCTGTGGCGTGAAAGGCATTCGTCCGGGCGCGTACTGTGAACCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGTGATGAACTGa                                                                                                                  <  1:415303/140‑1 (MQ=255)                 cgGCTTCTCGCTGGCGCAAAAAATGGTAGGCCGTGCCTGTGGCGTGAAAGGCATTCGTCCGGGCGCGTACTGTGAACCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGTGATGAACTGAAAg                                                                                                               >  2:669736/1‑140 (MQ=255)                       ctcGCTGGCGCAAAAAATGGTAGGCCGTGCCTGTGGCGTGAAAGGCATTCGTCCGGGCGCGTACTGTGAACCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGTGATGAACTGAAAGACCTgg                                                                                                         <  1:888780/140‑1 (MQ=255)                              gcgcAAAAAATGGTAGGCCGTGCCTGTGGCGTGAAAGGCATTCGTCCGGGCGCGTACTGTGAACCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCt                                                                                                  >  2:776157/1‑140 (MQ=255)                                    aaaaTGGTAGGCCGTGCCTGTGGCGTGAAAGGCATTCGTCCGGGCGCGTACTGTGAACCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCtt                                                                                            <  1:402873/140‑1 (MQ=255)                                             ggCCGTGCCTGTGGCGTGAAAGGCATTCGTCCGGGCGCGTACTGTGAACCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGa                                                                                   <  1:889190/140‑1 (MQ=255)                                                cGTGCCTGTGGCGTGAAAGGCATTCGTCCGGGCGCGTACTGTGAACCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCt                                                                                >  1:1100012/1‑140 (MQ=255)                                                  tGCCTGTGGCGTGAAAGGCATTCGTCCGGGCGCGTACTGTGAACCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTgg                                                                              >  1:1009578/1‑140 (MQ=255)                                                        tgGCGTGAAAGGCATTCGTCCGGGCGCGTACTGTGAACCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGc                                                                        >  1:64852/1‑140 (MQ=255)                                                        tgGCGTGAAAGGCATTCGTCCGGGCGCGTACTGTGAACCGAAAATGACTTCTGGAGGTTCCCAGGACACCACCGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGc                                                                        >  2:845604/1‑140 (MQ=255)                                                             tGAAAGGCATTCGTCCGGGCGCGTACTGTGAACCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGCAGTCt                                                                   >  1:691819/1‑140 (MQ=255)                                                             tGAAAGGCATTCGTCCGGGCGCGTACTGTGAACCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGCAGTCt                                                                   <  1:556302/140‑1 (MQ=255)                                                             tGAAAGGCATTCGTCCGGGCGCGTACTGTGAACCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGCAGTCt                                                                   >  1:877984/1‑140 (MQ=255)                                                             tGAAAGGCATTCGTCCGGGCGCGTACTGTGAACCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGCAGTCt                                                                   <  1:192915/140‑1 (MQ=255)                                                             tGAAAGGCATTCGTCCGGGCGCGTACTGTGAACCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGCAGTCt                                                                   >  2:370982/1‑140 (MQ=255)                                                             tGAAAGGCATTCGTCCGGGCGCGTACTGTGAACCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGCAGTCt                                                                   <  2:139243/140‑1 (MQ=255)                                                             tGAAAGGCATTCGTCCGGGCGCGTACTGTGAACCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGCAGTCt                                                                   <  2:139244/140‑1 (MQ=255)                                                                aaGGCATTCGTCCGGGCGCGTACTGTGAACCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGCAGTCTTTc                                                                <  1:284236/140‑1 (MQ=255)                                                                     aTTCGTCCGGGCGCGTACTGTGAACCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGCAGTCTTTCTGCca                                                           >  1:796622/1‑140 (MQ=255)                                                                          tCCGGGCGCGTACTGTGAACCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGCAGTCTTTCTGCCACACcg                                                      >  2:561624/1‑140 (MQ=255)                                                                          tCCGGGCGCGTACTGTGAACCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGCAGTCTTTCTGCCACACcg                                                      <  2:1009578/140‑1 (MQ=255)                                                                                         tgAACCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGCAGTCTTTCTGCCACACCGCGGCGTATCCGAAGc                                       >  2:907583/1‑140 (MQ=255)                                                                                         tgAACCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGCGATGCAGTCTTTCTGCCACACCGCGGCGTATCCGAAGc                                       >  2:988602/1‑140 (MQ=255)                                                                                            aCCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGCAGTCTTTCTGCCACACcg                                                      >  2:50377/1‑122 (MQ=255)                                                                                            aCCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGCAGTCTTTCTGCCACACcg                                                      <  1:50377/122‑1 (MQ=255)                                                                                              cGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGCAGTCTTTCTGCCACACCGCGGCGTATCCGAAGCCAGtt                                  <  2:4708/140‑1 (MQ=255)                                                                                              cGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGCAGTCTTTCTGCCACACCGCGGCGTATCCGAAGCCAGtt                                  <  1:669736/140‑1 (MQ=255)                                                                                                aaaaTGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGCAGTCTTTCTGCCACACCGCGGCGTATCCGAAGCCAGTTGa                                >  1:827806/1‑140 (MQ=255)                                                                                                 aaaTGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGACCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGCAGTCTTTCTGCCACACCGCGGCGTATCCGAAGCCAGTTGAc                               >  2:159047/1‑140 (MQ=255)                                                                                                  aaTGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGCAGTCTTTCTGCCACACCGCGGCGTATCCGAAGCCAGTTGACg                              <  2:523501/140‑1 (MQ=255)                                                                                                   aTGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGCAGTCTTTCTGCCACACCGCGGCGTATCCGAAGCCAGTTGACGt                             >  2:96619/1‑140 (MQ=255)                                                                                                     gACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGCAGTCTTTCTGCCACACCGCGGCCTATCCGAAGCCAGTTGACGTGa                           <  1:370982/140‑1 (MQ=255)                                                                                                        ttCTGTAGGTTCCCAGGACACCACCGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGCAGTCTTTCTGCCACACCGCGGCGTATCCGAAGCCAGTTGACGTGAaca                        <  2:1100012/140‑1 (MQ=255)                                                                                                               ggTTCCCAGGACACCACCGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGCAGTCTTTCTGCCACACCGCGGCGTATCCGAAGCCAGTTGACGTGAACACGcacca                 <  1:845604/140‑1 (MQ=255)                                                                                                                  tCCCAGGACACCACCGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGCAGTCTTTCTGCCACACCGCGGCGTATCCGAAGCCAGTTGACGTGAACACGCACCacac              >  2:734065/1‑140 (MQ=255)                                                                                                                     cAGGACACCACCGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGCAGTCTTTCTGCCACACCGCGGCGTATCCGAAGCCAGTTGACGTGAACACGCACCACACGCt           >  1:245962/1‑140 (MQ=255)                                                                                                                       ggACACCACCGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGCAGTCTTTCTGCCACACCGCGGCGTATCCGAAG‑CAGTTGACGTGAACACGCACCACACGCTGcc        >  2:299790/1‑140 (MQ=255)                                                                                                                         acacCACCGGCCCGATGACCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGCAGTCTTTCTGCCACACCGCGGCGTATCCGAAGCCAGTTGACGTGAACACGCACCACACGCTGCCg       <  1:159047/140‑1 (MQ=255)                                                                                                                           accaccGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGCAGTCTTTCTGCCACACCGCGGCGTATCCGAAGCCAGTTGACGTGAACACGCACCACACGCTGCCGGa     <  1:96619/140‑1 (MQ=255)                                                                                                                             caccGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGCAGTCTTTCTGCCACACCGCGGCGTATCCGAAGCCAGTTGACGTGAACACGCACCACACGCTGCCGGACt   <  2:64852/140‑1 (MQ=255)                                                                                                                             caccGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGCAGTCTTTCTGCCACACCGCGGCGTATCCGAAGCCAGTTGACGTGAACACGCACCACACGCTGCCGGACt   <  2:877984/140‑1 (MQ=255)                                                                                                                             caccGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGCAGTCTTTCTGCCACACCGCGGCGTATCCGAAGCCAGTTGACATGAACACGCACCACACGCTGCCGGACt   >  1:178794/1‑140 (MQ=255)                                                                                                                              accGGCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGCAGTCTTTCTGCCACACCGCGGCGTATCCGAAGCCAGTTGACGTGAACACGCACCACACGCTGCCGGACtt  >  2:752761/1‑140 (MQ=255)                                                                                                                                  gCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGCAGTCTTTCTGCCACACCGCGGCGTATCCGAAGCCAGTTGACGTGAACACGCACCACACGCTGcc        >  2:1157831/1‑130 (MQ=255)                                                                                                                                  gCCCGATGGCCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGCAGTCTTTCTGCCACACCGCGGCGTATCCGAAGCCAGTTGACGTGAACACGCACCACACGCTGcc        <  1:1157831/130‑1 (MQ=255)                                                                                                                                          |                                                                                                                                TCGCTGAGAGCGATCGCGGCTTCTCGCTGGCGCAAAAAATGGTAGGCCGTGCCTGTGGCGTGAAAGGCATTCGTCCGGGCGCGTACTGTGAACCGAAAATGACTTCTGTAGGTTCCCAGGACACCACCGGCCCGATGACCCGTGATGAACTGAAAGACCTGGCGTGCCTGGGCTTCTCGGCTGACCTGGTGATGCAGTCTTTCTGCCACACCGCGGCGTATCCGAAGCCAGTTGACGTGAACACGCACCACACGCTGCCGGACTT  >  NC_000913/132738‑133002

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑