Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 210,388 G→A 10.2% K5298K (AAG→AAA PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451210,3880GA10.2% 189.6 / 9.2 68K5298K (AAG→AAAPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (35/25);  new base A (5/2);  total (41/27)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTC  >  AE015451/210161‑210613
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                 
aTCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGtcat                                                                                                                                                                                                                           >  1:765838/1‑236 (MQ=255)
 tCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGtcat                                                                                                                                                                                                                           >  1:992518/1‑235 (MQ=255)
      tCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGt                                                                                                                                                                                                                     >  2:886847/1‑236 (MQ=19)
        ggcggTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGACTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGTCCGTCATCAcc                                                                                                                                                                                                                       <  1:1206565/232‑1 (MQ=255)
         gcggTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGcc                                                                                                                                                                                                                  <  1:129651/236‑1 (MQ=255)
          cggTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGATCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCg                                                                                                                                                                                                                 >  1:845671/1‑236 (MQ=255)
               gACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCa                                                                                                                                                                                                             <  1:1152629/235‑1 (MQ=255)
                          tGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGTGGCGAGATCACTTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGTGCCGGCAAAACCGtcag                                                                                                                                                                                                                           >  1:355995/1‑209 (MQ=255)
                                 ttCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGt                                                                                                                                                                                                   <  1:8465/227‑1 (MQ=255)
                                     agaACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGc                                                                                                                                                                                                               <  2:52685/211‑1 (MQ=255)
                                      gaACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCg                                                                                                                                                                                     >  2:996949/1‑236 (MQ=16)
                                          cTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTg                                                                                                                                                                                                                    <  2:189458/201‑1 (MQ=255)
                                           tGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACt                                                                                                                                                                                <  1:534898/236‑1 (MQ=255)
                                            ggAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCgg                                                                                                                                                                                                                >  1:389847/1‑203 (MQ=255)
                                            ggAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGc                                                                                                                                                                                                                   >  2:777971/1‑200 (MQ=255)
                                                   ggCGGTGCCGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCg                                                                                                                                                                          <  1:7955/234‑1 (MQ=255)
                                                      ggTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGCCATCACCg                                                                                                                                                                                                                      >  2:242521/1‑187 (MQ=255)
                                                       gTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGATCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTGACTGCCCCGGACAAc                                                                                                                                                                    >  1:12450/1‑236 (MQ=37)
                                                                                caccacCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGAccc                                                                                                                                             >  2:1069936/1‑234 (MQ=255)
                                                                                      cgacgaAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACg                                                                                                                                                                                                          >  1:470653/1‑167 (MQ=38)
                                                                                             gtggtgGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGc                                                                                                                                                                                          <  2:379581/176‑1 (MQ=38)
                                                                                             gtggtgGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGa                                                                                                                                                                                  >  2:1271879/1‑184 (MQ=255)
                                                                                               ggtggCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGa                                                                                                                               <  2:992513/233‑1 (MQ=255)
                                                                                                      aaGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCg                                                                                                                                                                                     >  2:526559/1‑172 (MQ=40)
                                                                                                              tGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCCCCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTa                                                                                                                                                               >  2:122288/1‑186 (MQ=255)
                                                                                                                     ccGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACg                                                                                                                                                 >  1:70346/1‑193 (MQ=40)
                                                                                                                     ccGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACg                                                                                                                                                 <  2:70345/193‑1 (MQ=40)
                                                                                                                     ccGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCaa                                                                                                                                                                                                            >  2:1251179/1‑134 (MQ=18)
                                                                                                                      cGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTc                                                                                                     >  1:773319/1‑236 (MQ=255)
                                                                                                                                cGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCATAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGtt                                                                                             <  2:674741/234‑1 (MQ=255)
                                                                                                                                                 ccTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTgg                                                                            <  2:1215776/234‑1 (MQ=255)
                                                                                                                                                        cATCACCCTGAACAACAACGGCGGCCTGCCGAGCAACAACCACAGTCCATTGAGCTTCACGCTGAGCTACGGCAATACCGTCATCACCGTGCCGGCCAACGGCAcc                                                                                                                                                                                                     >  2:899609/1‑106 (MQ=11)
                                                                                                                                                         aTCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGc                                                                                        <  2:560440/214‑1 (MQ=255)
                                                                                                                                                           cACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGt                                                                                                                     >  1:686810/1‑183 (MQ=40)
                                                                                                                                                             cccTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCAcc                                                                                                                                                                                                     <  1:734547/101‑1 (MQ=35)
                                                                                                                                                             cccTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCg                                                                                                                                                                                     <  2:146962/117‑1 (MQ=38)
                                                                                                                                                             cccTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCg                                                                                                                                                                                     >  1:146965/1‑117 (MQ=38)
                                                                                                                                                             cccTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGt                                                               >  1:664271/1‑235 (MQ=255)
                                                                                                                                                             cccTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCCTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCAcc                                                                                                                                                                                                     >  2:734542/1‑101 (MQ=25)
                                                                                                                                                               cTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAg                                                            >  1:893287/1‑236 (MQ=255)
                                                                                                                                                               cTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACg                                                                                                                                                 >  1:604545/1‑151 (MQ=25)
                                                                                                                                                                tGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGCGGGCACTGCAACCGTGACTGCCCCGGACAACGCCTACGTCGGCAGCACCGCCCCTGTCgcg                                                                                                                                      >  2:201954/1‑159 (MQ=255)
                                                                                                                                                                   ccAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTc                                                                                                                                                           >  2:621909/1‑137 (MQ=37)
                                                                                                                                                                     aacaaAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATGACCGTGCCGGCCAACgg                                                                                                                                                                                                         >  2:307483/1‑89 (MQ=18)
                                                                                                                                                                     aacaaAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTgg                                                                                                                  >  2:488633/1‑176 (MQ=255)
                                                                                                                                                                               gCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGCCCCTGTCGTGATGt                                                                                                                                  >  1:41139/1‑150 (MQ=255)
                                                                                                                                                                               gCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACcgac                                            <  1:501689/236‑1 (MQ=255)
                                                                                                                                                                                     cGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGAc                                                                                                                                                                                                                                 <  1:508834/49‑1 (MQ=38)
                                                                                                                                                                                                cacaGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAgg                           <  1:242697/236‑1 (MQ=255)
                                                                                                                                                                                                 agaGTGCATTGACGTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCAGTCGTGCTGTCGATCGCTACCGTTGGCGGTGCGGACgc                                                                                                      <  2:679629/158‑2 (MQ=255)
                                                                                                                                                                                                 acaGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCg                                                                                                           <  2:66873/155‑1 (MQ=255)
                                                                                                                                                                                                      gCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTg                                                                                                                                                                                                                    <  2:783790/45‑1 (MQ=38)
                                                                                                                                                                                                       cATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACgg                                                                                                                                                                                                         <  2:1264082/55‑1 (MQ=35)
                                                                                                                                                                                                       cATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACGACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAAc                       >  2:810619/1‑233 (MQ=255)
                                                                                                                                                                                                       cATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAg                                     >  1:109625/1‑219 (MQ=255)
                                                                                                                                                                                                       cATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGa                     >  1:491229/1‑235 (MQ=255)
                                                                                                                                                                                                       cATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGa                                                                                  >  1:556983/1‑174 (MQ=38)
                                                                                                                                                                                                         ttGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCaaa                                                                                                                                                                                                           >  2:400628/1‑50 (MQ=255)
                                                                                                                                                                                                            aCCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACCAGACGCCGGt                                                               >  2:437538/1‑188 (MQ=255)
                                                                                                                                                                                                                tCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCCACCGTTGGCGGTGCGGACGTCGGTAAGCTCGAg                                                                                         >  2:202499/1‑158 (MQ=255)
                                                                                                                                                                                                                   cGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTc        <  2:664266/236‑1 (MQ=255)
                                                                                                                                                                                                                    gCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCg                                                                                                                                                                                     >  1:445378/1‑62 (MQ=35)
                                                                                                                                                                                                                    gCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCg                                                                                                                                                                                     <  2:445373/62‑1 (MQ=35)
                                                                                                                                                                                                                    gCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGcccc                                                                                                                                                                           <  1:1061766/72‑1 (MQ=34)
                                                                                                                                                                                                                    gCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCg                                                                                                                >  2:587358/1‑131 (MQ=255)
                                                                                                                                                                                                                    gCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGc                                        >  2:942410/1‑203 (MQ=16)
                                                                                                                                                                                                                    gCTGAGCGACGGCAAGACCGGCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGc                                        <  2:815794/203‑1 (MQ=255)
                                                                                                                                                                                                                         gCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTc  >  1:130768/1‑236 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                 
ATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAGACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGCAGCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGGCGGTGCGGACGTCGGTAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTC  >  AE015451/210161‑210613

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