Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,948,369	C→T	S95F (TCC→TTC)	maoP →	macrodomain Ori protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,948,369	0	C	T	100.0%	217.4 / NA	61	S95F (TCC→TTC)	maoP	macrodomain Ori protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/34);  total (27/34)

CGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCAGTAGTTTCTGCAGATGGATGAGCAAACGGTCAATCGCTCGATAGC  >  NC_000913/3948233‑3948501                                                                                                                                         |                                                                                                                                     cgcgAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTcgg                                                                                                                                   <  2:858530/140‑1 (MQ=255)     aGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGt                                                                                                                               >  1:497897/1‑140 (MQ=255)          gTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAAcc                                                                                                                          >  2:648819/1‑140 (MQ=255)          gTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAAcc                                                                                                                          >  1:516756/1‑140 (MQ=255)           ttACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCg                                                                                                                         <  2:497897/140‑1 (MQ=255)           ttACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCg                                                                                                                         <  2:420131/140‑1 (MQ=255)           ttACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCg                                                                                                                         <  2:213261/140‑1 (MQ=255)            tACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTcgg                                                                                                                                   <  1:480134/129‑1 (MQ=255)            tACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTcgg                                                                                                                                   >  2:480134/1‑129 (MQ=255)                 aGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGtt                                                                                                                   >  1:41888/1‑140 (MQ=255)                     agagAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAg                                                                                                               <  2:761560/140‑1 (MQ=255)                      gagaAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAgg                                                                                                              <  1:884726/140‑1 (MQ=255)                        gaAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTg                                                                                                            <  2:840133/140‑1 (MQ=255)                          aaCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTgtt                                                                                                          <  1:1005339/140‑3 (MQ=255)                          aaCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCt                                                                                                          <  2:41888/140‑1 (MQ=255)                           aCTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTg                                                                                                         <  1:883888/140‑1 (MQ=255)                             tctTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGaa                                                                                                       <  1:964378/140‑1 (MQ=255)                              ctTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAg                                                                                                      >  2:203694/1‑140 (MQ=255)                                tCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGAc                                                                                                    <  2:980516/140‑1 (MQ=255)                                 cGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACt                                                                                                   >  1:678946/1‑140 (MQ=255)                                    aGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTaca                                                                                                >  1:746544/1‑140 (MQ=255)                                    aGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAACCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTaca                                                                                                <  2:742460/140‑1 (MQ=255)                                     gCAGTATGCCGTGGCGAACGTGAGCCAGTGATAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTacac                                                                                               <  2:27919/140‑1 (MQ=255)                                     gCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCt                                                                                                          <  2:977577/129‑1 (MQ=255)                                     gCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCt                                                                                                          >  1:977577/1‑129 (MQ=255)                                           tGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTc                                                                                         >  1:511949/1‑140 (MQ=255)                                              cGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGa                                                                                      >  2:484934/1‑140 (MQ=255)                                              cGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGa                                                                                      >  1:381361/1‑140 (MQ=255)                                                tGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACg                                                                                    >  2:209916/1‑140 (MQ=255)                                                 ggCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGa                                                                                   <  1:209916/140‑1 (MQ=255)                                                   cGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGAtt                                                                                 <  1:484934/140‑1 (MQ=255)                                                   cGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGAtt                                                                                 >  2:903137/1‑140 (MQ=255)                                                         tGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATaaa                                                                           >  1:194111/1‑140 (MQ=255)                                                                 tGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGAc                                                                                     <  2:720340/122‑1 (MQ=255)                                                                 tGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGAc                                                                                     >  1:720340/1‑122 (MQ=255)                                                                       aaGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGcc                                                             <  1:984018/140‑1 (MQ=255)                                                                       aaGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGcc                                                             <  1:82120/140‑1 (MQ=255)                                                                       aaGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCCGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGcc                                                             <  2:746544/140‑1 (MQ=255)                                                                           aGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCttt                                                         >  1:553254/1‑140 (MQ=255)                                                                             aaCGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCttttt                                                       >  2:533861/1‑140 (MQ=255)                                                                              aCGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACCCCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCtttttt                                                      <  1:1000038/140‑1 (MQ=255)                                                                               cgcgTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGG‑GCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTa                                                    <  2:516756/140‑1 (MQ=255)                                                                                  gtgtGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATg                                                  <  2:194111/140‑1 (MQ=255)                                                                                       gTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCagt                                             >  1:862972/1‑140 (MQ=255)                                                                                                  tGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCAGTAGTTTCt                                      >  1:698655/1‑136 (MQ=255)                                                                                                  tGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCAGTAGTTTCt                                      <  2:698655/136‑1 (MQ=255)                                                                                                  tGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCAGTAGTTTCTGCAg                                  >  1:41307/1‑140 (MQ=255)                                                                                                  tGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCAGTAGTTTCTGCAg                                  >  2:225343/1‑140 (MQ=255)                                                                                                     cgcgTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGGGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCAGTAGTTTCTGCAgatg                               >  2:135552/1‑140 (MQ=255)                                                                                                       cgTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCAGTAGTTTCTGCAgatgga                             >  2:648060/1‑140 (MQ=255)                                                                                                       cgTATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCAGTAGTTTCTGCAgatgga                             <  1:648060/140‑1 (MQ=255)                                                                                                         tATTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCAGTAGTTTCTGCAgatggatg                           <  2:897063/140‑1 (MQ=255)                                                                                                          aTTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCagta                                            <  2:889458/122‑1 (MQ=255)                                                                                                          aTTAAGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCagta                                            >  1:889458/1‑122 (MQ=255)                                                                                                              aGCGTCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCAGTAGTTTCTGCAGATGGATGAGCaa                      <  2:678946/140‑1 (MQ=255)                                                                                                                  tCCAAAACGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCAGTAGTTTCTGCAGATGGATGAGCAAACgg                  >  2:455132/1‑140 (MQ=255)                                                                                                                     aaaaCGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCAGTAGTTTCTGCAGATGGATGAGCAAACGGTCa               >  2:658481/1‑140 (MQ=255)                                                                                                                       aaCGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCAGTAGTTTCTGCAGATGGATGAGCAAACGGTCAAt             <  1:455132/140‑1 (MQ=255)                                                                                                                         cGTTTTCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCAGTAGTTTCTGCAGATGGATGAGCAAACGGTCAATCg           <  1:658481/140‑1 (MQ=255)                                                                                                                             ttCACACCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCAGTAGTTTCTGCAGATGGATGAGCAAACGGTCAATCGCTCg       <  2:90919/140‑1 (MQ=255)                                                                                                                                  acCCTTTTCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCAGTAGTTTCTGCAGATGGATGAGCAAACGGTCAATCGCTCGATAGc  <  2:511949/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                     CGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACCCTTTCCGGCGGTAAACCGCAGGTTGAAGGTGCTGAAGACTACACCGATTCTGACGATTAATAAAAAGGGCGAAATGCCCTTTTTTTATGTCAGTAGTTTCTGCAGATGGATGAGCAAACGGTCAATCGCTCGATAGC  >  NC_000913/3948233‑3948501

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