Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 201,340 C→T 5.3% T2282T (ACC→ACT PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451201,3400CT 5.3% 303.0 / 5.6 93T2282T (ACC→ACTPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (37/51);  new base T (2/3);  total (39/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.23e-01

CCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTT  >  AE015451/201114‑201561
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                              
ccAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGc                                                                                                                                                                                                                        <  1:345127/235‑1 (MQ=255)
 cAACAGCACCACGGACTCGGCCACTGAAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGAACCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAg                                                                                                                                                                                                                         >  2:925611/1‑233 (MQ=255)
     aGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGaa                                                                                                                                                                                                                  <  1:949879/236‑1 (MQ=12)
        accacGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACg                                                                                                                                                                                                                <  2:764300/235‑1 (MQ=255)
            cGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCaa                                                                                                                                                                                                            <  1:824107/235‑1 (MQ=255)
              tgCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAAcc                                                                                                                                                                                                         <  2:1037401/235‑1 (MQ=255)
              ggCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAAcc                                                                                                                                                                                                         <  1:854561/236‑1 (MQ=255)
               gCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTc                                                                                                                                                                                                              <  1:315583/230‑1 (MQ=255)
                   ggCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAGCCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCt                                                                                                                                                                                                    <  1:33209/236‑1 (MQ=255)
                   ggCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCt                                                                                                                                                                                                    <  1:390068/236‑1 (MQ=255)
                    gCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCtt                                                                                                                                                                                                   >  1:546061/1‑236 (MQ=255)
                     ccACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCt                                                                                                                                                                                                    <  2:1127987/234‑1 (MQ=255)
                       aCTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGc                                                                                                                                                                                                                        <  1:956588/212‑1 (MQ=39)
                       aCTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGc                                                                                                                                                                                                                        >  2:956586/1‑212 (MQ=255)
                       aCTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAGCCAACCTTCAc                                                                                                                                                                                                >  1:627861/1‑236 (MQ=255)
                           cAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTc                                                                                                                                                                                                              <  1:389360/218‑1 (MQ=40)
                                     aGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTc                                                                                                                                                                                                              <  2:1417551/208‑1 (MQ=18)
                                             tCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGc                                                                                                                                                                               <  1:1156109/231‑1 (MQ=255)
                                             ccTAGGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGc                                                                                                                                                                               <  1:1010541/230‑1 (MQ=255)
                                              cTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGt                                                                                                                                                                                                               >  2:1222020/1‑198 (MQ=38)
                                               tATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGAcg                                                                                                                                                                                                     >  2:972961/1‑206 (MQ=40)
                                                    ctgCACTAATGCGCCGGTGGTCAATTCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGc                                                                                                                                                                       <  1:105491/229‑1 (MQ=255)
                                                         cTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCAc                                                                                                                                                                                                >  2:571726/1‑202 (MQ=40)
                                                         cTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACtg                                                                                                                                                                                              >  2:1438711/1‑204 (MQ=40)
                                                            aTGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTc                                                                                                                                                                                                              <  1:72850/185‑1 (MQ=34)
                                                               cgcCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACGAAGGTCACCGACGAGCCAGGCACGCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTc                                                                                                                                                                                      <  1:1479381/206‑1 (MQ=255)
                                                               cgcCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACtgt                                                                                                                                                                                             <  1:465557/199‑1 (MQ=38)
                                                                    gtggtCAAGGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCAGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCt                                                                                                                                                  <  1:654106/236‑1 (MQ=255)
                                                                     tggtCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACAGCCGACCAGACCTCGGTGGCTGAGAACGTCAAAcc                                                                                                                                                                                                         >  2:112256/1‑181 (MQ=11)
                                                                       gtCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGACGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCgag                                                                                                                                                                                                                    >  2:8355/1‑168 (MQ=18)
                                                                            cgCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCAc                                                                                                                                                                   <  2:409978/210‑1 (MQ=255)
                                                                                   gACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTc                                                                                                                                                                                                              <  2:909274/162‑1 (MQ=33)
                                                                                   gACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCa                                                                                                                                                                                 >  2:650800/1‑191 (MQ=40)
                                                                                      gCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGcc                                                                                                                                                                             <  1:650802/191‑1 (MQ=40)
                                                                                       ccGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGa                                                                                                                                                >  2:1123608/1‑219 (MQ=255)
                                                                                         gTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCa                                                                                                                                                                                                 <  1:1488510/169‑1 (MQ=37)
                                                                                         gTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCa                                                                                                                                                                                                 >  2:1488507/1‑169 (MQ=37)
                                                                                         gTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACCGTACATATCAACACCGCCTTGGCCCCCGACCTGGCCGTGACCCTGAGCAACCACGCCCac                                                                                                                                  >  2:407420/1‑230 (MQ=11)
                                                                                            aGCGGCGCAGATGCGTGGAAGTTCGAAAACCGGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACt                                                                                                                                                                                               <  1:483678/168‑1 (MQ=37)
                                                                                                   cAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACt                                                                                                                                                                                               <  2:1418912/161‑1 (MQ=38)
                                                                                                         cGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCa                                                                                                                           <  2:894570/222‑1 (MQ=40)
                                                                                                                 ttCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAgcg                                                                                                      <  1:972081/235‑1 (MQ=255)
                                                                                                                   cGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAgcgcgc                                                                                                   >  1:1252240/1‑236 (MQ=255)
                                                                                                                   cGAAAACCTGAACATGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTAATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGcc                                                                                                                   >  1:921952/1‑220 (MQ=255)
                                                                                                                     aaaaCCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCCTCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTg                                                                                                                                                       <  2:1006651/182‑1 (MQ=255)
                                                                                                                     aaaaCCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTg                                                                                                                                                       >  1:1006653/1‑182 (MQ=255)
                                                                                                                      aaaCCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAgcgcg                                                                                                    <  2:776826/232‑1 (MQ=255)
                                                                                                                             aacctgGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATc                                                                                                                         >  2:915874/1‑204 (MQ=255)
                                                                                                                                ctgGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCacga                                                                                      >  1:1094508/1‑236 (MQ=255)
                                                                                                                                      aaGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCaacaac                                                                                                                                        >  2:305361/1‑180 (MQ=255)
                                                                                                                                      aaGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCaacaac                                                                                                                                        <  1:305363/180‑1 (MQ=255)
                                                                                                                                           cccGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTg                                                                                                                                                       <  1:1124117/160‑1 (MQ=255)
                                                                                                                                                  aGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAgg                                                                             <  1:657167/227‑1 (MQ=255)
                                                                                                                                                  aGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCAACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGt                                                                    >  1:1397994/1‑236 (MQ=255)
                                                                                                                                                     aCTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTa                                                                 <  1:944570/236‑1 (MQ=255)
                                                                                                                                                     aCTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTa                                                                 <  1:378237/236‑1 (MQ=255)
                                                                                                                                                       tACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGACACCAGCGCGcc                                                                                                  >  2:717575/1‑201 (MQ=255)
                                                                                                                                                        aCCGTCACTGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGACCGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTg                                                                                                                                                                         <  2:347652/129‑1 (MQ=255)
                                                                                                                                                           gTCACCGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATg                                                                        <  1:363624/223‑1 (MQ=255)
                                                                                                                                                            tCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGAcc                                                                                                                                                                                                                                >  2:100285/1‑71 (MQ=31)
                                                                                                                                                               cTGACGAGCCAGGCACTCCAAGCAAGGAAAGCGACATCGTCAAGGTCACCAACAGGGCAGACCAGACCTCGGTAGCCGAGAACGACAAACCGACCTTCAATATGCACGTCAATCAGACGCTGGGCCACGACATGCTCGt                                                                                                                                                        >  2:1443779/1‑139 (MQ=14)
                                                                                                                                                                cgACGAGCCAGGCACACCAGGTAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTg                                                                                                                                                       <  2:807281/138‑1 (MQ=255)
                                                                                                                                                                  aCGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGc                                                    <  1:400711/236‑1 (MQ=255)
                                                                                                                                                                   cGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCacgacg                                                                                    <  1:915876/203‑1 (MQ=255)
                                                                                                                                                                          ggCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTc                                                                                                                                                                                                              >  2:48410/1‑75 (MQ=34)
                                                                                                                                                                          ggCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAg                                            >  1:12395/1‑236 (MQ=255)
                                                                                                                                                                              cTCCTGTCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGa                                                                                                                                                                 <  2:1294034/115‑1 (MQ=255)
                                                                                                                                                                                 cAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGAt                                                                         <  2:1252237/200‑1 (MQ=255)
                                                                                                                                                                                            ggCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCAcc                                                                                                                                                                                               >  2:648080/1‑71 (MQ=38)
                                                                                                                                                                                                  aTCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGa                    <  1:194687/236‑1 (MQ=255)
                                                                                                                                                                                                   tCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTgg                     >  2:504434/1‑234 (MQ=255)
                                                                                                                                                                                                     gTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCACGTCTGCTGTGGatg                  <  1:994978/235‑3 (MQ=255)
                                                                                                                                                                                                          ggTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGcc                                                   >  2:719267/1‑197 (MQ=255)
                                                                                                                                                                                                          ggTCACCATCACGGCCGACCAGACCGCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGa                                                                          <  2:265462/174‑1 (MQ=255)
                                                                                                                                                                                                            tCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGt          <  1:694307/236‑1 (MQ=255)
                                                                                                                                                                                                            gcACCATCACCGCCGACCAGGCCTCGGTAGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGc                                                    <  1:31551/193‑1 (MQ=255)
                                                                                                                                                                                                               ccATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGa                                                                                                                                                                 >  1:611700/1‑82 (MQ=255)
                                                                                                                                                                                                               ccATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGa                                                                                                                                                                 <  2:611698/82‑1 (MQ=255)
                                                                                                                                                                                                               ccATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTcg                            <  2:546059/215‑1 (MQ=255)
                                                                                                                                                                                                                caccACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAgg                                                                             <  2:158235/162‑1 (MQ=255)
                                                                                                                                                                                                                 aTCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTcggcgg                        >  2:1485487/1‑217 (MQ=255)
                                                                                                                                                                                                                  tCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCg                                                                                                                                                         >  1:1090639/1‑87 (MQ=255)
                                                                                                                                                                                                                      ggCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTac                                                                                              <  1:353194/142‑1 (MQ=255)
                                                                                                                                                                                                                       gCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGtt  <  2:502427/233‑1 (MQ=255)
                                                                                                                                                                                                                        ccGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCa               >  1:846656/1‑219 (MQ=255)
                                                                                                                                                                                                                        ccGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGa                                                                                                                                                            >  1:1371036/1‑78 (MQ=255)
                                                                                                                                                                                                                        ccGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGa                                                                                                                                                            <  2:1371033/78‑1 (MQ=255)
                                                                                                                                                                                                                        ccGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCa                                                                                         >  2:730525/1‑145 (MQ=255)
                                                                                                                                                                                                                        ccGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGtt  >  1:488464/1‑232 (MQ=255)
                                                                                                                                                                                                                         cGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATct                                                                                                                        >  1:1261905/1‑112 (MQ=255)
                                                                                                                                                                                                                           aCCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGc                                                                                                                                                                               <  1:24534/57‑1 (MQ=255)
                                                                                                                                                                                                                           aCCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGc                                                                                                                                                                               >  2:24534/1‑57 (MQ=255)
                                                                                                                                                                                                                              aGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACAcc                                                                                           >  2:622577/1‑137 (MQ=255)
                                                                                                                                                                                                                               gACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGGGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGcc                                                   >  1:1084824/1‑176 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                              
CCAACAGCACCACGGGCTCGGCCACTGCAACCGCACCAGACAACGTCTATGTCGGCACTAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCTGCGTT  >  AE015451/201114‑201561

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