Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	163,469	C→T	G455G (GGC→GGT)	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,469	0	C	T	100.0%	172.2 / NA	50	G455G (GGC→GGT)	hrpB	putative ATP‑dependent RNA helicase HrpB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (36/14);  total (36/14)

GACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCG  >  NC_000913/163392‑163545                                                                              |                                                                             gacgacgaAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGggtgtg                                                                           <  1:38893/81‑1 (MQ=255)  acgacgaAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGggtgtgg                                                                          >  2:38967/1‑81 (MQ=255)   cgacgaAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGggtgtggc                                                                         <  1:38896/81‑1 (MQ=255)   cgacgaAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGggtgtggc                                                                         <  2:38894/81‑1 (MQ=255)   cgacgaAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGggtgtggc                                                                         <  2:38895/81‑1 (MQ=255)     acgaAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGggtgtggcgt                                                                       >  1:38969/1‑81 (MQ=255)       gaAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGCTGGGCAATAGTGACCTGGGTGTGGCGttt                                                                     >  2:38947/1‑81 (MQ=255)        aaGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGTGTGGCGtttt                                                                    >  2:38972/1‑81 (MQ=255)        aaGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGTGTGGCGtttt                                                                    >  1:38975/1‑81 (MQ=255)        aaGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGTGTGGCGtttt                                                                    >  1:38959/1‑81 (MQ=255)           ctgctACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGTGTGGCGTTTTcgc                                                                 <  1:38897/81‑1 (MQ=255)                aCCGCGGCAAAGATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAAt                                                            >  2:38989/1‑81 (MQ=255)                  cgcgGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATc                                                           >  1:38964/1‑80 (MQ=255)                   gcgGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATCaa                                                         >  1:38953/1‑81 (MQ=255)                    cgGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAAc                                                        >  1:38950/1‑81 (MQ=255)                     ggaaaaaaTTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAAc                                                        >  2:38963/4‑80 (MQ=255)                     ggCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAAc                                                        >  2:38977/1‑80 (MQ=255)                      gCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAAc                                                        >  2:38986/1‑79 (MQ=255)                      gCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAACCa                                                      <  2:38898/81‑1 (MQ=255)                        aaaaaTTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAACCAGc                                                    <  2:38899/81‑1 (MQ=255)                        aaaaaTTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAACCAGc                                                    >  2:38954/1‑81 (MQ=255)                          aaaTTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAACCAGCCt                                                  <  2:38900/81‑1 (MQ=255)                           aaTTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAACCAGCCTg                                                 >  1:38958/1‑81 (MQ=255)                           aaTTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAACCAGCCTg                                                 >  1:38993/1‑81 (MQ=255)                            aTTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAACCAGCCTgg                                                >  1:39003/1‑81 (MQ=255)                                   ccATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAAc                                         <  1:38901/81‑1 (MQ=255)                                    cATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACg                                        >  1:38960/1‑81 (MQ=255)                                     aTTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGt                                       <  2:38902/81‑1 (MQ=255)                                      ttCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACg                                        >  2:38982/1‑79 (MQ=255)                                       tctcGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAg                                     >  2:38976/1‑81 (MQ=255)                                       tctcGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAg                                     >  1:38957/1‑81 (MQ=255)                                       tctcGAAGAGCCGCCACGGATCGGCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAg                                     >  1:39005/1‑81 (MQ=255)                                          cGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGt                                    >  2:39004/1‑79 (MQ=255)                                           gAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAg                                 >  1:38971/1‑81 (MQ=255)                                           gAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAg                                 >  1:38974/1‑81 (MQ=255)                                            aaGAGCCGCCACGGATGGGCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAg                                 >  1:39009/1‑80 (MQ=255)                                            aaGAGCCGCCACGGATGGGCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGc                                <  1:38903/81‑1 (MQ=255)                                                  cgccACGGATGGGCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGt                          >  2:38962/1‑81 (MQ=255)                                                   gccACGGATGGGCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGt                          >  2:38985/1‑80 (MQ=255)                                                    ccACGGATGGGCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGtt                         >  2:38988/1‑80 (MQ=255)                                                       cGGATGGGCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTaaaa                     <  1:38904/81‑1 (MQ=255)                                                        ggATGGGCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAAc                    <  2:38905/81‑1 (MQ=255)                                                        ggATGGGCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAAc                    >  2:38956/1‑81 (MQ=255)                                                        ggATGGGCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAAc                    >  1:38996/1‑81 (MQ=255)                                                         gATGGGCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACg                   >  1:38965/1‑81 (MQ=255)                                                         gATGGGCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACg                   <  2:38906/81‑1 (MQ=255)                                                           tGGGCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCt                 <  1:38907/81‑1 (MQ=255)                                                            gggCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCt                 >  1:38978/1‑80 (MQ=255)                                                              gCAATAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACATAGTCAGCAACTGTTAAAACGCTTa               >  1:38995/1‑80 (MQ=255)                                                                 aTAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAacg           >  1:38980/1‑81 (MQ=255)                                                                 aTAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAacg           >  1:38981/1‑81 (MQ=255)                                                                  tAGTGACCTGGGTGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAacgt          <  1:38908/81‑1 (MQ=255)                                                                     tGACCTGGGTGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAacgtacg       >  2:38984/1‑81 (MQ=255)                                                                         cTGGGTGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGc   <  1:38909/81‑1 (MQ=255)                                                                         cTGGGTGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGc   <  1:38910/81‑1 (MQ=255)                                                                          tGGGTGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTACACGTACGTGGCg  >  2:38999/1‑81 (MQ=255)                                                                              |                                                                             GACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCG  >  NC_000913/163392‑163545

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑