Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,204,890 C→T G348S (GGT→AGT)  purD ← phosphoribosylamine‑‑glycine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,204,8900CT97.5% 277.6 / NA 83G348S (GGT→AGT) purDphosphoribosylamine‑‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  major base T (52/27);  minor base G (0/2);  total (53/30)
Fisher's exact test for biased strand distribution p-value = 1.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTTCTGCCACGGTATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAAATTCAATAACCTTCGGATTGCCCTGTTTGTCGATCATCAGGCCCGCGTAGAGAAAACCGGTGT  >  NC_000913/4204650‑4205135
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
gCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTTCTGCCACGGTATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTa                                                                                                                                                                                                                                              >  2:546064/1‑250 (MQ=255)
    gcgTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTTCTGCCACGGTATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCCCGTCCCTGGTGCGGTAATCAc                                                                                                                                                                                                                                         >  2:265743/1‑251 (MQ=255)
     gtttgaacgccgcgcCAGCCATTATCTTTCCGCCAGAAGGAGCCGTCCCAGGGGATATCTGTCAGTAAGCCAGGGGCGTGTTCCTGCGCTTCTCCCGCGGAAGGACCCAGCGCGGTTCCGCACAGGCCGCCCCCGCCGTTGGTCACTGCCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTTCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCAcc                                                                                                                                                                                                                                        <  1:78868/241‑1 (MQ=255)
     cgTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCGGTCACTAAGGCATAGGCGCGTTTCTGCGCTTCTGCCACGGAATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTAACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCCGGCCGTGGATCACGTCACTGGTGCGGTAATCAcc                                                                                                                                                                                                                                        >  1:587517/1‑251 (MQ=255)
        gcgagaGCGCCCCAGCCGATATCTTTCCGGCAGAAGCAGTCGCCCCAGTGGATATCGGTCATAAAGGCATAGGCGCGTTTCTGCGCTTCTGCCACGTTATGCCCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCg                                                                                                                                                                                                                                      <  1:202820/245‑1 (MQ=255)
         cGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTTCTGCCACGGTATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCCGGCCGTGGATCACGTCACTGGGGCGGTAATCACCCgg                                                                                                                                                                                                                                     >  2:218481/1‑250 (MQ=255)
         cGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTTCTGCCACGGTATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCgg                                                                                                                                                                                                                                     >  1:14635/1‑250 (MQ=255)
          gATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTTCTGCCACGGTATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCCCTGGTGCGGTAATCACCCGGat                                                                                                                                                                                                                                   >  1:349139/1‑251 (MQ=255)
          gATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTTCTGCCACGGTATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGat                                                                                                                                                                                                                                   >  1:308576/1‑251 (MQ=255)
          aagagcggcccAGCCTATATCTTTCCGGAAGAAGCGTTCGCCCCGGGGGATAGCGTTAATAAAGGCATAGGGGCTTTTCGGCGCTTCTGCCGCGGTAGGCCCCAGCGGGGTGACGCACAGGCGGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCTCCCGGat                                                                                                                                                                                                                                   <  1:56701/242‑1 (MQ=255)
                 cgcCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCACTGGCTATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTTCTGCCACGGTATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCAcc                                                                                                                                                                                                                            >  1:513508/1‑251 (MQ=255)
                   ccAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTTCTGCCACGGTATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCg                                                                                                                                                                                                                          >  1:568221/1‑251 (MQ=255)
                    cAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTTCTGCCACGGTATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGGAATCACCCGGATATCCACCCGc                                                                                                                                                                                                                         >  1:2045/1‑251 (MQ=255)
                            atCTTTCCGGCAGCAGCAGTCGTCCCAGTGGATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTTCTGCCACGGTATGACCCAGCGCGGTGCCGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGCTATCCACCCGcctccatcc                                                                                                                                                                                                                 >  2:5233/1‑243 (MQ=255)
                            atCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTTCTGCCACGGTATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCCCCCGGATATCCACCCGCAGCCATcc                                                                                                                                                                                                                 >  1:555254/1‑250 (MQ=255)
                                          aGCAGTCGTCCCAGTGAATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTTCTGCCACGGTCTGACCCAGAGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCgagaga                                                                                                                                                                                                    >  1:366320/1‑250 (MQ=255)
                                            cagcctccccagggaatagcgttCATTAAGGCATGGGCGCTTTCCGGCCCTTCGCCCACGGTATGCCCCGGGGGGGTGACCACAAGACCGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAgc                                                                                                                                                                                                 <  1:245408/229‑1 (MQ=255)
                                                            tatCGGTCCTTACTGCATACGCGCGTTTCTGCGCTACTGCCACGGTATGACCAAGCGCGGTGAAGCAAAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCAGCGAGGAACACTTTGCCCCCTACCACTTCTTCCAGCGGCAGGCCGTGGATCCCGTCACTGGTGCGGTCATCACCCGGACATCCACCCGCAGCCATCACCCCGCCGAGAGACGCGACTTCATccacctcc                                                                                                                                                                                 >  2:587657/1‑245 (MQ=255)
                                                            tatCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTTCTGCCACGGTATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCg                                                                                                                                                                                 >  1:527427/1‑251 (MQ=255)
                                                                   cATTAAGGCATCGGCGCGTTTCTGCGCTTCTGCCACGGTATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCCCTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCCCTCGGACGttt                                                                                                                                                                          >  2:84900/1‑251 (MQ=255)
                                                                    aTTAAGGCATAGGCGCGTTTCTGCGCTTCTGCCACGGTATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCCCTGGTGCGGTAATCAACCGGATATCCCCCCGCAGCCATCACCCCGCCGAGAGAAGCCCGTTCATACCCCTCCGACGtttt                                                                                                                                                                         >  2:12479/1‑251 (MQ=255)
                                                                      tAAGGCATAGGCGCGTTTCTGCGCTTCTGCCACGGTATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTtct                                                                                                                                                                       >  2:264050/1‑251 (MQ=255)
                                                                       aaGGCATAGGCGCGTTTCTGCGCTTCTGCCACGGTATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTtctc                                                                                                                                                                      >  1:325217/1‑251 (MQ=255)
                                                                       aaGGCATAGGCGCGTTTCTGCGCTTCTGCCACGGTATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTtctc                                                                                                                                                                      >  1:356573/1‑251 (MQ=255)
                                                                            aTAGGCGCGTTTCTGCGCTTCTGCCACGGTATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTcc                                                                                                                                                                  >  1:463161/1‑250 (MQ=255)
                                                                              aGGCGCGTTTCTGCGCTTATGCCACGGTATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGCGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGCAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGt                                                                                                                                                               >  2:551753/1‑251 (MQ=255)
                                                                                  gcgTTTCTGCGCTTCTGCCACGGTATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGCACGTTTTCTCGTCCAGTTTGc                                                                                                                                                           >  1:327593/1‑251 (MQ=255)
                                                                                        cTGCGCTTCTGCCACGGTATGACCCAGCGCGGTGACGCCCAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGCACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTcac                                                                                                                                                     >  1:563356/1‑251 (MQ=255)
                                                                                            gcTTCTGCCACGGAATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGgc                                                                                                                                                 >  1:490940/1‑251 (MQ=255)
                                                                                             cTTCTGCCACGGTATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTTCTCGTCATCCGCCAGTTTCGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGAGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGGAAGAAAAGCATTACTCCACCTCGGCCGGTTTCTCATCCAGTTTGCTGTCACAgcac                                                                                                                                                >  2:345276/1‑248 (MQ=255)
                                                                                                cTGCCACGGTATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGgccgcc                                                                                                                                             >  1:556303/1‑251 (MQ=255)
                                                                                                         tATGACCCAGCGCCGCGACGCCCAATACGCGCCCGCCGTTGGTAAATACATGCTCCTCATCGGCAAGCTTTACTCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCCGCGGCCGTCCGTGGATCCCTTCACTTGGGCGGTAATCAACCGGGTATCCACCCGGACCCATCCCCCCGCCGAGAGAAGCGCCTTCATCACCCTACGACGTTTTCTCGTACAACTTGCTTTCACAGGGCGCCAGGCagcac                                                                                                                                    >  2:313033/1‑248 (MQ=255)
                                                                                                          aTGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCagcgcc                                                                                                                                   >  2:69545/1‑247 (MQ=255)
                                                                                                          aTGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCt                                                                                                                                   >  1:417917/1‑251 (MQ=255)
                                                                                                             cacccggggtggggccccactgcaCGCCCCCCCCGTGGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTc                                                                                                                                  <  2:109305/227‑1 (MQ=255)
                                                                                                                  ggcggggggacaccaggtccgcgCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAg                                                                                                                            <  1:182703/232‑1 (MQ=255)
                                                                                                                        tGACGCACAGTACCCGCCCGCCGTTAGTCACTACCTGATCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACAATTTGCCGCCTCCCACTTCTTCCAGCGGCAGGCAGTGGATCACGTCACCGGAGCGGTAATAACCCGGATATACACCCGCAGCCATAACACCGCCGAGAGAAGCGCGTTCATGCCACTGGGACGTCTTCTCGTCCAGCTTGCTTTCACAGGCCGCCAGGCAGAGCTAAACCAGATCGGAc                                                                                                                     >  1:310121/1‑251 (MQ=255)
                                                                                                                          cccccacagcccgcccccgccGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTGTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACt                                                                                                                    <  1:265743/236‑1 (MQ=255)
                                                                                                                               caGTCCGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATgc                                                                                                              <  1:163977/251‑1 (MQ=255)
                                                                                                                               caGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCCTCACCACGCCGAGAGAAGCGCGTTCATCCCCCTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCCGGCAGAGCTCAACCAGATCGGACTTCATgc                                                                                                              >  2:386389/1‑251 (MQ=255)
                                                                                                                                  gtccgcccccgccGTTGGTCACTCCCTGCCCGTCATCCGCCAGTTTGGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCGCCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATgcgc                                                                                                            <  1:419969/244‑1 (MQ=255)
                                                                                                                                       cccccgccGTTGGTCCCTCCCTGCTCGTCTCCCGCCAGTTTTGTGCCCGCGTGGAACCCTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCAt                                                                                                       <  2:14635/249‑1 (MQ=255)
                                                                                                                                        gcccgccGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCAt                                                                                                       >  1:107536/1‑250 (MQ=255)
                                                                                                                                          ccgccgTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCATGCAGAGCTCAACCAGATCGGACTGCATGCGCAGCATAAt                                                                                                    <  1:92487/251‑1 (MQ=255)
                                                                                                                                               gTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCt                                                                                               >  2:427928/1‑251 (MQ=255)
                                                                                                                                               gTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCt                                                                                               >  1:490417/1‑251 (MQ=255)
                                                                                                                                               gTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCt                                                                                               >  1:198594/1‑251 (MQ=255)
                                                                                                                                                ttGGTCACTACCTGCTCGAAATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCCCCCGCAGCCATCACCACGCCGAAAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGACAGGAAGAGCTCAACCAGCTCGGACTTCATGCGCAGCATCATCGGCTg                                                                                              >  2:555228/1‑251 (MQ=255)
                                                                                                                                                      aCTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTc                                                                                        >  1:433608/1‑251 (MQ=255)
                                                                                                                                                         cccTGCTCGTCCTCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCgg                                                                                     <  2:98753/250‑1 (MQ=255)
                                                                                                                                                         cccTGCCCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACGGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAAGCGGCTGGGTTTcc                                                                                       <  2:312049/248‑1 (MQ=255)
                                                                                                                                                                    aTCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGc                                                                          >  2:40004/1‑251 (MQ=255)
                                                                                                                                                                    aTCCGCCAGTTTTGTCCCCGCGTGGAACACTTTGCCGCCTGCCCCTTCTTCCAGCGGCAGGCCGTGGATCACGACACTGGCGCGGTAATAACCCGGATATCCACCCGCAGCCATCACCACCCCGCGAGAGGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGGTTGCTTTCACACGCCGCCAGGCAGAGCTCAAACAGCTCGGACTTCATGCGCAGCATAATCGACTGTGTTTCCTGATCGcccaa                                                                            >  2:297601/1‑246 (MQ=255)
                                                                                                                                                                     tCCGCCATTTTTGTGCCCGCGCGGATCCCTTTGCCGCCTGACCCTGCTTCCAGCGTCAGGCCGTGGGTCACGGCACTGGTGCGGTAATCGCCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTGTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAATCAGATCGGTCTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGc                                                                          <  1:609630/250‑1 (MQ=255)
                                                                                                                                                                        gCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGCTCACGTCACTGGGGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACGCGGACGTTTTCTCGTCCAGTTTGCTTTCTCAGGCCGCCGGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCGAAGCGGcc                                                                      >  2:238394/1‑250 (MQ=255)
                                                                                                                                                                          ccagtttGTGCCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAg                                                                     <  2:2045/247‑1 (MQ=255)
                                                                                                                                                                            gTTTGGTGCCCGCGGGGAACCCTTCGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGGGGATCACTTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACGTCATGCGCAGCATAATCGGCGGGGTTTCCGGATCGCCAAAGCGGCAGTTa                                                                  <  2:240775/251‑1 (MQ=255)
                                                                                                                                                                               ttGTGCCCGCGTGGAACCCTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAAAt                                                               >  2:281932/1‑251 (MQ=255)
                                                                                                                                                                                gagtgcccgccgggcaccccttctcgccTTCCACTTCTTACCGCGGCAGGCCGTGGATCACGGCCCTGGTGCGGTAATCACCCTGATATCCACCCGCAGCCATCACCACGCCAAGAGAAGCGCGTTCATCCCTCTCCGACGTCTTCTCGTCCAGTTCGCTTTCACAGGCCGCCAGGCGTAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTaaac                                                               <  2:321859/227‑2 (MQ=255)
                                                                                                                                                                                            ggatcccTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGCCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAAATTCAATAACCTTc                                                   <  2:194034/245‑1 (MQ=255)
                                                                                                                                                                                                 acTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAAATTCAATAACCTTCGGAtt                                              <  1:427928/251‑1 (MQ=255)
                                                                                                                                                                                                  cTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAAATTCAATAACCTTCGGATTg                                             >  1:594715/1‑251 (MQ=255)
                                                                                                                                                                                                     tGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCGCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAAATTCAATAACCTGCGGATTGcc                                           <  1:118533/250‑1 (MQ=255)
                                                                                                                                                                                                      gccgccTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCCCTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAAATTCAATAACCTTCGGATTGCCCt                                         >  2:494054/1‑251 (MQ=255)
                                                                                                                                                                                                        cgccTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAAATTCAATAACCTTCGGATTGCCCTGt                                       <  2:308576/251‑1 (MQ=255)
                                                                                                                                                                                                          ccTGCCACGTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCCCCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATTCGCAGCATAATCGGCTGGGTTGCCGGATCGCCAAAGCGGCAGTTAAATTCAATAACCTTCGGATTGCCCTGttt                                     <  2:513508/251‑1 (MQ=255)
                                                                                                                                                                                                                   tcttcCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAAATTCAATAACCTTCGGATTGCCCTGTTTGTCGatca                             <  1:142755/250‑1 (MQ=255)
                                                                                                                                                                                                                    attaacaccggaGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAAATTCAATAACCTTCGGATTGCCCTGTGTGTCGatcatc                           <  2:30442/240‑1 (MQ=255)
                                                                                                                                                                                                                      tcctgCGGCCGGCAGTGGATCACGTCACGGGTGCGGATATCGCCCGCATCACCACCCGCATCCATCACCACGCCGAGAGAAGCGCGTTCAACCCACTCCGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAAACGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAAATTCAATAACCGTCGGATGTCACTGTTTGTctgtcat                            <  2:366320/244‑7 (MQ=255)
                                                                                                                                                                                                                       ccAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAAATTCAATAACCTTCGGATTGCCCTGTTTGTCGATCATCAgg                        >  2:174784/1‑251 (MQ=255)
                                                                                                                                                                                                                       ccAGCGGCAGGCCGTGGACCCCGATACTGGCGCGGAAATCTCCCGGATATCCACCCGCCGCCCTCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGCGCTCAACCAGATCGGACTTCATGCGCAGCATCATCGGCTGGGTTTCCGGATCGCCCAAGCGGCAGTTAAATTCCATAACCCTCGGATTGCCCTGTTTGTCGattctccgg                        >  2:76484/1‑244 (MQ=255)
                                                                                                                                                                                                                         aGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCAGCCCACTCGGACGTTTTCTCGTCCAGTTTGCTGTCACAGGCCGCCAGGCAGAGCGCATCCAGATCGGACTGCAGGCGCAGCAGAAGCGGCGGGGTTTCCGGATc                                                                                  <  2:86046/191‑1 (MQ=255)
                                                                                                                                                                                                                         aGCGGCAGGCCGTGGATAACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGGTCGCCAAAGCGGCAGTTAAATTCCATAAACTTCGGATTGCCCTTTTTGTCGATCATCAGGcc                      >  2:372837/1‑251 (MQ=255)
                                                                                                                                                                                                                              cAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATCTCCACCCGCAGCCATCACCACGCCGCGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACCGGCCGCCAGGCAGCGCTCAACCAGATCGGACTTCATGCGCAGCATCATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAACTTCAATAACCTTCGGATTGCCCTGTTTGTCGATCCACAGGCCCGCGt                 >  2:35554/1‑251 (MQ=255)
                                                                                                                                                                                                                              cAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAAATTCAATAACCTTCGGATTGCCCTGTTTGTCGATCATCAGGCCCGCGt                 <  1:546064/251‑1 (MQ=255)
                                                                                                                                                                                                                                ggCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAAATTCAATAACCTTCGGATTGCCCTGTTTGTCGATCATCAGGCCCGCGt                 >  1:505153/1‑249 (MQ=255)
                                                                                                                                                                                                                                ggCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAAATTCAATAACCTTCGGATTGCCCTGTTTGTCGATCATCAGGCCCGCGTag               <  1:251766/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAAATTCAATAACCTTCGGATTGCCCTGTTTGTCGATCATCAGGCCCGCGTaga              <  2:278341/251‑1 (MQ=255)
                                                                                                                                                                                                                                    gTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGCGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAAAtt                                                              >  1:50273/1‑200 (MQ=255)
                                                                                                                                                                                                                                    gTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAAAtt                                                              <  2:50273/200‑1 (MQ=255)
                                                                                                                                                                                                                                        aTCACGTCACTGGTGCGGTAATCACCAGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGAGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGCGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAAATTCAATAACCTTCGGATTGCCCTGTTTGTCGATCATCAGGCCCGCGTAGAGAAAAcc       >  1:231713/1‑251 (MQ=255)
                                                                                                                                                                                                                                            cGTCACTGGTGCGGTAATCACACGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAAATTCAATAAACTTCGGATTGCCCTGTTTGTCGGTTATCAGGCCCGAGTAGAGAAAACCgggg   >  2:479375/1‑249 (MQ=255)
                                                                                                                                                                                                                                             gTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAAATTCAATAACCTTCGGATTGCCCTGTTTGTCGATCATCAGGCCCGCGTAGAGAAAACCGgtgt  >  1:552436/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
GCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTTCTGCCACGGTATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAAATTCAATAACCTTCGGATTGCCCTGTTTGTCGATCATCAGGCCCGCGTAGAGAAAACCGGTGT  >  NC_000913/4204650‑4205135

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