Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	739,617	T→G	A37A (GCT→GCG)	phr →	deoxyribodipyrimidine photolyase (photoreactivation)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	739,617	0	T	G	100.0%	210.4 / NA	57	A37A (GCT→GCG)	phr	deoxyribodipyrimidine photolyase (photoreactivation)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (30/27);  total (30/27)

GCGCCATTCAGGAGTTTTATGACTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGAAAAAGGTATTCCTTTATTGTTCCGTGAAGTGGATGACTTTGTC  >  NC_000913/739489‑739752                                                                                                                                 |                                                                                                                                        gcgcCATTCAGGAGTTTTATGACTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCa                                                                                                                              >  1:123299/1‑140 (MQ=255) gcgcCATTCAGGAGTTTTATGACTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCa                                                                                                                              >  2:248906/1‑140 (MQ=255)  cgcCATTCAGGAGTTTTATGACTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAg                                                                                                                             >  1:699797/1‑140 (MQ=255)  cgcCATTCAGGAGTTTTATGACTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAg                                                                                                                             >  2:284622/1‑140 (MQ=255)   gcCATTCAGGAGTTTTATGACTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGt                                                                                                                            >  2:296608/1‑140 (MQ=255)      aTTCAGGAGTTTTATGACTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTggg                                                                                                                         <  1:585513/140‑1 (MQ=255)      aTTCAGGAGTTTTATGACTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTACTGGCGTTGTATATCGCGACACCACGCCAGTggg                                                                                                                         <  2:310337/140‑1 (MQ=255)        tCAGGAGTTTTATGACTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCg                                                                                                                       <  1:784253/140‑1 (MQ=255)                 ttATGACTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGa                                                                                                                      >  2:405808/1‑132 (MQ=255)                 ttATGACTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGa                                                                                                                      <  1:405808/132‑1 (MQ=255)                 ttATGACTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAAc                                                                                                              >  1:163657/1‑140 (MQ=255)                 ttATGACTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAAc                                                                                                              >  1:20499/1‑140 (MQ=255)                         aCCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTcgcc                                                                                                      >  2:204723/1‑140 (MQ=255)                         aCCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTcgcc                                                                                                      >  2:346534/1‑140 (MQ=255)                             aTCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGt                                                                                                  <  2:93649/140‑1 (MQ=255)                             aTCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGt                                                                                                  <  2:712911/140‑1 (MQ=255)                                 ggtctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAgg                                                                                              <  1:22474/140‑1 (MQ=255)                                 ggtctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAgg                                                                                              <  2:201406/140‑1 (MQ=255)                                  gtctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGc                                                                                             >  1:381371/1‑140 (MQ=255)                                   tctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCt                                                                                            <  2:123299/140‑1 (MQ=255)                                    ctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTg                                                                                           <  1:49150/140‑1 (MQ=255)                                     tggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGa                                                                                          >  2:194584/1‑140 (MQ=255)                                     tggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGa                                                                                          >  1:549421/1‑140 (MQ=255)                                       gcgaaaccAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAAc                                                                                        >  2:361314/7‑140 (MQ=255)                                         ttCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCcgcg                                                                                                   >  2:142925/1‑127 (MQ=255)                                         ttCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCcgcg                                                                                                   <  1:142925/127‑1 (MQ=255)                                             ccAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCa                                                                                  <  1:481579/140‑1 (MQ=255)                                                ggATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATg                                                                               <  2:153974/140‑1 (MQ=255)                                                      aCGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAAc                                                                         >  2:295584/1‑140 (MQ=255)                                                            gCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAAt                                                                    <  1:671557/139‑1 (MQ=255)                                                            gCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAAt                                                                    >  2:671557/1‑139 (MQ=255)                                                            gCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATg                                                                   >  2:597844/1‑140 (MQ=255)                                                             cACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATgg                                                                  >  1:33422/1‑140 (MQ=255)                                                               cGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGc                                                                <  1:684646/140‑1 (MQ=255)                                                                 aTAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTa                                                              <  2:425058/140‑1 (MQ=255)                                                                  tAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTAc                                                             >  2:92823/1‑140 (MQ=255)                                                                    aTCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAAt                                                                    <  1:466746/131‑1 (MQ=255)                                                                    aTCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAAt                                                                    >  2:466746/1‑131 (MQ=255)                                                                       tcGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATa                                                        >  2:59945/1‑140 (MQ=255)                                                                       tcGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATa                                                        >  1:668173/1‑140 (MQ=255)                                                                          cACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAgcg                                                     >  2:269986/1‑140 (MQ=255)                                                                          cACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAgcg                                                     >  2:357926/1‑140 (MQ=255)                                                                          cACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAgcg                                                     <  2:475035/140‑1 (MQ=255)                                                                              ggCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTg                                                 <  1:685823/140‑1 (MQ=255)                                                                              ggCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTg                                                 <  1:269986/140‑1 (MQ=255)                                                                               gCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGc                                                >  1:407213/1‑140 (MQ=255)                                                                               gCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGc                                                >  1:668168/1‑140 (MQ=255)                                                                                   ccgcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGaa                                            <  1:597844/140‑1 (MQ=255)                                                                                     gcctgccgcAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGaaaa                                          >  1:552755/1‑140 (MQ=255)                                                                                         gccgcAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGAAAAAGGt                                      >  2:599937/1‑140 (MQ=255)                                                                                          ccgcAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGAAAAAGGTa                                     <  2:407213/140‑1 (MQ=255)                                                                                          ccgcAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGAAAAAGGTa                                     <  2:549421/140‑1 (MQ=255)                                                                                           cgcAATTCGACTGCACGCGTGCTGGCGTTGTATAGCGCGACACCACGCCAGTCGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGAAAAAGGTAt                                    <  2:363035/140‑1 (MQ=255)                                                                                             cAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGAAAAAGGTATTc                                  <  2:366911/140‑1 (MQ=255)                                                                                                          cgcgTGCTGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGAAAAAGGTATTCCTTTATTGTTCCg                     <  2:476712/140‑1 (MQ=255)                                                                                                                 tGGCGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGAAAAAGGTATTCCTTTATTGTTCCGTGAAGTg              <  1:296608/140‑1 (MQ=255)                                                                                                                    cGTTGTATATCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGG‑AAAAGGTATTCCTTTATTGTTCCGTGAAGTGGATg          <  1:417882/140‑1 (MQ=255)                                                                                                                             tCGCGACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGAAAAAGGTATTCCTTTATTGTTCCGTGAAGTGGATGACTTTGTc  >  1:306115/1‑140 (MQ=255)                                                                                                                                 |                                                                                                                                        GCGCCATTCAGGAGTTTTATGACTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTGCCGCAATTCGTCTGCACGCGTGCTGGCGTTGTATATCGCTACACCACGCCAGTGGGCGACGCATAACATGTCGCCGCGTCAGGCTGAACTTATCAATGCTCAACTGAATGGGCTACAAATAGCGCTTGCGGAAAAAGGTATTCCTTTATTGTTCCGTGAAGTGGATGACTTTGTC  >  NC_000913/739489‑739752

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