Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,080,042 T→C 9.0% E531G (GAG→GGG)  tktA ← transketolase 1, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,080,0420TC 9.0% 359.9 / 2.3 130E531G (GAG→GGG) tktAtransketolase 1, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base T (53/62);  new base C (9/5);  total (63/67)
Fisher's exact test for biased strand distribution p-value = 2.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 132 total aligned reads displayed.

AGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGT  >  NC_000913/3079804‑3080278
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                              
ggcatctggTGCAGTACCCGTTTTCGGCAGTTCGGTATCCGGGTAAGCAGGCTCCTTCTTGGACAGTGGGTCGGTAGACGGCATGGAACCCACGCGCGCTTTCCCGCCTTCGGCAGTCAGTTTTGCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGct                                                                                                                                                                                                                                                <  2:674509/231‑1 (MQ=255)
   acggggTGCATTAACCGCTTTCGGCATTACGGATTCCCGGTAAGCAGCCTCCTGCTTGTCAAAGGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCtt                                                                                                                                                                                                                                             <  2:353070/234‑1 (MQ=255)
    aCGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAACTCAGTTCCGGCTGACCGGCGCAGCCTTTCAGCACATAACCACCGCCCGCGATGTTTGCCAGTGGCtttt                                                                                                                                                                                                                                             >  1:329759/1‑238 (MQ=255)
    aCGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGGCCGATGTTTGCCAGTTGCCCTTc                                                                                                                                                                                                                                            >  2:173642/1‑239 (MQ=255)
    aCGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCCCTTc                                                                                                                                                                                                                                            >  1:893758/1‑239 (MQ=255)
            aGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGGCTTTCAGCCCATAACCCCCGCGCGCGATGTTTGCCAGTTGCTCTTCAgttcgt                                                                                                                                                                                                                                     >  1:547152/1‑238 (MQ=255)
                    cTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCtgctgc                                                                                                                                                                                                                            >  1:48432/1‑239 (MQ=255)
                      ttCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGGTGTTTGCCAGTTGGCTTTCAggtcggtctggccgggc                                                                                                                                                                                                                          >  2:465628/1‑230 (MQ=255)
                       tcggcgtttcGGATTCCCGGTAAGCAGCCTCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACCCCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGcc                                                                                                                                                                                                                         <  2:361843/230‑1 (MQ=255)
                              caCGGATTCACGGTAAGCAGCATCCTGATTGTCAAATGCGTCGGTAGACGGCATGGACCCCACGCGCGCTTTCACGCCTTCGGCAGTCAGATTTTCGTAGCCAGCACCAGCCAGTTCAACTTCTGAACCGGTAGCGCTGAAAATCAGTTCCGGCTGACCCGCGCAGCCTTTCAGCACATCACCACCCCCCCCCCTGTTTGCAGGTTCCCCTTCCgttcgttctcgctcccccctgtttc                                                                                                                                                                                                                  >  2:383148/2‑222 (MQ=255)
                                     tcccggtaagcgcctctccgtttGTCAATTGGGTGGGTGGAGGGGCTGGACCCCACGCGCCTTTGCAGGCCTTCGGCAGTCAGTTTTTCGTAGGACGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCGCCGCGCGCGATGTTTGCCAGGTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGga                                                                                                                                                                                                           <  2:11912/218‑1 (MQ=255)
                                     tcccGGGAAGCCGCACCCTGTTTGTCAAATGGGTGGGTAGACGGCATGGCCCCCACGCGCGCTTTCACGCCTTCGACAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGga                                                                                                                                                                                                           <  2:706762/236‑1 (MQ=255)
                                     tCACGGTAAGCAGCATCCTGGTTGGCAAATGGGTCGGTAGGCGGAATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGga                                                                                                                                                                                                           <  1:406678/239‑1 (MQ=255)
                                     tCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGga                                                                                                                                                                                                           >  2:18497/1‑239 (MQ=255)
                                        cGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGggag                                                                                                                                                                                                        >  1:981254/1‑239 (MQ=255)
                                             agcagcATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACAACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGGTCTGACGGGGGAGGAtag                                                                                                                                                                                                   >  1:738697/1‑237 (MQ=255)
                                              ccctcctCCTTCTTTTCAAAGGCGTGGGGAGCCGGCAGGGACACCACGGGCGCTTTCACCCCTTGGGCAGTCAGTTTTGCGTAGGCCGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAg                                                                                                                                                                                                  <  2:893758/233‑1 (MQ=255)
                                                agctttcTGGTTTTCAAATGGGTGGGGTGGCGGACTGGACCCCACGGGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAg                                                                                                                                                                                                  <  2:226829/231‑1 (MQ=255)
                                                     ccTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACAAAACACCCGCGCGCGCTGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGCCGGGGGGGGATCAgggggggg                                                                                                                                                                                           >  2:1083429/1‑232 (MQ=255)
                                                       tGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTCCTCCGCCAGGTTCTGACGGGAGAGGATCAGTCCGGTCgg                                                                                                                                                                                         >  1:1002749/1‑239 (MQ=255)
                                                        tcTTGTCAAATGGGCGGGGAGGCGGGATGACCCCCCCCGCGCTTTTCCCGCCTGCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCggg                                                                                                                                                                                        <  2:615760/238‑1 (MQ=255)
                                                        gCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTGCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCggg                                                                                                                                                                                        <  1:761912/239‑1 (MQ=255)
                                                         cTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGAGCAGTGGgggcgggc                                                                                                                                                                                       >  2:65632/1‑239 (MQ=255)
                                                         cTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGCCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCCCCAGGTTCTGCCGGGAGAGGATCAGTGCGGTCggg                                                                                                                                                                                        >  1:480708/1‑238 (MQ=255)
                                                         cTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGAGGAAAATCAGTTCCGGCTGACCGGCGCAGGCTTTCAGCACATAACCACCGCGCGCGCGGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCCCCAGGGTCTTACGGGGGAGGGTCAGGGGGGGCggg                                                                                                                                                                                        >  2:573149/1‑238 (MQ=255)
                                                            gTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCACCAGCCAGGTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCCGTTCTTTCCTGCGCCCCCAGGTTCTGATGGGGGAGGAGCCGTGGGGTCGGGCCg                                                                                                                                                                                     >  2:1020352/1‑238 (MQ=255)
                                                            gTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGGAAGGATCAGGCCgggcggggcg                                                                                                                                                                                     >  2:309544/1‑238 (MQ=255)
                                                                atgcggtcggtagcggggcgggACCCCACGCGCTCTTTCACCCCTTCGGGAGTCAGTTTTTGGTAGGCAGCAACGGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTg                                                                                                                                                                                <  2:83521/219‑1 (MQ=255)
                                                                actGCGTCGGTAGGCGGTATGGACACCACGCTCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTg                                                                                                                                                                                <  1:269729/237‑1 (MQ=255)
                                                                aaTGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTg                                                                                                                                                                                <  1:465628/239‑1 (MQ=255)
                                                                 aTGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCCCGCGCGATGTTTGCCAGTTGCCCTTCATTTCGTTCCTGCTGCGCCAGGTTCTGCCGGGAGGGGATCAGGGGGGGCGGGCCGTCCTGa                                                                                                                                                                               >  2:1053943/1‑239 (MQ=255)
                                                                      tCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCTCGCGCGATGTTTGCCAGTTGCGCTTCAGTTCGTTCCTCCTGCGCCAGGTTCTGACGGGAGAGGAGCAGGGCGGTGGGCCCGTCCGGAcccct                                                                                                                                                                          >  1:430456/1‑235 (MQ=255)
                                                                      tCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTc                                                                                                                                                                          >  1:909704/1‑239 (MQ=255)
                                                                        ggtaagcGGCATGGCACCCACGCGCGCTTTCACGGCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCaa                                                                                                                                                                        <  1:65632/233‑1 (MQ=255)
                                                                        ggTAGACGGCATGGACACCACGCGCGCTTCCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCaa                                                                                                                                                                        <  1:68469/239‑1 (MQ=255)
                                                                        cgttgccGGCATGGCCCCCCCGCGCGTTTTCACGCCCGCGGCAGTCAGTTTTTCGGAGGCAGAAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCCCCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCaa                                                                                                                                                                        <  2:1040271/233‑1 (MQ=255)
                                                                           aggcGGCATGGACCCCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAaca                                                                                                                                                                      <  1:929053/235‑1 (MQ=255)
                                                                             aCGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCACAACCg                                                                                                                                                                  >  1:360545/1‑239 (MQ=255)
                                                                                   tGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGCGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTAttt                                                                                                                                                             >  1:385861/1‑238 (MQ=255)
                                                                                   tGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTGGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTCCTGACGGGAGGGGCTCAGGGCGGTCGGGCGGGCCGGACGGCCAACACCGCATGTc                                                                                                                                                            >  2:461314/1‑239 (MQ=255)
                                                                                         ccccccgctcTTTCACGCCTTCGGCAGTCAGTTTTTCTTAGGCAGCACCAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCCGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCAcgcg                                                                                                                                                      <  2:360545/234‑1 (MQ=255)
                                                                                          cacGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGTTCGGGCCGTCCTGACGCTCCACACCGCATTTCCCcgcg                                                                                                                                                      >  1:1058988/1‑238 (MQ=255)
                                                                                           acGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGAGGTTTGCCAGTTGCCCTTCAGTTCGTTCCTCCTGCGCCAGGTTCTGACGGGAGGGGATCGGTGGGGTCGGGCCGTCCTGACGCCCAAAACCGtcctctttaaacca                                                                                                                                                    >  2:1054359/1‑226 (MQ=255)
                                                                                            cgcgcgCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCCGGTTCTGACGGGAGAGGATCAGTGCGGTCTGGCCGTCCTGACGCTCAACCCCGCATTTCCACGCGAc                                                                                                                                                    >  1:192421/1‑238 (MQ=255)
                                                                                             gcgcgcTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGCCCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACcg                                                                                                                                                  <  1:1066926/239‑1 (MQ=255)
                                                                                             gcgcgcTTTCACGCCTTCGGCAGTCAGTTTTCCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACcg                                                                                                                                                  <  1:1036690/239‑1 (MQ=255)
                                                                                                  cTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCACCACCGTATTTCCACGCGACCGCGGAt                                                                                                                                             >  1:376036/1‑239 (MQ=255)
                                                                                                     ttCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCTCGCGCGATGTTTGCCAGTTGCTCTTCCGTTCGTTCCTGCTGCGCCCGGTTCTGCCGGAGGAGGACCAGGCCGGCCGGGCAGTCCCGACCCCACCACCCGTATTTCCCCCCGACcccgagtcc                                                                                                                                           >  2:943294/1‑234 (MQ=255)
                                                                                                            cTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTACTGCTGCGCCAGGTTCTGACGGAGGAGAATCAGGACGGTCGGGCCGCCCTGCCGCCACACCCCGGATTTCCAcccaccgggggtgtaaacctgg                                                                                                                                    >  2:979641/1‑218 (MQ=255)
                                                                                                              tCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGTTCAGTGGGGTCGGGCCGTCCTGACGCCCAACACCGTCTTTCCCACCGACCCCGGATTCAACCTggc                                                                                                                                  >  2:755273/1‑237 (MQ=255)
                                                                                                              tCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGGGCGGTCGGGCCGTCCTGACGCCCAACACCGGATTTCCACGCGACCGCGGATTCAACCTggcc                                                                                                                                 >  2:32656/1‑237 (MQ=255)
                                                                                                                 ccagctagtTTTTCTTAGGCAGCAACAGCCCGTTCAACTTCTGAACCGGTAGCCATGAAAATCAGTTCCGGCTGACCGGCCCAGTCTTTCAGCACATAACCACCGCGCGGGATGTTTGCCAGTTGCTCTTCGGTTCGTTCCTGCTGCGCCAGGTTCTGCCGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTcaca                                                                                                                              <  2:430456/233‑1 (MQ=255)
                                                                                                                  cagtcagtTTTCCGTTGGCAGCAACCGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTcaca                                                                                                                              <  1:32656/238‑1 (MQ=255)
                                                                                                                      cgttttttGGGAGGCCCCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGGTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGAc                                                                                                                         <  2:583909/236‑1 (MQ=255)
                                                                                                                      cagtTTTCCGCAGGGACCAACAGCCATGTCAACTTCTGAACCGGTAGCGATGAAAATCAGGTCCGGCTGACCGGCTCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGAc                                                                                                                         <  1:1010903/239‑1 (MQ=255)
                                                                                                                       tttttttGGTGGGCAACAAAAGCCCGTTCAACTTCTGAACCGGTAGCGATCAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACg                                                                                                                        <  2:981254/237‑1 (MQ=255)
                                                                                                                       tttttttCGTGGGCGGCGACAGCCATTTCCACTTCTTAACCGGTAGCGATGAAATACCGTTCCGGCTGACCGGCGCAGTCTTTCAGCCCATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACg                                                                                                                        <  2:859555/237‑1 (MQ=255)
                                                                                                                       tgtTTTTCGTTGGCCGTAACCGCCAGGTCAACCTCTGAACCGGTAGGCATGAAAATCCAGTCCGGCTGACCGGCGGAGTCTTTCAGCACATACCCACCGCGCGCGATGGTTGCCAGTTGCTCTGCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACg                                                                                                                        <  2:222362/238‑1 (MQ=255)
                                                                                                                         ttttcggtgggCAGCAACCGCCATTTCAACTTCTGAACCGGTAGCGATGAAAAGCGGTTCCGGCTGACCGGCGCAGTCTTTCAGCCCCTAGCCACCGCGCGCGATGGTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGcc                                                                                                                      <  2:485952/230‑1 (MQ=255)
                                                                                                                          ttttCGTTGGCAGCCACCGCCAGTTCAACTTCTGAACCGGTAGGGATGAAAAGCAGTTCCGGCTGCCCGGCGCAGTCTTTCAGCAACTAACCCCCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGcc                                                                                                                      <  2:109008/238‑1 (MQ=255)
                                                                                                                          ttatCGGAGGCAGGAACCGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACg                                                                                                                        <  1:157583/233‑1 (MQ=255)
                                                                                                                                       caACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCTGCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATGCAACCTGTGTCACAGGACGCCATGACGCCATGtta                                                                                                        >  1:280330/1‑238 (MQ=255)
                                                                                                                                              cAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCCAAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGtt                                                                                                  <  1:309544/238‑1 (MQ=255)
                                                                                                                                                 ttCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCCTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGCCCCACGGACGCCATGTAGACATGTTCGGGGTTAcg                                                                                               >  1:521928/1‑238 (MQ=255)
                                                                                                                                                    aaCTTCCGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCCCCGCGCGCGATGGTGGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCa                                                                                           <  2:779371/239‑1 (MQ=255)
                                                                                                                                                      cTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGGCACACGGGAGCCATTGACACATGTTCGGGGTTACGCGCagg                                                                                         >  2:500476/1‑238 (MQ=255)
                                                                                                                                                           gAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCCACCTGGTCACACGGCCGCCATGTAGACACGTTCGGGGCTACGCGCAGaaaaca                                                                                    >  1:904596/1‑237 (MQ=255)
                                                                                                                                                           gAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGCAGACATGTCCGGGGTTACgcgccggaaagc                                                                                    >  2:278802/1‑231 (MQ=255)
                                                                                                                                                             aCCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCg                                                                                   <  1:88403/238‑1 (MQ=255)
                                                                                                                                                              ccGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGAc                                                                                 <  1:140501/239‑1 (MQ=255)
                                                                                                                                                                      gttGAAAATCAGGTCCGGCTGACCGGCGCAGTCTTTCAGCACATACCCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCa                                                                          <  1:340112/236‑1 (MQ=255)
                                                                                                                                                                             aTCCGTTCCGGCTGGCCGGCGCCTTCTTTCAGCACAAAACCCCCGGGCGAGATGGTTGCCAGTTGCTCTTCAGTGCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGGGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTg                                                                  <  1:452386/239‑1 (MQ=255)
                                                                                                                                                                                          gACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCACCACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGGGAAGCGACCTGCTCAACCGGCTGGGGAGGCGGGcgg                                                     >  1:1083304/1‑237 (MQ=255)
                                                                                                                                                                                                cgcAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGCCCGCGGATTCAACCTGGCCACCCGGACGCCATGTAGACATTTTCGGGGTTACGCGCAGAGAAGCGACCTGCCCAACCGCCGGGGGAGTCGGGGCGCCTTc                                                >  2:59503/1‑238 (MQ=255)
                                                                                                                                                                                                       tttCAGCCCATAACCCCCGGGCGCGATGGTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGAc                                        <  1:348859/239‑1 (MQ=255)
                                                                                                                                                                                                       tttCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCCCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGAc                                        <  1:917617/239‑1 (MQ=255)
                                                                                                                                                                                                           cgcccaTAACCCCCGGGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGAt                                    <  1:1053943/235‑1 (MQ=255)
                                                                                                                                                                                                               caTAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGCCCGATggcg                                >  1:740483/1‑237 (MQ=255)
                                                                                                                                                                                                               caTAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAg                                >  1:684922/1‑239 (MQ=255)
                                                                                                                                                                                                               caTAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTAAGCGCAGAGAAGCGACCTGCTCAACCGGCTGGGGAGTCGGGCCGGCTTCGCCCAAACCggagggg                                >  2:806085/1‑233 (MQ=255)
                                                                                                                                                                                                               caTAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGCAGACATGTTCGGGGTTACGCCCAGAGAACCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGCGGGAg                                >  1:594799/1‑239 (MQ=255)
                                                                                                                                                                                                                 tAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGGGGGTCGGGCCGTCTTCGCCCAGACCGAGGGAGTc                              >  1:593861/1‑239 (MQ=255)
                                                                                                                                                                                                                   cccccccccccgggtGTTTGCCATTTTATCGTCAGTTGGTTCCTGGTGCCCCAGGTTCTGACGGGAGAGGGTCAGGGCGGTCGGTCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTTAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGt                            <  2:738697/225‑1 (MQ=255)
                                                                                                                                                                                                                   acccccgggggggATGTTTGCCGGTTGCTCTTCAGTTCGTTCCTGCTGCGCAAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGt                            <  2:687756/227‑1 (MQ=255)
                                                                                                                                                                                                                   accaccGCGGGGGATGTTTGCAGGTTGCCCTTCAGTTCGTTCCTGCTGCGCCAGGGTCGGGCGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGGATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGt                            <  2:480708/239‑1 (MQ=255)
                                                                                                                                                                                                                    cccccGCGGGGGGTGTTTGCCCGTTGCTTTTCCGTTCGGTCCTGGCGCCCCGGGTTCCGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTg                           <  2:1034265/236‑1 (MQ=255)
                                                                                                                                                                                                                       ccGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCGGCTCAACCGGCTGGGGAGTCGGGCCGTCTTCCCCCAGACCCAGGGAGGCGggggg                        >  2:889405/1‑238 (MQ=255)
                                                                                                                                                                                                                         cccgctgttGTTTGCAGTTCCTTTTCAGTTCGTTCCTGCTGCGCCAGTTTCTGCCGGGAGAGGGTCAGTGCGGTCGGGCCGTGCTGACGGTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTa                     <  2:535148/231‑1 (MQ=255)
                                                                                                                                                                                                                          cgcgcgGTGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGCCCTGCTCAACCGGACGGTGAGTCGGGCCGGCTTCCCCCAAACCgcagggggcgggggggtc                     >  2:89445/1‑222 (MQ=255)
                                                                                                                                                                                                                          cgcgcgATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGACGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCCGCCCGCTGGAGTCTTGGGGGTa                     >  1:704904/1‑239 (MQ=255)
                                                                                                                                                                                                                          cgcgcgATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGCCGTGGGGGTa                     >  1:514664/1‑239 (MQ=255)
                                                                                                                                                                                                                          cgcgcgATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGGGAGTCGGGCCGGCTTCGCCCCGACCGAgggggaggggggggga                     >  2:706416/1‑226 (MQ=255)
                                                                                                                                                                                                                          cgcgcgATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCCGGTGAGTCGGGCCGTCTTCGCCCCGACCGATGGGGCCGggggggga                     >  2:1061452/1‑235 (MQ=255)
                                                                                                                                                                                                                            cgcgATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACCCGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGAGACCTGCGCAACCGGCTGGGGAGTCGGGCCGCCTTCCCCCAGCCCGTTgggggggggggggaaaa                   >  2:699522/1‑224 (MQ=255)
                                                                                                                                                                                                                             cggttGTTTGCCAGTTGCCCTTTAGTTCGTTCCTGCTGCGGCGGGTTCTTGCGGGAGAGGGTCAGTGCGGTCGGGCCGACCTGACGCTCAACCCCGTATTTCCACGCGCCCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAAc                  <  2:192421/235‑1 (MQ=255)
                                                                                                                                                                                                                                    ttGCCAGTTGCTCTTCAGTTCGTCCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCAcc           <  1:278802/239‑1 (MQ=255)
                                                                                                                                                                                                                                     tGCCAGTTTGTCTTCAGTTCGTTCCTGGTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCt          <  1:225611/239‑1 (MQ=255)
                                                                                                                                                                                                                                        cttttGCTCTTTAGTTCGTTCCCGTGGCGCCAGTTTCTGGCGGGAGGGGATCAGGGCGGTCGGGCCGTCCTGCCGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGTTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGAc       <  2:418131/236‑1 (MQ=255)
                                                                                                                                                                                                                                           ttGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCGGGTGAGTCGGGCCGGCTTCGCCAAGACCCATGGAGCATTGGGTGTAAAACATCACCTGAc       >  1:164802/1‑236 (MQ=255)
                                                                                                                                                                                                                                             gCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCGGGGGAGTCGGGCCGTCTTCGCCCAGACCGATGGGGCCGGGGGGGTAAAACACCCCCCGAcccggt  >  2:225611/1‑234 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                              
AGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACCGGCGCAGTCTTTCAGCACATAACCACCGCGCGCGATGTTTGCCAGTTGCTCTTCAGTTCGTTCCTGCTGCGCCAGGTTCTGACGGGAGAGGATCAGTGCGGTCGGGCCGTCCTGACGCTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGT  >  NC_000913/3079804‑3080278

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