Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 197,866 G→C 20.3% A1124A (GCG→GCC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451197,8660GC20.3% 82.8 / 11.3 34A1124A (GCG→GCCPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (21/6);  new base C (5/2);  total (26/8)
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.89e-01

TGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTC  >  AE015451/197632‑198084
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                          
tGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCgg                                                                                                                                                                                                                           <  1:356433/236‑1 (MQ=255)
  tCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGt                                                                                                                                                                                                                         >  1:374371/1‑236 (MQ=255)
                        cacaGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACAGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGAtt                                                                                                                                                                                                   >  1:408478/1‑236 (MQ=255)
                                                             aCGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCgg                                                                                                                                                                                       <  2:301052/211‑1 (MQ=255)
                                                                                      aCTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTGGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAgc                                                                                                                                                                                             >  1:173483/1‑179 (MQ=255)
                                                                                                 accacTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCa                                                                                                                                 >  1:489647/1‑229 (MQ=255)
                                                                                                    acaaccgCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATc                                                                                                                                           >  2:298228/7‑216 (MQ=255)
                                                                                                        ctGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTc                                                                                                                   >  1:391242/1‑236 (MQ=255)
                                                                                                         tGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGAtt                                                                                                                                                                                                   >  1:375379/1‑155 (MQ=255)
                                                                                                                   aCCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCg                                                                                                                                                             >  1:255129/1‑183 (MQ=255)
                                                                                                                   aCCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGGAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCgg                                                                                                                                                                                       >  2:129402/1‑157 (MQ=255)
                                                                                                                     cGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCTGCCaa                                                                                                                                                                                   >  1:267810/1‑159 (MQ=255)
                                                                                                                      gCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGACGCggt                                                                                                                                                                                                                          >  1:139431/1‑118 (MQ=255)
                                                                                                                      gCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGc                                                                                                                                                    >  1:31945/1‑189 (MQ=255)
                                                                                                                           gCCAAGCGTGACCGAAGGCTGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCgg                                                                                                                                                                                                                           <  1:408480/113‑1 (MQ=255)
                                                                                                                             aaaGCGTGACCGAAGTCGGCGTGATCACCTACACCGTGACCCTGAGCATTCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTATAACCCAGGGCAGCGTCTATTTCCAGACCCCCGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCg                                                                                                                                          <  2:267809/191‑1 (MQ=255)
                                                                                                                                     aCCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTacaa                                                                                                                                                                       <  2:375378/155‑1 (MQ=35)
                                                                                                                                                                 tGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCACTGACGTCTACAACAACGGCTCGACCGTCAGCGTCCCCATCGAGAGCGCCACGGGCGGCAACTTCGAACAt                                                                                                             >  2:44266/1‑185 (MQ=255)
                                                                                                                                                                   aCCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATAc                                                        <  2:391241/236‑1 (MQ=255)
                                                                                                                                                                    cccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCgg                                                                                                                                                                                                                           <  2:491599/72‑1 (MQ=39)
                                                                                                                                                                      cTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCa                                                                                                                                                      <  2:270681/139‑1 (MQ=255)
                                                                                                                                                                      cTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCa                                                                                                                                                      >  1:270682/1‑139 (MQ=255)
                                                                                                                                                                      cTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCCGACTCCGAATCCGACCCCAGCCTCGACTg                                                                              >  1:47356/1‑211 (MQ=255)
                                                                                                                                                                      cTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGCTCGACTGTCAGCGTCACCATCGAGAACGCCACGGGCGGCAACTTCGAACATCTGACTCCGAACCCGACCCCAg                                                                                       >  2:483864/1‑202 (MQ=11)
                                                                                                                                                                                 tGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGACCGCCACGGGCGGCaagt                                                                                                                     >  2:268237/1‑159 (MQ=11)
                                                                                                                                                                                           ccGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCCGCGTCACCATCGAGAgcgc                                                                                                                                   >  2:416577/1‑137 (MQ=38)
                                                                                                                                                                                           ccGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGtcatca                                                                        >  1:222632/1‑196 (MQ=255)
                                                                                                                                                                                                gacagTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGa                                                                    >  2:599031/5‑195 (MQ=255)
                                                                                                                                                                                                aaCCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGa                                                                                                         >  1:220344/1‑158 (MQ=40)
                                                                                                                                                                                                  cagTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAgc                                  <  1:599032/225‑1 (MQ=255)
                                                                                                                                                                                                              ccAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGt                                                     <  2:222632/196‑1 (MQ=255)
                                                                                                                                                                                                               cAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTg                                       >  1:22051/1‑209 (MQ=255)
                                                                                                                                                                                                                          ccATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCCATGGCAATGTGACCGAAGACCAGCAGccttc  >  1:152029/1‑235 (MQ=255)
                                                                                                                                                                                                                          ccATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGa                                                                                 >  2:536284/1‑156 (MQ=255)
                                                                                                                                                                                                                                      aaGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGt                                                     >  2:38892/1‑172 (MQ=255)
                                                                                                                                                                                                                                      aaGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGAc                                                                                            <  2:643916/133‑1 (MQ=255)
                                                                                                                                                                                                                                      aaGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACAGGCGGCAACTTCGAACATCTGACTCCGAATCCGAc                                                                                            >  1:643917/1‑133 (MQ=255)
                                                                                                                                                                                                                                        gCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCCAGAGCGCCACGggc                                                                                                                            >  1:88074/1‑99 (MQ=18)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                          
TGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTC  >  AE015451/197632‑198084

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