Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 568,007 G→A 15.3% G29D (GGC→GAC)  ybcN → DLP12 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273568,0070GA15.3% 57.8 / 2.7 26G29D (GGC→GAC) ybcNDLP12 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (10/12);  new base A (2/2);  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 = 9.86e-01

ATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGCTTAAACCGTGGCGTGAGAAACGCAGTCTTTCCCAGAATGCACTCAGCCACATGTGGTACAGCGAAATCAGTGAATACCTCATCAGCAGGGGTAAAACGTTCGCCACTCCAGCTTGGGTAAAA  >  CP009273/567863‑568149
                                                                                                                                                |                                                                                                                                              
aTGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCaaa                                                                                                                                            <  2:234055/149‑1 (MQ=255)
    aGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGAAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAgggga                                                                                                                                              >  1:222721/1‑141 (MQ=255)
       aaggCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTc                                                                                                                                     <  2:137753/146‑1 (MQ=255)
       aaTGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTc                                                                                                                                     <  2:373302/149‑1 (MQ=255)
                          gCATCTAGCGTTTGACCTCCGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGCTTAAACCGTg                                                                                                                  <  2:54035/149‑1 (MQ=255)
                                 gCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGCTTAAACCGTGGCGTgag                                                                                                           >  2:50987/1‑149 (MQ=255)
                                                 ggggggACGCGCGAACCTCCCACAAGACGGCACTCAATTACCTCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATTTGGAGGAGGGCACGTAAAATCCCAGCCTGCTTAAACCGTGGCGTGAGAAACGCAGTc                                                                                                 <  1:159300/141‑1 (MQ=255)
                                                   ggggAAGCGGGAACCCCTCACAAGACGGCATCAAATTACGTCGCGGCATCTTCACCGCTATCGGTCGGCAGATTCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGCTTAAACCGTGGCGTGAGAAACGCAGTCtt                                                                                               <  1:335919/143‑1 (MQ=255)
                                                    gggAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGCTTAAACCGTGGCGTGAGAAACGCAGTCTTTCCCAGa                                                                                        >  2:38343/1‑149 (MQ=255)
                                                      gAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGCTTAAACCGTGGGGTGAGAAACGCAGTCTTTCCCAGAAt                                                                                      >  2:123069/1‑149 (MQ=255)
                                                      gAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGCTTAAACCGTGGCGTGAGAAACGCAGTCTTTcc                                                                                            >  1:126679/1‑143 (MQ=255)
                                                             gAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATACTTTCGCATGGTGCTTAAACCGTGGCGTGAGAAACGCAGTCTTTCCCAGAATg                                                                                     <  1:203299/143‑1 (MQ=255)
                                                                        aaGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGACAAATGCTTTCGCATGGTGCTTAAACCGTGGCGTGAGAAACGCCGTCTTTCCCCGAATGCACCCAGCcac                                                                          >  1:154617/1‑143 (MQ=255)
                                                                        aaGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGACAAATGCTTTCGCATGGTGCTTAAACCGTGGCGTGAGAAACGCAGTCTTTCCCAGAATGCACTCAGCcac                                                                          >  1:154618/1‑143 (MQ=255)
                                                                             ggCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGTCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGCTTACACCGTGGCGTGAGAAACGCAGTCTTTCCCAGAATGCACTCAGCCACAtgtg                                                                     <  1:253341/143‑1 (MQ=255)
                                                                                       tACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGCTTAAACCGTGGCGTGAGAAACGCAGTCTTTCCCAGAATGCACTCAGCCAAATGTGGTACAGCGAAATCaat                                                     >  2:154791/1‑147 (MQ=255)
                                                                                            cgcgGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGCTTAAACCGTGGCGTGAGAAACGCAGTCTTTCCCAGAATGCACTCAGCCACATGTGGTACAGCGAAATCAg                                                      <  1:302020/143‑1 (MQ=255)
                                                                                             gcgGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGCTTAAACCGTGGCGTGAGAAACGCAGTCTTTCCCAGAATGCACTCAGCCACATCTGGTACAGCGAAAGCAGt                                                     >  1:76662/1‑143 (MQ=255)
                                                                                                     ttCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGCTTAAACCGTGGCGTGAGAAACGCAGTCTTTCCCAGAATGCACTCAGCCACATGTGGTACAGCGAACTCAGTGAATCCCt                                             >  1:323718/1‑143 (MQ=255)
                                                                                                          cGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGCTTAAACCGTGGCGTGAGAAACGCAGTCTTTCCCAGAATGCACTCAACCACATGTGGTACAGCGAAATCAGTGAAGACCtcctca                                        >  1:90747/1‑143 (MQ=255)
                                                                                                                    cggCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGCTTAAACCGTGGCGTGAGAAACGCAGTCTTTCCCAGAATGCACTCAGCCACATGTGGTACAGCGAAATCAGTGAATACCTCATCAGCAGGGGTAAAACGtt                        <  2:310268/149‑1 (MQ=255)
                                                                                                                             cAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGCTAAAACCGTGGCGTGAGAAACGCAGTCTTTCCCAGAATGCACTCAGCCACATGTGGTACAGCGAAATCAGTGAATACCTCATCAGCAGGGGTAAAACGTTCGCCACTcc               >  2:312137/1‑149 (MQ=255)
                                                                                                                              aGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGCTTAAACCGTGGCGTGAGAAACGCAGTCTTTCCCAGAATGCACTCAGCCACATGTGGTACAGCGAAATCGGTGAGTACCGCATCAGCAGGGGCAAAACGTTCGCCCCTccg              >  2:180866/1‑148 (MQ=255)
                                                                                                                                    aTCTGGAGGAGGGCAAATGCTTTCGCATGGTGCTTAAACCGTGGCGTGAGAAACGCAGTCTTCCCCAGAATGCACTCAGCCACATGCGGTACAGCGAAATCAGTGAATACCTCATCAGCAGGGGTAAAACGTTCGCCACTCCAGCTTgg        <  2:293650/149‑1 (MQ=255)
                                                                                                                                         gaggagGGCAAATGCTTTCGCATGGTGCTTAAACCGTGGCGTGAGAAACGCAGTCTTTCCCAGAATGCACTCAGCCACATGTGGTACAGCGAAATCAGTGAATACCTCATCAGCAGGGGTAAAACGTTCGCCACTCCAGCTTg         <  1:361294/143‑1 (MQ=255)
                                                                                                                                          aggaggACAAATGCTTTCGCATGGTGCTTAAACCGTGGCGTGAGAAACGCAGTCTTTCCCAGAATGCACTCAGCCACATGTGGTACAGCGAAATCAGTGAATACCTCATCAGCAGGGGTAAAACGTTCGCCACTCCAGCTTGGGTaaaa  <  2:176660/149‑1 (MQ=255)
                                                                                                                                          aggaggACAAATGCTTTCGCATGGTGCTTAAACCGTGGCGTGAGAAACGCAGTCTTTCCCAGAATGCACTCAGCCACATGTGGTACAGCGAAATCAGTGAATACCTCATCAGCAGGGGTAAAACGTTCGCCACTCCAGCTTGGGTaaaa  <  2:176659/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                              
ATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATCGCGGCAACTTCACCGCTATCGGTCGGCAGATCCAGCCTTATCTGGAGGAGGGCAAATGCTTTCGCATGGTGCTTAAACCGTGGCGTGAGAAACGCAGTCTTTCCCAGAATGCACTCAGCCACATGTGGTACAGCGAAATCAGTGAATACCTCATCAGCAGGGGTAAAACGTTCGCCACTCCAGCTTGGGTAAAA  >  CP009273/567863‑568149

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