Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,113,851 G→A H208Y (CAT→TAT)  fpr ← flavodoxin/ferredoxin‑NADP(+) reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,113,8510GA93.0% 244.2 / 6.4 86H208Y (CAT→TAT) fprflavodoxin/ferredoxin‑NADP(+) reductase
Reads supporting (aligned to +/- strand):  ref base G (3/3);  new base A (44/36);  total (47/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.09e-01

ACCTCCTTGTCTGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATGGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGACCAGATTTTTGAAGCGATCTAAATCTTTACCTAGTT  >  NC_000913/4113606‑4114096
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
aCCTCCTGGACTGTGGTCACCATCACTCATTTATCAACAGTTTGATGCCGGCATATCGCGGGATGATCAATCAGGTGAGTGCGGATTTGCTCATCGTGACGTTTATCGATACGCAACCGCTTCTCAGTAATGCCCCGCTGTCATATGGCCCGGTCGGCGACGTAAATATTTCGTCATCTGCAGGGTCTCTTTCCGCCACTGTTGTGCAGCGCGCACCATCTGTGGAATGCCGCACAGACTCACctggctgg                                                                                                                                                                                                                                                  >  2:409407/1‑251 (MQ=255)
  ctcctTGTCTGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGtt                                                                                                                                                                                                                                                >  2:192923/1‑251 (MQ=255)
    cctTGTCTGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTc                                                                                                                                                                                                                                              <  1:457345/251‑1 (MQ=255)
        gTCTGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGGGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTTGTTTCTTTa                                                                                                                                                                                                                                          <  1:170149/251‑1 (MQ=255)
            gtgGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTc                                                                                                                                                                                                                                       >  1:17691/1‑250 (MQ=255)
            gtgGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGAGCACCATCTGTGGATTGCCGCACAGCATAACATAGCTGGTTTCTTTATTc                                                                                                                                                                                                                                       >  1:486739/1‑250 (MQ=255)
               gTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCg                                                                                                                                                                                                                                   <  2:480643/251‑1 (MQ=255)
                   tcatcaCTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGccg                                                                                                                                                                                                                               >  2:379111/1‑249 (MQ=255)
                   ccatcaCTCATTTATCAACCGTGTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGGTTATCGATCAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCGTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTGTATTCATCGGCAg                                                                                                                                                                                                                               <  2:479840/250‑1 (MQ=255)
                     atcaCTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCGCAACGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGTCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGc                                                                                                                                                                                                                             <  2:861128/251‑1 (MQ=255)
                       caCTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCGTCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATGGCTGGTTTCTTTATTCATCGGCAGGcc                                                                                                                                                                                                                            >  1:424386/1‑250 (MQ=255)
                                 tCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCttt                                                                                                                                                                                                                 >  1:246678/1‑251 (MQ=255)
                                       gTTTGATGACGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACCTAGATGGTTTCTTTATTCTTCGGCAGGCCAATCGTGATTTCCAGtt                                                                                                                                                                                                           >  1:664762/1‑251 (MQ=255)
                                                cGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACttt                                                                                                                                                                                                  <  1:737086/251‑1 (MQ=255)
                                                 ggCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTc                                                                                                                                                                                                 >  1:334983/1‑251 (MQ=255)
                                                   cATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCa                                                                                                                                                                                                >  1:503506/1‑250 (MQ=255)
                                                      gtCGCGGGATGATCAATAAGATGAGTGCGTATTTGCTCATCGTTCCGTTTATCGATAAGTAGCCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGTGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATGGCTGGTTTCTTTATTCATCGGCATGCCAATCGTGCTTTCCAGTTCCCCTCTTTCAATTa                                                                                                                                                                                            <  2:399256/250‑1 (MQ=255)
                                                      atCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAAtt                                                                                                                                                                                             >  1:114616/1‑250 (MQ=255)
                                                        cgcgGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTAGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCGGTGGATTGCCGCACAGCATCACATAGCTGGGTTCTTTATTCATCGGCAGGCCCATCGTGCTTGCCAGTTCCCCACTTTCAATTAAt                                                                                                                                                                                          >  2:489227/1‑251 (MQ=255)
                                                            ggATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCg                                                                                                                                                                                      <  1:17295/251‑1 (MQ=255)
                                                            ggATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTAATCGGCAGGCAAACCTGGCTTTCCAGTTCCCCACTTTCAATTAATGCCg                                                                                                                                                                                      >  1:315612/1‑251 (MQ=255)
                                                                gatCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTAt                                                                                                                                                                                  >  1:434892/1‑251 (MQ=255)
                                                                     aTAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGCTAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGACGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCTTCGGCAGGCCAATCGTGCTTTACAGCTCCCCACTTTCAATTAATGCCGGAAtccgtc                                                                                                                                                                             >  2:630143/1‑251 (MQ=255)
                                                                        aGATGAGTGCGGCCTTGCTCATCGTACCGTTTATCGAACAGTAACCGCTTCCCAGTAATGCTCCGCTGTCATATGGCCCAGCCGGAGACGTCAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCCTCACCTAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGCGCTTTCCACTTCCCCACTTTCAATTTATGCCGGTATCCTTCCgg                                                                                                                                                                          >  2:408912/1‑251 (MQ=255)
                                                                              gTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCCATCGGGCTTTCCATTTTCCCCCTTTCAATTAATGCCGGTATCCGTCCGGGGAGCg                                                                                                                                                                    >  2:530422/1‑251 (MQ=255)
                                                                                   gATTTGCTCATCGGGCCGTTTATCGATAAGTAACCGATTACCAGTAATGCTCCACTGTCATCTGGCCCGGTCGGCGCCGCAAATATTTAGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGAAAAGCATAACATAGCTGGTTTCTTTCTTCATCGGCAAGCCACTCGGTATTTACAGTTCCCCTGTTTCAATTAATGCCAGTATCCGTACGGGGGGCGAcccc                                                                                                                                                               >  2:566807/1‑251 (MQ=255)
                                                                                     tttGCTCATCGTGCCGGTTATCGATAAGTAACCGCTTACCAGTAATTCTCCGCGTTCATATGGCCCGGTCGGCGTCGTAAATGGTTCGTCATCTGCCGGGTCCCTTGCAGCACCTGTTGTGTATCGCGCTCCATCTGTGGATTTCCTCACATCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCCACCCCg                                                                                                                                                              <  2:322384/250‑1 (MQ=255)
                                                                                       tGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTg                                                                                                                                                           >  1:936028/1‑251 (MQ=255)
                                                                                           caccGTGCCGTTTAGCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGt                                                                                                                                                       <  2:33350/248‑1 (MQ=255)
                                                                                              cGTGCCTTTTATCGATAAGTTACCGCTTACCAGTAATGCTTCGCTGTCATATGGCCCGGTCGGCTACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCTACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGTCCCCGCTGCCGttt                                                                                                                                                     <  2:850575/250‑1 (MQ=255)
                                                                                               gggccgttgtacggtAAGTAACCGATTTCCAGTATGGTTCCGCTGTCATATGGCCCGGTCGGCGGCGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTAGTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGGGAGCGACCCCGCTGCCGttgcc                                                                                                                                                   <  2:732624/237‑4 (MQ=255)
                                                                                                tGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCAGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGTGAGACCCCGCTGCCGtttacc                                                                                                                                                  >  2:362904/1‑248 (MQ=255)
                                                                                                 gCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCg                                                                                                                                                 >  1:796246/1‑251 (MQ=255)
                                                                                                       tATCGATAAGTAACCGCTTACCGGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGa                                                                                                                                            <  1:35138/250‑1 (MQ=255)
                                                                                                        aTCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCtt                                                                                                                                                                                                                                            >  1:929050/1‑153 (MQ=255)
                                                                                                        aTCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTAGCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCtt                                                                                                                                                                                                                                            <  2:929050/153‑1 (MQ=255)
                                                                                                          cGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGaccac                                                                                                                                        >  1:636151/1‑251 (MQ=255)
                                                                                                          cGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGATTCCCGACTGaccac                                                                                                                                        >  1:152512/1‑251 (MQ=255)
                                                                                                           gATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGaccacc                                                                                                                                       <  1:875173/251‑1 (MQ=255)
                                                                                                              aaGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCCAAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATCACGGTATCCGTCCGGTGAACGACCCCGCTGCCGTTTCCCGACTGCCCACCggc                                                                                                                                    >  2:386595/1‑249 (MQ=255)
                                                                                                               agtatcctcTTACCAGAAATGCTCCGCTGTCATATGGCCCGGTTGGCGACGGAAATGTTTCGTCATCTGCCGGGTCACTTTCAGCAACTTTTGTTTATCGCTCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAAGGGTGCTTTCCAGTTCGCCACTTTCAATTAATTCCGCTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTGCCGACTGACCACCGTCt                                                                                                                                   <  2:432163/243‑1 (MQ=255)
                                                                                                               aGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCt                                                                                                                                   <  1:387087/251‑1 (MQ=255)
                                                                                                                  aaCCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCCACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTACTGCCGGTATCCGTCCGGTGCGCGACCCCGCTGCCGTTTTCCGACTGACCACCGTATGa                                                                                                                                 >  2:227318/1‑250 (MQ=255)
                                                                                                                      gCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATgcg                                                                                                                            <  1:544286/251‑1 (MQ=255)
                                                                                                                         tACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATGAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCa                                                                                                                          <  2:417309/250‑1 (MQ=255)
                                                                                                                         gaacAGTAATGCCCCGCTGGCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCGGCCGTTTCCCGACTGACCACCGTCTGAATGCGCa                                                                                                                          <  2:308156/247‑1 (MQ=255)
                                                                                                                          ccaaggaaTGCTCCGCGGCTATATGGGCCGGTCGGCGACGTAAATGTTTCGTCATCTGACGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCCCAGCATCACATAGCTGTTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCATTTCCCCACTTTCAATGAATGCCGGTATCCGTCCGGTGAGCGATCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATTCGCAg                                                                                                                         <  2:61420/244‑1 (MQ=255)
                                                                                                                          aCCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTTCATTAATGCCTGTATCCGTCCGGTGGGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGactgtgccgt                                                                                                                        >  2:842039/1‑242 (MQ=255)
                                                                                                                          aCCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGt                                                                                                                        <  1:699368/251‑1 (MQ=255)
                                                                                                                              gTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGAGCTTTCCAGTTCACCCCTTTCAATTAATTCCGGTATCCGTCCGGTGAACGCCCCCGCCGCCGTTTCCCGACTGACCACCGTCTGACGGCGCAGttttt                                                                                                                    >  2:815791/1‑250 (MQ=255)
                                                                                                                                      ccGCTTTCATATGGCCCGGTCAGCGACGTAAATGTTTCGTCATCTGCCGGGACTCTATCAGCAACTGGTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCGTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTGTCAATTAATGCCGGTATCCGGCCGGTGAGCGACCCCGCGGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCGTTCGTat                                                                                                            <  2:477175/251‑2 (MQ=255)
                                                                                                                                      ccGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGGCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAg                                                                                                            <  1:227318/251‑1 (MQ=255)
                                                                                                                                          tGTCATATGGCCCGGTCGGCGCCGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGt                                                                                                         <  2:819457/250‑1 (MQ=255)
                                                                                                                                            tCATATGGCCCGGTCGGCGACGTCAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTGATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGCCCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGtttt                                                                                                      >  2:202775/1‑251 (MQ=255)
                                                                                                                                                     ccGGTCGGCGACTTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTcc                                                                                             <  2:552995/251‑1 (MQ=255)
                                                                                                                                                     ccGGTCGGCGACGTTAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTcc                                                                                             <  2:425739/251‑1 (MQ=255)
                                                                                                                                                           ggcaaCGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGAGGATTGCCGCACAGAATCACATAGATGGTTTCTTTATTCATAGACAGGCAAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATACGTCCGGTACGCGACCCCGCTGCAGTTTCCAGACTGACCACCGTCTGAATGCGCCGTTTTCATTAGCAGAGTTTTTCCAGTTCCTGCATc                                                                                       >  2:545720/5‑251 (MQ=255)
                                                                                                                                                           ggCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATc                                                                                       >  2:234179/1‑251 (MQ=255)
                                                                                                                                                             cGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATc                                                                                       >  1:539886/1‑249 (MQ=255)
                                                                                                                                                               aCGTAAATGTTTCGTCATCTGCCGGGTCTCTTTAAGCAACGGTTGTGTATCGCGCACCATCTGTGGATTGCTGCACTGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTGCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGATCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTg                                                                                   <  2:901959/251‑1 (MQ=255)
                                                                                                                                                                cGTAACTGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTGCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGGGAGCGACCCCGCTGCCGTTTCCCGACCGACCCCCGTCTGACTGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCCTCAgggg                                                                                  >  2:442576/1‑248 (MQ=255)
                                                                                                                                                                         ttCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGCGGATTGCCGCACAGCCTCACCTAGCTGGTTTATTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCCACGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCt                                                                         >  1:924443/1‑251 (MQ=255)
                                                                                                                                                                            gTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCTGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTaa                                                                      <  1:9038/251‑1 (MQ=255)
                                                                                                                                                                             tCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGATTTCCAGTTCCCCCCTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCAGACTGACCACCGTCTTAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCaaaacgctta                                                                       >  2:501448/1‑242 (MQ=255)
                                                                                                                                                                                 cTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGACGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCTAGTTCCTGCATCAGTGGCAAATAGCTTATGTcgg                                                                 >  1:264909/1‑251 (MQ=255)
                                                                                                                                                                                  tGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGAGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTcggc                                                                >  1:454098/1‑251 (MQ=255)
                                                                                                                                                                                  tGCCGGGTCTCTTTCAGAAACTGTTGTGTATCGCGCACCATCTGTCGATTGCCGCACAGCATAACCTAGATGGTTTCTTTATTCATCGGCAGGCCCCTCGCGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGCGCGCTACCCCGCTAACGTTTCCCAACTGACCAACGTCTGAATGCATAGTTTTCCTTCGGAGCGTTTTTCCAGTTCCCGCATCACTGGCAATTAGCTTAAGTaggc                                                                >  2:55445/1‑251 (MQ=255)
                                                                                                                                                                                   gccggtcctctTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTTGCAAATAGCTTAAGTcggcg                                                               <  2:973457/244‑1 (MQ=255)
                                                                                                                                                                                      gggTCTCTTTCAGCAACTGTTGAGTATCGCGCACAATCTGTGGATTGCCGCACAGCATCACATAGCTGGATTCTTTATTCATCGGCAGGCCAATCGTGCTTTACAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGACAAGCGACCCCGCTGACGTTTCCCGACTGAACACCGTCTGACTGCGCAGTTTTCCTTCGTAGCGGTTTTCCCGTTCCTGCATCAGTGGCCAATAGCttcattctgcggct                                                            >  2:374579/1‑239 (MQ=255)
                                                                                                                                                                                                cAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTgcgg                                                  >  1:106342/1‑251 (MQ=255)
                                                                                                                                                                                                cAGCAACTGTTGTGAATAGCGCAACATCTGTGGCTTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCATGCTTTCCAGTTCCCAACTTTCACTTAATGCCGGTATCCGTTCGGTGAGCGACCCCACTGACGTTTCCCCCCTGAACACCGTCCGCATGCGGAGTATTCCTTCGTCGAGCTTTTCCAGTTTCTGCCTCCGTAGCAACTATCTTAACTCGGCGGCATAACgaacgg                                                  >  2:630309/1‑246 (MQ=255)
                                                                                                                                                                                                 aGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCCCCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTAGGCGGCATAACGTgcggc                                                 >  2:786594/1‑251 (MQ=255)
                                                                                                                                                                                                     aCTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATGGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGc                                                                          <  2:459498/222‑1 (MQ=255)
                                                                                                                                                                                                     aCTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATGGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGc                                                                          >  1:129708/1‑222 (MQ=255)
                                                                                                                                                                                                     aCTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATGGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGc                                                                          >  1:459498/1‑222 (MQ=255)
                                                                                                                                                                                                     aCTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATGGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGc                                                                          <  2:129708/222‑1 (MQ=255)
                                                                                                                                                                                                      cttttgtGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGATTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTACCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGGCTGACCACCGTCTGAATACGCAGGTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGAATCACGTGCGGCGTGcg                                            >  1:928427/4‑250 (MQ=255)
                                                                                                                                                                                                           gtttttCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCTGCAGGCCACTCGTGCTGTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCTGTGAGCGCCCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAgg                                       <  2:709257/246‑1 (MQ=255)
                                                                                                                                                                                                                 cgcgCACCAGCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTAGTCAGCGGCAGGCCAATCGTGCGTTCCAGTGCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCTCTTCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTTGCAAAGAGCTTAACGGGGCGGCATAACGGGCTCCGTGCACCAGGACCAg                                  <  2:486739/250‑1 (MQ=255)
                                                                                                                                                                                                                    gcACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGACCAGAttt                              >  2:212487/1‑251 (MQ=255)
                                                                                                                                                                                                                       ccATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGACCAGATTTTTg                           <  1:234179/251‑1 (MQ=255)
                                                                                                                                                                                                                        cATCTGTGGATTGCCGCACAGCATCCCATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCCGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGGCGCGTTTTTCCCGTTCCTGCATCAGTGGCAAATAGCTTATGTCGGCGGCATAACGTGCGGCGTGCACCAGGAACAGATTTTTgg                          >  2:395980/1‑250 (MQ=255)
                                                                                                                                                                                                                             gtgGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGACCAGATTTTTGAAGCGa                     <  1:799038/251‑1 (MQ=255)
                                                                                                                                                                                                                                gaatGCCACACAGCAGCTCATAGCTGGTATCTTTTTTCATCGGCAGGCCCATCGCGCTTTCCAGTGCCCATCTACCAATTAATGCCTGTATCTGTCCTGCGAGCCACCCCGCTGCCGTTTCCCGACGGACCACCGTCTGAATGCGCAGTTTACCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGTCGGCATTACGTGCGGCATGTACCAGGAGCAGATTTTTGAAGCGATCt                  <  2:602956/248‑1 (MQ=255)
                                                                                                                                                                                                                                gATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGACCAGATTTTTGAAGCGATCt                  <  2:270856/251‑1 (MQ=255)
                                                                                                                                                                                                                                         actgCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTGCAAGTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGACCAGATTTTTGAAGCGATCTAAATCTTTa         <  2:434892/248‑1 (MQ=255)
                                                                                                                                                                                                                                                cacaTAGCTGGTTTCTTTATTCATTGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGACCAGCTTTTTGAAGCGATCTAAATCTTTACCTAGtt  >  2:295634/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
ACCTCCTTGTCTGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATGGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGACCAGATTTTTGAAGCGATCTAAATCTTTACCTAGTT  >  NC_000913/4113606‑4114096

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