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%	256.3 / NA	72	C290Y (TGT→TAT)	hrpB	putative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (32/40);  total (32/40)

CCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGCAGCGTGCCGGGCG  >  NC_000913/162878‑163072                                                                                                |                                                                                                    cctccCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGcc                                                                                                 <  2:1330036/100‑1 (MQ=255) cctccCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGcc                                                                                                 >  1:1260986/1‑100 (MQ=255)  ctccCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGccc                                                                                                >  2:260270/1‑100 (MQ=255)  ctccCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGccc                                                                                                >  1:1739858/1‑100 (MQ=255)   tccCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAGCCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGccc                                                                                                <  2:1960166/99‑1 (MQ=255)   tccCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCa                                                                                               >  1:109091/1‑100 (MQ=255)   tccCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCa                                                                                               <  2:1200799/100‑1 (MQ=255)    cccGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAg                                                                                              <  1:126852/100‑1 (MQ=255)    cccGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAg                                                                                              >  1:449241/1‑100 (MQ=255)    cccGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAg                                                                                              <  1:1817800/100‑1 (MQ=255)     ccGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAgg                                                                                             >  2:354054/1‑100 (MQ=255)     ccGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTAGTGGTGGATTATGCCCAgg                                                                                             <  2:2332868/100‑1 (MQ=255)       ggCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAg                                                                                           <  1:656387/100‑1 (MQ=255)            cGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCgtgt                                                                                      >  1:2065161/1‑100 (MQ=255)             gCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCgtgt                                                                                      <  1:1470478/99‑1 (MQ=255)             gCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGtgtg                                                                                     >  2:962397/1‑100 (MQ=255)               aaGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGgc                                                                                   <  1:390610/100‑1 (MQ=255)                aGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGgcg                                                                                  >  2:1204442/1‑100 (MQ=255)                     tGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGttt                                                                             >  1:1136624/1‑100 (MQ=255)                      gcgcAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGtttt                                                                            <  2:459442/100‑1 (MQ=255)                        gcAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGa                                                                          <  2:2088452/100‑1 (MQ=255)                        gcAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGa                                                                          <  2:109899/100‑1 (MQ=255)                          aaaGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTTGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATc                                                                        <  2:1253122/100‑1 (MQ=255)                          aaaGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATc                                                                        <  1:354056/100‑1 (MQ=255)                          aaaGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATc                                                                        <  1:1112011/100‑1 (MQ=255)                           aaGTGGTGCTGGCGACCCATATTGCTGATCCCAGTTTATCCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCTGGAGCGTGTGGCGCGTTTTGATc                                                                        <  2:872394/99‑1 (MQ=255)                            aGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCcg                                                                      <  2:511763/100‑1 (MQ=255)                            aGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCcg                                                                      <  1:184052/100‑1 (MQ=255)                              tggtgCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCcgcg                                                                    <  2:576020/100‑1 (MQ=255)                                gtgCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCcgcg                                                                    >  2:540926/1‑98 (MQ=255)                                  gCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACg                                                                >  2:728022/1‑100 (MQ=255)                                      gcaaCCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGc                                                            <  2:1823531/97‑1 (MQ=255)                                      gCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGc                                                            <  2:1136622/100‑1 (MQ=255)                                          ccAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTAc                                                        >  2:819414/1‑100 (MQ=255)                                             atatTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTAcgcg                                                     <  2:449239/100‑1 (MQ=255)                                               atTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGAc                                                   >  2:1268383/1‑100 (MQ=255)                                                  gCTGACACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGAAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTg                                                 >  2:1252848/1‑99 (MQ=255)                                                  gCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACt                                                  <  2:247068/98‑1 (MQ=255)                                                  gCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGa                                                >  2:644531/1‑100 (MQ=255)                                                    tGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGAtt                                              <  1:342397/100‑1 (MQ=255)                                                       aaCCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACt                                           <  1:1542616/100‑1 (MQ=255)                                                        atctgTTTAATCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTTATACGCGCACGGGGCTTACGCGACTGATTACTc                                          <  2:1539711/96‑1 (MQ=255)                                                        aCCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTc                                          <  1:1114348/100‑1 (MQ=255)                                                         ccAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCa                                         <  1:1204445/100‑1 (MQ=255)                                                           aGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATCATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAAc                                       >  2:493471/1‑100 (MQ=255)                                                            gTTTAACCATTGCAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCagat                                      >  1:1539715/1‑97 (MQ=255)                                                               tAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGt                                   <  2:2086558/100‑1 (MQ=255)                                                                   cATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGc                               <  2:1739853/100‑1 (MQ=255)                                                                   cATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGc                               >  2:515289/1‑100 (MQ=255)                                                                    aTTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGcc                              >  1:944720/1‑100 (MQ=255)                                                                      tcaaGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAg                            <  2:1416977/98‑1 (MQ=255)                                                                      tGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAg                            <  2:2065154/100‑1 (MQ=255)                                                                      tGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAg                            <  2:109089/100‑1 (MQ=255)                                                                       gAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAgg                           >  1:21609/1‑100 (MQ=255)                                                                       gAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAgg                           <  1:644532/100‑1 (MQ=255)                                                                       gAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAga                           >  1:1416981/1‑99 (MQ=255)                                                                           gTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATc                       >  1:1047160/1‑100 (MQ=255)                                                                             aTTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCa                     <  2:944718/100‑1 (MQ=255)                                                                               tCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATg                   >  1:1831325/1‑100 (MQ=255)                                                                               tCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATg                   <  2:848523/100‑1 (MQ=255)                                                                                cGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGa                  >  1:372249/1‑100 (MQ=255)                                                                                 gTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGAc                 >  1:900039/1‑100 (MQ=255)                                                                                   cTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGc               >  2:2288927/1‑100 (MQ=255)                                                                                      gtggtggATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGTCAGGCATCCATGACGCAGc            <  2:2329540/100‑1 (MQ=255)                                                                                       tggtggATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGCAGCg           <  1:515291/100‑1 (MQ=255)                                                                                       tggtggATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGCAGCg           <  2:372247/100‑1 (MQ=255)                                                                                        ggtggATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGCAGCGt          <  1:260272/100‑1 (MQ=255)                                                                                        ggtggATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGCAGCGt          >  2:85586/1‑100 (MQ=255)                                                                                        ggtggATTATGCACAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGCAGCGt          >  1:1535487/1‑100 (MQ=255)                                                                                           ggATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGCAGCGTGcc       >  1:2078554/1‑100 (MQ=255)                                                                                             aTTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAGCGCGTTAGCCAGGCATCCATGACGCAGCGTGCCgg     >  1:1191778/1‑100 (MQ=255)                                                                                             aTTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGCAGCGTGCCgg     >  2:596556/1‑100 (MQ=255)                                                                                                gtgCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGCAGCGTGCCGGgcg  <  1:1293249/100‑1 (MQ=255)                                                                                                |                                                                                                    CCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGCAGCGTGCCGGGCG  >  NC_000913/162878‑163072

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