Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 199,318 T→G 11.0% T1608T (ACT→ACG PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451199,3180TG11.0% 205.5 / 7.2 80T1608T (ACT→ACGPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (30/39);  new base G (3/6);  total (35/45)
Fisher's exact test for biased strand distribution p-value = 7.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

TCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGA  >  AE015451/199083‑199521
                                                                                                                                                                                                                                           |                                                                                                                                                                                                           
tCACTAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACt                                                                                                                                                                                                             <  1:661786/236‑1 (MQ=14)
                       ttCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCCGGCAACGAAGGCGACATCGt                                                                                                                                                                                      >  1:263650/1‑236 (MQ=18)
                       ttCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGt                                                                                                                                                                                      >  1:611411/1‑236 (MQ=35)
                                   gATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAg                                                                                                                                                                                  >  2:685352/1‑228 (MQ=34)
                                   gATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAg                                                                                                                                                                                  <  1:685353/228‑1 (MQ=34)
                                           tactGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCt                                                                                                                                                                   <  1:781466/235‑1 (MQ=34)
                                                     gtGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCaa                                                                                                                                                                                                     >  2:858840/1‑191 (MQ=33)
                                                     gtGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGCCACCATCACGGCTGACCAGACTTc                                                                                                                                                        >  1:242219/1‑236 (MQ=25)
                                                      tGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACGTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACtt                                                                                                                                                         <  2:611410/234‑1 (MQ=37)
                                                        gCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCAcc                                                                                                                                                                            >  1:1014300/1‑213 (MQ=34)
                                                          aGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACagc                                                                                                                                                                                        >  1:364259/1‑197 (MQ=21)
                                                            ccAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGc                                                                                                                                                  <  1:85853/235‑1 (MQ=35)
                                                                      caccGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTc                                                                                                                                        <  1:957756/235‑1 (MQ=37)
                                                                           gTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTc                                                                                                                                        <  1:926480/230‑1 (MQ=37)
                                                                              cGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGa                                                                                                                                                                                                  >  2:424365/1‑169 (MQ=33)
                                                                              cGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCAt                                                                                                                                                                          >  2:46862/1‑193 (MQ=34)
                                                                               ggCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCt                                                                                                                              <  1:354851/236‑1 (MQ=37)
                                                                               ggCCACCGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTc                                                                                                                                        <  1:915434/226‑1 (MQ=38)
                                                                                   aCTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCa                                                                                                                                                                        >  2:736309/1‑190 (MQ=35)
                                                                                   aCTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCg                                                                                                                                  >  2:639325/1‑228 (MQ=14)
                                                                                   aCTGCAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCg                                                                                                                                  <  1:639326/228‑1 (MQ=9)
                                                                                     ggTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAgg                                                                                                                                                                                               <  2:193651/164‑1 (MQ=18)
                                                                                           cgccCCGGACAACGTCTATGGCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTc                                                                                                                                        <  2:524166/214‑1 (MQ=38)
                                                                                             ccccGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGcc                                                                                                                                   >  2:1044144/1‑217 (MQ=40)
                                                                                             ccccGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGcc                                                                                                                                   <  1:1044145/217‑1 (MQ=39)
                                                                                                  gACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCATCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAAcc                                                                                                            >  1:628438/1‑235 (MQ=255)
                                                                                                  gACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGAccc                                                                                                                                                              >  1:177382/1‑184 (MQ=37)
                                                                                                  gACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACCCCAGGCaa                                                                                                                                                                                                     >  1:1064393/1‑146 (MQ=17)
                                                                                                  gACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGAACt                                                                                                                                                                                                             >  2:767804/1‑138 (MQ=11)
                                                                                                    caccGTCTATGTCGGCACCAACCAGCCGGTAGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGCGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGACGAGAACGTCAAGCCGACCTTCACCGTGCACGTCa                                                                                                               <  1:872624/227‑1 (MQ=255)
                                                                                                          cTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGcg                                                                                                                                                                   >  1:178739/1‑171 (MQ=25)
                                                                                                           tATGTCGGCACCACCCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTc                                                                                                                <  1:385998/222‑1 (MQ=255)
                                                                                                           tATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGTAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGc                                                                                                  <  1:100960/236‑1 (MQ=255)
                                                                                                               tCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGcc                                                                                                                                                 <  2:818202/185‑1 (MQ=39)
                                                                                                               tCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCa                                                                                                                           <  1:371259/207‑1 (MQ=38)
                                                                                                               tCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCg                                                                                                       >  1:249827/1‑227 (MQ=255)
                                                                                                                cGGCACCAACCAACCGGTGGTCAATTCCATCGACGAGGTCAGCGGTACAGATGCATGGAAGTTCGAAAACCTGAACCTGCACAAAACCCCTGTTAACACCGATGTCACATACGAGCCAGGCCCTCCAGGCAACGAAGGCCACATc                                                                                                                                                                                        <  2:481296/145‑1 (MQ=11)
                                                                                                                   cACCAACCTGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGACCCTGGACAAGACCCCTGTCAGCACCGAGGTCACCTACGAGCCAGGCACTCCACGCAACGAAGGCGACTTCGTCAAGGTCAccacca                                                                                                                                                                        <  2:936344/158‑1 (MQ=11)
                                                                                                                    accaaccaGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCg                                                                                                                                                >  2:1020729/1‑181 (MQ=37)
                                                                                                                     ccaaccaGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGc                                                                                                         >  2:837157/1‑219 (MQ=18)
                                                                                                                     ccaaccaGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCga                                                                                                                                               >  1:146693/1‑181 (MQ=18)
                                                                                                                     ccaaccaGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACGCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGc                                                                                                         <  1:837158/219‑1 (MQ=255)
                                                                                                                      caaccaGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACt                                                                                                                                                                                                             <  2:898554/118‑1 (MQ=18)
                                                                                                                      caaccaGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACt                                                                                                                                                                                                             >  1:898555/1‑118 (MQ=33)
                                                                                                                          cggCCGGTGGTGAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCa                                                                                                                           <  2:568938/194‑1 (MQ=255)
                                                                                                                          caGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGACACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCAcc                                                                                                                         >  2:579874/1‑198 (MQ=39)
                                                                                                                            gCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACtt                                                                                                                                                         <  2:5829/164‑1 (MQ=18)
                                                                                                                                gtggtCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGAccc                                                                              >  1:14077/1‑235 (MQ=38)
                                                                                                                                 tcgtCAACGCCATCGGCGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCAc                                                                                                                                                                             <  2:177382/137‑1 (MQ=255)
                                                                                                                                    tCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGTACGTCAAGCCGACCTTCACCCTGCACGTCAACCAGc                                                                                                         <  2:952103/204‑1 (MQ=255)
                                                                                                                                           cATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCAACGAGGTCACCGACGAGACAGGCACACCAGGCAACGAAGGCGAc                                                                                                                                                                                           >  1:76457/1‑115 (MQ=255)
                                                                                                                                                aCGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGc                                                                                                                                    >  2:894246/1‑165 (MQ=39)
                                                                                                                                                   cgGTCAGCGGCGCAGAGGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTg                                                                                  <  2:529671/212‑1 (MQ=39)
                                                                                                                                                    ggTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACg                                                                                             <  1:249538/200‑1 (MQ=38)
                                                                                                                                                       cAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCg                                                                                                                                                       >  2:168500/1‑139 (MQ=37)
                                                                                                                                                       cAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCa                                                                                               >  1:64130/1‑195 (MQ=38)
                                                                                                                                                          cggcgCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACGCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCa                                                   <  1:933196/236‑1 (MQ=255)
                                                                                                                                                           ggcgCAGAGGCAGGGAAGTCCGAAAGCCTGACCCTGGACAAGACCCCGGTCAGCGCCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCAcc                                                                                                                         <  2:241644/165‑1 (MQ=14)
                                                                                                                                                             cgcAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGcaaca                                                                     >  2:933195/1‑215 (MQ=255)
                                                                                                                                                                  aTGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGc                                                                                                      <  2:1064392/177‑1 (MQ=39)
                                                                                                                                                                          aaGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCCACAACGCCCAGGTCACCATCCAGGCTGGTga                                         >  1:32387/1‑230 (MQ=255)
                                                                                                                                                                          aaGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACGCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTc                                                                                                                                        <  1:517607/135‑1 (MQ=21)
                                                                                                                                                                           aGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTgg                                                                                       <  1:1105387/183‑1 (MQ=255)
                                                                                                                                                                               cGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAgcgc                                <  1:120824/234‑1 (MQ=40)
                                                                                                                                                                                gaaactctgaacctgGAAAAGACCCCGGTCACCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACCTCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCa                                                                                                                           <  2:1016008/136‑1 (MQ=11)
                                                                                                                                                                                gAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGc                                                                 >  2:463921/1‑200 (MQ=37)
                                                                                                                                                                                 aaaaCCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCa                                                   >  1:785776/1‑213 (MQ=38)
                                                                                                                                                                                 aaaaCCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTgaga                                       <  1:696965/225‑1 (MQ=40)
                                                                                                                                                                                  aaaCCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGca                                                                        >  1:404290/1‑191 (MQ=37)
                                                                                                                                                                                  aaaCCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGca                                                                        <  2:404289/191‑1 (MQ=37)
                                                                                                                                                                                  aaaCCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGACAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCga                                                                                                                                               >  2:676818/1‑120 (MQ=25)
                                                                                                                                                                                       tgaacctgGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGcaaca                                                                     <  2:412757/189‑1 (MQ=35)
                                                                                                                                                                                           cctgGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTg                                                                            >  2:519790/1‑178 (MQ=35)
                                                                                                                                                                                                          ggTCAGCACCGAGGTCATCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCTCGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGcaacaat                                                                   <  2:584413/172‑2 (MQ=21)
                                                                                                                                                                                                          ggTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCaacaac                                                                   <  2:1107459/172‑1 (MQ=35)
                                                                                                                                                                                                            tCAGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCCGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGa  >  1:547247/1‑235 (MQ=255)
                                                                                                                                                                                                             cAGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGa                                                                                            <  2:842363/144‑1 (MQ=255)
                                                                                                                                                                                                              aGCGCCGAGGTCACCGACGCGCCAGGCACGCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATAACGGCTgg                                                                                                                                                                 <  2:239956/74‑2 (MQ=11)
                                                                                                                                                                                                                cACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTga                                         <  2:628437/192‑1 (MQ=255)
                                                                                                                                                                                                                acccGAGGTCACCGACGGGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCAc                                                                                                                          <  2:297040/109‑1 (MQ=255)
                                                                                                                                                                                                                    ggggTCACCGACGAGCCAGGTACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGAt    <  2:547246/223‑1 (MQ=255)
                                                                                                                                                                                                                          aCCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTc                                                                                                                                        <  2:288877/87‑1 (MQ=38)
                                                                                                                                                                                                                             gacgaGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAATGCTGGTGAGACCAGCGCGCCGTACACCCacga                 >  2:812639/1‑203 (MQ=255)
                                                                                                                                                                                                                              acgaGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCCCCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAg                                   <  2:420175/184‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                           
TCAATAACCACTCGGAGCTGTACTTCAAGCTGACCGATGGCACTACTGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGTAGCCGCCCCGGACAACGTCTATGTCGGCACCAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCTGACCAGACTTCGGTAGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGTGAGACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGA  >  AE015451/199083‑199521

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: