Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,193,177 C→T E134K (GAG→AAG)  thiF ← sulfur carrier protein ThiS adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,193,1770CT98.4% 425.0 / NA 122E134K (GAG→AAG) thiFsulfur carrier protein ThiS adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (57/63);  minor base G (0/2);  total (57/65)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 122 total aligned reads displayed.

GCAGGTGTCTCTATACCGCTTAATAACTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCC  >  NC_000913/4192930‑4193420
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
gCAGGTGTCTCTATACCGCTTAATAACTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATACCCCCGCCCCCCGGGCCACCCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTg                                                                                                                                                                                                                                                    <  2:290383/251‑1 (MQ=255)
      gtCTCTATCCCGTTTAATACCTTAATGCCTTCAAGTGCCGGCAAAGTTCCCACCACCCCGACCCCCGGGCCAACCGCGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCGGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCGCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTGCTGGCGAGt                                                                                                                                                                                                                                              <  2:257483/251‑1 (MQ=255)
          ctATACCGCTTAATAACTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGTTAAGCGGCGGGTTGAGTGCCCCGCAGGCGGCATTAATCTTCTGGCGAGTCGcc                                                                                                                                                                                                                                          >  1:527080/1‑251 (MQ=255)
           tataCCGCTTAATAACTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGAAAAGCGGCGTGTTGGGTGCCACGCCGGCGGCACTAATCTTCTGGCGAGTCGccc                                                                                                                                                                                                                                         >  2:25225/1‑250 (MQ=255)
               ccGCTTAATAACTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATAtt                                                                                                                                                                                                                                     <  1:394522/251‑1 (MQ=255)
                cGCTTAATAACTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCCTTAATCTTCTGGCGAGTCGCCATATTg                                                                                                                                                                                                                                    >  1:615950/1‑251 (MQ=255)
                 gCTTAATAACTAAAGGGCTCCCAGTGCCTGCAAAGTGCCCAGACCCCCGCCCACCGGGCCACCCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGt                                                                                                                                                                                                                                   <  2:229046/251‑1 (MQ=255)
                    taataaCTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCg                                                                                                                                                                                                                                 <  1:181189/250‑1 (MQ=255)
                    taataaCTTAAGGGCTTCCAGTGCCTGCAAAGTGCCCAAACCCCCGCCCACCGGGCCAACCCCGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCg                                                                                                                                                                                                                                 <  2:136962/250‑1 (MQ=255)
                             aaTGGCTTCCAGTGCCTGCAAAGAGCCCATAACCACGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGCTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCCTCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCCTTAATCTTCTGGCGCGTCGCCCAATTTTCGGTACAGTccg                                                                                                                                                                                                                       >  2:316329/1‑249 (MQ=255)
                              aTGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGGGTCGCCCTATTGTCGGGACAGTCg                                                                                                                                                                                                                       >  2:312339/1‑249 (MQ=255)
                              aTGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGGGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCCTATTGTCGGTACAGTCg                                                                                                                                                                                                                       >  1:713008/1‑249 (MQ=255)
                               tGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGGGCCACGCAGGCGGCATTCATCTTCTGGCGAGTCGCCATATTTTCGGTACCGTCGAGc                                                                                                                                                                                                                    >  2:386707/1‑251 (MQ=255)
                                 gCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTTTCGGTACAGTCGAGcac                                                                                                                                                                                                                  <  2:718536/251‑1 (MQ=255)
                                     ccAGTGCCTGCAAAGTGCCCATAACCCCGCCCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACcaca                                                                                                                                                                                                              <  1:312064/251‑1 (MQ=255)
                                     ccAGTGCCTGCAAAGTGCCCATAACCCCGACCCCCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGcaccac                                                                                                                                                                                                               <  1:25225/250‑1 (MQ=255)
                                     ccAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCCTTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGcaccacc                                                                                                                                                                                                              >  2:626233/1‑250 (MQ=255)
                                             tGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCg                                                                                                                                                                                                     >  2:132688/1‑251 (MQ=255)
                                              gCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCCCGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCgg                                                                                                                                                                                                    >  2:475513/1‑250 (MQ=255)
                                               cAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACAGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTTCCACGCAGGCGGCATTAATCTTCTGGGGAGGCGACCTATTGTCTGTACAGTCTGACACCACCTCGGCCCggg                                                                                                                                                                                                   >  2:567024/1‑249 (MQ=255)
                                                    ggCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAAcc                                                                                                                                                                                              <  1:262430/250‑1 (MQ=255)
                                                    cgCCCATACCCCCGACCCCCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAAcc                                                                                                                                                                                              <  1:53253/250‑1 (MQ=255)
                                                     ccccATACCCCCGACCCCCGGTCCACCCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCg                                                                                                                                                                                             <  2:86812/250‑1 (MQ=255)
                                                         aTAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCCTTAATCTTCTGGCGAGTCGCCCTATTGTCGGGTCAGTCGAGCACCCCATCGGCCCCTGCAACCGCATc                                                                                                                                                                                         >  1:607989/1‑251 (MQ=255)
                                                           aaCCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGGGCACCACATCGGCCCgggccaacgcctctt                                                                                                                                                                                       >  2:169038/1‑237 (MQ=255)
                                                                    caccGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGccc                                                                                                                                                                              >  2:530871/1‑250 (MQ=255)
                                                                    caccGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCt                                                                                                                                                                              <  1:132688/251‑1 (MQ=255)
                                                                     accGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTTTCGGTACCGTCGCGCCCCACATCGGCCCGTGGAACCCCACCtttttaagcccc                                                                                                                                                                             >  2:567217/1‑243 (MQ=255)
                                                                       cGGGCCAACCACGCCCGAAGTGAGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAACCGGCGTGTTGAGGGCCACGCAGGCGGCATTAATCTTTTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGACTccc                                                                                                                                                                           >  2:504103/1‑249 (MQ=255)
                                                                        gggCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCCTCTTTTAACGCCTCAcc                                                                                                                                                                          >  1:672006/1‑251 (MQ=255)
                                                                           ccAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGGACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGt                                                                                                                                                                       >  1:190280/1‑251 (MQ=255)
                                                                             caCCACGCCCGCCGTGCGGCAGTTGCTTTCTGGCTCCTGGTTATCTGTCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTCCCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTa                                                                                                                                                                     <  2:456557/250‑1 (MQ=255)
                                                                                  cgcccgccGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCAGTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAAcc                                                                                                                                                                  <  1:98025/249‑1 (MQ=255)
                                                                                     ccgccgTGCGGCAGTTGCGTTCTGGCTCCTGTGTATCGGGCCACAGGCAGCGGTAACACCCCGGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATGAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCTGCCCGTGCAACCGCAGCTTTTAACGCCTCACCCGTTAACCgttgt                                                                                                                                                             <  2:295221/251‑1 (MQ=255)
                                                                                       gccgTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGttgttg                                                                                                                                                           <  1:512143/251‑1 (MQ=255)
                                                                                         cttGCGGCAGATGCGTTCTGGCTCCTGGTTATCTGGCCACCGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGAACGCCAAATCCGACCGCGCTGGCGGCGATAACCGGAGCGTTGAGTGCCACGCCGGCGGCATTAATCCTCTGGCGAAGAGCCATTTTTGCGGTACCGGCGAGCACCACATCGGCCCCGTCAACCGCATCTTTTAAACCCTCCCCCCTTAAcccgtgatgt                                                                                                                                                          >  2:276723/3‑242 (MQ=255)
                                                                                          gTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTa                                                                                                                                                         >  1:87982/1‑250 (MQ=255)
                                                                                          gTGCGGCAGGTGCGTTCGGGCTCCGGGTTATCGGGCCACAGGCAGCGGTAACACCCCTGCCCCCAGGGCGGCGTCAGTACCAGCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTaa                                                                                                                                                        <  2:589319/251‑1 (MQ=255)
                                                                                            gCGGCAGTGGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCGCCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAAt                                                                                                                                                       <  2:286833/250‑1 (MQ=255)
                                                                                             cGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGc                                                                                                                                                     <  1:687276/251‑1 (MQ=255)
                                                                                              ggCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCTGCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGATTGCCTCGCAGGCGGCATTAAGCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCt                                                                                                                                                    <  2:19932/251‑1 (MQ=255)
                                                                                               gCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCCACCGCATCTTTTAACTCCCCACCCGTTAACCGTTGTTGTAATGCtt                                                                                                                                                   >  2:486553/1‑250 (MQ=255)
                                                                                                       gTTCTTGCTCCCGTTTATCGGGCCACCGGCAGCGGTAACACCACTGCTCCCAGGGCGGCGTCCTTACCACCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTATCCGTTGTTGTAATGCTGTCAGTTGa                                                                                                                                           <  2:406927/251‑1 (MQ=255)
                                                                                                         tCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGaa                                                                                                                                          <  2:610814/250‑1 (MQ=255)
                                                                                                               tCCTGGGTATCTGCCCACAGGCAGCGGTACCTCCCCTGCTCCCAGGGCGGCGTCAGTTCCACCAACTGCCCGCCAAATCCGACCGCGCTGGCGGTTATAAGCGGCGTGTTGAGTGCCCCGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCAGCACATCGGCCCGTGCAACCGCATCTTTGAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCgg                                                                                                                                    <  2:399039/250‑1 (MQ=255)
                                                                                                                ccTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCggg                                                                                                                                   <  1:133530/250‑1 (MQ=255)
                                                                                                                  tggatATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAAATGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGAGCCCCGCAGGCGGCATTAATCTTCTGGCGAATCCCCATATTGGCGGTAAAGTCGAGCACCACATCGGCCCGGGCAACCGCATCTTTTAACGCCCCACCCGTAAACCGTTGGTGTACTGCTGGCACTTGGATGTCggggtt                                                                                                                                >  2:497850/5‑248 (MQ=255)
                                                                                                                      tagcTGGCAACAGGCAGCGGTAACACCCCGGCTCCGAGGGCGGGGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGGCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCGTTTAACGCCTCACCCGGTAACCGGTGTTTTAATGCTGTCAGTTGAATGTCGGGATTCaa                                                                                                                             <  2:466750/247‑1 (MQ=255)
                                                                                                                      tATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCaa                                                                                                                             <  2:634662/250‑1 (MQ=255)
                                                                                                                      tATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGGTGTAATGCTGTCAGTTGAATGTCGGGATTCaa                                                                                                                             <  2:430817/250‑1 (MQ=255)
                                                                                                                           ggCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAGCCGCGTCGTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACtgt                                                                                                                         <  2:657336/249‑1 (MQ=255)
                                                                                                                           ggCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACtgt                                                                                                                         <  1:530871/249‑1 (MQ=255)
                                                                                                                           ggCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACt                                                                                                                           <  1:647101/247‑1 (MQ=255)
                                                                                                                            gCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGGACAGTCGAGCACCACATCGGCCCGTGCAACCGCCTCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTc                                                                                                                      >  2:213645/1‑251 (MQ=255)
                                                                                                                             ccACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATTCTGTCAGTTGAATGTCGGGATTCCACTGTGTcc                                                                                                                     >  2:228618/1‑250 (MQ=255)
                                                                                                                             ccACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAAGCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCa                                                                                                                     <  2:175497/251‑1 (MQ=255)
                                                                                                                                caGGCAGCGGTAACCCCCCGTCATCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAAGCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGGGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCGCCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCTACTGTGTCAGTc                                                                                                                  <  2:641450/251‑1 (MQ=255)
                                                                                                                                      gCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAAGGTCGGGATTCAACTGTGTCAGTCGCTGtt                                                                                                            <  1:214254/251‑1 (MQ=255)
                                                                                                                                        ggTAACGCCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTGCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGGTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGGTGTTgg                                                                                                          <  1:83841/251‑1 (MQ=255)
                                                                                                                                          aaaccccccTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCAGCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCt                                                                                                        <  2:146834/246‑1 (MQ=255)
                                                                                                                                               ccccTGCTCCCAGGGAGGCGTCAGTACCATCAACTGACCGCCAACTCAACCCGCGCTGGCGGTGACACGCGGCGTGTCGAGTGCCCCGCAGGGGGAGTTAATCTTCTGGCGAGACCCCATATTTCCGGAACAGTCAAGCCCCACATCGGCCCGCGCAACCGAATCTTTTAACGACTCACCAGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGCTTCCACTGTGTCAGTCGCTGTTGGCTGACCt                                                                                                   >  1:421923/1‑251 (MQ=255)
                                                                                                                                                cccTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTg                                                                                                  <  1:626233/251‑1 (MQ=255)
                                                                                                                                                   tGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCg                                                                                                >  1:290140/1‑250 (MQ=255)
                                                                                                                                                    gCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACAGCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGa                                                                                               <  2:719733/250‑1 (MQ=255)
                                                                                                                                                     cTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGAtt                                                                                             >  1:9603/1‑251 (MQ=255)
                                                                                                                                                          aGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCCCGCAGGAGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGCTTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCggg                                                                                        >  1:737177/1‑251 (MQ=255)
                                                                                                                                                            ggcggcGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTGGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCGCCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGCTGTAATGCTGTCAGTTGAATGTCGGGATTCAACGGTGTCAGTCGCTGTTGGCTGACCGGCGATTTCGGGcg                                                                                      <  2:200798/251‑1 (MQ=255)
                                                                                                                                                                 cgTCAGTACCATCCACAGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGGGCACCACATCGGCCCGTGCCACCGCCTCTTTTAAAGCCTCACCCGGTAACCCGTGGTGTAATTCTGTCAGGTTAAAGTCGGGGTTCAACTGTGTAAGTCGCTTTTGGCTGACCTGAGATTTCGGGcgattgg                                                                                 >  2:142840/1‑248 (MQ=255)
                                                                                                                                                                  gTCAGTACCATCAACTGGCCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCGGTCAGTTGAATGTCGGGATGCAACTGTGTCAGTCGCTGTTGGCTGACCGGCGATTTCGGGcgatcgat                                                                                <  1:259752/251‑1 (MQ=255)
                                                                                                                                                                     aGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCCGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCTGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATc                                                                             >  1:544264/1‑251 (MQ=255)
                                                                                                                                                                      gTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGCTAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTCACCGTTGTTGTAATGCTGTCAGTTGTATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCt                                                                            >  2:493902/1‑251 (MQ=255)
                                                                                                                                                                        aCCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCtt                                                                           <  1:66712/250‑1 (MQ=255)
                                                                                                                                                                            tCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGGTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTg                                                                      >  1:225451/1‑251 (MQ=255)
                                                                                                                                                                                            cccGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTc                                                      <  1:332517/250‑1 (MQ=255)
                                                                                                                                                                                             ccGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCGTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTc                                                      <  2:343122/250‑1 (MQ=255)
                                                                                                                                                                                               gACCGCGCTTGAGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAGCTTCGGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACAGCGGCCCGTGCAACCGCATCTTTTAACGCCGCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGtt                                                   <  1:358230/251‑1 (MQ=255)
                                                                                                                                                                                               gACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTTAGTCGCTGTTGGCTGACCTGCGGTTTCGGGCGATAGATATCTTCCGTGGTAAAGAGGAtttagttcg                                                   >  2:289670/1‑245 (MQ=255)
                                                                                                                                                                                               gACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGAttt                                                         <  1:289670/245‑1 (MQ=255)
                                                                                                                                                                                                aCCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTg                                                  <  1:386707/251‑1 (MQ=255)
                                                                                                                                                                                                aCCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTg                                                  >  2:239054/1‑251 (MQ=255)
                                                                                                                                                                                                aCCGCGCTGGCGGGGATAAGCGGCGAGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCCCCACATCGGCCCGCGCCACCGCCTCTTTTAACGCCTCCCCCGTTAACCGTTGTTTTAATGCTGTCCGTTGAATGTCGGGATTCCAATGTGTCAGTCGCTGTTGGCTGACCTGCGCTTTCGGGCGATCGATATCTTCCGGGGTAAAGTGGATTTGTCGTTg                                                  >  2:525893/1‑251 (MQ=255)
                                                                                                                                                                                                  cgcgCTGGCGGTGATAAGCGGCGGGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAAACGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAAATGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGTGCGATATCTTCAGGGGGAAAGAGGTTTTGTCGTTGcg                                                >  2:543605/1‑250 (MQ=255)
                                                                                                                                                                                                        ggCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTg                                           <  1:643739/250‑1 (MQ=255)
                                                                                                                                                                                                        ggCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTg                                           <  1:493902/250‑1 (MQ=255)
                                                                                                                                                                                                                 aaGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGc                                  <  2:300207/250‑1 (MQ=255)
                                                                                                                                                                                                                      gcgTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGACATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATCTCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACAtcg                            >  2:369788/1‑251 (MQ=255)
                                                                                                                                                                                                                              aGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATc                    <  1:654713/251‑1 (MQ=255)
                                                                                                                                                                                                                              aGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATc                    <  1:504103/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCCCCACATCGGCCCGTGCAATCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTg                  <  1:497850/250‑1 (MQ=255)
                                                                                                                                                                                                                                  ccACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGcc                <  1:312339/251‑1 (MQ=255)
                                                                                                                                                                                                                                   cACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCa               <  1:653607/251‑1 (MQ=255)
                                                                                                                                                                                                                                     cGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGt             >  1:540017/1‑251 (MQ=255)
                                                                                                                                                                                                                                      gCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGGCGAGCACCACATCGGCCCGTGCACACGCATCTTTTAACGCCTAACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTAACCTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGt             >  2:478062/1‑250 (MQ=255)
                                                                                                                                                                                                                                           gggcATTAACCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGAACGTGCAACCGCCTCTTTTACCGCCTCACCCGTTAACCGGTGTTGTAATGTTGTCAGTTGAATGTCGGGATTCATCTGTGTCAGTCGCTGTTGGCTGATCTGCGATTTCGGTCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATAGTCGTCATCTGCCAGTACCAGc       <  1:565337/250‑1 (MQ=255)
                                                                                                                                                                                                                                           cggcATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCCCCCGTTAAACGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGc       >  2:662776/1‑251 (MQ=255)
                                                                                                                                                                                                                                           cggcATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGAGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCACCTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGc       >  2:55700/1‑251 (MQ=255)
                                                                                                                                                                                                                                            ggcATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCg      <  2:530133/251‑1 (MQ=255)
                                                                                                                                                                                                                                             gcATTAATCTTCTGGCGAGTCGCCATATTTTCGGTACAGTCGAGCACCACATCGGCCCGTGCATCCGCATCTTGTAACGCCTCGCCCGTTAACCGTTGTTGTAATGCTGTCAGTGGAATGTCGGGATTCAACTGTGTCAGTGGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCGGAGTGCTTAAAGGCACAGCGTTTTCAGCGGGCAGTACCAGCGt     <  2:732120/251‑1 (MQ=255)
                                                                                                                                                                                                                                               aTTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTcc   >  1:597824/1‑251 (MQ=255)
                                                                                                                                                                                                                                                ttAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTccc  >  1:172178/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
GCAGGTGTCTCTATACCGCTTAATAACTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCC  >  NC_000913/4192930‑4193420

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: