Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 197,506 T→C 14.1% N1004N (AAT→AAC)  PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451197,5060TC14.1% 101.1 / 9.4 48N1004N (AAT→AAC) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (19/22);  new base C (5/2);  total (24/24)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

CCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGT  >  AE015451/197300‑197724
                                                                                                                                                                                                              |                                                                                                                                                                                                                          
ccGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCa                                                                                                                                                                                                                    >  1:214986/1‑215 (MQ=21)
      aaTGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGcc                                                                                                                                                                                                              <  2:214982/215‑1 (MQ=17)
                                                                    cgtCGAACAGCTGCCCCCGCACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCCCCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTg                                                                                                                                                                                                    <  1:57889/161‑1 (MQ=255)
                                                                      gcGACCAGCTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTg                                                                                                                                                                                                    <  1:68432/160‑1 (MQ=25)
                                                                                        aCCCGACTCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTg                                                                                                                                                                                                    <  1:10406/143‑1 (MQ=25)
                                                                                         cccGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCg                                                                                                                         <  2:182346/217‑1 (MQ=21)
                                                                                           cGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTccca                                                                                                                                                                                         >  1:99324/1‑149 (MQ=25)
                                                                                           cGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                                                                                               <  2:422339/177‑1 (MQ=25)
                                                                                                  ggCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGAGATCCTGAGCAATCCTGccc                                                                                                                                                                                                                     <  1:21921/116‑1 (MQ=17)
                                                                                                                               aTGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGcc                                                                                                                                                                                                              <  1:436819/94‑1 (MQ=33)
                                                                                                                                  ccaccaCCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCgg                                                              >  1:225088/1‑235 (MQ=35)
                                                                                                                                         cgcgAGCCTGACGGCGAGCCCGTCGGTGACCGGAGGAGGCGTGATCGCCTACACCGTGATCCTGAGCAATCCTGcc                                                                                                                                                                                                                      <  1:153002/76‑1 (MQ=255)
                                                                                                                                         cgcgACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAAc                                                                                                                                                                                        <  1:107918/106‑1 (MQ=35)
                                                                                                                                             aCCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTATCCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGAATTCCAGACCCCgg                                                                                                                           >  2:36581/1‑163 (MQ=25)
                                                                                                                                                cTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAAt                                                                                                                      >  2:259007/1‑165 (MQ=37)
                                                                                                                                                cTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACg                                                                                                      >  2:408190/1‑181 (MQ=25)
                                                                                                                                                      gCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGc                                                                                                                                                                                                               <  2:247356/70‑1 (MQ=35)
                                                                                                                                                      gCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCCACACCGCGATCCTGAGCCATCCTGCCCAGACGCCGGTACCCGTGACCCTGTCCCACGGCAAAACCATTAc                                                                                                                                                                          >  2:107914/1‑107 (MQ=14)
                                                                                                                                                         aGCCCGCCGGTCACCGAAGGCAGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGcc                                                                                                                                                                                                              <  2:173652/68‑1 (MQ=255)
                                                                                                                                                         aGCCAGACGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCa                                                                                                                                                                                   <  1:210920/95‑1 (MQ=255)
                                                                                                                                                              gTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGGCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGa                                 <  1:448286/236‑1 (MQ=255)
                                                                                                                                                                  gTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGcc                                                                                                                                                                                                              <  2:496064/59‑1 (MQ=35)
                                                                                                                                                                          aGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGa                                                                                                                        <  2:183943/137‑1 (MQ=12)
                                                                                                                                                                               gcgTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCg                                                                                       <  1:259011/165‑1 (MQ=18)
                                                                                                                                                                                cgTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGa                                 >  2:299226/1‑218 (MQ=14)
                                                                                                                                                                                cgTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGAACCTGTc                                                                                                                                                                                            >  2:514518/1‑63 (MQ=255)
                                                                                                                                                                                     ccgccTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGAttt                                                                                                                                      <  1:517752/109‑1 (MQ=9)
                                                                                                                                                                                             acCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGt  <  1:292033/236‑1 (MQ=21)
                                                                                                                                                                                               cGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATCACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGa                                             <  2:222626/191‑1 (MQ=18)
                                                                                                                                                                                                  gacccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCACACCATTACCGTTGAAGCCg                                                                                                                                                               >  2:213463/4‑74 (MQ=255)
                                                                                                                                                                                                  gATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                                                                                               <  1:213467/74‑1 (MQ=255)
                                                                                                                                                                                                    tCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTg                                                                                                                                                                                              <  1:69774/41‑1 (MQ=32)
                                                                                                                                                                                                    tCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTTCAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTaa                                                           >  2:351441/1‑172 (MQ=17)
                                                                                                                                                                                                    tCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGATCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGCAAttt                                                        >  1:501903/1‑175 (MQ=25)
                                                                                                                                                                                                    cccTGAGCAATCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGa                                 >  1:150613/2‑198 (MQ=17)
                                                                                                                                                                                                    cccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAggg                                                                                                                                                  >  2:405759/2‑85 (MQ=255)
                                                                                                                                                                                                    cccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                                                                                               <  1:340863/71‑1 (MQ=25)
                                                                                                                                                                                                    cccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGc                                                                                                                                              >  2:62674/2‑89 (MQ=37)
                                                                                                                                                                                                    cccTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCGACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCACGAACCCGAccc                              >  1:316627/2‑199 (MQ=11)
                                                                                                                                                                                                    cccTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGt                                                             >  1:205827/2‑170 (MQ=255)
                                                                                                                                                                                                      cTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                                                                                               >  2:517748/1‑70 (MQ=18)
                                                                                                                                                                                                      cTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATCACCGtg                                                                                                                                                                      >  2:45733/1‑62 (MQ=255)
                                                                                                                                                                                                      cTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGGCCAACTGTCAAAACCTTACCGTTCAAGCCGGCaaccc                                                                                                                                                        >  2:243461/1‑74 (MQ=14)
                                                                                                                                                                                                      cTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAAtt                                                         >  1:97532/1‑172 (MQ=255)
                                                                                                                                                                                                      cTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGATCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGGAAt                                                          >  1:1510/1‑171 (MQ=11)
                                                                                                                                                                                                      cTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                                                                                               >  1:361951/1‑70 (MQ=255)
                                                                                                                                                                                                          gCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                                                                                               >  1:342226/1‑66 (MQ=39)
                                                                                                                                                                                                          gCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                                                                                               <  2:342222/66‑1 (MQ=40)
                                                                                                                                                                                                              tCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACGATGACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGcca                                                                                                                        <  1:384465/101‑3 (MQ=11)
                                                                                                                                                                                                              |                                                                                                                                                                                                                          
CCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGT  >  AE015451/197300‑197724

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