Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,584,934	(T)7→8	coding (1553/2280 nt)	ydeP ←	putative oxidoreductase YdeP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,584,927	1	.	T	97.9%	152.8 / ‑6.6	48	N520K (AAC→AAA)	ydeP	putative oxidoreductase YdeP
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base T (19/28);  minor base C (1/0);  total (20/28)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCG‑TTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCT  >  NC_000913/1584790‑1585061                                                                                                                                           |                                                                                                                                       ctctGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCG‑tt                                                                                                                                      <  2:418327/140‑1 (MQ=255)  tctGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCG‑ttt                                                                                                                                     >  1:709813/1‑140 (MQ=255)      gggTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCG‑ttttttt                                                                                                                                 >  1:567088/1‑140 (MQ=255)           gTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCG‑ttttttttgcat                                                                                                                            >  2:425527/1‑135 (MQ=255)            tgctgcCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTTGCATg                                                                                                                            >  1:151815/1‑140 (MQ=255)            tgctgcCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGGTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCCGCGCACCGTTTTTTTTGCATg                                                                                                                            >  2:434255/1‑140 (MQ=255)              ctgcctgcGCGATTCGCGCGACGACTGCATACTCTGAATTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTTGCATGTc                                                                                                                          <  1:378687/140‑1 (MQ=255)                   tgcGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGCTTTTTTTTCATGTCAAtat                                                                                                                     >  2:755067/1‑138 (MQ=255)                    gcgcGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTTTCATGTCAAtctc                                                                                                                    >  2:84352/1‑140 (MQ=255)                      gcgATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCAACGCTTTTTTTGCATGTCAATCTCGc                                                                                                                  >  1:266740/1‑140 (MQ=255)                        gATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTTGCATGTCAATCTCGCTa                                                                                                                <  1:74930/140‑1 (MQ=255)                        gATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTTGCATGTCAATCTCGCTa                                                                                                                <  2:197682/140‑1 (MQ=255)                        gATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTTGCATGTCAATCTCGCTa                                                                                                                <  2:567088/140‑1 (MQ=255)                        gATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTTGCATGTCAATCTCGCTa                                                                                                                <  2:531284/140‑1 (MQ=255)                           tCCCGCGGCCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTTGCATGTCAATCTCGCTACGt                                                                                                             >  2:716691/1‑140 (MQ=255)                                  aCCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGa                                                                                                      <  1:505385/140‑1 (MQ=255)                                   ccACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGAc                                                                                                     <  1:84352/140‑1 (MQ=255)                                   ccACTGCACACTCTGATCTCAGGATTACACCGGCGGGTTTTAACACGCTACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGAc                                                                                                     <  2:346541/140‑1 (MQ=255)                                      cTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCgg                                                                                                  <  1:112923/140‑1 (MQ=255)                                      cTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCgg                                                                                                  <  1:716691/140‑1 (MQ=255)                                      cTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCgg                                                                                                  <  2:520531/140‑1 (MQ=255)                                      cTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCgg                                                                                                  <  1:198036/140‑1 (MQ=255)                                       tGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGc                                                                                                 <  2:112493/140‑1 (MQ=255)                                       tGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGc                                                                                                 <  2:158811/140‑1 (MQ=255)                                       tGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGc                                                                                                 <  1:425527/140‑1 (MQ=255)                                        gCACACTCTGATTTCAGAATTAAACCGGCGGTTTTAAACCCGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCa                                                                                                <  2:287136/140‑1 (MQ=255)                                           cacTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGaa                                                                                             <  2:151815/140‑1 (MQ=255)                                             ctctGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCCCAACGGTTACCGCCTGCGCACCGTTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAata                                                                                           <  2:470965/140‑1 (MQ=255)                                             ctctGATTTCAGCATTACACCGGCGGGGTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAata                                                                                           >  1:140746/1‑140 (MQ=255)                                              tctGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTTGAATGTCAATCTCACTGCGTCCCAGGACCGGCAGAatat                                                                                          >  1:213339/1‑140 (MQ=255)                                                        gCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCg                                                                                >  1:708222/1‑140 (MQ=255)                                                         cATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGt                                                                               <  1:755067/140‑1 (MQ=255)                                                              cacCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTTTCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGt                                                                          >  2:45158/1‑140 (MQ=255)                                                                 cggcggGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCaa                                                                       <  1:427056/140‑1 (MQ=255)                                                                   gcggGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACa                                                                     >  2:640565/1‑140 (MQ=255)                                                                   gcggGTTTTAACACGCCACGCGAGGCATGAATCATCGACAATGAATCCTCAACGGGTACCGCCTGCGCACCGTTTTTTTTGCCTGTCAATCTCGCCACCTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACa                                                                     >  1:433940/1‑140 (MQ=255)                                                                          ttAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATgaga                                                              <  2:709813/140‑1 (MQ=255)                                                                                        gAGGCATAAATCATCGACATTGAATCCTCAAGGGTTACCTCCTGCGCACCGTTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAg                                                <  2:343696/140‑1 (MQ=255)                                                                                                      cGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTAcc                                  >  2:768679/1‑140 (MQ=255)                                                                                                         cATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGcc                               <  2:431984/140‑1 (MQ=255)                                                                                                         cATTGAATCCTCAACGGGTACCGCCTGCGCATCGTTTTTTTTGCATGTGAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGcc                               <  2:266740/140‑1 (MQ=255)                                                                                                          aTTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATA‑‑‑‑‑GCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCt                         >  2:598005/1‑140 (MQ=255)                                                                                                            tGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCaaa                            >  1:618096/1‑140 (MQ=255)                                                                                                            tGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCaaa                            >  2:566277/1‑140 (MQ=255)                                                                                                                   tCAACGGTTACCGCCTGCGCACCGTTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACCGATGAGAGGGGTTAAGGTTAGTGGCTACGTGTACCGCCAAATCTAAtt                     >  2:318946/1‑140 (MQ=255)                                                                                                                   tCAACGGTTACCGCCTGCGCACCGTTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAAtt                     <  2:3839/140‑1 (MQ=255)                                                                                                                     aaCGGTTACCGCCTGCGCACCGTTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGc                   >  2:815721/1‑140 (MQ=255)                                                                                                                               gCCTGCGCACCGTTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTa         <  2:140746/140‑1 (MQ=255)                                                                                                                                  tGCGCACCGTTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTATGCTTAGTGACTACGTGCACCGCCAAGTCTAATTGCGTTAACGGTACAg      >  2:283957/1‑140 (MQ=255)                                                                                                                                  tGCGCACCGTTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAg      >  1:659354/1‑140 (MQ=255)                                                                                                                                   gcgcACCGTTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGc                   >  1:18337/1‑126 (MQ=255)                                                                                                                                   gcgcACCGTTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGc                   <  2:18337/126‑1 (MQ=255)                                                                                                                                    cgcACCGTTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAgcg    >  2:699524/1‑140 (MQ=255)                                                                                                                                      cACCGTTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGcgct  <  2:522362/140‑1 (MQ=255)                                                                                                                                           |                                                                                                                                       CTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCG‑TTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCT  >  NC_000913/1584790‑1585061

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