Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 195,706 C→T 9.2% N404N (AAC→AAT PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451195,7060CT 9.2% 227.1 / 13.2 87N404N (AAC→AATPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (40/38);  new base T (5/3);  total (46/41)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

CGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCAC  >  AE015451/195487‑195936
                                                                                                                                                                                                                           |                                                                                                                                                                                                                                      
cGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGACCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCACCCCTGCCCAGACGCCgg                                                                                                                                                                                                                        >  1:1041348/1‑236 (MQ=255)
                  gaaTGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAg                                                                                                                                                                                                                               >  1:71188/2‑211 (MQ=21)
                              caacaaCGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGcc                                                                                                                                                                                                                          <  1:1428219/204‑1 (MQ=25)
                              caacaaCGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGcc                                                                                                                                                                                                                          >  2:1428216/1‑204 (MQ=37)
                                    cGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACg                                                                                                                                                                                                   >  2:672143/1‑221 (MQ=12)
                                       ttCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTg                                                                                                                                                                                 >  1:39820/1‑236 (MQ=16)
                                        tCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGa                                                                                                                                                                                <  1:996755/236‑1 (MQ=255)
                                                       aCCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTg                                                                                                                                                                                                          >  1:737002/1‑195 (MQ=255)
                                                                        aggcggTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTcc                                                                                                                                                >  1:795760/2‑236 (MQ=255)
                                                                             gTAATTTCGCGCAACTGACCCCGAATCCGACACCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGcc                                                                                                                                                                                                                          <  1:166284/157‑1 (MQ=9)
                                                                                aTTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGccg                                                                                                                                                                           <  1:991877/201‑3 (MQ=2)
                                                                                  tgcGAGCAACTGACCCCGACTCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACCATTAcc                                                                                                                                                                                     <  1:661016/187‑1 (MQ=255)
                                                                                   tCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGtga                                                                                                                                                                                                                     <  2:1298854/156‑3 (MQ=17)
                                                                                             gACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACACGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGcc                                                                                                                                                                                                                                  <  1:1168743/133‑1 (MQ=255)
                                                                                             gACCCCGAAGCCGACCCCGGCTCAGACCACGGTCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGcc                                                                                                                                                                                                                          <  1:1372992/141‑1 (MQ=255)
                                                                                                ccggAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGccg                                                                                                                                                                           <  1:124073/182‑3 (MQ=2)
                                                                                                   gaacccGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGa                                                                                                                                                                                >  1:631477/5‑177 (MQ=255)
                                                                                                   gAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGa                                                                                                                                                                                <  2:631475/177‑1 (MQ=255)
                                                                                                        cGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                           <  1:543885/177‑1 (MQ=40)
                                                                                                         gACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCGCTTACACCGTGACCCTGAGCAACCCCGCCCACACGCCgg                                                                                                                                                                                                                        >  1:1280170/1‑131 (MQ=11)
                                                                                                             ccGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGcc                                                                                                                                                                                                                          >  1:1101981/1‑125 (MQ=35)
                                                                                                              cGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAAt                                                                                                                                  <  2:599781/212‑1 (MQ=25)
                                                                                                               ggCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGtt                                                                                                                                                                                  <  2:1238793/163‑1 (MQ=12)
                                                                                                                   cAGCCCACGACCACCGACTCGGTCGACGCCACCACCGCGACCCTCACGGCGAGCCCGTCGGTCGCCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACgg                                                                                                                                                                                                  <  2:457636/143‑1 (MQ=14)
                                                                                                                              cAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAAc                                                                                                                                                                                                    <  1:1236845/130‑1 (MQ=255)
                                                                                                                              cAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAAc                                                                                                                                                                                                    >  2:1236842/1‑130 (MQ=255)
                                                                                                                                aCGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACggc                                                                                                                                                                                                                          <  2:369280/106‑3 (MQ=25)
                                                                                                                                        gTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGcc                                                                                  >  1:491548/1‑234 (MQ=255)
                                                                                                                                        gTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGcc                                                                                                                                                                                                                                  >  1:552864/1‑90 (MQ=38)
                                                                                                                                          cgatgccaccaccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                           >  1:983427/6‑143 (MQ=255)
                                                                                                                                          cgatgccaccaccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                           <  2:983425/138‑1 (MQ=255)
                                                                                                                                          cGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGCAACCGTGAACCTGTCCAACGGc                                                                                                                                                                                                 >  1:273702/1‑121 (MQ=255)
                                                                                                                                          cGACACCACTACCGCGACCCTGACGACGAGCCCGTCGGTCACCGAAGGAGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCa                                                                                 >  1:839929/1‑233 (MQ=255)
                                                                                                                                          cGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGcc                                                                                                                                                                                                >  1:519775/1‑122 (MQ=25)
                                                                                                                                                      cgcgACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACgg                                                                                                                 >  1:1253799/1‑189 (MQ=37)
                                                                                                                                                      cgcgACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTc                                                                   >  1:1409269/1‑235 (MQ=255)
                                                                                                                                                             cTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGcc                                                                                                                                                                                                                          <  2:1179412/77‑1 (MQ=35)
                                                                                                                                                             cTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGcc                                                                                                                                                                                                                          <  1:1084815/77‑1 (MQ=34)
                                                                                                                                                             cTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAACCGACTAcctc                                                                                                                                                                                   >  2:493501/1‑114 (MQ=255)
                                                                                                                                                             cTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGa                                                                                                                                >  1:55301/1‑167 (MQ=38)
                                                                                                                                                             cTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGc                                                            >  1:787387/1‑235 (MQ=255)
                                                                                                                                                                aCGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGcc                                                                                                                                                                                                                          >  1:875439/1‑74 (MQ=34)
                                                                                                                                                                   gCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGAcc                                                                                                                                           >  1:389806/1‑150 (MQ=39)
                                                                                                                                                                        cccGTCGGTCACCGAAGGCGGCGTGATCACTTACACCTGTACCCTAAACACATCTGCCCACCCACCAGTAACCGTGACACTGTCCAACCCCCAAACGATTACCGTAGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACggtt                                                                                                               >  1:705607/1‑171 (MQ=11)
                                                                                                                                                                          cctCGGTCACCGAAATCTGCGTAATCACCTACACTGTGACCCTGAGCAACCCTGCCCAGACg                                                                                                                                                                                                                            <  1:785565/60‑1 (MQ=11)
                                                                                                                                                                          cGTCGGTGACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                           <  2:797562/111‑1 (MQ=255)
                                                                                                                                                                          cGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGccc                                                                                                                                                                                                                                 <  1:555799/57‑1 (MQ=38)
                                                                                                                                                                             ccgTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGcc                                                                                                                                                                                                                          <  1:905914/59‑1 (MQ=37)
                                                                                                                                                                             cGGTCACCGAAGGGGGCGTGATCACTTACCCCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCa                                                                                                                                    <  2:313529/147‑1 (MQ=255)
                                                                                                                                                                               gtgaCCGAAGGTGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGt                                                                                                                             <  2:175697/149‑1 (MQ=16)
                                                                                                                                                                               gtgaCCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACa                                                              <  1:64729/212‑1 (MQ=9)
                                                                                                                                                                                tCACCGAAGGTGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGt                                                                                                                                                        <  2:772713/124‑1 (MQ=255)
                                                                                                                                                                                tCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAGCGATCACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCa                                                                                                      <  2:378212/174‑1 (MQ=12)
                                                                                                                                                                                     gAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTg                                                                                                                                                                                                                >  1:798660/1‑63 (MQ=255)
                                                                                                                                                                                             cgTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAAc                                                                                                                                                                                            <  2:272921/75‑1 (MQ=255)
                                                                                                                                                                                               tGATCACTTACACCGTGATCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                           <  1:1479188/90‑1 (MQ=255)
                                                                                                                                                                                               tGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGGCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGAcc                                                                                                          >  1:464821/1‑155 (MQ=255)
                                                                                                                                                                                               tGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGCCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTc                                                                                                              >  1:223860/1‑151 (MQ=255)
                                                                                                                                                                                                gATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATc                         >  1:739644/1‑235 (MQ=255)
                                                                                                                                                                                                 aTCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGac                                                                >  1:1139372/1‑194 (MQ=37)
                                                                                                                                                                                                    aCTTACACCGTGACCCTGAGCAACCCTGCCCAGACg                                                                                                                                                                                                                            >  2:518272/1‑36 (MQ=40)
                                                                                                                                                                                                            cGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCg                                                                                                                                                                            <  2:1041621/76‑1 (MQ=255)
                                                                                                                                                                                                            cGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGTCCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCg                                              <  1:220796/202‑1 (MQ=255)
                                                                                                                                                                                                               gACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGa                                                                 >  2:487983/1‑180 (MQ=9)
                                                                                                                                                                                                               gACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGtt                                                                                                                                                                                  >  2:265211/1‑67 (MQ=255)
                                                                                                                                                                                                               gACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGtt                                                                                                                                                                                  <  1:265213/67‑1 (MQ=255)
                                                                                                                                                                                                                 cccTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCgagacc                                                                                     >  2:40758/1‑156 (MQ=38)
                                                                                                                                                                                                                   cTGAGCAATCCTGCCCAGACGCCGGTAACAGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGCAAAACCCCGGGCAGCGTCGATTTCAAGACCCCGGCCAATGACGTCTACAACAACGGCTCg                                                                                                             >  2:723497/1‑132 (MQ=11)
                                                                                                                                                                                                                   cTGAGCAACGCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGAcaccac     <  1:1406306/236‑1 (MQ=255)
                                                                                                                                                                                                                   cTGAGCAACCCTGCCCAGACGCCGGTGACCGTGAACCTGTCCAACGGCCAAACCATTAc                                                                                                                                                                                      >  2:1082816/1‑59 (MQ=11)
                                                                                                                                                                                                                   cTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACgg                                                                                                                                                                                                  >  1:976053/1‑47 (MQ=255)
                                                                                                                                                                                                                   cTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAAc                                                                                                                                                                                            <  2:1039354/53‑1 (MQ=255)
                                                                                                                                                                                                                   cTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAg                                                                                                                                              >  1:1007246/1‑99 (MQ=255)
                                                                                                                                                                                                                   cTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGaatc                                                >  2:488177/1‑191 (MQ=2)
                                                                                                                                                                                                                   cTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCCGCGTCCCCAt                                                                                            >  2:465546/1‑149 (MQ=255)
                                                                                                                                                                                                                   cTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCGCGACCGTCAGCGTCACCATCgagaac                                                                                     >  2:420081/1‑154 (MQ=255)
                                                                                                                                                                                                                      aGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGAcaccaccac  <  1:487985/236‑1 (MQ=255)
                                                                                                                                                                                                                      aGCAACCCTGCCCAGACGCCGGTCACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAg                                                                                                                                              <  2:1444158/96‑1 (MQ=12)
                                                                                                                                                                                                                      aGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAGACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGt                                                                                                                                                        <  2:1422964/86‑1 (MQ=255)
                                                                                                                                                                                                                      aGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAAc                                                                                                                                                                                                    <  1:933246/42‑1 (MQ=255)
                                                                                                                                                                                                                      aGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAAc                                                                                                                                                                                                    >  2:933244/1‑42 (MQ=255)
                                                                                                                                                                                                                      aGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTgg                                                                             <  2:510574/161‑1 (MQ=40)
                                                                                                                                                                                                                       gCAACCCTGCCCAGACGCCGGTAACTGTGAACCTGTCCAACg                                                                                                                                                                                                   >  1:1194313/1‑42 (MQ=255)
                                                                                                                                                                                                                       gCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCgg                                                                                                                                                                           <  2:334991/66‑1 (MQ=255)
                                                                                                                                                                                                                       gCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCCTc                                                                                           >  2:391955/1‑146 (MQ=255)
                                                                                                                                                                                                                       gCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGAcaccaccac  >  1:1404849/1‑235 (MQ=255)
                                                                                                                                                                                                                         aaCCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAAc                                                                                                                                                                                                    >  2:433815/1‑39 (MQ=255)
                                                                                                                                                                                                                           |                                                                                                                                                                                                                                      
CGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCAC  >  AE015451/195487‑195936

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