Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,428,087 C→T E214K (GAG→AAG)  insH5 ← IS5 family transposase and trans‑activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,428,0870CT98.1% 281.8 / NA 108E214K (GAG→AAG) insH5IS5 family transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (63/42);  minor base G (0/2);  total (64/44)
Fisher's exact test for biased strand distribution p-value = 1.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 110 total aligned reads displayed.

GCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTG  >  NC_000913/1427841‑1428332
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
gcCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTctct                                                                                                                                                                                                                                                   >  1:267918/1‑251 (MQ=255)
      tGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCCGCCAGTCCACCTCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATg                                                                                                                                                                                                                                              >  2:290670/1‑250 (MQ=255)
      tGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATg                                                                                                                                                                                                                                              >  1:756191/1‑250 (MQ=255)
       gTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGa                                                                                                                                                                                                                                             <  1:111379/250‑1 (MQ=255)
        tATCTGGCTTCCAAGAGGCTGAACGGCCGCTTGATGGTGCGAATCGTGGGCTCCACCCTGGCACGGATTCTGGCTGTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCTGCGATCAGCCAGTCCACATCCTCCTCGGCCAGCTACTCGCGCTGTTGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGaa                                                                                                                                                                                                                                            <  2:662292/250‑1 (MQ=14)
         aTCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTTGTAGCCGGCATCGGCTGGGACAAATTGCTTCTCTCCTTGaggc                                                                                                                                                                                                                                          >  2:18507/1‑248 (MQ=255)
         aTCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAg                                                                                                                                                                                                                                           <  1:539494/250‑1 (MQ=255)
            tGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAgc                                                                                                                                                                                                                                       >  1:9922/1‑250 (MQ=255)
              gCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGTCCGGATGCTGGCTTTCATGGATTCGATGGTGATGGCCGTTTTGTTCTTGCGCGTATTCTGCTTCAAGGTTTTAACCTTGCCGGGCCGCTCGGCGATCATCCAGTCCACATCCACCTCGGCCAGCTCCTCGCTCTGTGTCTCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTTCTTCTCTCCATGAAGCAGAtt                                                                                                                                                                                                                                     <  2:562222/251‑1 (MQ=14)
                    aCGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTAccc                                                                                                                                                                                                                                 >  1:408983/1‑249 (MQ=255)
                      gAAGCCGACCTGCCGCTTGGTGATGCGCAACGGGTGCTCCACCCTGACACGGATGCTGGCTTTCCTGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGGTTCTGCTTCCAGGTTTTTACCTTGCCGGGACGCTCGACGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGATCCTTGGTAGCCGGCATCGGCTGAGCCAAAGTGCTTCTCTCCATGAAGCAgcttcaccatct                                                                                                                                                                                                                             >  2:69947/1‑240 (MQ=17)
                      gAAGCCGAACTGCCGCTTGATGATGCGAATCGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCt                                                                                                                                                                                                                             <  1:548391/251‑1 (MQ=255)
                      gAAGCCGAACTGACGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCCTCGGCTGGGACAAATTGCTTCTCTTCATTAAGCAGATTACCCAGCt                                                                                                                                                                                                                             >  1:184556/1‑251 (MQ=255)
                          ccGAACTGCCGCTTGATGAGGCGAAACGGGTGCGCCGCCCTGGCACGGATGCTGGGTTTCATGTATTCGATGTTGTTGGCCGTTTTGTGCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCAGGAAGCAGATTACGCAGCTGAtt                                                                                                                                                                                                                         <  2:358833/251‑1 (MQ=17)
                              aCTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAgg                                                                                                                                                                                                                     >  1:183497/1‑251 (MQ=255)
                                  ccgcttgcttaagataaaCGGGTGCTCCTCCCTGGCCTTGATTCTGCCTTTCATGTATTCTATGTTGATGTCCGTGTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGTTCAGCCAGTCCACAGCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCGTCTCTCCATGAAGCAGATTACCCAGCTGAGTGAGGTCag                                                                                                                                                                                                                 <  2:781175/236‑2 (MQ=14)
                                      ttGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCTCTCCTTGGTAGCCGGCATCGGCTGAGCCAAATTGCTTCTCTCCCTGAAGCAGATTACCCAGCTGTTTGAGGTCATGCTc                                                                                                                                                                                                             >  2:669863/1‑251 (MQ=255)
                                        gatgctGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCTGCGCTCAGCCAGTCCCCATCCACCTCGGCCAGCTCCTCACGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGGCAAATAGCTTCTCTCCATGACGCCGATTACCCAGCTGTTTgcggcaatgctctt                                                                                                                                                                                                           >  2:210503/1‑238 (MQ=14)
                                                gAAACGGGTGCTCCACCCTGGCCCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGATTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGcc                                                                                                                                                                                                      >  2:838084/1‑248 (MQ=255)
                                                gAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGcc                                                                                                                                                                                                      <  1:838084/248‑1 (MQ=255)
                                                  aaCGGGTGCTCCACCCTGGCACGGATTCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTTTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCTGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCggt                                                                                                                                                                                                 <  1:745479/251‑1 (MQ=255)
                                                    cGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGAAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGACGCAGATTACCCAGCTGATGGAGGTCATGCTCGTTGGCCGCggggg                                                                                                                                                                                               >  1:164043/1‑248 (MQ=255)
                                                     gggTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCCTGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCggtggt                                                                                                                                                                                              >  2:725544/1‑251 (MQ=255)
                                                     gggTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCggtggt                                                                                                                                                                                              <  1:602605/251‑1 (MQ=255)
                                                                      aCGGTTGCTGGCTTTAAGGTAGTCGATTTTGATGGCCGTTTGGTTCTTGCGCGGATTCTGCCTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACAGCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCAGCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCa                                                                                                                                                                             <  2:494549/251‑1 (MQ=14)
                                                                      aCGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTc                                                                                                                                                                              >  1:722800/1‑250 (MQ=255)
                                                                       cGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAg                                                                                                                                                                            <  1:669863/251‑1 (MQ=255)
                                                                          aTGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGcc                                                                                                                                                                         <  1:574253/251‑1 (MQ=255)
                                                                                  tttCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTgg                                                                                                                                                                 >  1:240609/1‑251 (MQ=255)
                                                                                   ttCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTCCCTTGCCCGGACGCTCGGCGATCAGCCAGTCTACATTCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGCCCCATTGCTTCTCTCCATGAAGCGGATTACCAAGCTGATTGAAGACCTGCTTGTTGGCCGCAGTGGTGCACAGGCTGTGGGTCAGGCCACTCTTGGc                                                                                                                                                                >  2:37137/1‑251 (MQ=14)
                                                                                       tGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCg                                                                                                                                                            >  1:634785/1‑251 (MQ=255)
                                                                                          aTTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGac                                                                                                                                                          >  1:668397/1‑250 (MQ=255)
                                                                                            tCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCGCCAGGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGacac                                                                                                                                                        <  1:568986/250‑1 (MQ=255)
                                                                                             cGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGcttctc                                                                                                                                                                                                                                                    <  1:672435/157‑1 (MQ=255)
                                                                                                 gTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAAtgt                                                                                                                                                  >  1:501704/1‑251 (MQ=255)
                                                                                                 gTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAAtgg                                                                                                                                                  >  1:171256/1‑250 (MQ=255)
                                                                                                  ttGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTCCCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCGCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGa                                                                                                                                                                                                                                       <  2:160688/165‑1 (MQ=255)
                                                                                                  ttGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGa                                                                                                                                                                                                                                       >  1:160688/1‑165 (MQ=255)
                                                                                                   tGATGGCCGTTTTGTTCTTGCGCGGATTCGGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGCTCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAACAGATTACCCAGCTGCTTGAGGACATTCTCGTTGTCCGGGGTGGTGACCAGGCTGTGGTTCAGCACACTCTTGCCATCGACACCAATGTgg                                                                                                                                                >  2:326251/1‑251 (MQ=14)
                                                                                                       ggCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCt                                                                                                                                            <  1:481666/251‑1 (MQ=255)
                                                                                                       ggCCGTTGTGTGCTTGCGGGGATTCTGCTTAATGGTTTTGACCTTGCCGGGACCCTCGGCGATCAGCCTGTCCACATCCACCTAGGCCATCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGTCCGCGGTGGTGACCAGGCTGTGGGTCAGGCAACTCTTGGCATCGACACCAATGTGGGCCt                                                                                                                                            <  2:723847/251‑1 (MQ=14)
                                                                                                       ggCCGTTATGTTATTGCGCGGAATCTGCTTCAAGGATCTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCGGCGGCGCTCCTTGGTAGCCGGCATCGGCTGAGGCAAATTGCTTCTCTCCATGAAGCAGGTCACCCAGCTGCTTGGGGTCATGCTCGTTGGCCGTGGTGGTGACCCTGCTGGGGGTCAGGCAACTCTTGGCATAGACTCCACTGTGGGCCt                                                                                                                                            >  2:516765/1‑251 (MQ=14)
                                                                                                        gacGTTTTGTTCTTGCGCGGTTTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGAGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCTACACCAATGTGGGCCtt                                                                                                                                           >  2:635050/3‑251 (MQ=255)
                                                                                                           gTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCAt                                                                                                                                        >  1:219131/1‑251 (MQ=255)
                                                                                                                   cTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGt                                                                                                                                >  1:24465/1‑251 (MQ=255)
                                                                                                                     tGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGc                                                                                                                              >  1:44770/1‑251 (MQ=255)
                                                                                                                     tGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGGGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACCCCCATGTGGGCCTTCAGGCCAAAgggc                                                                                                                              >  2:22337/1‑248 (MQ=255)
                                                                                                                        gcgGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCCGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCAc                                                                                                                           >  2:788376/1‑251 (MQ=255)
                                                                                                                        gcgGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAACGTGCCAc                                                                                                                           >  1:34657/1‑251 (MQ=255)
                                                                                                                        gcgGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCAc                                                                                                                           >  1:389007/1‑251 (MQ=255)
                                                                                                                        gcgGATTCGGCTTCAAGGTTTTTTCCTGGCCGGGACGCTCGGCGAGCAGCCAGTCCACATCCTCCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTTCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCAc                                                                                                                           <  2:211613/251‑1 (MQ=14)
                                                                                                                         cgGATTCTGCTTCAAGGTTTTTACCGTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGATGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACt                                                                                                                          <  1:75051/251‑1 (MQ=255)
                                                                                                                          ggATTCTGCTTTAAGGTTTTGCCCTTGCCGTGCCGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCGTGGTAGCCGGCATCGGTTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCTCTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTg                                                                                                                         <  2:505192/251‑1 (MQ=11)
                                                                                                                          ggATTCTGCGTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTACACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCCCTCTTGGCATCGACACCAATTGGGGCCTTCATGCCAAAGTGCCACTg                                                                                                                         >  2:239200/1‑251 (MQ=255)
                                                                                                                           gATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGa                                                                                                                        <  1:153301/251‑1 (MQ=255)
                                                                                                                                tGCTTAAAGGTATTTACCTTGCAGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCAAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCCGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGTTTgac                                                                                                                   >  2:70039/1‑249 (MQ=11)
                                                                                                                                      aaGGTTTTTACCTTGCCGGGCCGCTGGGCGTCTAGCCAGTCCACTTCCACCTCGGCCATCTCCTCGCGCGGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTGCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGTCACCAATGTGGGCATTCATGCCAAAGTGCCACTGATTGCCTTTCtt                                                                                                             <  2:219451/251‑1 (MQ=11)
                                                                                                                                      aaGGTTTCTACTTTGCCGGGACGCTCGGCGAGCGGCCGGTCCACTACCACCTCGGCCATCTCCTCGGGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTTGCAGCGACACCAATGTGGGCCTTCAGGCCAAAGTGCCACTGATTGCCTTTcgt                                                                                                             <  2:864551/251‑3 (MQ=11)
                                                                                                                                         gTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGt                                                                                                          >  1:13599/1‑251 (MQ=37)
                                                                                                                                         gTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGt                                                                                                          >  1:863760/1‑251 (MQ=37)
                                                                                                                                         gTTTTTAACTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTACGGCCACTCTTGGCATCGACACCACTGTGGGCCTTCATGCCAAAGGGCCACTGATTGCCTTTCTTGGt                                                                                                          >  2:355607/1‑251 (MQ=255)
                                                                                                                                                 cTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGCCCAGGCTGTGGGTCAGGCCACTCTTGGCAACGACACCAATGTGGGCCTTCATGCCAAAGAGCCCCTGAGTGCCTTTCTTGGTCTGATGc                                                                                                   >  2:133248/1‑250 (MQ=255)
                                                                                                                                                    gCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCAtct                                                                                               >  1:448521/1‑251 (MQ=21)
                                                                                                                                                     ccGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCAtctc                                                                                              >  1:799608/1‑251 (MQ=21)
                                                                                                                                                     ccGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCAtctc                                                                                              <  1:758690/251‑1 (MQ=21)
                                                                                                                                                      cGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTcc                                                                                             <  1:504475/251‑1 (MQ=21)
                                                                                                                                                      cGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCGCCTTGGTAGCCGGCATCGGCTGAGACAAAGTGCTTCTCTCCGTGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCATGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCAtctac                                                                                             <  1:492328/251‑3 (MQ=255)
                                                                                                                                                            ccTCGGCGAGCAGCCAGTCCACATCCCCCTCGGCCAGCGCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCGTCTCTCCAGGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGAGGCATCTCCGGATcg                                                                                       <  2:543015/250‑1 (MQ=11)
                                                                                                                                                                   gATCAGCCAGTCCACATCCACCGCGGCCAGCTAATCGCGCTGTGGCGCTCCTTAGTAGCCGGCATCGGCTGAGACACATTGCTTCTCTCCATGAATCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGAGGTGACCAGGCTGAGGGTCCGGCCACTCTTTGCATCGGAACCAATGTGGGCCTTCCTGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGCTGCATCTCCGGATcgtggtgat                                                                                >  2:274469/1‑244 (MQ=11)
                                                                                                                                                                    aTCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTTGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAACTTGCTTCTCTCCATGAAGCAGATTCCCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGTCCTTCATGCCAAAGGGCCACTGATTGCCTTGCTTGGTCTGATGCATCTCCGGATCGCGTtgctg                                                                               >  2:457732/1‑251 (MQ=11)
                                                                                                                                                                      cAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTgctgct                                                                             <  1:833454/251‑1 (MQ=18)
                                                                                                                                                                        agctaGTCCACATCCTCCCCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCCCTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGctc                                                                            <  1:428648/247‑1 (MQ=255)
                                                                                                                                                                         gCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGctct                                                                           >  1:659972/1‑251 (MQ=18)
                                                                                                                                                                          ccAGTCCACATCCCCCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCAGGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCtt                                                                          <  2:448521/251‑1 (MQ=18)
                                                                                                                                                                               ccACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGtt                                                                      <  1:234726/250‑1 (MQ=18)
                                                                                                                                                                                  caTCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGATTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGGCATGCTCGTTGGCAGCGGAGGTGACCAGGCTGTGGGTCAGGCCACTATTGGCATCGACACCAATGTGAGCCTTCATGCCAACGTGCCACTGCTTGCATTTCTTGGTCTGATGCATCTCCGGATCGCGTAGCTGCTCTTTGTTCTTg                                                                  >  1:102033/1‑251 (MQ=11)
                                                                                                                                                                                   aTCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTgg                                                                 >  1:761418/1‑251 (MQ=18)
                                                                                                                                                                                      ccccccgcccAGCTCCCCGCGCTGTGGCCCGCCTTGGTCGCCGGCATCGGCTGAGACAAATTGCTGCTCTCCAGGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGTCTGTGGGTCAGGCCACTCTTGGCATCGACGCCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTag                                                              <  2:44770/243‑1 (MQ=11)
                                                                                                                                                                                      cACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGt                                                                >  2:674407/1‑249 (MQ=18)
                                                                                                                                                                                          ccgtccAGCTCCTAGCGCTGTGGCGCTCCTTGGTTTCAGGCAGCGGCTGAGACAATTGGCTTTTCGCCATGAAGCATATTACCCAGCTGATGGAGGTCAGGCTCGTTGGCCGCGGTGGTGTCCAGGCTGTGGGTCAGGCCTCTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGACTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCGCTTTGTTCTTGGTAGAGCt                                                          <  2:15505/247‑1 (MQ=11)
                                                                                                                                                                                           cGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTCGAGCTg                                                         <  1:752618/251‑1 (MQ=14)
                                                                                                                                                                                           cGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTg                                                         >  2:222368/1‑251 (MQ=21)
                                                                                                                                                                                             gCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGACGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGCTTGAGGTAATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTggg                                                       >  1:320679/1‑251 (MQ=255)
                                                                                                                                                                                                aGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCTGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACACAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGAGCCCAGGATGTGGGTCAGGCCACTCTTTGCATCGACCCCAATGTGGGCCTTCATGCAAAAGTGCCACTGCTTGCATTTCTTTGTCTGATGCATCTCCGGATCGCGTTGCTTCTCTTTGTTCTTTGGAGAGCTGGGTGc                                                    >  2:33203/1‑251 (MQ=11)
                                                                                                                                                                                                    cctcTCGCTGTGGCGCTCCGTGTTAGCCGGCATCGGCTGAGACAAATTGCTTCTCGCCATGAAGCAGATTACCCAGCCGATTGAGGGCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATTTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCGTGGTAGAGCTGGGTGCCt                                                  <  2:34657/249‑1 (MQ=11)
                                                                                                                                                                                                    cctcGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAACTTGCTTCTCTCCATGAAGCAGATTACCCAGCTGCTTGAGGTCATGATCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCACTGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCAGTTCGCGTTGCGGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTc                                                 >  2:513365/1‑250 (MQ=11)
                                                                                                                                                                                                       cgcgCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAatg                                             <  1:725544/251‑1 (MQ=21)
                                                                                                                                                                                                                 cgcTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCAt                                   <  1:215882/251‑1 (MQ=21)
                                                                                                                                                                                                                     ccTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTCGAGCTGGTTGCCTCAATGATGTTGTCATccac                               >  2:378725/1‑251 (MQ=9)
                                                                                                                                                                                                                     ccTTGGTAGCCGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATccac                               <  1:18507/251‑1 (MQ=21)
                                                                                                                                                                                                                          gTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGCAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTTGTAGAGCTGGGGGCCTCAATGATGGTGGCATCCACCaa                            >  2:815411/1‑249 (MQ=14)
                                                                                                                                                                                                                              ccGGCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGACAGGCCACTCTTGGCATCGACACCAATGAGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGCTCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCTCCAAAGTGcc                      >  2:634152/1‑251 (MQ=255)
                                                                                                                                                                                                                                ggCATCGGCTGAGACAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCTGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCtt                    <  1:69947/251‑1 (MQ=37)
                                                                                                                                                                                                                                  cATCGGCTGACACAAATTGCTTCTCTCCTTGAAGATTAATACCCAGCTGCTTGAGCGCATGCTCCTTGGCCGGGGTGGTGACCAGGCTGTGGGTCAGTACACTCTTGGCATGGACACCAATGTGGGCCTTCATGCCAAAGTGTCACTGATTGCCTTGCTTTGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTTCCTTgg                  <  2:189986/251‑1 (MQ=11)
                                                                                                                                                                                                                                          gagaCAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTcatcat          >  1:325489/1‑251 (MQ=38)
                                                                                                                                                                                                                                              cAAATTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTgg                  >  1:822089/1‑239 (MQ=38)
                                                                                                                                                                                                                                                aaTTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGcc    >  1:505627/1‑251 (MQ=38)
                                                                                                                                                                                                                                                 aTTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCt   >  1:242829/1‑251 (MQ=38)
                                                                                                                                                                                                                                                 aTTGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCt   >  1:757011/1‑251 (MQ=38)
                                                                                                                                                                                                                                                   tGCTTCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTg  <  1:88086/250‑1 (MQ=39)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
GCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTACCCAGCTGATTGAGGTCATGCTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTG  >  NC_000913/1427841‑1428332

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