Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 213,106 T→C 5.1% N6204N (AAT→AAC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451213,1060TC 5.1% 325.8 / 2.6 95N6204N (AAT→AACPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (52/38);  new base C (3/2);  total (55/40)
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.95e-01

GTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCA  >  AE015451/212871‑213332
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                  
gTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAAt                                                                                                                                                                                                                                    <  1:260566/236‑1 (MQ=40)
gTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGCCACTGCCCCGGACAAt                                                                                                                                                                                                                                    >  1:305370/1‑236 (MQ=255)
    cgcgAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGt                                                                                                                                                                                                                                  >  1:304739/1‑234 (MQ=255)
                 gTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCaa                                                                                                                                                                                                                   <  1:1068621/236‑1 (MQ=255)
                     ccGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACg                                                                                                                                                                                                                 <  1:203860/234‑1 (MQ=38)
                                  cgaAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAAt                                                                                                                                                                                                                                    >  2:1025725/1‑202 (MQ=255)
                                       tCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATc                                                                                                                                                                                             <  1:379246/236‑1 (MQ=39)
                                         gTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCACCGACCCTGTCGTGATGt                                                                                                                                                                                                  >  2:1389119/1‑229 (MQ=255)
                                                     aCCGCCCCCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAg                                                                                                                                                                                <  1:190906/235‑1 (MQ=255)
                                                     aCCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAg                                                                                                                                                                                <  2:678932/235‑1 (MQ=255)
                                                      ccGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTc                                                                                                                                                                                                 <  1:1211870/217‑1 (MQ=39)
                                                      ccGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGAt                                                                                                                                                                                              >  2:1361338/1‑220 (MQ=255)
                                                         ccACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGgcg                                                                                                                                                                             >  1:545341/1‑234 (MQ=40)
                                                           aCCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGAt                                                                                                                                                                                              >  1:1075835/1‑215 (MQ=39)
                                                               cGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGAc                                                                                                                                                                                                               >  1:220121/1‑194 (MQ=38)
                                                                  cGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACg                                                                                                                                                                                                                 >  2:1201941/1‑189 (MQ=40)
                                                                  cGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACg                                                                                                                                                                                                                 <  1:1201944/189‑1 (MQ=255)
                                                                    gTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCa                                                                                                                                                                 <  1:741001/235‑1 (MQ=255)
                                                                    gTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACCGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGCTCGCAACCGTTGAAGGCGCGGACGCTCGCac                                                                                                                                                                >  2:489884/1‑235 (MQ=255)
                                                                         cGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGACACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACg                                                                                                                                                                                                                 <  2:5657/182‑1 (MQ=37)
                                                                          gAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTAcc                                                                                                                                                                                                                     >  2:225590/1‑177 (MQ=38)
                                                                          gAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGATGTCGGCAAGTTc                                                                                                                                                            <  2:1028421/234‑1 (MQ=255)
                                                                                     gaTTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACg                                                                                                                                                >  1:694442/1‑235 (MQ=40)
                                                                                              cacCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAg                                                                                                                                       <  1:934984/235‑1 (MQ=255)
                                                                                                   tCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGTACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGAt                                                                                                                                                                                                    >  2:769430/1‑169 (MQ=18)
                                                                                                   tCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACa                                                                                                                                         >  2:599879/1‑228 (MQ=255)
                                                                                                      cGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACaa                                                                                                                                        >  2:1160243/1‑226 (MQ=255)
                                                                                                      cGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGAcc                                                                                                                                                                                                              >  2:324390/1‑156 (MQ=14)
                                                                                                      cGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCTACCGTTGAAg                                                                                                                                                                                >  1:1207508/1‑186 (MQ=14)
                                                                                                       gCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCCACCGtt                                                                                                                                                                                    >  2:1165981/1‑181 (MQ=37)
                                                                                                       gCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGtt                                                                                                                                                                     >  1:370447/1‑196 (MQ=40)
                                                                                                       gCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGtt                                                                                                                                                                     >  2:1092288/1‑196 (MQ=255)
                                                                                                       gCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCa                                                                                                                             <  1:1224129/236‑1 (MQ=255)
                                                                                                          gACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAg                                                                                                                                       <  1:106880/223‑1 (MQ=255)
                                                                                                             caacaaAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACACTCGACAAg                                                                                                                                       <  2:1475283/220‑1 (MQ=255)
                                                                                                              aacaaAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGt                                                                                                                                                              >  2:1368617/1‑196 (MQ=255)
                                                                                                              aacaaAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGAc                                                                                                                                     >  1:92415/1‑221 (MQ=255)
                                                                                                              aacaaAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTc                                                                                                                      >  1:1025651/1‑236 (MQ=255)
                                                                                                              aacaaAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTc                                                                                                                      >  1:1275065/1‑236 (MQ=255)
                                                                                                                     aCGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTg                                                                                                                                                                                                      <  1:749669/149‑1 (MQ=39)
                                                                                                                        gtcTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACcgac                                                                                                            <  1:565932/234‑1 (MQ=255)
                                                                                                                        gCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCa                                                                                                                          <  2:1038778/222‑1 (MQ=255)
                                                                                                                             ccGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCTCCGGACAATGGCTACGTCGGTACCAACGACCCTGTCGTGATGTCCATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGc                                                                                                        <  1:1436794/235‑1 (MQ=255)
                                                                                                                              cGATTAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACcgac                                                                                                            <  2:1286374/230‑1 (MQ=11)
                                                                                                                              cGATTAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACTGTCACTGCCCCGGATAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACcgac                                                                                                            >  1:1286377/1‑230 (MQ=11)
                                                                                                                                          acaGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTAGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTCGACAAGACGCCGGTCAGCACGTCGGTGa                                                                                                                 >  1:820115/1‑213 (MQ=255)
                                                                                                                                           caGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTc                                                                                                                              <  1:1106882/199‑1 (MQ=255)
                                                                                                                                                 aTTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGCGCAACTGACGCTGGACAAGACGCCGGTCAGCAc                                                                                                                         >  2:1197067/1‑198 (MQ=255)
                                                                                                                                                 aTTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAg                                                                                   >  1:419993/1‑236 (MQ=255)
                                                                                                                                                 aTTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAg                                                                                   >  1:187146/1‑236 (MQ=255)
                                                                                                                                                  ttGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGCTg                                                                                                                                                                                                   >  2:231742/1‑123 (MQ=38)
                                                                                                                                                   tGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACcgac                                                                                                            <  2:1305620/209‑1 (MQ=255)
                                                                                                                                                   tGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACcgac                                                                                                            >  1:1305623/1‑209 (MQ=255)
                                                                                                                                                    gACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAgg                                                                                  <  2:819862/234‑1 (MQ=255)
                                                                                                                                                     aCCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGAcac                                                                                                                                               >  1:901866/1‑170 (MQ=255)
                                                                                                                                                      ccTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCaa                                                                                                                                                                >  2:395361/1‑154 (MQ=255)
                                                                                                                                                      ccTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGAc                                                                              >  1:1181782/1‑236 (MQ=255)
                                                                                                                                                        ttCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGAcct                                                                            >  1:1397823/1‑234 (MQ=255)
                                                                                                                                                        ttCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACtt                                                                            >  1:1411885/1‑236 (MQ=255)
                                                                                                                                                        ttCACGCGGAGCGACGGCAAACCCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCa                                                                                                                                                                 <  2:124573/151‑1 (MQ=255)
                                                                                                                                                           acccTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAATGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTc                                                                                                                                                                                                                           >  2:614625/4‑90 (MQ=255)
                                                                                                                                                           acccTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTc                                                                                                                                                                                                                           <  1:614627/87‑1 (MQ=255)
                                                                                                                                                             gCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACg                                                                                                                                    >  1:1290278/1‑175 (MQ=255)
                                                                                                                                                             gCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTAc                                                                                                 >  2:705095/1‑210 (MQ=255)
                                                                                                                                                             gCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCa                                                                       >  1:399695/1‑236 (MQ=255)
                                                                                                                                                               tGAGCGACGGCAAAACCGTCATCACCGTGCCGGCGAATGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTc                                                                                                                              <  2:1051212/179‑1 (MQ=255)
                                                                                                                                                                  gCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGCTGTCGATCGCAAcc                                                                                                                                                                                       >  2:696286/1‑119 (MQ=255)
                                                                                                                                                                       ggCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACcgac                                                                                                            <  1:145594/189‑1 (MQ=255)
                                                                                                                                                                            aaCCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCa                                                           <  2:1028111/233‑1 (MQ=255)
                                                                                                                                                                              ccGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACACGACGCCGGTCAGCACGGCGGTGGCCGACGAGCCGGGGACCCCAGGc                                                                                          >  2:178779/1‑200 (MQ=255)
                                                                                                                                                                                gTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGCCGCCGGTCAGCACGTCGGTGAc                                                                                                                >  1:767229/1‑176 (MQ=255)
                                                                                                                                                                                      aCCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAg                                                                                                                                       <  2:454412/147‑1 (MQ=255)
                                                                                                                                                                                      aCCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAg                                                                                                                                       >  1:454414/1‑147 (MQ=255)
                                                                                                                                                                                        cGTGCCGGCCAACGGCACCGCGGGCCCTGCCCCCGTCACGGCCCCGGACAATGTCTACGTCGGTACCAACGCCCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCAc                                                                <  2:352656/216‑1 (MQ=255)
                                                                                                                                                                                              ggCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACcgac                                                                                                            >  2:1129405/1‑166 (MQ=255)
                                                                                                                                                                                              ggCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACcgac                                                                                                            <  1:1129408/166‑1 (MQ=255)
                                                                                                                                                                                               gCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCTGTAc                                                                                                                                                                                                                      >  2:289210/1‑59 (MQ=255)
                                                                                                                                                                                               gCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGGACCCCAGGCacc                                                                                       >  2:93596/1‑184 (MQ=255)
                                                                                                                                                                                                  aaCGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGa                                                                                                                                                                                               >  1:125752/1‑79 (MQ=255)
                                                                                                                                                                                                  aaCGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGa                                                                                                                                                                                               <  2:125751/79‑1 (MQ=21)
                                                                                                                                                                                                       cACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTc                                                                                                                              >  2:229046/1‑139 (MQ=255)
                                                                                                                                                                                                       cACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTc                                                                                                                              <  1:229047/139‑1 (MQ=255)
                                                                                                                                                                                                          cGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTc                           <  1:984823/235‑1 (MQ=255)
                                                                                                                                                                                                           gTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTc                           <  1:591360/234‑1 (MQ=255)
                                                                                                                                                                                                            cgggTACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGt                                                                                                                                                                                                          <  2:60189/57‑1 (MQ=255)
                                                                                                                                                                                                             gggCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACa                                                                                                                                         >  2:77285/1‑122 (MQ=255)
                                                                                                                                                                                                               gCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTg                                                                                                                                                                                                      >  2:1353659/1‑59 (MQ=255)
                                                                                                                                                                                                                   tGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACg                                                         <  1:79401/196‑1 (MQ=255)
                                                                                                                                                                                                                   ggCCACCGGCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACGTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTc                           <  2:32790/225‑1 (MQ=255)
                                                                                                                                                                                                                         cGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGcc                                    >  1:715146/1‑211 (MQ=255)
                                                                                                                                                                                                                          gTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCAcc            <  2:1181779/234‑1 (MQ=255)
                                                                                                                                                                                                                           tCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGAcca                                            >  2:212812/1‑200 (MQ=255)
                                                                                                                                                                                                                             acggCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCt                                            <  2:143879/196‑1 (MQ=255)
                                                                                                                                                                                                                              cTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACg                                                                                      <  2:865052/156‑1 (MQ=255)
                                                                                                                                                                                                                                gCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGc                                                                                 >  2:279441/1‑159 (MQ=255)
                                                                                                                                                                                                                                  cccGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCa  >  1:84717/1‑236 (MQ=255)
                                                                                                                                                                                                                                      gACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGcac      >  1:21813/1‑228 (MQ=255)
                                                                                                                                                                                                                                       aCAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCa                                                                                             >  2:289558/1‑140 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                  
GTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATTACCTACACCATCACGCTGACCAACAAAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTGAAGGCGCGGACGTTGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCA  >  AE015451/212871‑213332

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