Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 211,096 C→G 7.8% V5534V (GTC→GTG PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451211,0960CG 7.8% 196.5 / 6.7 62V5534V (GTC→GTGPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (31/26);  new base G (3/2);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

TGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGTTGGCCGAGAACGTCAA  >  AE015451/210864‑211329
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                         
tGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACt                                                                                                                                                                                                                                        <  1:876066/236‑1 (MQ=34)
            gcAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACaa                                                                                                                                                                                                                             >  1:1058616/1‑235 (MQ=34)
              aaCGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACg                                                                                                                                                                                                                           <  2:882287/235‑1 (MQ=33)
                    tCGGTTCAAGTGACCGACACCCCCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGt                                                                                                                                                                                                                    <  1:97725/236‑1 (MQ=18)
                          cAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCAc                                                                                                                                                                                                              <  1:1082897/236‑1 (MQ=35)
                                  cGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGAcc                                                                                                                                                                                                      <  1:1025144/236‑1 (MQ=37)
                                   gACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCAGGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGTCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAATGTCTACGTCGTCACCAACGa                                                                                                                                                                                                        <  2:562281/233‑1 (MQ=255)
                                     caccacCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCt                                                                                                                                                                                                    <  1:488288/235‑1 (MQ=37)
                                                  gtggtgGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTc                                                                                                                                                                                         <  2:1009198/233‑1 (MQ=40)
                                                     gtggCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGATCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATc                                                                                                                                                                                     <  2:561246/234‑1 (MQ=255)
                                                         ccAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGt                                                                                                                                                                                                                    <  1:17667/199‑1 (MQ=34)
                                                           aaGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGCGGGCACTGCTACCGTCACTGCCCCGGACAACGCCTACGTCGGCACCAACGAccc                                                                                                                                                                                                     >  1:333870/1‑212 (MQ=25)
                                                                gACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCg                                                                                                                                                                                                                                  >  1:768195/1‑178 (MQ=37)
                                                                  cTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGt                                                                                                                                                                      <  1:508879/236‑1 (MQ=40)
                                                                            ttCGGTAACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCACTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTGACTGCCCCGGACAACGTCTACGTCGGCGCTAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGt                                                                                                                                                              <  1:883757/233‑1 (MQ=255)
                                                                                     cGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACGACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAACACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACg                                                                                                                                                                                                         <  1:926527/182‑1 (MQ=25)
                                                                                            ggcgAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAAc                                                                                                                                             <  1:928055/235‑1 (MQ=38)
                                                                                                agaTCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTg                                                                                                                                           >  2:927439/1‑233 (MQ=37)
                                                                                                     aCCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTgg                                                                                                                                    <  1:488974/235‑1 (MQ=39)
                                                                                                                  cGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCaa                                                                                                                                              >  2:522777/1‑212 (MQ=35)
                                                                                                                  cGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACa                                                                                                                                 >  2:848173/1‑225 (MQ=39)
                                                                                                                  cGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGACACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGCCGTGATg                                                                                                                                                                                           >  1:174612/1‑167 (MQ=255)
                                                                                                                  cGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTg                                                                                                                                                                                                                                       <  2:822662/123‑1 (MQ=0)
                                                                                                                  cGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTCACTg                                                                                                                                                                                                                                       >  1:822665/1‑123 (MQ=17)
                                                                                                                   gCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGCGGTGCGGACGTCg                                                                                                                                                            >  2:939164/1‑197 (MQ=37)
                                                                                                                   gCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTg                                                                                                                                     >  2:697410/1‑220 (MQ=38)
                                                                                                                   gCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCa                                                                                                                     >  1:323005/1‑236 (MQ=37)
                                                                                                                       aCCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTg                                                                                                                                                                                              >  1:971357/1‑159 (MQ=25)
                                                                                                                          aacaaAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGt                                                                                                                                                          >  2:535370/1‑192 (MQ=35)
                                                                                                                          aacaaAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGCCTACGCCGGCAc                                                                                                                                                                                                              >  2:955409/1‑140 (MQ=21)
                                                                                                                               aGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCa                                                                                                         >  1:854274/1‑236 (MQ=39)
                                                                                                                                    gTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACcgac                                                                                                    >  1:108166/1‑236 (MQ=37)
                                                                                                                                       tGCCAATCAACAACCACAGTGCATTGAGCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGa                                                                                                                                                                         >  1:1025796/1‑164 (MQ=25)
                                                                                                                                                  aacCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTc                                                                                                                                                                                         >  2:589512/1‑137 (MQ=17)
                                                                                                                                                           gCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAAccc                                                                             >  1:737445/1‑236 (MQ=38)
                                                                                                                                                            cATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTc                                                                                                                                                             <  2:980384/155‑1 (MQ=32)
                                                                                                                                                            cATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTc                                                                                                                                                             >  1:980387/1‑155 (MQ=32)
                                                                                                                                                             aTTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTc                                                                                                                      <  1:802241/193‑1 (MQ=38)
                                                                                                                                                             aTTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTc                                                                                                                      >  2:802238/1‑193 (MQ=14)
                                                                                                                                                             aTTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGCGGTGCGGACGTCGGCAAGTTCGAGCAACTGACACTGGACAAGACCCCGGTCAGCAcga                                                                                                               >  1:921334/1‑198 (MQ=16)
                                                                                                                                                               tGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGTCACTGCTACCGTCACTACCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTg                                                                                                                                                                           >  2:322928/1‑138 (MQ=17)
                                                                                                                                                                  ccTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCgg                                                                                                           >  2:1073971/1‑198 (MQ=37)
                                                                                                                                                                  ccTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCAccc                                                                                       >  2:1078122/1‑219 (MQ=38)
                                                                                                                                                                       aCGCTGAGCGACTGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGCGGTGCGGATGTCGGCAAGTTCGAGCAACTGa                                                                                                                                          <  1:105956/163‑1 (MQ=255)
                                                                                                                                                                         gCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAACCCGGGCGGCAGCAACg                                                               >  1:471990/1‑236 (MQ=38)
                                                                                                                                                                         gCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGGCGGTAAGTTCGAGCaa                                                                                                                                              >  2:225586/1‑157 (MQ=14)
                                                                                                                                                                              gCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCAcc                                                                                                       >  2:499933/1‑191 (MQ=38)
                                                                                                                                                                                 aCGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTCCCGTCACcga                                                                                                     >  2:320091/1‑190 (MQ=37)
                                                                                                                                                                                    gCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAACCCGGGCGgcag                                                                    <  2:640349/220‑1 (MQ=39)
                                                                                                                                                                                     cAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTgg                                                   <  1:675090/236‑1 (MQ=38)
                                                                                                                                                                                                       gCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTgg                                                                                            <  2:511737/177‑1 (MQ=38)
                                                                                                                                                                                                       gCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAAc                                                                               >  2:265611/1‑190 (MQ=38)
                                                                                                                                                                                                        ccGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGc                                >  1:752491/1‑236 (MQ=38)
                                                                                                                                                                                                         cGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCa                                    <  2:640540/231‑1 (MQ=38)
                                                                                                                                                                                                         cGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGcc                               <  1:958156/236‑1 (MQ=37)
                                                                                                                                                                                                           gCCAACGGCACCGTGGGCACTGCTACCGTCACTTCCCCGGACAACGTCTCCGTCGGCACCAACGACCCTGTCGTGATGGCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGccat                                                                                                                >  2:347324/1‑150 (MQ=255)
                                                                                                                                                                                                                 ggCCCCGTGGGCACTGCTCCCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAACCCg                                                                            <  2:221108/183‑1 (MQ=255)
                                                                                                                                                                                                                  gCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCa                                       <  2:618167/219‑1 (MQ=38)
                                                                                                                                                                                                                   cACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACt                      <  2:854271/235‑1 (MQ=37)
                                                                                                                                                                                                                        ccggCACTGCCACCGTGACTGCCCCGGACAACGTCTAAGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGTGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGAc                            >  2:596275/3‑224 (MQ=12)
                                                                                                                                                                                                                                 cTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAAcc                                                                              >  1:611493/1‑165 (MQ=35)
                                                                                                                                                                                                                                      gTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGTTGGCCGAGAACGTCaa  <  1:908015/236‑1 (MQ=35)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                         
TGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGAGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGTTGGCCGAGAACGTCAA  >  AE015451/210864‑211329

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