Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,885,926 C→A 12.2% E523* (GAA→TAA)  yoaA ← putative ATP‑dependent helicase, DinG family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,885,9260CA12.2% 88.4 / 4.4 33E523* (GAA→TAA) yoaAputative ATP‑dependent helicase, DinG family
Reads supporting (aligned to +/- strand):  ref base C (14/15);  new base A (2/2);  total (16/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.00e-01

CGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGCCCTACCCCCTGTTTGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGATCATCCGGCGAGGTAAACGGCAATTTGTCGATAATTACCAATGACAATGTATCGCCACGCACGTCCACCCCTTCCCAGAAACTGCTGGTTGCCACAAGAAGCGCATTACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAG  >  CP009273/1885666‑1886193
                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                           
cgttgtgcCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGCCCTACCCCCNGTTCGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGTTCATCCGGCGAGGTAAACGGCAATTTGTCGATAATTACCAATGACAATGTATCGCCACGTACGTCCACCCCTTCCCAGAAACTGCTgg                                                                                                                                                                                                                                                               <  1:550664/269‑1 (MQ=255)
                  cATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGCCCTACCCCCTGTTTGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGATCATCCGGCGAGGTACACGGCAATTTGTCGATAATTACCCATGACAATGTATCGCCACGCACGTCCCCCCCTTCCCAGAAACTGCTGGTTGCCACAAGAACCGCAt                                                                                                                                                                                                                                             >  1:348993/1‑275 (MQ=255)
                    tCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGCCCTACCCCCTGTTTGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCACGCGCGCTTTTAACAGTGGATCATCCGGCGAGGTAAACGGCAATTTGTCGATAATTACCAATGACACTGTATCGCAACGCCCAACCACCCCTTCCCAAACACGGCTGGTTGCCAACAGACGCCCAtt                                                                                                                                                                                                                                            >  2:735991/1‑274 (MQ=255)
                                                          tcgtcgGCGTCGCGGATCAGTCGCCCTACCAACTGTTTGAGAGTAATGACGGCATCTGGTAGTTGCACTTAATCGAGCGGGTCGCCACCGAGCAAACGACAAGACTCCATGCGCGCTTTTAACAGTGGATCATCCGGCGCGGTAAACGGCAATTTGTCGATAATTACCAATGACAATGTATCGCCACGCACGTCCCCCCCTTCCCACACACTGCTGGTTGCCACAAGAAGCGCATTACAGGCACTGACCAAGAACTGCAACACNTACCCTTTGc                                                                                                                                                                                                      >  2:68880/1‑274 (MQ=255)
                                                               ggCGTCGCGAATCAGTCGCCCTACCCCCTGTTTGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGATCATCCGGCGAGGTAAACGGCAATTTGTCGATAATTACCAATGACAATGTATCGCCACGCACGTCCACCCCTTCCCAGAAActg                                                                                                                                                                                                                                                                   >  1:41102/1‑208 (MQ=255)
                                                                   tCGCGAATCAGTCGCCCTACCCCCTGTTTGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGATCATCCGGCGAGGTAAACGGCAATTTGTCGATAATTACCAATGACAATGTATCGCCACGCACGTCCACCCCTTCCCAGAAACTGCTgg                                                                                                                                                                                                                                                               <  2:41101/208‑1 (MQ=255)
                                                                                                    gTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGATCATCCGGCGAGGTAAACGGCAATTTGTCGATAATTACCAATGACAATGTATCGCCACGCACGTCCACCCCTTCCCAGAAACTGCTGGTTGCCACAAGAAGCGCATTACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGc                                                                                                                                                                                        <  1:416096/246‑1 (MQ=255)
                                                                                                                          gCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACCGTGGATCATCCGGCGAGGTAAACGGCAATTTGTCGATAATTACCAATGACAATGTATCGCAACGCACGTCCACCCCTTCCCAGAAACTGCTGGTTGCaac                                                                                                                                                                                                                                                        >  1:470220/1‑160 (MQ=255)
                                                                                                                                       cGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGATCATCCGGCGAGGTAAACGGCAATTTGTCGNTAATTACCNATGACAATGTATCGCCACGCACGTCCACCCCTTCCCAGAAACTGCTGGTTGCCACAAGAAGCGCATTACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGTAACTGCTCGGCCAGATCGCGCATCATGGCGTGCg                                                                                                                        <  1:487354/275‑1 (MQ=255)
                                                                                                                                       cGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGATCATCCGGCGAGGTAAACGGCAATTTGTCGATAATTACCAATGACAATGTATCGCCACGCACGTCCACCCCTTCCCAGAAACTGCTGGTTGCCACAAGAAGCGCATTACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCg                                                                                                                        <  1:736005/275‑1 (MQ=255)
                                                                                                                                                    cgcAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGATCATCCGGCGAGGTAAACGGCAATTTGTCGATAATTACCAATGACAATGTATCGCCACGCACGTCCACCCCTTACCAGAAACTGCTGGTTGCCACAAGAAGCGCATTACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTcccc                                                                                                                                                                                            >  2:640421/1‑194 (MQ=255)
                                                                                                                                                        aaCGACAATCTTCCATGCGCGCTTTTAACAGTGGATCATCCTGCGAGGTAAACGGCAATTTGTCGATAATTACCAATGACAATGTATCGCCACGCACGTCCACCCCTTACCAGAAACTGCTGGTTGCCACAAGAAGCGCATTACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGc                                                                                                                                                                                        <  1:640435/194‑1 (MQ=255)
                                                                                                                                                             cAATCTTCCATGCGCGCTTTTAACAGTGGATCATCCGGCGAGGTAAACGGCAATTTGTCGATAATTACCAATGACAATGTATCGCCACGCACGTCCACCCCTTCCCAGAAACTGCTGGTTGCCACAAGAAGCGCATTACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAgg                                                                                                                                                                             <  1:350953/200‑1 (MQ=255)
                                                                                                                                                                   tCCATGCGCGCTTTTAACAGTGGATCATCCGGCGAGGTAAACGGCAATTTGTCGATAATTACCAATGACAATGTATCGCCACGCACGTCCACCCCTTCCCAGAAACTGCTGGTTGCCACACGAAGCGCATTAcc                                                                                                                                                                                                                                         >  2:238193/1‑134 (MQ=255)
                                                                                                                                                                       tGCGCGCTTTTAACAGTGGATCATCCGGCGAGGTAAACGGCAATTTGTCGATAATTACCAATGACAATGTATCGCCACGCACGTCCACCCCTTCCCAGAAACTGCTGGTTGCCACAAGAAGCGCATTACCGgcg                                                                                                                                                                                                                                     <  1:238196/134‑1 (MQ=255)
                                                                                                                                                                                        ggATCATCCGGCGAGGTAAACGGCAATTTGTCGATAATTACCAATGACAATGTATCGCCACGCACGTCCACCCCTTCCCAGAAACTGCTGGt                                                                                                                                                                                                                                                              >  1:564589/1‑92 (MQ=255)
                                                                                                                                                                                        ggATCATCCGGCGAGGTAAACGGCAATTTGTCGATAATTACCAATGACAATGTATCGCCACGCACGTCCACCCCTTCCCAGAAACTGCTGGTTGCCACAAGAAGCGCATTACCGGCGCTGAc                                                                                                                                                                                                                                >  2:257387/1‑122 (MQ=255)
                                                                                                                                                                                            catcCGGCGAGGTAAACGGCAATTTGTCGATAATTACCAATGACAATGTATCGCCACGCACGTCCACCCCTTCCCAGAAACTGCTGGTTGcc                                                                                                                                                                                                                                                          <  2:564577/92‑1 (MQ=255)
                                                                                                                                                                                            catcCGGCGAGGTAAACGGCAATTTGTCGATAATTACCAATGACAATGTATCGCCACGCACGTCCACCCCTTCCCAGAAACTGCTGGTTGCCACAAGAAGCGCATTACCGGCGCTGACAAAt                                                                                                                                                                                                                            <  1:257390/122‑1 (MQ=255)
                                                                                                                                                                                            catcCGGCGAGGTAAACGGCAATTTGTCGATAATTACCAATGACAATGTATCGCCACGCACGTCCACCCCTTCCCAGAAACTGCTGGTTGCCACAAGAAGCGCATTACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAg                                                                                                                                                                          <  2:208859/172‑1 (MQ=255)
                                                                                                                                                                                                  gCGAGGTAAACGGCAATTTGTCGATAATTACCAATGACAATGTATCGCCACGCACGTCCACCCCTTc                                                                                                                                                                                                                                                                             <  2:267739/67‑1 (MQ=255)
                                                                                                                                                                                                            cGGCAATTTGTCGATAATTACCAATGACAATGTATCGCCACGCACGTCCACCCCTTCCCAGAAACTGCTGGTTGCCACAAGAAGCGCATTACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGCTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGGACAAAGCATAAAACAAAGACCGTGGTTAGATTCGCTGATCGGTGGCAGCATTGCCGCTAACtag                                                    >  2:292197/1‑272 (MQ=255)
                                                                                                                                                                                                                         aTAATTACCAATGACAATGTATCGCCACGCACGTCCACCCCTTCCCAGCAACTGCTGGTTGCCACAAGAAGCGCATTACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCTTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAaacaac                                                                                                   >  1:226789/1‑214 (MQ=255)
                                                                                                                                                                                                                                  aaggACAATGTATCGCCACGCACGGCCACCCCTTCCCAGAAACTGCTGGTTGCCACAAGAAGCGCATTACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCtggttggt                             <  1:292202/272‑1 (MQ=255)
                                                                                                                                                                                                                                   aTGACAATGTATCGCCACGCACGTCCACCCCTTCCCAGAAACTGCTGGTTGCCACAAGAAGCGCATTACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATAc                                                                                                                                                                                >  2:224224/1‑127 (MQ=255)
                                                                                                                                                                                                                                       cAATGTATCGCCACGCACGTCCACCCCTTCCCAGAAACTGCTGGTTGCCACAAGAAGCGCATTACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGa                                                                                                                                                                            <  1:224226/127‑1 (MQ=255)
                                                                                                                                                                                                                                                  cacgcacgTCCACCCCTTCCCAGAAACTGCTGGTTGCCACAAGAAGCGCATTACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAg                                                                                                              >  2:73512/1‑178 (MQ=255)
                                                                                                                                                                                                                                                     tcacgTCCACCCCTTCCCAGAAACTGCTGGTTGCCACAAGAAGCGCATTACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCgcagc          <  1:68881/274‑1 (MQ=255)
                                                                                                                                                                                                                                                      cacgTCCACCCCTTCCCAGAAACTGCTGGTTGCCACAAGAAGCGCATTACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATa                                                                                                          <  1:73513/178‑1 (MQ=255)
                                                                                                                                                                                                                                                        cgTCCACCCCTTCCCAGAAACTGCTGGTTGCCACAAGAAGCGCATTACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGcaca       <  1:460179/275‑1 (MQ=255)
                                                                                                                                                                                                                                                           ccACCCCTTACCAGAAACTGCTGGTTGCCACAAGAAGCGCATTACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCt                                                                                                                                                                                          >  2:465566/1‑93 (MQ=255)
                                                                                                                                                                                                                                                               cccTTACCAGAAACTGCTGGTTGCCACAAGAAGCGCATTACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGcaa                                                                                                                                                                                      <  1:465573/93‑1 (MQ=255)
                                                                                                                                                                                                                                                                   tCCCAGAAACTGCTGGTTGCCACAAGAAGCGCATTACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCTCGGAACTGCTCGGCCAGCTCGCGCATCCTGGCGTGCGAGGTACCAAGCATAACACAACGACCGTTCGTCACTTCCACGGACCGTCCCAGCATTGCCGCTACCTTGCGCGCAGACCGTGGCTGGCTGGTTTGCGGGNGAGTtcgccggnnnnacac  >  2:27411/1‑258 (MQ=255)
                                                                                                                                                                                                                                                                   tCCCAGAAACTGCTGGTTGCCACAAGAAGCGCATTACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCgc                                                                      >  2:342318/1‑201 (MQ=255)
                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                           
CGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGCCCTACCCCCTGTTTGAGAGTAATGACGGCATCTGGTAGTTGCACTTCATCGAACGGGTCGCCACCGCGCAAACGACAATCTTCCATGCGCGCTTTTAACAGTGGATCATCCGGCGAGGTAAACGGCAATTTGTCGATAATTACCAATGACAATGTATCGCCACGCACGTCCACCCCTTCCCAGAAACTGCTGGTTGCCACAAGAAGCGCATTACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAG  >  CP009273/1885666‑1886193

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