Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 210,465 G→C 14.0% S5324T (AGC→ACC)  PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451210,4650GC14.0% 103.7 / 5.0 34S5324T (AGC→ACC) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (20/9);  new base C (2/3);  total (22/12)
Fisher's exact test for biased strand distribution p-value = 3.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCAC  >  AE015451/210311‑210694
                                                                                                                                                          |                                                                                                                                                                                                                                     
acCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACa                                                                                                                                                          <  2:453583/232‑1 (MQ=18)
              caacaaAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCa                                                                                                                                           <  2:719700/233‑1 (MQ=255)
                  aaaGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGc                                                                                                                                            <  2:236756/228‑1 (MQ=255)
                    aGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGCGCAACTGACg                                                                                                                                                                 >  2:588540/1‑205 (MQ=255)
                    aGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAg                                                                                                                                             >  2:265807/1‑225 (MQ=255)
                    aGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTg                                                                                                                                                              >  2:667113/1‑208 (MQ=255)
                         gCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGcaa                                                                                                                                                                                                          >  2:563719/1‑157 (MQ=25)
                                               tGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAAc                                                                                                        >  1:528596/1‑235 (MQ=255)
                                                 cATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACTGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTGACTGCTCCGGACAACGTCTACGTCGGCGCTAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGAc                                                                                                                                                           >  1:167798/1‑182 (MQ=11)
                                                    tGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGc                                                                                                  >  1:528310/1‑236 (MQ=255)
                                                    tGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGACTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAg                                                                                                             >  1:526849/1‑225 (MQ=255)
                                                         ttCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGcc                                                                                                                                                   <  1:310727/182‑1 (MQ=255)
                                                           cACGCTAAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTg                                                                                            <  1:730350/235‑1 (MQ=255)
                                                            aCGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCAcc                                                                                                                                >  1:679475/1‑198 (MQ=255)
                                                            aCGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCAcc                                                                                                                                <  2:679474/198‑1 (MQ=255)
                                                              gCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCACGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTGCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAAc                                                                                                        >  1:62505/1‑220 (MQ=255)
                                                                tGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTAc                                                                                                                  >  2:80803/1‑208 (MQ=255)
                                                                                       aCCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCGTCACCGcagaccagac                                                               >  1:720247/1‑236 (MQ=255)
                                                                                              cGGCCAACGGCACCGTGGGCACTGCCACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCAcc                                                                          <  1:485630/218‑1 (MQ=255)
                                                                                              cGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTg                                                        >  1:341094/1‑236 (MQ=255)
                                                                                               ggCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGc                                                                                                                         >  2:311858/1‑170 (MQ=255)
                                                                                                               ggCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCaa                                          <  2:341093/233‑1 (MQ=255)
                                                                                                               ggCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAAcc                                       <  1:535319/236‑1 (MQ=255)
                                                                                                                gCACTGCCACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAAcc                                       <  1:509196/235‑1 (MQ=255)
                                                                                                                gCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGTCGGTGCGGACGTCGGTACGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGCCGGTCACCGACGAGcc                                                                                                                        >  2:202748/1‑154 (MQ=255)
                                                                                                                  aCTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCgc                                                                                                >  2:438298/1‑175 (MQ=255)
                                                                                                                         ccGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCa                               >  1:316902/1‑234 (MQ=255)
                                                                                                                                 gCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGt                     >  1:217215/1‑236 (MQ=255)
                                                                                                                                   cccGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGt                           >  2:596842/1‑228 (MQ=255)
                                                                                                                                   cccGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGACGGTCACCGACGAGCCAGGCACCCCAGGCaa                                                                                                         >  2:605577/1‑150 (MQ=255)
                                                                                                                                                tACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACTGACCTTCAc                              <  1:89238/212‑1 (MQ=255)
                                                                                                                                                    tCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCgg                                                          >  1:588380/1‑180 (MQ=255)
                                                                                                                                                    tCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGCCACCATCAcc                                                                          >  1:75273/1‑164 (MQ=255)
                                                                                                                                                    gcGGCACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCAc  <  1:211021/235‑1 (MQ=38)
                                                                                                                                                          |                                                                                                                                                                                                                                     
ACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCAC  >  AE015451/210311‑210694

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