Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,815,883	+C	coding (667/687 nt)	rph ←	truncated RNase PH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,815,879	1	.	C	100.0%	113.0 / NA	29	E224G (GAA→GGA)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (13/16);  total (13/16)

CAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGA  >  NC_000913/3815741‑3816009                                                                                                                                             |                                                                                                                                   cAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattcgatt‑c                                                                                                                                   <  2:641223/140‑1 (MQ=255) cAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattcgatt‑c                                                                                                                                   <  2:571415/140‑1 (MQ=255) cAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattcgatt‑c                                                                                                                                   <  2:27501/140‑1 (MQ=255)      aaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑ccccct                                                                                                                              >  1:277124/1‑138 (MQ=255)             gcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCCTCGGGCCa                                                                                                                        <  1:142760/140‑1 (MQ=255)                 atatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCCTCGGGCCagag                                                                                                                     <  2:831406/140‑1 (MQ=255)                 atatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCCTCGGGCCAGAGc                                                                                                                    >  1:50835/1‑140 (MQ=255)                      gTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCCTCGGGCCAGAGCcaaca                                                                                                               >  1:640512/1‑140 (MQ=255)                        ttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCCTCGGGCCAGAGCCAACAAg                                                                                                             <  2:708015/140‑1 (MQ=255)                           aTGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCCTCGGGCCAGAGCCAACAAGATg                                                                                                          >  1:834971/1‑140 (MQ=255)                             gCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCCCCGGGCCAAAGCCATCAAGATGAg                                                                                                        <  2:230280/140‑1 (MQ=255)                                ttCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                                                                                     <  2:106221/140‑1 (MQ=255)                                ttCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                                                                                     <  2:581462/140‑1 (MQ=255)                                     tcctcATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCtt                                                                                                <  1:357535/140‑1 (MQ=255)                                        tcATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCAt                                                                                             <  2:614333/140‑1 (MQ=255)                                                 ttttCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACg                                                                                    >  1:925229/1‑140 (MQ=255)                                                    tCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCt                                                                                 >  2:806269/1‑140 (MQ=255)                                                    tCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCt                                                                                 >  2:323867/1‑140 (MQ=255)                                                     cTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTc                                                                                >  1:484991/1‑140 (MQ=255)                                                       aCAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGc                                                                              >  1:208169/1‑140 (MQ=255)                                                           aC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTc                                                                          >  1:631507/1‑140 (MQ=255)                                                              aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTg                                                                        >  1:752548/1‑140 (MQ=255)                                                              aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTg                                                                        <  1:510726/140‑1 (MQ=255)                                                                aaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGcc                                                                      >  2:134620/1‑140 (MQ=255)                                                                aaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGcc                                                                      <  2:640512/140‑1 (MQ=255)                                                                   aaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTc                                                                   <  1:262662/140‑1 (MQ=255)                                                                   aaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTc                                                                   <  1:350133/140‑1 (MQ=255)                                                                          aCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATTGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCa                                                            <  2:277124/140‑1 (MQ=255)                                                                            tcatcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACt                                                          >  2:963311/1‑140 (MQ=255)                                                                                                   tttGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGt                                   <  2:208169/140‑1 (MQ=255)                                                                                                         aGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCAc                             <  1:584414/140‑1 (MQ=255)                                                                                                                        gTCGCTACAATGGATTCGATTCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGt              <  2:774863/140‑1 (MQ=255)                                                                                                                         tCGCTACAATGGATTCGATTCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGtc             >  1:828459/1‑140 (MQ=255)                                                                                                                                 aTGGATTCGATTCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGc     <  1:963311/140‑1 (MQ=255)                                                                                                                                  tGGATTCGATTCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTGCTGCCGTCCCCTGCGCTTCAATGGTGGGGCCGTCTTCGGTCATCACTGCGTTCATGTCGGGCTCTGCGGCa    <  2:925229/140‑1 (MQ=255)                                                                                                                                    gattcgattcCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCaga  >  1:114886/1‑140 (MQ=255)                                                                                                                                             |                                                                                                                                   CAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGA  >  NC_000913/3815741‑3816009

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