Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,485,031	T→C	G657G (GGT→GGC)	hrpA →	ATP‑dependent RNA helicase HrpA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,485,031	0	T	C	100.0%	108.8 / NA	33	G657G (GGT→GGC)	hrpA	ATP‑dependent RNA helicase HrpA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/16);  total (17/16)

GCGAACCGGCGGAGTATCGCGAAATTCACATTGCGTTGCTGACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACG  >  NC_000913/1484901‑1485165                                                                                                                                   |                                                                                                                                       gCGAACCGGCGGAGTATCGCGAAATTCACATTGCGTTGCTGACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCCCGTAACGTGCGTTTCTCCATCCTCCCCGGCTCTGGTTTa                                                                                                                               >  2:326263/1‑140 (MQ=255)    aaCCGGCGGAGTATCGCGAAATTCACATTGCGTTGCTGACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGCTCTGGTTTATTc                                                                                                                            >  1:506489/1‑140 (MQ=255)                  cgcgAAATTCACATTGCGTTGCTGACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGCTCTGGTTTATTCAAAAAACCGCCTaa                                                                                                              <  2:120060/140‑1 (MQ=255)                    cgAAATTCACATTGCGTTGCTGACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGCTCTGGTTTATTCAAAAAACCGCCTAAAt                                                                                                            >  1:604034/1‑140 (MQ=255)                     gAAATTCACATTGCGTTGCTGAACGCTTTTCTTTCCCTTATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGCTCTGGTTTATTCAAAAAACCGCCTAAATg                                                                                                           <  1:326263/140‑1 (MQ=255)                         ttCACATTGCGTTGCTGACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTGCCCGGCTGTGGTTTATACAAAAAACCGCCTAAATGGGTa                                                                                                       >  2:364782/1‑140 (MQ=255)                                tGCGTTGCTGACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGCTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAAtggtgg                                                                                                <  1:210037/140‑1 (MQ=255)                                 gCGTTGCTGACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGCTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGc                                                                                               >  2:124426/1‑140 (MQ=255)                                    ttGCTGACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGCTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGa                                                                                            <  1:166918/140‑1 (MQ=255)                                         gACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGCTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTgg                                                                                       <  1:514485/140‑1 (MQ=255)                                                   cTTTCCCATATCGGCGTGAAAGCTGTCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGCTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAg                                                                             <  1:364782/140‑1 (MQ=255)                                                        gtttatCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGCTCTGGTTTATTCAAATAACCGCCTTAATAGGTAATGGTGGCGGAACGGGTAGAAACCAgccgcc                                                                        >  2:248952/4‑140 (MQ=255)                                                        ccATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGCTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAgccgcc                                                                        >  1:669796/1‑140 (MQ=255)                                                         cATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGCTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCt                                                                       >  1:94867/1‑140 (MQ=255)                                                                              gATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGCTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGcgcg                                                  >  1:631331/1‑140 (MQ=255)                                                                               aTAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGCTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGt                                                 >  2:810982/1‑140 (MQ=255)                                                                                tAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGCTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTa                                                <  1:56272/140‑1 (MQ=255)                                                                                           ataCCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGCTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGaa                                     <  1:124426/140‑1 (MQ=255)                                                                                           ataCCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGCTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGaa                                     <  2:754894/140‑1 (MQ=255)                                                                                             aCCGGCGCACGTAACGCTCGTTTCTCCATCTTCCCGGGCCCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATg                                   <  2:308451/140‑1 (MQ=255)                                                                                               cGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGCTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATggg                                 >  2:523529/1‑140 (MQ=255)                                                                                               cGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGCTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATAGCTGCGCGTATCGACCCGGAATggg                                 >  1:308451/1‑140 (MQ=255)                                                                                                 gcgcACGTAACGCGCGTTTCTCCATCTTCCCCGGCTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTg                               >  1:37422/1‑140 (MQ=255)                                                                                                      cGTAACGCGCGTTTCTCCATCTTCCCCGGCTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGcc                          >  1:621309/1‑140 (MQ=255)                                                                                                             cgcgTTTCTCCATCTTCCCCGGCTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCt                   <  1:453361/140‑1 (MQ=255)                                                                                                             cgcgTTTCTCCATCTTCCCCGGCTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCt                   <  2:669796/140‑1 (MQ=255)                                                                                                               cgTTTCTCCATCTTCCCCGGCTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCa                 >  2:51396/1‑140 (MQ=255)                                                                                                                 tttCTCCATCTTCCCCGGCTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTc                  <  1:178286/137‑1 (MQ=255)                                                                                                                 tttCTCCATCTTCCCCGGCTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTc                  >  2:178286/1‑137 (MQ=255)                                                                                                                   tctcCATCTTCCCCGGCTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCAt             <  2:631331/140‑1 (MQ=255)                                                                                                                          cTTCCCCGGCTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTa      <  2:604034/140‑1 (MQ=255)                                                                                                                            tCCCCGGCTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTaaa    <  2:94867/140‑1 (MQ=255)                                                                                                                              cccGGCTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGAATTGCTGCGCGTATCGACCGGGAATCGATGGAGCCGTTTGCGCAGCATTTAATTAAACg  >  2:405306/1‑140 (MQ=255)                                                                                                                                   |                                                                                                                                       GCGAACCGGCGGAGTATCGCGAAATTCACATTGCGTTGCTGACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGGGGGCGCATTGCTGCGCGTATCGACCCGGAATGGGTGGAGCCAGTTGCTCAGCATTTGATTAAACG  >  NC_000913/1484901‑1485165

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