Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 162,916 C→T 55.6% A271V (GCG→GTG)  hrpB → predicted ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb162,9160CT55.6% 3.6 / 33.1 27A271V (GCG→GTG) hrpBpredicted ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base C (4/8);  new base T (7/8);  total (11/16)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGC  >  W3110S.gb/162849‑162984
                                                                   |                                                                    
ttGTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGTGAc                                                                   >  1:450161/1‑71 (MQ=255)
  gTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGAcc                                                                  >  1:1019947/1‑70 (MQ=255)
   tCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGTGACCaa                                                                <  1:1478038/71‑1 (MQ=255)
                  cGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGTGACCAATATTGCTGAAACCAg                                                 <  1:989875/71‑1 (MQ=255)
                    aaaaGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGtt                                               >  1:423264/1‑71 (MQ=255)
                        gCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAAc                                           <  1:265053/71‑1 (MQ=255)
                          gATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGTGACCAATATTGCTGAAACCAGTTTAACCa                                         <  1:1303586/71‑1 (MQ=255)
                             cctccCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGTGACCAata                                                              >  1:2626552/1‑47 (MQ=255)
                             cctccCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGTGACCAATATTGCTGAAACCAGTTTAACCAtt                                       >  1:158629/1‑70 (MQ=255)
                             cctccCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAata                                                              >  1:2535018/1‑47 (MQ=255)
                             cctccCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCAtt                                       >  1:915397/1‑70 (MQ=255)
                              ctccCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTg                                      <  1:906089/70‑1 (MQ=255)
                              ctccCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTg                                      <  1:862058/70‑1 (MQ=255)
                                      aCCGCAAGGGATGCGCAAAGTGGTGCTGGTGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTc                             <  1:1086278/71‑1 (MQ=255)
                                      aCCGCAAGGGATGCGCAAAGTGGTGCTGGTGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTc                             >  1:793710/1‑71 (MQ=255)
                                      aCCGCAAGGGATGCGCAAAGTGGTGCTGGTGACCAATATTGATGAAACCAGTTTAACCATTGAAGGTATTc                             <  1:1733434/71‑1 (MQ=255)
                                       ccGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCg                            <  1:811824/71‑1 (MQ=255)
                                          cAAGGGATGCGCAAAGTGGTGCTGGTGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTa                                <  1:33709/64‑1 (MQ=255)
                                            aGGGATGCGCAAAGTGGTGCTGGTGACCAATATTGCTGAAACCAGTTTAACCAt                                        >  1:874144/1‑54 (MQ=255)
                                            aGGGATGCGCAAAGTGGTGCTGGTGACCAATATTGCTGAAACCAGTTTAACCAt                                        >  1:307775/1‑54 (MQ=255)
                                            aGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCtgg                       <  1:1604781/71‑1 (MQ=255)
                                                 tGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCtggtggtg                  <  1:1607961/71‑1 (MQ=255)
                                                   cgcAAAGTGGTGCTGGTGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCtggtggtgg                 <  1:2123849/70‑1 (MQ=255)
                                                   cgcAAAGTGGTGCTGGTGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCtgg                       >  1:1038439/1‑64 (MQ=255)
                                                    gcAAAGTGGTGCTGGTGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGAt               <  1:463501/71‑1 (MQ=255)
                                                        aGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATtgtg           <  1:2718003/71‑1 (MQ=255)
                                                                  gCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGc  <  1:1806914/70‑1 (MQ=255)
                                                                   |                                                                    
TTGTCGCTGAACGATCAGCGAAAAGCGATCCTCCCGGCACCGCAAGGGATGCGCAAAGTGGTGCTGGCGACCAATATTGCTGAAACCAGTTTAACCATTGAAGGTATTCGTCTGGTGGTGGATTGTGCCCAGGAGC  >  W3110S.gb/162849‑162984

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