Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	163,281	T→C	L393L (TTA→CTA)	hrpB →	putative ATP‑dependent RNA helicase HrpB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	163,281	0	T	C	100.0%	145.2 / NA	40	L393L (TTA→CTA)	hrpB	putative ATP‑dependent RNA helicase HrpB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/19);  total (21/19)

GGAGATCTTACAAAGCGATCTTTCCGGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACC  >  NC_000913/163154‑163409                                                                                                                                |                                                                                                                                 ggAGATCTTACAAAGCGATCTTTCCGGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGCTACAAATGCTgg                                                                                                                      <  1:870359/140‑1 (MQ=255)   agttCTTACAAGGCGTTCTTTCCGGTTTGTGGATGAATTTCTTGAAGTGGGGAGGCAGCGACCGGGCGGAGATGGGCTGGTTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGCTACAAATGCTgggg                                                                                                                    <  2:485388/137‑1 (MQ=255)   agaTCTTACAAAGCGATCTTTCCGGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGCTACAAATGCTgggg                                                                                                                    >  1:332607/1‑140 (MQ=255)   agaTCTTACAAAGCGATCTTTCCGGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGCTACAAATGCTgg                                                                                                                      <  1:394684/138‑1 (MQ=255)   agaTCTTACAAAGCGATCTTTCCGGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGCTACAAATGCTgg                                                                                                                      >  2:394684/1‑138 (MQ=255)      ttttACAAAGCGTTTTTCCTGGTTGGTTGATGGAATTACGGGAGTGGGGTCGCAGCGACCCGGCGGGGATGGGCTGGTTGGAAGAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGCTACAAATGCTGGGGGCa                                                                                                                 <  2:447974/138‑1 (MQ=255)        ttAAAAAGGGATTTTTCCGGTTTGCTGAGGAAATTACTGAAAGGGGGATGCAGCGATCGGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGCTACAAATGCTGGGGGCACt                                                                                                               <  2:716972/140‑1 (MQ=255)             aaGCGATCTTTCCGGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGCTACAAATGCTGGGGGCACTGGAgg                                                                                                          >  2:504608/1‑140 (MQ=255)               gCGATCTTTCCGGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGCTACAAATGCTGGGGGCACTGGAGGGt                                                                                                        >  1:63923/1‑140 (MQ=255)                cGATCTTTCCGGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGCTACAAATGCTGGGGGCACTGGAgggct                                                                                                       >  1:485388/1‑138 (MQ=255)                 gATCTTTCCGGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGCTACAAATGCTGGGGGCACTGGAGGGTGa                                                                                                      <  1:94974/140‑1 (MQ=255)                 gATCTTTCCGGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGCTACAAATGCTGGGGGCACTGGAGGGTGa                                                                                                      >  1:453716/1‑140 (MQ=255)                 gATCTTTCCGGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGCTACAAATGCTGGGGGCACTGGAGGGTGa                                                                                                      >  2:94974/1‑140 (MQ=255)                       tCCGGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGCTACAAATGCTGGGGGCACTGGAGGGTGAACGGCt                                                                                                <  2:326813/140‑1 (MQ=255)                             ttGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGCTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTgc                                                                                          >  2:269951/1‑140 (MQ=255)                                     ggAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGCTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGc                                                                                  >  1:58449/1‑140 (MQ=255)                                          tACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGCTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCaaaaa                                                                             >  2:152617/1‑140 (MQ=255)                                               cAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGCTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTGGTGCGCAAGGGCAAAAAATGGc                                                                        >  1:742732/1‑140 (MQ=255)                                                 aTGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGCTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAg                                                                      >  1:938477/1‑140 (MQ=255)                                                 aTGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGCTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAg                                                                      >  2:909976/1‑140 (MQ=255)                                                   ggggATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGCTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAgcg                                                                    <  2:453716/140‑1 (MQ=255)                                                       aTGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGCTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAg                                                                      >  1:102341/1‑134 (MQ=255)                                                       aTGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGCTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAg                                                                      <  2:102341/134‑1 (MQ=255)                                                              gATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGCTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGa                                                         >  1:502806/1‑140 (MQ=255)                                                                       cAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGCTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGttt                                                <  1:152617/140‑1 (MQ=255)                                                                            gggcgggcgggATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGCTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCgg                                           <  2:742732/135‑1 (MQ=255)                                                                            gAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGCTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCgg                                           <  1:504608/140‑1 (MQ=255)                                                                            gAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGCTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCgg                                           <  1:909976/140‑1 (MQ=255)                                                                            gAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGCTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCgg                                           <  1:269951/140‑1 (MQ=255)                                                                               ctggctggATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGCTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCaa                                        >  1:603100/1‑140 (MQ=255)                                                                                 ggcgggATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGCTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATg                                      <  2:659178/140‑1 (MQ=255)                                                                                                cAGTAGTGAATCTACTGGCCGCGAAACGTCTGCTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCgaag                       <  2:63923/140‑1 (MQ=255)                                                                                                          tCTACTGGCCGCGAAACGTCTGCTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAgcg                            >  1:732083/1‑125 (MQ=255)                                                                                                          tCTACTGGCCGCGAAACGTCTGCTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAgcg                            <  2:732083/125‑1 (MQ=255)                                                                                                           cTACTGGCCGCGAAACGTCTGCTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAAcgacgacga            <  2:332607/140‑1 (MQ=255)                                                                                                             aCTGGCCGCGAAACGTCTGCTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAg          >  2:895030/1‑140 (MQ=255)                                                                                                             aCTGGCCGCGAAACGTCTGCTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAg          >  1:953042/1‑140 (MQ=255)                                                                                                                     cgAAACGTCTGCTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTAcc  >  1:139632/1‑140 (MQ=255)                                                                                                                           gTCTGCTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAAcgacgacga            >  2:41830/1‑124 (MQ=255)                                                                                                                           gTCTGCTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAAcgacgacga            <  1:41830/124‑1 (MQ=255)                                                                                                                                |                                                                                                                                 GGAGATCTTACAAAGCGATCTTTCCGGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACC  >  NC_000913/163154‑163409

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