Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,052,265 T→C 14.8% T10A (ACA→GCA)  pcrA ← ATP‑dependent DNA helicase PcrA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,052,2650TC14.8% 62.3 / 5.4 27T10A (ACA→GCA) pcrAATP‑dependent DNA helicase PcrA
Reads supporting (aligned to +/- strand):  ref base T (13/10);  new base C (2/2);  total (15/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

AAAACATTGTATGGTGAGACATCTTTTTCGTCTAATAAATAAGCAATTCTATGTGTTAAAACACGTGTCTTCCCTGAACCAGCACCTGCCATAATTAACAATGGTCCTTCTGTTGTCTTTACAGCTTCACTTTGCTCTGTATTCATATGATTTAATAACGCATTCATTTACTAGACTCCTTTATTTTTACTGTTTTTAAAGCTTTTTTAATATCTTTATAAATAATATCACCGACGA  >  CP000730/2052126‑2052362
                                                                                                                                           |                                                                                                 
aaaaCATTGTATGGTGAGACATCTTTTTCGTCTAATAAATAAGCAATTCTATGTGTTAAAACACGTGTCTTCCCTGAACCAGCACCTGCCATAATTAACAATGGTCCTTCTGTTGTCTTTACAGCTTCACTTTGCTCTGtc                                                                                                  <  2:246331/141‑2 (MQ=255)
   aCATTGTATGGTGAGACATCTTTTTCGTCTAATAAATAAGCAATTCTATGTGTTAAAACACGTGTCTTCCCTGAACCAGCACCTGCCATAATTAACAATGGTCCTTCTGTTGTCTTTACAGCTTCACTTTGCTCTGTATTc                                                                                               >  1:4509/1‑141 (MQ=255)
         tATGGTGAGACATCTTTTTTTTTTTTTTTATAAGCAATTCTATGTGTTAAAACACGTGTCTTCCCTGAACCAGCACCTGCCATAATTAACAATGGTCCTTCTGTTGTCTTTACAGCTTCACTTTGCTCTGTATTCATATGa                                                                                         >  2:146932/1‑141 (MQ=255)
         tATGGTGAGACATCTTTTTCGTCTAATAAATAAGCAATTCTATGTGTTAAAACACGTGTCTTCCCTGAACCAGCACCTGCCATAATTAACAATGGTCCTTCTGTTGTCTTTACAGCTTCACTTTGCTCTGTATTCATATGa                                                                                         >  1:29677/1‑141 (MQ=255)
          aTGGTGAGACATCTTTTTCGTCTAATAAATAAGCAATTCTATGTGTTAAAACACGTGTCTTCCCTGAACCAGCACCTGCCATAATTAACAATGGTCCTTCTGTTGTCTTTACAGCTTCACTTTGCTCTGTATTCATATGAt                                                                                        <  2:5743/141‑1 (MQ=255)
                 gaCATCTTTTTCGTCTAATAAATAAGCAATTCTATGTGTTAAAACACGTGTCTTCCCTGAACCAGCACCTGCCATAATTAACAATGGTCCTTCTGTTGTCTTTACAGCTTCACTTTGCTCTGTATTCATATGATTtaataa                                                                                 <  1:61328/141‑1 (MQ=255)
                             tttttttttataaGCAATTCTATGTGTTAAAACACGTGTCTTCCCTGAACCAGCACCTGCCATAATTAACAATGGTCCTTCTGTTGTCTTTACAGCTTCACTTTGCTCTGCATTCATATGATTTAATAACGCATTCATTTa                                                                     >  2:10112/10‑141 (MQ=255)
                             tttttttttataaGCAAATCTATGTGTTAAAACACGTGTCTTTCCTGAACCAGCACCTGCCATAATTAACAATGGTCCTTCTGTTGTCTTTACAGCTTCACTTTGCTCTGCATTCATATGATTTAATAACGCATTCATTTa                                                                     >  2:42761/10‑141 (MQ=255)
                                taataaATAAGCAATTCTATGTGTTAAAACACGTGTCTTCCCTGAACCAGCACCTGCCATAATTAACAATGGTCCTTCTGTTGTCTTTACAGCTTCACTTTGCTCTGTATTCATATGATTTAATAACGCATTCATTTACTa                                                                  >  1:171212/1‑141 (MQ=255)
                                                     tgtTAAAACACGTGTCTTCCCTGAACCAGCACCTGCCATAATTAACAATGGTCCTTCTGTTGTCTTTACAGCTTCACTTTGCTCTGTATTCATATGATTTAATAACGCATTCATTTACTAGACTCCTTTATTTTTACTGtt                                             >  2:205697/1‑141 (MQ=255)
                                                         aaaaCACGTGTCTTCCCTGAACCAGCACCTGCCATAATTAACAATGGTCCTTCTGTTGTCTTTACAGCTTCACTTTGCTCTGTATTCATATGATTTAATAACGCATTCATTTACTAGACTCCTTTATTTTTACTGTTTTTa                                         <  1:146932/141‑1 (MQ=255)
                                                         aaaaCACGTGTCTTCCCTGAACCAGCACCTGCCATAATTAACAATGGTCCTTCTGTTGTCTTTACAGCTTCACTTTGCTCTGTATTCATATGATTTAATAACGCATTCATTTACTAGACTCCTTTATTTTTACTGTTTTTa                                         <  2:29677/141‑1 (MQ=255)
                                                           aaCACGTGTCTTCCCTGAACCAGCACCTGCCATAATTAACAATGGTCCTTCTGTTGTCTTTACAGCTTCACTTTGCTCTGTATTCATATGATTTAATAACGCATTCATTTACTAGACTCCTTTATTTTTACTGTTTTTaaa                                       <  1:212464/141‑1 (MQ=255)
                                                           aaCACGTGTCTTCCCTGAACCAGCACCTGCCATAATTAACAATGGTCCTTCTGTTGTCTTTACAGCTTCACTTTGCTCTGTATTCATATGATTTAATAACGCATTCATTTACTAGACTCCTTTATTTTTACTGTTTTTaaa                                       <  2:251259/141‑1 (MQ=255)
                                                             cacGTGTCTTCCCTTAACCACCGCCTGCCATAATTAACAATGGTCCTTCTGTTGTCTTTACAGCTTCACTTTGCTCTGTATTCATATGATTTAATAACGCATTCATTTACTAGACTCCTTTATTTTTACTGTTTTTAAAGc                                     <  1:137763/141‑1 (MQ=255)
                                                                gtgtCTTCCCTGAACCAGCACCTGCCATAATTAACAATGGTCCTTCTGTTGTCTTTACAGCTTCACTTTGCTCTGTATTCATATGATTTAATAACGCATTCATTTACTAGACTCCTTTATTTTTACTGTTTTTAAAGCttt                                  >  2:110569/1‑141 (MQ=255)
                                                                 tgtCTTCCCTGAACCAGCACCTGCCATAATTAACAATGGTCCTTCTGTTGTCTTTACAGCTTCACTTTGCTCTGTATTCATATGATTTAATAACGCATTCATTTACTAGACTCCTTTATTTTTACTGTTTTTAAAGCtttt                                 <  1:99266/141‑1 (MQ=255)
                                                                      tCCCTGAACCAGCACCTGCCATAATTAACAATGGTCCTTCTGTTGTCTTTACAGCTTCACTTTGCTCTGTATTCATATGATTTAATAACGCATTCATTTACTAGACTCCTTTATTTTTACTGTTTTTAAAGCTTTTTTAat                            <  2:4509/141‑1 (MQ=255)
                                                                            aaCCAGCACCTGCCATAATTAACAATGGTCCTTCTGTTGTCTTTACAGCTTCACTTTGCTCTGTATCCATATGATTTAATAACGCATTCATTTACTAGACTCCTTTATTTTTACTGTTTTTAAAGCTTTTTTAATATCttt                      >  2:113130/1‑141 (MQ=255)
                                                                              ccAGCACCTGCCATAATTAACAATGGTCCTTCTGTTGTCTTTACAGCTTCACTTTGCTCTGTATTCATATGATTTAATAACGCATTCATTTACTAGACTCCTTTATTTTTATTTTTTTTTAAGCTTTTTTAATATCTTtat                    >  2:191267/1‑141 (MQ=255)
                                                                              ccAGCACCTGCCATAATTAACAATGGTCCTTCTGTTGTCTTTACAGCTTCACTTTGCTCTGTATTCATATGATTTAATAACGCATTCATTTACTAGACTCCTTTATTTTTACTGTTTTTAAAGCTTTTTTAATATCTTtat                    >  2:143693/1‑141 (MQ=255)
                                                                                 gCACCTGCCATAATTAACAATGGTCCTTCTGTTGTCTTTACAGCTTCACTTTGCTCTGTATTCATATGATTTAATAACGCATTCATTTACTAGACTCCTTTATTTTTACTGTTTTTAAAGCTTTTTTAATATCTTTataaa                 >  2:254253/1‑141 (MQ=255)
                                                                                 gCACCTGCCATAATTAACAATGGTCCTTCTGTTGTCTTTACAGCTTCACTTTGCTCTGTATTCATATGATTTAATAACGCATTCATTTACTAGACTCCTTTATTTTTACTGTTTTTAAAGCTTTTTTAATATCTTTataaa                 >  2:163485/1‑141 (MQ=255)
                                                                                 gCACCTGCCATAATTAACAATGGTCCTTCTGTTGTCTTTACAGCTTCACTTTGCTCTGTATTCATATGATTTAATAACGCATTCATTTACTAGACTCCTTTATTTTTACTGTTTTTAAAGCTTTTTTAATATCTTTataaa                 >  2:102170/1‑141 (MQ=255)
                                                                                 gCACCTGCCATAATTAACAATGGTCCTTCTGTTGTCTTTACAGCTTCACTTTGCTCTGCATTCATATGATTTAATAACGCATTCATTTACTAGACTCCTTTATTTTTACTGTTTTTAAAGCTTTTTTAATATCTTTataaa                 <  1:10112/141‑1 (MQ=255)
                                                                                 gCACCTGCCATAATTAACAATGGTCCTTCTGTTGTCTTTACAGCTTCACTTTGCTCTGCATTCATATGATTTAATAACGCATTCATTTACTAGACTCATTTATTTTTACTGTTTTTAAAGCTTTTTTAATATCTTTataac                 <  1:42761/141‑2 (MQ=255)
                                                                                  cACCGGCCATAATTAACAATGGTCCTTCTGTTGTCTTTACAGCTTCACTTTGCTCTGTATTCATATGATTTAATAACGCATTCATTTACTAGACTCCTTTATTTTTACTGTTTTTAAAGCTTTTTTAATATCTTTataaaa                >  2:210516/1‑140 (MQ=255)
                                                                                                aaCAATGGTCCTTCTGTTGTCTTTACAGCTTCACTTTGCTCTGTATTCATATGATTTAATAACGCATTCATTTACTAGACTCCTTTATTTTTACTGTTTTTAAAGCTTTTTTAATATCTTTATAAATAATATCACcgacga  <  1:8937/141‑1 (MQ=255)
                                                                                                                                           |                                                                                                 
AAAACATTGTATGGTGAGACATCTTTTTCGTCTAATAAATAAGCAATTCTATGTGTTAAAACACGTGTCTTCCCTGAACCAGCACCTGCCATAATTAACAATGGTCCTTCTGTTGTCTTTACAGCTTCACTTTGCTCTGTATTCATATGATTTAATAACGCATTCATTTACTAGACTCCTTTATTTTTACTGTTTTTAAAGCTTTTTTAATATCTTTATAAATAATATCACCGACGA  >  CP000730/2052126‑2052362

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