Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,618,717 T→G A163A (GCT→GCG deoA → thymidine phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,618,7170TG98.7% 209.7 / NA 77A163A (GCT→GCGdeoAthymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (45/31);  minor base C (1/0);  total (46/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

ATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGGCTGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGT  >  NC_000913/4618476‑4618962
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                         
aTAAACATTCCACGGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCGTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGt                                                                                                                                                                                                                                                  <  1:52124/251‑1 (MQ=255)
         ccACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATCTTCCGATGATCTCTGGTCGCGGCCTCGGGCATACTGGCGGTACGCCCGACAAACTGGAATCCCTCCCTGGCTTCGGCATTTTTCCGGGAGACAACCGTTTCCCCGAAATTATTAAAGACGCCGGCGGGGCGATTTTCGGGCCGTCCACATCACCGGCTCCGGCGGATAAACGTTTCTAcgcg                                                                                                                                                                                                                                         >  2:156873/1‑251 (MQ=255)
         ccACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTAcgc                                                                                                                                                                                                                                          <  1:545736/250‑1 (MQ=255)
             gggTGGGGTCGTCGATGTGACTTCGCTGATTTTGGGGCCGATGGTCGCAGCCTGCGGCGTCTAGATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCAGCCCTTGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAATACGTCTGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACTTTTCTACTCGAcct                                                                                                                                                                                                                                    <  2:258595/250‑2 (MQ=255)
              ggTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTAAGCTCGACAAAATGGGAACCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGAACAGCTCACTGGGCCCGGCGGAAAAACGTTTCTAAGCGCCCCg                                                                                                                                                                                                                                   >  1:591890/1‑251 (MQ=255)
              ggTGGCGTCGGCGATGTGACTGCGCTGATGTTGGGGCCGATGGTCGCATCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCTCTGGCTCCGGCGGATAAACGTTTCTACGCGAcc                                                                                                                                                                                                                                     <  1:223892/249‑1 (MQ=255)
                 ggCGTCGGCGATGGGACTTCGCTGATGTTGGGGCCGATGGTTGCAGCCTGCGGCGGCTATATTCCGAGGATCTCTGGTCGCGGCCGCGGTCATACTGGCGGTACGCTCGACAACCTGGAATCCATCCCTGGCTTCGACATTTTCCCTGATGAGAACCGTTTCCGCGAAATTAGTAAAGGCGTCGGCGTGGCGATTATCGGTCAGCCCAGCTCACTGGCTCGGGCGGATAAACGTTTTTACTCGGCCCGTg                                                                                                                                                                                                                                 <  2:410833/250‑1 (MQ=255)
                  gCGTCGGCGATGTGACTTCGCTGATGTTGGGGTCGATGGTCGCCGCCTGCGGCGGCTAGATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAAGCCATCCCTGGCTTCGCCATTGTCCCGGATGACAACCGTTTCGGCGAAATTATTAAAGTCGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGTGat                                                                                                                                                                                                                              <  1:102959/251‑1 (MQ=255)
                   cGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGCGAccc                                                                                                                                                                                                                                    <  1:16951/245‑1 (MQ=255)
                    gTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACCAACTGGAATCCATCCCTGGCTTCGACATTTTTCCGGATGACAACCGTTTCCGCGAAATTTTTAAAGACGTCGGGGTGGCGATTATCGGTCCGACCAGCTCACTGGCTACGGCGGATAAACGGTTCTTCGCGGCCCggggtat                                                                                                                                                                                                                             >  1:288583/1‑245 (MQ=255)
                         cGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAAAGTTTCTACGGGACCCGTGATATTAccc                                                                                                                                                                                                                      >  1:216157/1‑250 (MQ=255)
                         cGATGTGACTTCGCTGATGTTGGGGCCGATGGGCGCAGCCTGCGGCGGCTATATTCCGAGGATCGCTTGTCGAGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTGTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGCGCTCCGTGATATTACCg                                                                                                                                                                                                                      <  2:272316/251‑1 (MQ=255)
                           aTGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCCGGCTTCCACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTCTCGGTCAGACCCGCTCACTGGATCCGGCGGGTAAAAGTTTTTACGCGAACCCGGATAtttacccc                                                                                                                                                                                                                    >  2:6125/1‑245 (MQ=255)
                                        gATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGGCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATACCTGGCTTCGCCCTTTTCCCGGATGAAAACCGTTTCCGCGAAATTATTTAAGACGTCGGCGTGGGGATTATCGGGCCGACCCGCTCACTGGCCCCCGCGGGTAAACGTTTTTACACGGCCCGTggcatttccgcacacgtgggctccc                                                                                                                                                                                                         >  2:462830/1‑227 (MQ=255)
                                           gTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATc                                                                                                                                                                                                     <  1:264980/250‑1 (MQ=255)
                                             tGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGCCAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCCTTTCCGCGAAATTATTAAAGACGTCCGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGCGACACGTGATATTTCCGCACCCGTGGACTCCATcccc                                                                                                                                                                                                  >  1:48214/1‑250 (MQ=255)
                                              ggggCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATccct                                                                                                                                                                                                  <  2:544011/250‑2 (MQ=255)
                                              ggggACGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACGGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCg                                                                                                                                                                                                  <  1:355238/250‑1 (MQ=255)
                                                         tCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACCAACTGGAATCCCTCCCTGGCTTCGACCTTTTCCCGGGTGGCCACCGTTTTCGCGAAATTATTAAAGACGGCCGCGGGGCGATTCTCGGTCAGAACCGCTCACTGGCTCCGGCGGAGAAACGTTTTTACGCGCCCCGGACTTTTACCCCCACCGCGGGCTCCATCCCGATGaaccacgcc                                                                                                                                                                                     >  2:368005/1‑243 (MQ=255)
                                                           gCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTAGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTCCCGCAACCGTGGACTCCATCCCGCTGATCACCGCctc                                                                                                                                                                                   <  1:227531/251‑1 (MQ=255)
                                                           gCAGCCTGCGGCGGCTATATTCCGATGATATCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGGCAAACTGGAATCCATCCCTCGCTTCGACATTTTCCCGGAGGACAACCGTTTCCGCGAAAATATTTAACACGTTGGCGTGGCGGTTATCGGTCAGACCAGGTCACTGGCCACGGCGGATAAACGTTTCTACGCGGCCCGGGGTAATTCCCGAACCGCGGAATCCATCCCGccggtccccccctc                                                                                                                                                                                   >  2:326361/1‑238 (MQ=255)
                                                            cAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAAACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACCCCCTCCCCCTGATCAccccccct                                                                                                                                                                                  >  1:339791/1‑245 (MQ=255)
                                                              gCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTTCGCGAAAATATTAAAGACGGCGGCGTGGGGATTATCGGTCCGACCAGCTCACTGGCTCCGGCGGGTAAACGTTTCCACGCGGCCCGGGATATTTCCGCCACCGTGGACTCCCTCCCgatgatttcaccctcttt                                                                                                                                                                                >  1:3506/1‑239 (MQ=255)
                                                                   cggcggtattatTCCGATGACTTCTGGTCGGGGCCTCGGTCATACGGGCGGTCGGCCCGAAAACCGTGAACCCACCCCGGGCTTCGACTTTTTCCCGGATGACAACCTTTTCCGCGAAATTATTAAAGCCGCCGGCTTGGCGATTATCGGTCAGACCAGCTCACTGGCTCGGGCGGATAAACGTTTCTACGCGCCCCGTGATATTCCCGCACCCGTGGACGCCATCCCGCTGATCACCGCCTCTATTCTgg                                                                                                                                                                           <  2:71654/242‑1 (MQ=255)
                                                                    ggcggcTCTATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACACACTGGAAACCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTTTACGCGCCCCGGGATATTACCGCAACCCTGGACTCCAACCCGCTGAACACCGCCTCTTTTCTGGc                                                                                                                                                                          >  2:419566/1‑251 (MQ=255)
                                                                       ggcTATATCCCGATGTTCTCTGGTCGCGGCCCCGGTCATCCGGGCGGTACGCTCGACAAACTGGAACCCATCCCTGGCTTCGACATTTTCCCGGAGGACAACCGTTTCCGCGAAATTATAAAAGACGTCGGCGTGGCGATTATCGGTCAGCCCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACCCCACCCCGCTGATCACCGCCTCTATTCTGGCgaa                                                                                                                                                                       <  1:8066/251‑1 (MQ=255)
                                                                             atTCCGATGATCTCTGGTCGCGGCCTCGGTCATCCGGGCGGTACGCTCGACAATCTGGAACCCATCCCTGGCTGCGAAATTTTCCCGGATGACACCCGTTTCCGCAAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACGGGCTCCGGCGGATAAACGTTTCTACGCGCCCCGTGATATTACCGAACCCGTGGACTCCACCCCGCGGATCACCGCCTCTATTCTGGCGAAGAAAc                                                                                                                                                                  <  1:160915/250‑1 (MQ=255)
                                                                              ttCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGCCAAACTGGAAACCATCCCTGGCTTCGACCTTTTCCCGGATGACAACCATTTACGCGAAATTATTAAAGAACTCGGCGTGTAGCATATCGGTCAGACCAGATCACAGGACCCCGCGGATACACGTTTCTACGCGAACCGGGATCTTACAGCAACCGTGGGATCCCTCCCGCTGATCAACGCCTCTATTCTGGAGAAGAAACtg                                                                                                                                                                >  1:368932/1‑250 (MQ=255)
                                                                                     gacttcgggTCGCGGCCTCGGTCATACTGGCTGTACGCCCAACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATCACAACCTTTCCCGCGAAATTATTAAAGCCTTCGGCGTGGCTATTACCGGTCAGCCCAGCTCACTGGCTCCGGCGGAAAACCGTTTCTAGGCGACCCGTGATATTACCGCACCCGGGGACTCCACCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAg                                                                                                                                                         <  1:233998/244‑1 (MQ=255)
                                                                                      atCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACCAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAAACGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTAAGACCAACTCACTGGCTCCGGCGGAAAAACGGTTCTACGCGACACGCGATATTAACGCAAACGTGGACTCCATCCCGCTGATCACCGCCTCTATTACGGCGCAGAAACTGGCGGCAgg                                                                                                                                                        >  1:123000/1‑251 (MQ=255)
                                                                                             gccgcgGCCTCGGTCATACGGGCGGTCCGCTCGAAAAACGGGAATCCACCCCTGGCTTCGACATTTCCCCGGATGACACCCGTTTCCGCGAAATTATAAAAGCCGTCGGCTTGGCGATTACCGGTCAGACCAGCTCACTGGCTCCGGCGAATAAACTTTTCTACGCGACCCGTGATATTCCCGCACCCGTGGACTCCACCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGAc                                                                                                                                                 <  1:286996/249‑1 (MQ=255)
                                                                                                gcgGCCTCGGTCATACTGGCGGTACGCTCCCCAAACTGGAGTCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGCCAGACCAGCTCACCGGCCCCGGCGGGTAAACGTTTCTACGCGACCCGCGCTCTTACCGCACCCGTGGACTCCATCCCGACGATCCCCGCCTCTATTCTGGCGGAGAAACtggacggggggttggagacg                                                                                                                                              >  2:589056/1‑232 (MQ=255)
                                                                                                            atccTGGCGGTACGCTCAAAAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTCCCGCGAAATTATTAAAGACGGCGGGGTGGCGATTACCGGTCAGACCAGCTCACGGGCTCCGGCGGATAAACGTTTCTCCGCGCCCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGAc                                                                                                                                  <  1:214960/248‑1 (MQ=255)
                                                                                                                    ggTACGCCCGACAAACTGGAAGCCATCCCTGGCTTCGATATTTCCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGGGGGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAAGGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGt                                                                                                                          <  1:434005/251‑1 (MQ=255)
                                                                                                                      aaCGCTGGACAAACTGGAATCCATCCGGGGCTTCGAAATTTTCCCGGAGTACAACCGTTTCCGCGAAATGATTAAAGCCGTGGGCGTGGCGATAATCGGGCAGACCAGCTCACGGGCTCCGGCGGATAAACGTTTCGACGCGCCCCGGTATATTCCGGCAACCGGGGACTCCTCCCCGGTGATCCCCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTgg                                                                                                                        <  1:589056/250‑1 (MQ=255)
                                                                                                                        cGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGCCTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGCCGCGCTGGTGATGGACGTGAAAGTgggg                                                                                                                      >  2:190327/1‑250 (MQ=255)
                                                                                                                         gCTCGACACCCTGGACTCCATCCCTGGATTCGCCATTTTCCCGGATGACAACCGTTTCCGCGAACCTATTAACGACGACGCCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCTAGCACTCCATCCCGCTGCTCCCCGCCTCTATTCTGGCTAAGAAACTGGCGGAAGGTCTGGACGCGCAGGTGATGGACGGGAAAGAGGGtc                                                                                                                     >  2:359188/1‑250 (MQ=255)
                                                                                                                            cGACAACCTGGAACCCACCCCTGGCTCCGTCTTTTCCCCGGATGACACCCGTTTCCGCAAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAACCGTTTCTCCGCGCCCCGTGATATTCCCCCACCCGTGGACTCCTCCCCGCTGATCTCCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAg                                                                                                                    <  1:62263/249‑1 (MQ=255)
                                                                                                                            cGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAgcg                                                                                                                  >  1:135312/1‑251 (MQ=255)
                                                                                                                              aCAAAATGGATTCCACCCCTGGCTTCGCTTTTTTCCCGAATGACAACCGTTCCCGCAAAATAATAAAAGACGTCGGCGTGGGGATTATCGGTCAGACCAGCCCTCGGGCCCCGGCGGATAAACGTTTCTACGCGCCCGGTGATATTACCGCAACCGTGGCTTCCATCCCGCGGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAgcggc                                                                                                                <  1:279164/251‑1 (MQ=255)
                                                                                                                                  aCTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGAATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGt                                                                                                            >  1:19679/1‑251 (MQ=255)
                                                                                                                                      gaaaTCCATCCCTGCCTTCTAAATTTCCCCGGATGACAACCGTTCCCCCAAAATTATTAAACCCGCCGCCGGGGCGATTTCGGGTCAGCCCACCTCCTGGGCTCCGGCGGATAACCGTTTCTACGGGTCCCGGGATTTTACCTCACCGGGGGCCTCCAACCCCCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTa                                                                                                         <  1:204550/249‑1 (MQ=255)
                                                                                                                                         accccactctgtttccGACATTTCCCCGGATGACACCCTTTGCCGCAAAATTATTAAAGCCGGCGGCGTGGAGTTTATGGGTCAGACCAGCCCCCGGGCCCCGGCGGATAACGGTTTCTCCGCGCCCTGTGATATTCCCCCAACCGTGGACTCCAGCCCGTGGACCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGc                                                                                                      <  2:556260/236‑1 (MQ=255)
                                                                                                                                          gccatccCTGTCTGCGACATTTTCCCGTATGACATCCGTTTCCGCGAAATTATTAAAGACGCCGGCTTGGCGATTAGCGTTCAGCCCAGCTCACTGTCCCGGGCGGATAAACTTTTCTACGCGACCCGGGATTTTTCCGCAACCGGGGACCCCACCCGTGGGACCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGcc                                                                                                     <  1:61376/250‑1 (MQ=255)
                                                                                                                                            catccTTGGCTACGACATTTTCCCGGATCCCAACCGTTTCCGCGCAATTATTAAAGCCGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCgc                                                                                                   >  2:618924/1‑250 (MQ=255)
                                                                                                                                                ccTGGCTTCGACATTTTGCCGGATGACAAACGTTTCCGCGAAATTATTCAAGACGTCGGAGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTAAGCGACCCGTGATATTACCGCAACCGTGGACCCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGGGATGGACGCGGAAGAGGGTAGCGGCGCGGTTATGCCGCCCt                                                                                                >  2:307707/1‑250 (MQ=255)
                                                                                                                                                  tGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATAACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGGCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACg                                                                                             >  1:513870/1‑251 (MQ=255)
                                                                                                                                                   ggCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGAGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGAAAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTAc                                                                                              <  1:179503/249‑1 (MQ=255)
                                                                                                                                                        cGACTTTTTCCCGGATAAAAACCTTTGCCGCGAAATTATAAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACGGGCTCCGGCGGATAAACGTTTCTACGCGACCCGTGATTTTCCCGCACCCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAActct                                                                                       <  1:340089/251‑1 (MQ=255)
                                                                                                                                                        cGACATTTTCCAGGCTGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTCTCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCCACCGTGGACTCCTTCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGCCTGGACGCGCTGGTGATGGACGTGGAAGTGGGTAGCGGCGCCTTTATGCCGCCCTACGAAccct                                                                                       >  2:557999/1‑248 (MQ=255)
                                                                                                                                                                  acGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAAACGTGGAATCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGGCTGGACGCGCTGGTGATGGACGGGAAAAGGGGTAGCGGCGCGGTTCGTCCGGCCTAGGGACTCTCTGAGGCCCt                                                                             >  2:162722/2‑251 (MQ=255)
                                                                                                                                                                      aTGACAACCGTTTCAGCGAACTTATTAAAGACGTCGACGTGGCGCTTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAACCGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGGGATGGACGTGAAAGTGGGTAGGGGCGCGTTTCAGCCGACCTACGAACTCTCTGAAGCCCTTGcc                                                                         >  2:244263/1‑251 (MQ=255)
                                                                                                                                                                       tGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGCAAGACCAGCTCACTGGCTCCGGCCGAAACACGTTTCTACGCGACCCGGGATATTACCGCAAACGTGGACTCCATCCCGCCGATCCCCGCCTCAATTCTGGCGAAGAAAATGGCGGCAGGACTGGACGCGATGGTGATGGACGTGACAGAGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCCTGCCg                                                                        >  1:372847/1‑251 (MQ=255)
                                                                                                                                                                           aaCCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGACACTGGCGGAAGGTCTGGACGCGCTGGGGATGGACGTGAACGTGGGTAGCGGCGCGGTTATGCCGACCTACGAACTCTCTGCAGCCCTTGCCGAAGc                                                                    >  2:459372/1‑251 (MQ=255)
                                                                                                                                                                           aaCCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTATGGCCGCGCTGGTGATGGCCGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGc                                                                    >  1:627215/1‑251 (MQ=255)
                                                                                                                                                                                  tCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGACGGTCTGGGAGCGCTGGTGCTGGACGGGGAAGTGGGTAGGGGCGCGTTTATGCCGGCCTACGAACTCTCTTAAGCCCTTGCCGAAGCGgttgtt                                                             >  2:332425/1‑251 (MQ=255)
                                                                                                                                                                                       aaaaTTATAAAAGCCGGCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGCCTCCGGCGGAAAAACGTTCCTACGCGCCCCGTGATATTCCCGCACCCGTGGACCCCACCCCGCGGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGt                                                        <  1:360611/250‑1 (MQ=255)
                                                                                                                                                                                        aaaTTATTAAAGGCGCCGGCGGGGCGATTACCGGCCAGCCCGGCTCACTGGCTCCGGCGGATAAACGTTTCTCCGCGTCCCGTGATATTCCCGCAACCGTGGCCTCCATCCGGCTGATCACCGCCTCTAGTCGGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTg                                                       <  1:414425/251‑1 (MQ=255)
                                                                                                                                                                                            tattaAAACCGTCGGCGTGGCGATTGTCGGTAAGCCCGGCCCGCTGGCCCCGGCGGATAAACTTTTCTCCTCGACCCGTTATATCCCCGCACCCGTGGACCCCCCCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCt                                                    <  1:307707/250‑1 (MQ=255)
                                                                                                                                                                                              ttaAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGGTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAAc                                                 >  1:430499/1‑251 (MQ=255)
                                                                                                                                                                                                 aaGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCCGGGCTAACGgc                                              >  1:434436/1‑251 (MQ=255)
                                                                                                                                                                                                 aaGACGTCGGCGTGAAGCATATCGGTCAGACCAGCTAACTGGCTCCGGGGGATAAACGTTTCTACGCGACCCGGGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGGGATGGACGGGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCGGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACggg                                              >  2:391794/1‑250 (MQ=255)
                                                                                                                                                                                                       tCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGCCCTACGGACTCTCTGAAGCCCTTGCCGGAGCGATTGTTGGCGTGGCCAAAGGCGCTGGc                                        >  2:174130/1‑251 (MQ=255)
                                                                                                                                                                                                        cGGCGTGGCGATTATCGGTCCGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGGACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCg                                       >  2:109441/1‑251 (MQ=255)
                                                                                                                                                                                                           cGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTACCGGCGCTGGCGTgc                                    >  2:593540/1‑251 (MQ=255)
                                                                                                                                                                                                           cGTGGCGATTATCGGTCAGACCACCTCACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGCAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGCAATCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGAGCTGGCGTgg                                    >  1:153408/1‑250 (MQ=255)
                                                                                                                                                                                                            gTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGCAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTgcg                                   >  1:492749/1‑251 (MQ=255)
                                                                                                                                                                                                            gTGGCGATTATCGGTAAGACCAGCTAACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCAGCTGATCACCGCATCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTgcg                                   >  2:368600/1‑251 (MQ=255)
                                                                                                                                                                                                                       cGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAAAGGCGCTGGCGTGCGCACCACCgcgc                        >  2:509016/1‑251 (MQ=255)
                                                                                                                                                                                                                         gTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGATGGCGTGCGCACCACCGCgctg                      >  2:347686/1‑251 (MQ=255)
                                                                                                                                                                                                                             aGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGGGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGACGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGGTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCCGCTc                   >  2:356066/1‑251 (MQ=255)
                                                                                                                                                                                                                               acccagctcacggccTCCGGCGGATAAACTTGTCTACGCGACCCGTGATATTCCCACACCCGTGGACTCCATCCCGCTAATCGCCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCa                  <  1:2032/237‑1 (MQ=255)
                                                                                                                                                                                                                                ccAGCTCACTGGCTCCGGCGGATAAACGTTTCTACGCGCCCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGGGGGTAGCGGCGCGTTTATGCCGCCCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCAcc                >  2:548592/1‑251 (MQ=255)
                                                                                                                                                                                                                                     tCACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGCTCCCCGACTCTATTCTGGCGAAGACACTGGCGGACGGTCAGGACGCGCTGGTGTTGGACGCGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGCACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCAACACCGCGGTGCTCAccaccat           >  1:320680/1‑246 (MQ=255)
                                                                                                                                                                                                                                     tCACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACAt           >  1:96556/1‑251 (MQ=255)
                                                                                                                                                                                                                                     tCACTGGCTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACAt           >  1:335731/1‑251 (MQ=255)
                                                                                                                                                                                                                                          ggCTCCGGCGGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATc      >  1:435560/1‑251 (MQ=255)
                                                                                                                                                                                                                                                cggcggATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCCCCGACATGAATCAGGt  >  1:95857/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                         
ATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGGCTGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGT  >  NC_000913/4618476‑4618962

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