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,9730GA92.5% 201.4 / 6.2 67C290Y (TGT→TAT) hrpBputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (1/4);  new base A (34/28);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 1.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

CTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGCAGCGTGCCG  >  NC_000913/162879‑163068
                                                                                              |                                                                                               
ctccCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGccc                                                                                            >  1:857071/1‑100 (MQ=255)
  cccGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAg                                                                                          >  2:1035337/1‑100 (MQ=255)
    cGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGa                                                                                        >  1:1179786/1‑100 (MQ=255)
     ggCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAg                                                                                       >  2:1730543/1‑100 (MQ=255)
      gCACCGCAAGGGAGGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGc                                                                                      <  1:1845205/100‑1 (MQ=255)
       cACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCg                                                                                     <  2:991037/100‑1 (MQ=255)
         ccGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGACTATGCCCAGGAGCgtg                                                                                   >  2:1560538/1‑100 (MQ=255)
          cGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCgtgt                                                                                  >  1:1565044/1‑100 (MQ=255)
             aaGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGgc                                                                               <  2:1269196/100‑1 (MQ=255)
                ggATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGcgcg                                                                            <  1:630339/100‑1 (MQ=255)
                   tGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGttt                                                                         <  2:633381/100‑1 (MQ=255)
                    gcgcAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCTTGTGGCGCGtttt                                                                        <  1:1200606/100‑1 (MQ=255)
                    gcgcAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGtttt                                                                        <  2:22215/100‑1 (MQ=255)
                     cgcAAAGTGGTGCTGGTGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTg                                                                       <  1:1139822/100‑1 (MQ=255)
                      gcAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGa                                                                      <  1:672788/100‑1 (MQ=255)
                      gcAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGa                                                                      <  2:1455843/100‑1 (MQ=255)
                        acaGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATc                                                                    <  1:1184069/98‑1 (MQ=255)
                        aaaGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATc                                                                    <  2:177718/100‑1 (MQ=255)
                             ggtgCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCgcgc                                                               >  1:1861515/1‑100 (MQ=255)
                                gCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACg                                                            <  2:1530877/100‑1 (MQ=255)
                                     ngACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCt                                                       >  2:681096/2‑100 (MQ=255)
                                     ngACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCt                                                       >  2:681080/2‑100 (MQ=255)
                                      gACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCtt                                                      >  2:725763/1‑100 (MQ=255)
                                        ccAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTAc                                                    >  1:1509476/1‑100 (MQ=255)
                                          aaTATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTAcgc                                                  >  1:178650/1‑100 (MQ=255)
                                          aaTATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTAcgc                                                  >  2:61771/1‑100 (MQ=255)
                                           atatTGCTTAAACCAGTTTATCCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCTAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTAcgcg                                                 <  1:1918296/100‑1 (MQ=255)
                                            tatTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTAcgc                                                  <  2:148691/98‑1 (MQ=255)
                                             atTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGAc                                               >  2:1032073/1‑100 (MQ=255)
                                               tGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTg                                             <  1:470111/100‑1 (MQ=255)
                                               tGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTg                                             <  1:472457/100‑1 (MQ=255)
                                                 cTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGAt                                           >  1:1434872/1‑100 (MQ=255)
                                                  tGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGAtt                                          <  2:969345/100‑1 (MQ=255)
                                                    aaaCCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTAc                                        >  1:1093360/1‑100 (MQ=255)
                                                         aGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAAc                                   >  2:725764/1‑100 (MQ=255)
                                                          gTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAAcg                                  >  2:1303720/1‑100 (MQ=255)
                                                            ttAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAAcgcg                                <  2:37721/100‑1 (MQ=255)
                                                             tAACCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGTCTGATTACTCAACGCGt                               >  1:384307/1‑100 (MQ=255)
                                                              aaCCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCTCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGtt                              <  1:1615029/100‑1 (MQ=255)
                                                              aaCCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGtt                              <  1:70758/100‑1 (MQ=255)
                                                               aCCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTa                             >  1:380178/1‑100 (MQ=255)
                                                               aCCATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTa                             <  1:512975/100‑1 (MQ=255)
                                                                ccATTGAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAg                            >  1:1127483/1‑100 (MQ=255)
                                                                  aTTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGcc                          <  2:412118/100‑1 (MQ=255)
                                                                    tGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAg                        <  2:624475/100‑1 (MQ=255)
                                                                     gAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAgg                       >  1:203018/1‑100 (MQ=255)
                                                                     gAAGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAgg                       >  1:1029592/1‑100 (MQ=255)
                                                                      aaGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGc                      >  2:928035/1‑100 (MQ=255)
                                                                      aaGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGc                      <  2:1509484/100‑1 (MQ=255)
                                                                      aaGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGc                      >  1:1294193/1‑100 (MQ=255)
                                                                       aGGTATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCa                     <  1:1288861/100‑1 (MQ=255)
                                                                          tATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATcc                  >  1:29415/1‑100 (MQ=255)
                                                                          tATTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATcc                  <  2:1135867/100‑1 (MQ=255)
                                                                           aTTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCa                 >  1:226885/1‑100 (MQ=255)
                                                                           aTTCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCa                 <  2:652942/100‑1 (MQ=255)
                                                                             tCGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATg               <  1:512721/100‑1 (MQ=255)
                                                                              cGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGa              <  2:720088/100‑1 (MQ=255)
                                                                              cGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATgn              <  1:681079/100‑2 (MQ=255)
                                                                              cGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATgn              <  1:681095/100‑2 (MQ=255)
                                                                              cGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGa              >  1:841822/1‑100 (MQ=255)
                                                                              cGTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGa              >  1:1241370/1‑100 (MQ=255)
                                                                               gTCTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGAc             >  1:1860725/1‑100 (MQ=255)
                                                                                 cTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGc           >  1:1093209/1‑100 (MQ=255)
                                                                                 cTGGTGGTGGATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGc           >  1:151096/1‑100 (MQ=255)
                                                                                    gtggtggATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGCAGc        >  1:1697257/1‑100 (MQ=255)
                                                                                         ggATTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGCAGCGTGcc   >  1:928357/1‑100 (MQ=255)
                                                                                           aTTATGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTGCTCAACGCGTTAGCAAGGCATCCATGACGCAGCGTGCCg  >  1:1063215/1‑99 (MQ=255)
                                                                                              |                                                                                               
CTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGATTACTCAACGCGTTAGCCAGGCATCCATGACGCAGCGTGCCG  >  NC_000913/162879‑163068

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