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_000913162,9730GA98.1% 180.9 / ‑4.0 53C290Y (TGT→TAT) hrpBputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (28/24);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

ATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGCAGCGTGCCGG  >  NC_000913/162876‑163069
                                                                                                 |                                                                                                
aTCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATtatg                                                                                                <  1:1463376/100‑1 (MQ=255)
  cctccCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGcc                                                                                              >  2:939996/1‑100 (MQ=255)
   ctccCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGccc                                                                                             <  1:563874/100‑1 (MQ=255)
     cccGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGAGGTGGATTATGCCCAg                                                                                           >  2:415378/1‑100 (MQ=255)
              gCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGtgtg                                                                                  <  1:532563/100‑1 (MQ=255)
              gCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGtgtg                                                                                  <  1:1559662/100‑1 (MQ=255)
                    gATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGt                                                                            >  2:330210/1‑100 (MQ=255)
                      tGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGttt                                                                          >  1:586404/1‑100 (MQ=255)
                      tGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGttt                                                                          >  1:235338/1‑100 (MQ=255)
                      tGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGttt                                                                          >  2:373136/1‑100 (MQ=255)
                       gcgcAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGtttt                                                                         <  1:490103/100‑1 (MQ=255)
                       gcgcAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGtttt                                                                         <  2:849625/100‑1 (MQ=255)
                        cgcAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTg                                                                        <  2:1243008/100‑1 (MQ=255)
                        cgcAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTg                                                                        >  2:1484984/1‑100 (MQ=255)
                        cgcAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTg                                                                        >  1:770289/1‑100 (MQ=255)
                                    cTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACgg                                                            <  1:253499/100‑1 (MQ=255)
                                     tGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACggg                                                           >  2:980440/1‑100 (MQ=255)
                                      ggCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACgggg                                                          <  1:305091/100‑1 (MQ=255)
                                       gCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGc                                                         <  2:1349448/100‑1 (MQ=255)
                                           ccAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTAc                                                     >  2:751941/1‑100 (MQ=255)
                                            cAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTAcg                                                    >  2:1196596/1‑100 (MQ=255)
                                              atatTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTAcgcg                                                  <  2:1238563/100‑1 (MQ=255)
                                                atTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGAc                                                <  2:834384/100‑1 (MQ=255)
                                                 ttGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTAcgcg                                                  <  2:78282/97‑1 (MQ=255)
                                                   gCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGa                                             >  1:1224992/1‑100 (MQ=255)
                                                    cTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGAt                                            <  1:1337128/100‑1 (MQ=255)
                                                    cTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGAt                                            >  2:187901/1‑100 (MQ=255)
                                                      gAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCGGCGCACGGGGCTTACGCGACTGATTa                                          >  2:1340469/1‑100 (MQ=255)
                                                       aaaCCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTAc                                         >  2:1109290/1‑100 (MQ=255)
                                                       aaaCCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTAc                                         <  1:1439069/100‑1 (MQ=255)
                                                        aaCCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACt                                        >  2:1280233/1‑100 (MQ=255)
                                                        aaCCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACt                                        >  2:206846/1‑100 (MQ=255)
                                                         aCCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACTGGGCTTACGCGACTGATTACTc                                       <  2:912012/100‑1 (MQ=255)
                                                         aCCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTc                                       <  1:1482737/100‑1 (MQ=255)
                                                          ccAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCa                                      <  2:1129166/100‑1 (MQ=255)
                                                                 aaCCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGtt                               >  2:553515/1‑100 (MQ=255)
                                                                   ccATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCATGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAg                             <  1:81401/100‑1 (MQ=255)
                                                                     aTTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGcc                           >  1:154320/1‑100 (MQ=255)
                                                                        gAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAgg                        <  2:182541/100‑1 (MQ=255)
                                                                         aaGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGc                       >  2:846707/1‑100 (MQ=255)
                                                                         aaGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGc                       >  2:401572/1‑100 (MQ=255)
                                                                          aGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCa                      <  2:622725/100‑1 (MQ=255)
                                                                           ggTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCa                      <  2:622822/99‑1 (MQ=255)
                                                                           ggTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGATTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCAt                     >  1:1167786/1‑100 (MQ=255)
                                                                             tATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATcc                   >  2:1247324/1‑100 (MQ=255)
                                                                             tATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATcc                   <  1:631393/100‑1 (MQ=255)
                                                                              aTTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGACTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCa                  <  2:502059/100‑1 (MQ=255)
                                                                                   tCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACg             >  1:14051/1‑100 (MQ=255)
                                                                                      ggtggtggATTATGCCCAGGAGCATGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGCAg          >  1:1281631/1‑100 (MQ=255)
                                                                                         ggtggATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGCAGCGt       <  2:1079488/100‑1 (MQ=255)
                                                                                         ggtggATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGCAGCGt       <  2:1078951/100‑1 (MQ=255)
                                                                                           tggATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGCAGCGTGc     >  2:404329/1‑100 (MQ=255)
                                                                                           tggATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGCAGCGTGc     >  2:404284/1‑100 (MQ=255)
                                                                                              aTTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGCAGCGTGCCgg  >  1:76101/1‑100 (MQ=255)
                                                                                                 |                                                                                                
ATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGCAGCGTGCCGG  >  NC_000913/162876‑163069

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: