Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 207,343 T→C 17.0% A4283A (GCT→GCC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451207,3430TC17.0% 91.0 / 2.2 39A4283A (GCT→GCCPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (12/20);  new base C (2/5);  total (14/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACGTCGGTCACCGACGAGCCGGGCACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTG  >  AE015451/207218‑207571
                                                                                                                             |                                                                                                                                                                                                                                    
aaGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACGTCGGTCACCGACGAGCCGGGCACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAg                                                                                                                        >  1:588400/1‑236 (MQ=32)
                        aaGACCCCGGTCAGCACGTCGGTCACCGACGAGCCGGGCACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCATACACCCACGCTgcgc                                                                                                 >  1:62050/1‑235 (MQ=33)
                             cccGGTCAGCACGTCGGTCACCGACGAGCCGGGCACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCa                                                                                                                         <  2:639369/206‑1 (MQ=33)
                                  tCAGCACGTCGGTCACCGACGAGCCGGGCACCCCAGGCAAGGAAGGCGACTTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAgcgcggc                                                                                                                  >  1:350619/1‑206 (MQ=17)
                                  tCAGCACGTCGGTCACCGACGAGCCGGGCACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGa                                                                                       >  1:164031/1‑235 (MQ=37)
                                   cAGCACGTCGGTCACCGACGAGCCGGGCACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAgc                                                                                                                       >  2:831811/1‑202 (MQ=32)
                                   cAGCACGTCGGTCACCGACGAGCCGGGCACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACAGCCGACCAGACATCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCggc                                                                                                                                >  1:140737/1‑192 (MQ=14)
                                         gTCGGTCACCGACGAGCCGGGCACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCACGGCCGGTGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTc                                                                                   >  1:217231/1‑232 (MQ=37)
                                         gTCGGTCACCGACGAGCCGGGCACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTg                                                                                                                               >  1:61938/1‑188 (MQ=32)
                                                   gacgaGCCGGGCACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCATa                                                                                                               <  2:394042/194‑1 (MQ=35)
                                                       aGCCGGGCACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCa                                                                  <  1:854267/235‑1 (MQ=255)
                                                          cGGGCACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGGCACCATCAAGGCCGGTGaaa                                                                                                                            <  1:601978/174‑1 (MQ=18)
                                                           gggCACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCATACAccc                                                                                                          >  1:657836/1‑191 (MQ=35)
                                                           gggCACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCATACAccc                                                                                                          <  2:657833/191‑1 (MQ=35)
                                                              ccccccAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGAt                                                                                         <  1:209889/203‑1 (MQ=38)
                                                                    gggCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGcc                                                                   <  2:486240/220‑1 (MQ=21)
                                                                       cAACGAAGGCGACTTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCg                                                                                           <  1:417607/194‑1 (MQ=21)
                                                                        aaCGAAGGCGACTTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCAc                                                <  1:240068/236‑1 (MQ=37)
                                                                            aaGGCGACTTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTcc                                            <  1:206171/236‑1 (MQ=37)
                                                                              ggCGACTTGGTCAAGGTCACCATCACAGCCGCCCAGACCTCGGTAGCCGAGACCGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTa                                                                                                               <  2:116669/167‑1 (MQ=11)
                                                                              ggCGACATGGCCACGGTCACCACCCCAGCCGCCCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTAACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCATACACCCACGCTGCGCAAg                                                                                              <  2:189929/184‑1 (MQ=11)
                                                                                cGACTTGGTCAAGGTCACCAGCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCtg                                        <  1:273125/236‑1 (MQ=37)
                                                                                  aCTTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTgcgc                                                                                                 >  2:86675/1‑177 (MQ=14)
                                                                                     tGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGa                                    >  1:858308/1‑235 (MQ=39)
                                                                                      ggtcacggtcacCATCACAGCCGACCAGACCTCGGTGGCTGAGACCGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGAt                                                                                         <  1:89816/181‑1 (MQ=255)
                                                                                             gTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGAc                             <  2:588397/234‑1 (MQ=37)
                                                                                                  cATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGAt                                                                                         <  2:164030/169‑1 (MQ=40)
                                                                                                  cATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCg                        <  2:356526/234‑1 (MQ=37)
                                                                                                    tCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAg                                                                                                                        >  2:689887/1‑136 (MQ=34)
                                                                                                     cacaGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCtg                                        >  2:644782/1‑215 (MQ=39)
                                                                                                         gCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCtt                  <  1:26406/233‑1 (MQ=35)
                                                                                                          ccGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTcc                                            <  1:9228/206‑1 (MQ=35)
                                                                                                          ccGACCAGACCTCGGGGGCCGAGAACGGCAAACCGACCTTCAACGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAgc                                                                                                                       <  2:73067/131‑1 (MQ=255)
                                                                                                                aGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCATACAccc                                                                                                          >  2:652984/1‑138 (MQ=255)
                                                                                                                aGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCATACAccc                                                                                                          <  1:652987/138‑1 (MQ=255)
                                                                                                                aGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCa                                                                                                                         <  2:597453/123‑1 (MQ=35)
                                                                                                                       ggTGGCTGAGAACGTCAAACCGGCCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTg  <  2:413648/235‑1 (MQ=25)
                                                                                                                       ggTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTg  <  1:831814/235‑1 (MQ=35)
                                                                                                                        gTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTg  <  2:565944/234‑1 (MQ=35)
                                                                                                                             |                                                                                                                                                                                                                                    
AAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACGTCGGTCACCGACGAGCCGGGCACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTCGCACCTTCGAAAACCTGCAACTG  >  AE015451/207218‑207571

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