Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 199,915 T→C 18.0% F1807F (TTT→TTC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451199,9150TC18.0% 108.2 / 8.6 43F1807F (TTT→TTCPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (21/14);  new base C (2/6);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 1.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

GATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAAC  >  AE015451/199696‑200142
                                                                                                                                                                                                                           |                                                                                                                                                                                                                                   
gaTCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGttt                                                                                                                                                                                                                                     >  1:446791/1‑220 (MQ=38)
                 cccTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGc                                                                                                                                                                                                    >  1:42305/1‑236 (MQ=40)
                  ccTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTgg                                                                                                                                                                                                                     >  2:421974/1‑218 (MQ=255)
                       ccAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGcccc                                                                                                                                                                                                 >  1:421116/1‑233 (MQ=255)
                                   gTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACcgacga                                                                                                                                                                                   >  1:438619/1‑235 (MQ=255)
                                                     aCTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAg                                                                                                                                                                                                                <  2:250564/188‑1 (MQ=255)
                                                         ggAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAg                                                                                                                                                                                              >  1:203820/1‑202 (MQ=25)
                                                                 aCTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACAt                                                                                                                                                     >  1:251488/1‑235 (MQ=255)
                                                                   ttCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGAc                                                                                                                                                                                                                   >  2:411826/1‑171 (MQ=255)
                                                                        gCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGc                                                                                                                                                                            <  2:273173/205‑1 (MQ=14)
                                                                        gCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGc                                                                                                                                                                            >  1:273177/1‑205 (MQ=21)
                                                                                aCGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGAc                                                                                                                                                       >  1:402031/1‑218 (MQ=255)
                                                                                      ccacTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAg                                                                                                                                              <  2:471762/221‑1 (MQ=255)
                                                                                         cTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGcc                                                                                                                                                                                <  2:282316/184‑1 (MQ=2)
                                                                                         cTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGcc                                                                                                                                                                                >  1:282320/1‑184 (MQ=255)
                                                                                          tGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGcc                                                                                                                                                                                <  2:134719/183‑1 (MQ=255)
                                                                                                   ggCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACg                                                                                                                                  <  1:380340/220‑1 (MQ=255)
                                                                                                           tcAGCCCCCCGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGc                                                                                                                                                                                                    <  2:89064/145‑1 (MQ=255)
                                                                                                             aGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGc                                                                                                                                                                                 >  2:240301/1‑163 (MQ=255)
                                                                                                             aGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGc                                                                                                                                                                                 <  1:240305/163‑1 (MQ=255)
                                                                                                                accacGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAaacctga                                                                                                                                                                                                                           <  2:150409/118‑1 (MQ=9)
                                                                                                                accacGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAaacctga                                                                                                                                                                                                                           >  1:150413/1‑118 (MQ=255)
                                                                                                                    cGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCaaa                                                                                                 >  1:489785/1‑236 (MQ=255)
                                                                                                                       gCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGaaa                                                                                                                                                                                                                                 >  2:139321/1‑105 (MQ=255)
                                                                                                                                        gCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTg                                                                                                                >  1:468043/1‑201 (MQ=255)
                                                                                                                                                 aaCGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCa                                                                       >  2:436693/1‑233 (MQ=255)
                                                                                                                                                    gTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAg                                                                                                                                              <  2:243813/159‑1 (MQ=14)
                                                                                                                                                     tCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTg                                                                <  1:436697/236‑1 (MQ=255)
                                                                                                                                                      cTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTgg                                                                                                                                                                                                                     >  1:256618/1‑86 (MQ=38)
                                                                                                                                                       tATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTgg                                                               <  1:216949/235‑1 (MQ=255)
                                                                                                                                                       tATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGGTTGAAAACCTGAACATGGAAAAGACCCCg                                                                                                                                                                                                          >  2:214025/1‑96 (MQ=255)
                                                                                                                                                              gCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGa                                                        <  1:421978/235‑1 (MQ=255)
                                                                                                                                                                aCCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGa                                                        <  2:446787/233‑1 (MQ=255)
                                                                                                                                                                aCCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCaa                                                                          >  2:433213/1‑215 (MQ=255)
                                                                                                                                                                 ccAACCAGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACcgac                                                                                                                                                                                     <  1:255652/107‑1 (MQ=12)
                                                                                                                                                                    aCGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCa                                                                                                                          <  2:42302/163‑1 (MQ=40)
                                                                                                                                                                         ccGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTc                                                                                                    <  1:358100/180‑1 (MQ=25)
                                                                                                                                                                                  aaTGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTc                                                                                                    >  2:189478/1‑171 (MQ=38)
                                                                                                                                                                                  aaTGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTc                                                                                                    <  1:189482/171‑1 (MQ=38)
                                                                                                                                                                                       cATCGACGCGGTCAGGGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCTGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAAt                               <  1:483092/235‑1 (MQ=255)
                                                                                                                                                                                                 gTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGCACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGACGGCCACATCGTCa                                                                                                                                                >  2:101254/1‑112 (MQ=255)
                                                                                                                                                                                                     gcggcgCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGcc                                                                    <  2:417637/184‑1 (MQ=14)
                                                                                                                                                                                                              aTGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACACCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGAc                                                       >  1:289786/1‑188 (MQ=17)
                                                                                                                                                                                                                      aaGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAAc  >  2:324811/1‑233 (MQ=255)
                                                                                                                                                                                                                           |                                                                                                                                                                                                                                   
GATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACTGTCGTCGTGGCAGCCAACAGCACCACGGGCTCGGCGACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGATCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTTGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTGGCGAAAC  >  AE015451/199696‑200142

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