Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	162,973	G→A	C290Y (TGT→TAT)	hrpB →	putative ATP‑dependent helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	162,973	0	G	A	100.0%	275.8 / NA	76	C290Y (TGT→TAT)	hrpB	putative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/48);  total (28/48)

TCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGCAGCGTGCC  >  NC_000913/162877‑163067                                                                                                 |                                                                                               tcctccCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGc                                                                                             <  1:1056416/100‑1 (MQ=255) tcctccCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGc                                                                                             <  2:2007916/100‑1 (MQ=255)   ctccCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTCGTGGATTATGccc                                                                                           <  2:1266282/100‑1 (MQ=255)   ctccCGGCACCGCAAGGGAAGCGCAAAGTGGTGCTGGCGCCCAATATTGCGGAAACCAGTTTAACCATTGAAGGTATTCGTCTGACGGTGGATTATGcac                                                                                           >  2:1287737/1‑98 (MQ=255)    tccCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCa                                                                                          <  1:325138/100‑1 (MQ=255)     cccGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAg                                                                                         >  2:35550/1‑100 (MQ=255)      ccGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAgg                                                                                        >  2:1410566/1‑100 (MQ=255)        ggCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAg                                                                                      >  2:729730/1‑100 (MQ=255)        ggCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAg                                                                                      >  2:1747105/1‑100 (MQ=255)            ccGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCgtg                                                                                  <  2:1901197/100‑1 (MQ=255)            ccGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCgtg                                                                                  <  1:65327/100‑1 (MQ=255)            ccGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCgtg                                                                                  <  1:65297/100‑1 (MQ=255)            ccGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCgtg                                                                                  <  1:65273/100‑1 (MQ=255)            ccGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCgtg                                                                                  <  1:65248/100‑1 (MQ=255)             cGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCgtgt                                                                                 >  1:1863658/1‑100 (MQ=255)              gCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATGTTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGtgtg                                                                                <  2:1082545/100‑1 (MQ=255)                 aGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGgcg                                                                             <  1:729733/100‑1 (MQ=255)                 aGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGgcg                                                                             >  2:407080/1‑100 (MQ=255)                     aTGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGtt                                                                         >  2:2022850/1‑100 (MQ=255)                       gcgcAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGtttt                                                                       <  1:29238/100‑1 (MQ=255)                       gcgcAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGtttt                                                                       <  2:1991189/100‑1 (MQ=255)                       gcgcAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGtttt                                                                       <  1:1597036/100‑1 (MQ=255)                         gcAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGa                                                                     <  1:1533018/100‑1 (MQ=255)                         gcAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGa                                                                     <  2:1674296/100‑1 (MQ=255)                            aaGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATcc                                                                  >  2:1619128/1‑100 (MQ=255)                             aGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCcg                                                                 <  2:722144/100‑1 (MQ=255)                             aGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCcg                                                                 <  1:1247768/100‑1 (MQ=255)                              gtggtgCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCcgc                                                                >  2:188747/1‑100 (MQ=255)                               tggtgCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCcgcg                                                               <  2:291367/100‑1 (MQ=255)                                ggtgCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCgcgc                                                              >  2:1592436/1‑100 (MQ=255)                                 gtgCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCa                                                             <  1:188748/100‑1 (MQ=255)                                  tgCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCAc                                                            >  1:620975/1‑100 (MQ=255)                                     tGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACggg                                                         >  1:1981590/1‑100 (MQ=255)                                      ggCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACgggg                                                        >  2:1020793/1‑100 (MQ=255)                                      ggCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACgggg                                                        <  1:1981828/100‑1 (MQ=255)                                       gCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGc                                                       <  2:1981572/100‑1 (MQ=255)                                       gCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGc                                                       <  1:2022868/100‑1 (MQ=255)                                       gCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGc                                                       <  1:1410575/100‑1 (MQ=255)                                          aCCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTa                                                    <  2:1320141/100‑1 (MQ=255)                                           ccAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTAc                                                   <  2:1740179/100‑1 (MQ=255)                                               tatTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGa                                               <  1:1992191/100‑1 (MQ=255)                                                atTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGAc                                              >  2:1055149/1‑100 (MQ=255)                                                 ttGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACt                                             <  2:653424/100‑1 (MQ=255)                                                  tGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTg                                            <  2:2045696/100‑1 (MQ=255)                                                  tGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTg                                            <  2:439486/100‑1 (MQ=255)                                                  tGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTg                                            <  2:960475/100‑1 (MQ=255)                                                  tGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTg                                            <  1:478586/100‑1 (MQ=255)                                                    cTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGAt                                          >  1:1286340/1‑100 (MQ=255)                                                     tGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGAtt                                         <  1:1747119/100‑1 (MQ=255)                                                       aaaCCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTAc                                       <  1:1748230/100‑1 (MQ=255)                                                        aaCCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACt                                      >  2:712929/1‑100 (MQ=255)                                                        aaCCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACt                                      <  1:408869/100‑1 (MQ=255)                                                         aCCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTc                                     <  2:624078/100‑1 (MQ=255)                                                          ccAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCa                                    >  2:1694858/1‑100 (MQ=255)                                                           cAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCaa                                   >  1:632016/1‑100 (MQ=255)                                                           cAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCaa                                   <  2:1045684/100‑1 (MQ=255)                                                            aGTTTAACCATTGAAGGTATTCGTCTGGTGGTGCATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAAc                                  <  2:1702780/100‑1 (MQ=255)                                                             gTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAAcg                                 <  2:234575/100‑1 (MQ=255)                                                             gTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAAcg                                 <  1:1158307/100‑1 (MQ=255)                                                               ttAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAAcgcg                               >  1:1958494/1‑100 (MQ=255)                                                                 aaCCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGtt                             >  1:1066073/1‑100 (MQ=255)                                                                   ccATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAg                           >  2:280470/1‑100 (MQ=255)                                                                     aTTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGcc                         <  2:1440388/100‑1 (MQ=255)                                                                        gAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTGAACGCGTTAGCCAgg                      >  1:2008243/1‑100 (MQ=255)                                                                        gAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAgg                      <  2:1286334/100‑1 (MQ=255)                                                                             tATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATcc                 >  2:770011/1‑100 (MQ=255)                                                                              aTTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCa                >  2:1184232/1‑100 (MQ=255)                                                                               ttCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCAt               <  1:1592448/100‑1 (MQ=255)                                                                                 cGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGa             <  2:892335/100‑1 (MQ=255)                                                                                  gTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGAc            <  1:1619141/100‑1 (MQ=255)                                                                                  gTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGAc            <  1:1088169/100‑1 (MQ=255)                                                                                     tggtggtggATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGCa         <  2:1721108/100‑1 (MQ=255)                                                                                       gtggtggATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGCAGc       >  1:823994/1‑100 (MQ=255)                                                                                       gtggtggATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGCAGc       <  1:1020798/100‑1 (MQ=255)                                                                                         ggtggATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGCAGCGt     >  2:1064780/1‑100 (MQ=255)                                                                                            ggATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGCAGCGTGcc  >  2:1509465/1‑100 (MQ=255)                                                                                                 |                                                                                               TCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGCAGCGTGCC  >  NC_000913/162877‑163067

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