Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 162,927 G→A A275T (GCT→ACT)  hrpB → putative ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913162,9270GA96.2% 75.7 / ‑6.2 26A275T (GCT→ACT) hrpBputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (12/13);  total (12/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGGCTTCGCGCATCGGCAGTGATGTATTGCTCTGCCCGCTGTATGGCGCGTTGTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGAT  >  NC_000913/162797‑163027
                                                                                                                                  |                                                                                                    
acTGGCTTCGCGCATCGGCAGTGATGTATTGCTCTGCCCGCTGTATGGCGCGTTGTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAAcc                                                                                              >  2:678395/1‑139 (MQ=255)
          cgcATCGGCAGTGATGTATTGCTCTGCCCGCTGTATGGCGCGTTGTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCa                                                                                    <  1:230717/139‑1 (MQ=255)
                         gTATTGCTCTGCCCGCTGTATGGCGCGTTGTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTc                                                                     >  2:123995/1‑139 (MQ=255)
                         gTATTGCTCTGCCCGCTGTATGGCGCGTTGTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTc                                                                     <  1:214105/139‑1 (MQ=255)
                            ttGCTCTGCCCGCTGTATGGCGCGTTGTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGaaa                                                                                                <  2:420555/109‑1 (MQ=255)
                            ttGCTCTGCCCGCTGTATGGCGCGTTGTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGaaa                                                                                                >  1:420555/1‑109 (MQ=255)
                                    cccGCTGTATGGCGCGTTGTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTc                                                                     >  2:1464618/1‑128 (MQ=255)
                                    cccGCTGTATGGCGCGTTGTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTc                                                                     <  1:1464618/128‑1 (MQ=255)
                                     ccGCTGTATGGCGCGTTGTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGAtt                                                         <  1:678395/139‑1 (MQ=255)
                                     ccGCTGTATGGCGCGTTGTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGAtt                                                         <  1:1297813/139‑1 (MQ=255)
                                              ggCGCGTTGTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAgg                                                <  1:123995/139‑1 (MQ=255)
                                                cgcgTTGTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAg                                              >  2:1088865/1‑139 (MQ=255)
                                                   gTTGTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGAtt                                                         >  2:517030/1‑125 (MQ=255)
                                                   gTTGTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGAtt                                                         <  1:517030/125‑1 (MQ=255)
                                                     tGTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGAtt                                                         >  1:1115070/1‑123 (MQ=255)
                                                     tGTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGAtt                                                         <  2:1115070/123‑1 (MQ=255)
                                                      gTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCTACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGtgtg                                        <  1:378794/139‑1 (MQ=255)
                                                                  tCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGAtt                                                         <  2:1099268/110‑1 (MQ=255)
                                                                  tCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGAtt                                                         >  1:1099268/1‑110 (MQ=255)
                                                                                            cGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCATGGGGCTTACGCGACTGAt  >  2:146193/1‑139 (MQ=255)
                                                                                              cAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCcg                         <  2:1003733/114‑1 (MQ=255)
                                                                                              cAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCcg                         >  1:1003733/1‑114 (MQ=255)
                                                                                                      gcgcAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCa                                                  >  2:514877/1‑81 (MQ=255)
                                                                                                      gcgcAAAGTGGTGCTGGCGACCAATATTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCa                                                  <  1:514877/81‑1 (MQ=255)
                                                                                                                             atatTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGAc      <  2:1047082/102‑1 (MQ=255)
                                                                                                                             atatTACTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGAc      >  1:1047082/1‑102 (MQ=255)
                                                                                                                                  |                                                                                                    
ACTGGCTTCGCGCATCGGCAGTGATGTATTGCTCTGCCCGCTGTATGGCGCGTTGTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGCGTGTGGCGCGTTTTGATCCGCGCACGGGGCTTACGCGACTGAT  >  NC_000913/162797‑163027

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