Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 5,426,415 A→T 7.7% L499* (TTG→TAG)  hrpB ← ATP‑dependent RNA helicase HrpB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,426,4150AT 7.7% 142.6 / 2.6 51L499* (TTG→TAG) hrpBATP‑dependent RNA helicase HrpB
Reads supporting (aligned to +/- strand):  ref base A (22/25);  new base T (2/2);  total (24/27)
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

CCAGCCGATACTCCGCCCCACCCGCCCGGCGCTGCTGGGCGACCCGGTCAGGGTAGGCCAGCGCCAGCAGGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACGAGCACGCTGCACACCACCCTGGCCGCCACGCGCTGCCTTGCTCTCACCGCTGACCAGCGCAAGCCGGCTGTGCAAGTCGGCACCACCGCCACGCTGGATGTCC  >  ACB122/5426287‑5426549
                                                                                                                                |                                                                                                                                      
ccAGCCGATACTCCGCCCCACCCGCCCGGCGCTGCTGGGCGACCCGGTCAGGGTAGGCCAGCGCCAGCAGGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAgg                                                                                                                                   <  1:197925/134‑1 (MQ=255)
ccAGCCGATACTCCGCCCCACCCGCCCGGCGCTGCTGGGCGACCCGGTCAGGGTAGGCCAGCGCCAGCAGGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAgg                                                                                                                                   <  6:97109/134‑1 (MQ=255)
      gATACTCCGCCCCCCCCGCCCGGCGCTGCTGGGCGACCCGGTCAGGGTAGGCCAGCTCCAGCAGGGCGCCCAGACAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGCCCTAGCCACGCAACAGGCCACGGt                                                                                                                            >  8:163348/1‑135 (MQ=255)
      gATACTCCGCCCCACCCGCCCGGCGCTGCTGGGCGACCCGGTCAGGGTAGGCCAGCGCCAGCAGGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGCCCGGCGCCCCGACCTTGCCACGCAACAGGCCACGGt                                                                                                                            >  1:270859/1‑135 (MQ=255)
      gATACTCCGCCCCACCCGCCCGGCGCTGCTGGGCGACCCGGTCAGGGTAGGCCAGCGCCAGCAGGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGt                                                                                                                            >  6:59538/1‑135 (MQ=255)
      gATACTCCGCCCCACCCGCCCGGCGCTGCTGGGCGACCCGGTCAGGGTAGGCCAGCGCCAGCAGGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGt                                                                                                                            >  4:198936/1‑135 (MQ=255)
       aTACTCCGCCCCACCCGCCCGGCGCTGCTGGGCGACCCGGTCAGGGTAGGCCAGCGCCAGCAGGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCTACAGGCCACGGTa                                                                                                                           >  4:32206/1‑135 (MQ=255)
       aTACTCCGCCCCACCCGCCCGGCGCTGCTGGGCGACCCGGTCAGGGTAGGCCAGCGCCAGCAGGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCTACAGGCCACGGTa                                                                                                                           >  4:28352/1‑135 (MQ=255)
                       gcccgGCGCTGCTGGGCGACCCGGTCAGGGTAGGCCAGCGCCAGCAGGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACa                                                                                                                                     >  4:69509/1‑109 (MQ=255)
                       gcccgGCGCTGCTGGGCGACCCGGTCAGGGTAGGCCAGCGCCAGCAGGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACa                                                                                                                                     <  3:69509/109‑1 (MQ=255)
                          cgGCGCTGCGGGGCGACCCGGTCAGGGTAGGCCAGCGCCAGCAGGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACGa                                                                                                        <  5:32182/135‑1 (MQ=255)
                            gcgcTGCTGGGCGACCCGGTCAGGGTAGGCCAGCGCCAGCAGGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCTACAGGCCACGGTATTGCCGGGCCAGCTGACGAGc                                                                                                      <  3:32206/135‑1 (MQ=255)
                            gcgcTGCTGGGCGACCCGGTCAGGGTAGGCCAGCGCCAGCAGGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCTACAGGCCACGGTATTGCCGGGCCAGCTGACGAGc                                                                                                      <  3:28352/135‑1 (MQ=255)
                            gcgcTGCGGGGCGACCCGGTCAGGGTAGGCCAGCGCCAGCAGGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACGAGc                                                                                                      <  8:32183/135‑1 (MQ=255)
                              gctgctGGGCGACCCGGTCAGGGTAGGCCAGCGCCAGCAGGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACGAGCAc                                                                                                    <  5:59538/135‑1 (MQ=255)
                               ctgctgGGCGACCCGGTCAGGGTAGGCCAGCGCCAGCAGGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACGAGCACg                                                                                                   >  7:178473/1‑135 (MQ=255)
                                        gACCCGGTCAGGGTAGGCCAGCGCCAGCAGGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACGAGCACGCTGCAcacc                                                                                          <  7:71795/135‑1 (MQ=255)
                                          cccGGTCAGGGTAGGCCAGCGCCAGCAGGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACGAGCACGCTGCAcaccac                                                                                        >  5:239358/1‑135 (MQ=255)
                                                    gTAGGCCAGCGCCAGCAGGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACGAGCACGCTGCACACCACCCTGgccgcc                                                                              <  6:1704/135‑1 (MQ=255)
                                                      aGGCCAGCGCCAGCAGGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACGAGCACGCTGCACACCACCCTGGCCGCCAc                                                                            <  3:147667/135‑1 (MQ=255)
                                                      aGGCCAGCGCCAGCAGGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACGAGCACGCTGCACACCACCCTGGCCGCCAc                                                                            <  7:238737/135‑1 (MQ=255)
                                                      aGGCCAGCGCCAGCAGGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACGAGCACGCTGCACACCACCCTGGCCGCCAc                                                                            <  5:189386/135‑1 (MQ=255)
                                                           agcgccagcAGGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACGAGCAc                                                                                                    <  3:232046/106‑1 (MQ=255)
                                                           agcgccagcAGGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACGAGCAc                                                                                                    >  4:232046/1‑106 (MQ=255)
                                                           agcgccagcAGGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACGAGCACGCTGCACACCACCCTGGCCGCCACGCGCt                                                                       >  7:303896/1‑135 (MQ=255)
                                                             cgccagcAGGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACGAGCACGCTGCACACCACCCTGGCCGCCACGCGCTGc                                                                     <  4:248964/135‑1 (MQ=255)
                                                                cagcagGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACg                                                                                                         >  1:131277/1‑96 (MQ=255)
                                                                cagcagGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACg                                                                                                         <  2:131277/96‑1 (MQ=255)
                                                                cagcagGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACGAGCACGCTGCACACCACCCTGGCCGCCACGCGCTGCCtt                                                                  >  1:91819/1‑135 (MQ=255)
                                                                cagcagGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACGAGCACGCTGCACACCACCCTGGCCGCCACGCGCTGCCtt                                                                  <  4:29528/135‑1 (MQ=255)
                                                                cagcagGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACGAGCACGCTGCACACCACCCTGGCCGCCACGCGCTGCCtt                                                                  >  7:257367/1‑135 (MQ=255)
                                                                cagcagGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACGAGCACGCTGCACACCACCCTGGCCGCCACGCGCTGCCtt                                                                  >  3:36710/1‑135 (MQ=255)
                                                                cagcagGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACGAGCACGCTGCACACCACCCTGGCCGCCACGCGCTGCCtt                                                                  >  1:22657/1‑135 (MQ=255)
                                                                cagcagGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACGAGCACGCTGCACACCACCCTGGCCGCCACGCGCTGCCtt                                                                  >  3:29528/1‑135 (MQ=255)
                                                                cagcagGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACGAGCACGCTGCACACCACCCTGGCCGCCACGCGCTGCCtt                                                                  >  3:206272/1‑135 (MQ=255)
                                                                  gcagGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACGAGCACGCTGCACACCACCCTGGCCGCCACGCGCTGCCTTGc                                                                <  5:147953/135‑1 (MQ=255)
                                                                  gcagGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACGAGCACGCTGCACACCACCCTGGCCGCCACGCGCTGCCTTGc                                                                <  5:147951/135‑1 (MQ=255)
                                                                  gcagGGCGCCCAGCCAGCGTGGGGGGTCCGGGGCGGCGACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACGAGCACGCTGCACACCACCCTGGCCGCCACGCGCTGCCTTGc                                                                <  5:162980/135‑1 (MQ=255)
                                                                          cccAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCCCGCAACAGGCCACGGTATTGCCGGGCCAGCTGCCGAGCACGCTGCACACCCCCCTGGCCGCCACGCGCTGCCTTGCTCTCCCCg                                                        >  2:125543/1‑135 (MQ=255)
                                                                            cagccagcGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACGAGCACGCTGCACACCACCCTGGCCGCCACGCGCTGCCTTGCTCTCACCGCt                                                      <  4:269789/135‑1 (MQ=255)
                                                                            cagccagcGGGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACGAGCACGCTGCACACCACCCTGGCCGCCACGCGCTGCCTTGCTCTCACCGCt                                                      <  2:270859/135‑1 (MQ=255)
                                                                             agccagcGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCGTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACGAGCACGCTGCACACCACCCTGGCCGCCACGCGCTGCCTTGCTCTCACCGCTg                                                     >  1:93008/1‑135 (MQ=255)
                                                                                aagcGTGGGTGGTCCGGGTCGGCGACCGCCGCACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACGAGCACGCTGCACACCACCCTGGCCGCCACGCGCTGCCTTGCTCTCACCGCTGa                                                    >  3:296571/2‑133 (MQ=255)
                                                                                         tggtCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACGAGCACGCTGCACACCACCCTGGCCGCCACGCGCTGCCTTGCTCTCACCGCTGACCAGCGCAAGc                                         <  6:239358/135‑1 (MQ=255)
                                                                                         tggtCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACGAGCACGCTGCACACCACCCTGGCCGCCACGCGCTGCCTTGCTCTCACCGCTGACCAGCGCAAGc                                         <  7:224887/135‑1 (MQ=255)
                                                                                          ggtCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACGAGCACGCTGCACACCACCCTGGCCGCCACGCGCTGCCTTGCTCTCACCGCTGACCAGCGCAAGcc                                        >  1:251641/1‑135 (MQ=255)
                                                                                          ggtCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACGAGCACGCTGCACACCACCCTGGCCGCCACGCGCTGCCTTGCTCTCACCGCTGACCAGCGCAAGcc                                        >  7:198648/1‑135 (MQ=255)
                                                                                                   cGGCGACCGGCGCACCGACCTTGCCACGCAACAgg                                                                                                                                   <  7:262278/35‑1 (MQ=255)
                                                                                                   cGGCGACCGGCGCACCGACCTTGCCACGCAACAgg                                                                                                                                   >  8:262278/1‑35 (MQ=255)
                                                                                                       gACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACGAGCACGCTGCACACCACCCTGGCCGCCACGCGCTGCCTTGCTCTCACCGCTGACCAGCGCAAGCCGGCTGTGCAAGTc                           <  4:112761/135‑1 (MQ=255)
                                                                                                              gcACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACGAGCACGCTGCACACCACCCTGGCCGCCACGCGCTGCCTTGCTCTCACCGCTGACCAGCGCAAGCCGGCTGTGCAAGTCGGcacca                    <  7:155569/135‑1 (MQ=255)
                                                                                                                                aaCAGGCCACGGTATTGCCGGGCCAGCTGACGAGCACGCTGCACACCACCCTGGCCGCCACGCGCTGCCTTGCTCTCACCGCTGACCAGCGCAAGCCGGCTGTGCAAGTCGGCACCACCGCCACGCTGGATGTcc  <  6:218631/135‑1 (MQ=255)
                                                                                                                                aaCAGGCCACGGTATTGCCGGGCCAGCTGACGAGCACGCTGCACACCACCCTGGCCGCCACGCGCTGCCTTGCTCTCACCGCTGACCAGCGCAAGCCGGCTGTGCAAGTCGGCACCACCGCCACGCTGGATGTcc  <  8:257367/135‑1 (MQ=255)
                                                                                                                                |                                                                                                                                      
CCAGCCGATACTCCGCCCCACCCGCCCGGCGCTGCTGGGCGACCCGGTCAGGGTAGGCCAGCGCCAGCAGGGCGCCCAGCCAGCGTGGGTGGTCCGGGTCGGCGACCGGCGCACCGACCTTGCCACGCAACAGGCCACGGTATTGCCGGGCCAGCTGACGAGCACGCTGCACACCACCCTGGCCGCCACGCGCTGCCTTGCTCTCACCGCTGACCAGCGCAAGCCGGCTGTGCAAGTCGGCACCACCGCCACGCTGGATGTCC  >  ACB122/5426287‑5426549

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