Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	162,927	G→A	A275T (GCT→ACT)	hrpB →	putative ATP‑dependent helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	162,927	0	G	A	100.0%	119.9 / NA	38	A275T (GCT→ACT)	hrpB	putative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/20);  total (18/20)

CAACTGGCTTCGCGCATCGGCAGTGATGTATTGCTCTGCCCGCTGTATGGCGCGTTGTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATG  >  NC_000913/162795‑163055                                                                                                                                     |                                                                                                                                 caacTGGCTTCGCGCATCGGCAGTGATGTATTGCTCTGCCCGCTGTATGGCGCGTTGTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGaaa                                                                                                                            <  1:107201/139‑1 (MQ=255)             cgcATCGGCAGTGATGTATTGCTCTGCCCGCTGTATGGCGCGTTGTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCa                                                                                                                <  2:1371996/139‑1 (MQ=255)                             tATTGCTCTGCCCGCTGTATGGCGCGTTGTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTc                                                                                                 <  2:966700/138‑1 (MQ=255)                             tATTGCTCTGCCCGCTGTATGGCGCGTTGTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTc                                                                                                 >  1:966700/1‑138 (MQ=255)                                                   cgcgTTGTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTaa                                                                                                                   >  2:1591524/1‑98 (MQ=255)                                                   cgcgTTGTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTaa                                                                                                                   <  1:1591524/98‑1 (MQ=255)                                                   cgcgTTGTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCAtt                                                                                                              <  1:726957/103‑1 (MQ=255)                                                   cgcgTTGTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCAtt                                                                                                              >  2:726957/1‑103 (MQ=255)                                                      gTTGTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGAtt                                                                                     <  1:1270914/125‑1 (MQ=255)                                                      gTTGTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGAtt                                                                                     >  2:1270914/1‑125 (MQ=255)                                                        tGTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTa                                                                                                                    <  2:1454214/92‑1 (MQ=255)                                                        tGTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTa                                                                                                                    >  1:1454214/1‑92 (MQ=255)                                                                    aTCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGa                                                         >  1:1493256/1‑139 (MQ=255)                                                                            aaaGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGt                                                                                                  >  1:1351949/1‑90 (MQ=255)                                                                            aaaGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGt                                                                                                  <  2:1351949/90‑1 (MQ=255)                                                                             aaGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCAc                                                <  2:936195/139‑1 (MQ=255)                                                                                   tcctccCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCt                                          >  1:1369702/1‑139 (MQ=255)                                                                                                 cAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCtggtg                                                                                            >  1:1124361/1‑75 (MQ=255)                                                                                                 cAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCtggtg                                                                                            <  2:1124361/75‑1 (MQ=255)                                                                                                 cAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGgc                                                                  <  1:471429/101‑1 (MQ=255)                                                                                                 cAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGgc                                                                  >  2:471429/1‑101 (MQ=255)                                                                                                     ggATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCa                        <  2:1369702/139‑1 (MQ=255)                                                                                                         gcgcAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTc                                                                                                 >  1:1312844/1‑62 (MQ=255)                                                                                                         gcgcAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTc                                                                                                 <  2:1312844/62‑1 (MQ=255)                                                                                                         gcgcAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGAt                                                        <  2:570728/103‑1 (MQ=255)                                                                                                         gcgcAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGAt                                                        >  1:570728/1‑103 (MQ=255)                                                                                                               aGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACTCGTTAGc              <  1:1032584/139‑1 (MQ=255)                                                                                                                   gtgCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCAc                                                >  2:1519024/1‑101 (MQ=255)                                                                                                                   gtgCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCAc                                                <  1:1519024/101‑1 (MQ=255)                                                                                                                    tgCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATtgtg                                                                                  <  2:1256068/66‑1 (MQ=255)                                                                                                                    tgCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATtgtg                                                                                  >  1:1256068/1‑66 (MQ=255)                                                                                                                    tgCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGAt                              >  1:1088333/1‑118 (MQ=255)                                                                                                                    tgCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGAt                              <  2:1088333/118‑1 (MQ=255)                                                                                                                         gCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCa                        <  2:1645419/119‑1 (MQ=255)                                                                                                                         gCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCa                        >  1:1645419/1‑119 (MQ=255)                                                                                                                         gCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAAcgc                    <  2:463714/123‑1 (MQ=255)                                                                                                                         gCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAAcgc                    >  1:463714/1‑123 (MQ=255)                                                                                                                           gACCAATATTACTGAAACCGGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATg  >  2:775152/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                                 CAACTGGCTTCGCGCATCGGCAGTGATGTATTGCTCTGCCCGCTGTATGGCGCGTTGTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATG  >  NC_000913/162795‑163055

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