Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	162,973	G→A	100%	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%	205.7 / ‑6.7	63	C290Y (TGT→TAT)	hrpB	putative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (29/33);  total (30/33)

ATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGCAGCGTGC  >  NC_000913/162876‑163066                                                                                                  |                                                                                              aTCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATtatg                                                                                             <  2:464978/100‑1 (MQ=255)  tcctccCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGc                                                                                            <  2:562905/100‑1 (MQ=255)   cctccCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGcc                                                                                           <  1:2056144/100‑1 (MQ=255)   cctccCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTCCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGcc                                                                                           <  1:2405088/100‑1 (MQ=255)     tccCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAATGTATTCGTCTGGTGGTGGATTATGccc                                                                                          >  1:1144417/1‑99 (MQ=255)      cccGTCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGccc                                                                                          <  1:664039/98‑1 (MQ=255)      cccGGCACCGCAAGGGATGCGCAAAGTGGGGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCa                                                                                         <  1:138990/99‑1 (MQ=255)       ccGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAgg                                                                                       <  1:2102121/100‑1 (MQ=255)         ggCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAg                                                                                     >  2:1264526/1‑100 (MQ=255)           cACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCg                                                                                   <  1:963427/100‑1 (MQ=255)             ccGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCgtg                                                                                 <  1:1401694/100‑1 (MQ=255)              cacaatgtaTGCTCAGAGTGTTTCTTGCGACTAATATTTCTCAAACCAGTTTAACCAGTGAATGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCgt                                                                                  <  1:2247524/90‑1 (MQ=255)              cGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCgtgt                                                                                >  1:1499316/1‑100 (MQ=255)               gtaaGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGtgtg                                                                               >  1:560961/3‑100 (MQ=255)                cAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTgg                                                                              <  1:208106/100‑1 (MQ=255)                cAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTgg                                                                              <  2:933344/100‑1 (MQ=255)                       tGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAAGCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGttt                                                                       <  2:333390/100‑1 (MQ=255)                        gcgaAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGtttt                                                                      <  1:690976/100‑1 (MQ=255)                         cgcAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGtttt                                                                      >  1:2001639/1‑99 (MQ=255)                           cAAAGTGGTGCTGGCTACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGAt                                                                   <  1:1908414/100‑1 (MQ=255)                               gtggtgCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCcgc                                                               >  2:663046/1‑100 (MQ=255)                                tggtgTTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTCTGATCcg                                                                <  1:1255957/98‑1 (MQ=255)                                 ggtgCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCgcgc                                                             >  1:1175128/1‑100 (MQ=255)                                    gCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACg                                                          <  1:648154/100‑1 (MQ=255)                                     cTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCAc                                                           <  1:935246/98‑1 (MQ=255)                                     cTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGCATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACgg                                                         >  2:1302617/1‑100 (MQ=255)                                       ggCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACgggg                                                       <  2:1830067/100‑1 (MQ=255)                                         ccatcAATATTGCTGAAACCGGTTTAACCGTTGCAGGTATTCGTCTGGGGGTGGAGTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCTCGCACggggg                                                      <  2:1579763/95‑2 (MQ=255)                                         cGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGc                                                      <  1:291777/99‑1 (MQ=255)                                           aCCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTa                                                   >  1:265894/1‑100 (MQ=255)                                            ccAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTAc                                                  >  2:977588/1‑100 (MQ=255)                                                 atTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGAc                                             >  1:1395185/1‑100 (MQ=255)                                                 atTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGAc                                             <  2:1818437/100‑1 (MQ=255)                                                    gCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGa                                          >  2:2338664/1‑100 (MQ=255)                                                    gCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGa                                          >  2:2122196/1‑100 (MQ=255)                                                    gCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGa                                          >  1:2171435/1‑100 (MQ=255)                                                    gCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGa                                          >  1:1242461/1‑100 (MQ=255)                                                    gCTGAAACCAGGTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTAcgcg                                               <  1:509545/95‑1 (MQ=255)                                                       gAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTa                                       <  1:738869/100‑1 (MQ=255)                                                          aCCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTc                                    >  2:1564075/1‑100 (MQ=255)                                                          aCCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTc                                    >  2:925753/1‑100 (MQ=255)                                                          aCCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTAATc                                    <  2:706279/100‑1 (MQ=255)                                                              gTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAAcg                                <  1:2332001/100‑1 (MQ=255)                                                              gTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAAcg                                >  2:2331973/1‑100 (MQ=255)                                                                ttgaCCATTGACGGCATTCGTCTGGTGGTGTATTATGCCCAGGAGCGTGTGGCGCGTTTTGGTCCGCGGACGGGGCTTTCGCGACTGATTACTCAAcgcg                              <  2:1349945/97‑1 (MQ=255)                                                                     cATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGc                         >  1:1900983/1‑100 (MQ=255)                                                                      aTTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGcc                        <  2:177116/100‑1 (MQ=255)                                                                       ttGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCa                       <  2:281889/100‑1 (MQ=255)                                                                        tGAAGGTATTCGTCGGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAg                      <  1:873830/100‑1 (MQ=255)                                                                         gAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAgg                     >  2:1082379/1‑100 (MQ=255)                                                                         gAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAgg                     >  2:2117848/1‑100 (MQ=255)                                                                         gAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAgg                     >  1:864365/1‑100 (MQ=255)                                                                          aaGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGc                    <  2:1399691/100‑1 (MQ=255)                                                                             gTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATc                 >  2:519532/1‑100 (MQ=255)                                                                              tATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATcc                <  2:1134812/100‑1 (MQ=255)                                                                              tATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATcc                >  1:504496/1‑100 (MQ=255)                                                                                   gTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGAc           <  1:1355832/100‑1 (MQ=255)                                                                                     cTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGc         >  2:1065880/1‑100 (MQ=255)                                                                                       ggtggtggATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGAAg       >  2:1548073/1‑100 (MQ=255)                                                                                        gtggtggATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGCAGc      <  2:2276079/100‑1 (MQ=255)                                                                                         tggtggATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTCCGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGCAGCg     >  2:1559622/1‑100 (MQ=255)                                                                                           gtggATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGCAGCGTg   >  1:1737392/1‑100 (MQ=255)                                                                                           gtggATTATGCCCAGGAGCGCGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGCAGCGTg   >  2:1466869/1‑100 (MQ=255)                                                                                            tggATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGCAGCGTGc  >  2:764406/1‑100 (MQ=255)                                                                                                  |                                                                                              ATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGCAGCGTGC  >  NC_000913/162876‑163066

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