Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 66,784 A→T intergenic (‑234/+51) araD ← / ← araA L‑ribulose‑5‑phosphate 4‑epimerase AraD/L‑arabinose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091366,7840AT98.3% 312.4 / ‑1.5 120intergenic (‑234/+51)araD/araAL‑ribulose‑5‑phosphate 4‑epimerase AraD/L‑arabinose isomerase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base T (59/57);  total (61/59)
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.83e-01

Only 100 of 120 total aligned reads displayed.

GATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTGTCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGT  >  NC_000913/66538‑67029
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
gaTCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGCCGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGCATGCCGGATGCGCTTTGCTTATCCGGCCTCCAAAATCGCAGCGTGTCGGCCTGAt                                                                                                                                                                                                                                                     >  2:404430/1‑250 (MQ=18)
    ttctAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCACTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAggc                                                                                                                                                                                                                                                >  2:414524/1‑249 (MQ=255)
      ctAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGTCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCTCATCAGGCGTTGAATGACGGATGCGCTTTGCTTATCCGGCCTACCAAATCGGAACGCGTAGGCCTGATAAGACTCGCCAGCGTCGCATACGACGTTGCATGCCGGATGCGCTTTGCTTATCCGGCCTACGAAATCGCAGCGTGTAGGCCTGATATGAcgc                                                                                                                                                                                                                                              >  2:173395/1‑251 (MQ=255)
        aataTGTTGACTCCTTCGTGCCAGATGCGCTTTGCTTATCCGGCCTACAAAAGCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGAGGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGAcgc                                                                                                                                                                                                                                              <  1:438982/246‑1 (MQ=21)
         aCATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTACGCCTGCTAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGCTAAGACGCGTc                                                                                                                                                                                                                                           >  2:241907/1‑251 (MQ=255)
                ggcgccTTCGTGCCGGATGCGCGTTGCTCAACCCGGCTACAAAATCGCAGCGTGTAGGCCGGATAAGAAGCGCCAGCGTCGAAGCAGGCGTTGAAGGCCGGATGCGCTTTGCGTATCCGGACTACAAAATCGGAACGCGTAGGCCTGAGAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCg                                                                                                                                                                                                                                    <  2:196713/247‑1 (MQ=18)
                     cTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAGTCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACTCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCCGGCGTTAAATGCCGGATGCGCTTTGCTTATCCGGCCTACCAAATCGCAGCGCGTAGGCCTGATAAGACGCGGCAGCGgcgccac                                                                                                                                                                                                                                >  2:708766/1‑247 (MQ=18)
                       tCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCCGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACCAAATCGGAACGAGTAGGCCTGATAAGCAGCGCCAGCGTCGCACACGGCGTTGAATGCCGGATGCGCTTTGCTTATCCTGCCGACAACCTCGAAGCGCGTAGGCCTGATGAGACGCGTCAGCGTCGTATCAgg                                                                                                                                                                                                                             >  2:908744/1‑251 (MQ=17)
                       tCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGTCCTGATAAGACGCGTCAGCGTCGCATCa                                                                                                                                                                                                                               <  1:883853/249‑1 (MQ=255)
                        cGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGATTATCCGGCCTACAAAATCGCCGCGTGTGGGCCTGTTCAGACGCGTCAGCGTCGCATCAggg                                                                                                                                                                                                                            >  2:929328/1‑250 (MQ=18)
                           gCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGtt                                                                                                                                                                                                                         >  1:306699/1‑251 (MQ=255)
                           gCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGtt                                                                                                                                                                                                                         >  1:135556/1‑251 (MQ=255)
                             cGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTAc                                                                                                                                                                                                                       <  1:8746/251‑1 (MQ=255)
                               gATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACac                                                                                                                                                                                                                     <  1:556778/251‑2 (MQ=255)
                               gATGCGCTGTGCTCGTCCGGCCGACAAATGCGCAGTGTTGAGGCCTGATAAGACGCGCCCGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCGTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGGCGCGCCAGTGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACac                                                                                                                                                                                                                     <  2:524771/251‑2 (MQ=17)
                                    gcTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACAcctg                                                                                                                                                                                                                <  2:413114/251‑3 (MQ=255)
                                       ttGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGAt                                                                                                                                                                                                              <  1:176282/250‑1 (MQ=255)
                                       ttGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGAt                                                                                                                                                                                                              <  1:791142/250‑1 (MQ=255)
                                        tGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGc                                                                                                                                                                                                            >  1:313327/1‑251 (MQ=255)
                                             aTCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACTCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCAGACCTACAAAATCGCAGCGTGTAGGCCTGATCAGAAGCGTCCGCGTCGCATCAGGCGGTACACAcctgcttcgaccc                                                                                                                                                                                                       >  2:734619/1‑240 (MQ=18)
                                             aTCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCTGCCTACAAAATCGCAGCGTGTAGGCCTGCTAAGACGCGTCGGGGTCGCATCAGGCGTTACACACCTGATGCgact                                                                                                                                                                                                        >  2:47930/1‑247 (MQ=21)
                                             aTCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAACGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACTCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCgact                                                                                                                                                                                                        <  2:442209/250‑4 (MQ=21)
                                              tCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTAc                                                                                                                                                                                                      >  1:95149/1‑251 (MQ=255)
                                                 ggCCTACAAATGCGCAGCGTGTAGGCCTAAGAAGACGCGACAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCGTAGCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTa                                                                                                                                                                                                   <  2:512221/251‑1 (MQ=18)
                                                          aaTCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGaaa                                                                                                                                                                                          <  1:655068/251‑1 (MQ=255)
                                                           aTCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGaaat                                                                                                                                                                                         >  1:511647/1‑250 (MQ=255)
                                                                 gCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTa                                                                                                                                                                                    >  1:271955/1‑250 (MQ=18)
                                                                     gtAGGCCGGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGGCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATac                                                                                                                                                                                <  2:420781/250‑1 (MQ=18)
                                                                              aTAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGCCGAAATCCGTAATACACTTCGTTcc                                                                                                                                                                      >  1:612003/1‑251 (MQ=18)
                                                                               tAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGCTGCGCTTTGCTTATCCGGCCTACAAAATCGCACCGTGTATGCCTGCTAAGACGCGTCAGCGTCGCACCAGGCGTTACACACCTGATGCGCCTACTTAGCGACGAAATCCGTAATCCCCTTCTTTccc                                                                                                                                                                     >  2:693596/1‑250 (MQ=18)
                                                                                aaGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCca                                                                                                                                                                     <  1:970517/250‑1 (MQ=18)
                                                                                aaGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCca                                                                                                                                                                     <  1:524073/250‑1 (MQ=18)
                                                                                       gcccgCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCGTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCcagcgcagcg                                                                                                                                                             <  2:341194/247‑1 (MQ=18)
                                                                                            cGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCt                                                                                                                                                        >  1:987664/1‑251 (MQ=18)
                                                                                                gCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTa                                                                                                                                                     >  1:475848/1‑250 (MQ=18)
                                                                                                    cAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCTTGTAGGCCTGCTCAGACTCTTCAGCGTCGCATCAGGCGTTACACACCTGCTGCGACTACATATCGACGACATCCGTAATACCCTTCGTTCCAGCGCAGCGCGGCTTTACACGc                                                                                                                                                >  2:6759/1‑251 (MQ=17)
                                                                                                     aGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGCTAAGACGAGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATACGTAATACACTTCGTTCCCGCGCAGCGCGTCTTTAAACgat                                                                                                                                               >  2:575496/1‑249 (MQ=18)
                                                                                                     aGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCt                                                                                                                                               >  1:439771/1‑251 (MQ=255)
                                                                                                      ggCGTTGAATGTCGGCTGCGCTTTGCTTACCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCTCATCAGGCGTTGAAGGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTg                                                                                                                                              <  2:135323/251‑1 (MQ=18)
                                                                                                            gAATGCCGGTAGCTCTTTGCTTATCCGGCCAACCAAAACGGAACGCGGAGGCCTGAAAAGACGCGCCGGCGTCGCATCAGGTGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAAGCCGTAATACACTGCGTTCCAGCGCAGCGCGTTTTTAAACGCTGGCAGGc                                                                                                                                        <  2:571900/251‑1 (MQ=18)
                                                                                                            gAATGCCGGATTCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGGAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAgg                                                                                                                                         <  2:799294/250‑1 (MQ=18)
                                                                                                                   ggATGCGCTTTGCTTATCCGGCCTACAAGATCGGACCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGCATGCCGGATGCGCTTTGCTTATCCGGCCTACAAACTCGCAGCGTGTAGGCCTGATAAGACGCGTCAGGGTCGCATCAGGCATTACACCCCTTATGCGACTCCTTAGCGACACAATCCGAAATACGCTTGGTTCCAGCGCCTGGCGTCTTTAAACGCTAGCAGGCgtggtgc                                                                                                                                 >  2:475929/1‑247 (MQ=9)
                                                                                                                    gATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACCCCCCTGATGCGACTACTTAGCGTCGAAATCCGTAATAAACTTCGTTCCAGCGCACCGCGTCTTTAACCGCTGGCAGGCggttatcg                                                                                                                                >  2:223981/1‑244 (MQ=17)
                                                                                                                          cTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGCCGAAATCCGTAATACACTTCGTTCCAGAGCAGCGCGTCTTTAAACGCTGGCAGGCgggtattcttatc                                                                                                                           >  1:537598/1‑238 (MQ=18)
                                                                                                                          cTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTatc                                                                                                                           >  1:372730/1‑250 (MQ=18)
                                                                                                                            ttGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGTCGCGTCAGCGTCGCATCAGACGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGCTCCAGCGCAGCGCGTCTTTAAACGCTGACAGGCggatatcgtttacaat                                                                                                                        >  2:13493/1‑236 (MQ=18)
                                                                                                                              gCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGGCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTatcaatca                                                                                                                      <  1:177175/251‑1 (MQ=21)
                                                                                                                              gCTTAGCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTatcaatca                                                                                                                      <  1:600919/251‑1 (MQ=21)
                                                                                                                               cTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTCCACACCTGATGCGCCTACTTAGCGACGACATCCGTAATACCCTTCGTTCCAGCGCAGCGCGACTTTAAACGCTGGCCGGCggggatcgttgtcactcac                                                                                                                     >  2:691235/1‑233 (MQ=17)
                                                                                                                                ttATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGAACGTTATCAATCAcc                                                                                                                    >  1:158268/1‑251 (MQ=18)
                                                                                                                                 tcaccgtccagcAAAATCGGAACGCGTAGGCCTGCTTAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATTCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTGAGCGACGAAATCCGTAATTCACTTCGTTCCAGCGCCGCGCGTCTTTAAACGCTGTCAGGCGGGTATCGTTATCAATCAc                                                                                                                     <  2:67473/238‑1 (MQ=18)
                                                                                                                                      ggCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGCCTAATTAGCGACGAAATCCGTAATAAACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCggggatcggttttaatcaacctggtt                                                                                                              >  2:878655/1‑226 (MQ=18)
                                                                                                                                       cccTACAAAATCGTTACGCGTAGGCCTGATAATGCGCGGCGGCGTCGCAGCAGGCGTTGAATGCCGGATGCGATTTGCTTACACGGGCTACAAAATCGCATCGTGTAGGCCTGATAAGACGCGTCAGCGTAGCATCAGACGTTACACACCTTAGGCGTCTATTTAGCGACGTAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTTGCAGGCGGGTATCGTTATCAATCTCCGTTAttt                                                                                                             <  2:835880/250‑1 (MQ=255)
                                                                                                                                           aCAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCaa                                                                                                          <  1:575701/250‑1 (MQ=21)
                                                                                                                                                      acgcgTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCAtct                                                                                              <  1:738691/250‑1 (MQ=25)
                                                                                                                                                        gcgTAGGCCTGATAAGACGCGCCAGCTTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCCCCGTGATTTCAATGTCGTGCATCTCg                                                                                            <  1:215389/251‑1 (MQ=25)
                                                                                                                                                        gcgTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCg                                                                                            <  1:172253/251‑1 (MQ=25)
                                                                                                                                                                  gATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTACTACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCATTCACCGTGATTTCAATGTCGTGCATCTCGGTGAACTGgc                                                                                  >  1:515350/1‑251 (MQ=25)
                                                                                                                                                                   aTAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGgcg                                                                                 <  1:459363/251‑1 (MQ=25)
                                                                                                                                                                   aTAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGgcg                                                                                 <  1:621931/251‑1 (MQ=25)
                                                                                                                                                                   aTAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGgcg                                                                                 <  1:47930/251‑1 (MQ=25)
                                                                                                                                                                            gcCAGCGTCGCATCAGGCGTTGAATGCCGGCTGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGAGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCCGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTATCAATCACCGTACTTTTCATGGCGTGCATATCGGCGGAATGGCGCATCTCGtt                                                                        >  1:756626/1‑251 (MQ=18)
                                                                                                                                                                            gcCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGtt                                                                        <  1:693596/251‑1 (MQ=37)
                                                                                                                                                                               aGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGGTCATCTCGGCGAACTGGCGCATATccttt                                                                       >  1:300319/1‑245 (MQ=25)
                                                                                                                                                                                gCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAgg                                                                    <  1:414524/251‑1 (MQ=37)
                                                                                                                                                                                ctcttcgcctcctgCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGCCGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTTTCCATCCCCGTGATTTCCATGTCGTGCATCTCGGCGAACTGGCGCATATCCTTggg                                                                     >  1:168354/14‑249 (MQ=25)
                                                                                                                                                                                  gTCGCATCAGGCGTTGAATGCCGGATGCGCTTTTCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTATCAATCAACGTGATTTCAATGTCGTGCATCTCGGCGAAATGGCGCCAATCGTTGAGGtt                                                                  >  1:217463/1‑251 (MQ=25)
                                                                                                                                                                                      cATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGGTATAAATCACAGTGGTTTCAATGTAGTGCATCTCGGCGAAATGGGGCATATCGGTGAGGttaagg                                                              >  1:844684/1‑247 (MQ=21)
                                                                                                                                                                                      cATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGAACGTTATCAATCACCGGGATTTTTATATCGTGCATCTCGGCCGACTGGCGCATATTGTTGTGGttaggc                                                              >  2:175604/1‑247 (MQ=21)
                                                                                                                                                                                         cAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCa                                                           <  1:929328/251‑1 (MQ=37)
                                                                                                                                                                                         cAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTATCAATCCCCGGGATTTCAATGTCGTGGATCTCCGCGAACTGGCGCATATCGTTGAGGTTCAgcgc                                                            >  1:429794/1‑247 (MQ=25)
                                                                                                                                                                                               gTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTATCCATCACCGTGATTTCAATGTCGTACATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAgcgccttgct                                                      >  1:599797/1‑242 (MQ=37)
                                                                                                                                                                                                 ttgaggCCGGATGCGCTTTGCGTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCATGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTGATCAATCACCGTGATTTCAAGGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTg                                                    <  2:440771/245‑1 (MQ=37)
                                                                                                                                                                                                   aaTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGCCGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGTTAGTTATAAAACACCGGGGTTTCAATGTTGTGTATCTCGGCGCACTGGCGCCTATCCTTTAGGTGACGTGCATGGccgcc                                                  >  2:601392/1‑246 (MQ=16)
                                                                                                                                                                                                    aTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGGTCAGCGCATGGCTGa                                                   >  1:461821/1‑248 (MQ=37)
                                                                                                                                                                                                      gCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACg                                              <  2:313327/251‑1 (MQ=38)
                                                                                                                                                                                                        cggatggtcTTTGATTATCCGGCCTACAAAAACGCAGCGGGTAGGCCCGATAAGACGCGGCAGCGTCGCATCAGTCGTTACACATCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTAGCGTTATCAAGCACCGTGATTTTAATGTCGTGCATCTCGGCGAACTGTCGCATATCGTTGAGGTTCAGCGCATGGCTGAATACGGt                                            <  2:687803/243‑1 (MQ=255)
                                                                                                                                                                                                           attagcTTTGGTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTCCTTAGCGACGAAATTCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATg                                         <  2:327825/247‑1 (MQ=37)
                                                                                                                                                                                                             gcgcTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATgg                                        <  1:137351/250‑1 (MQ=38)
                                                                                                                                                                                                              cgcTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTGCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTg                                      <  1:218024/251‑1 (MQ=38)
                                                                                                                                                                                                               gcTTGGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTgc                                     <  2:372730/251‑1 (MQ=38)
                                                                                                                                                                                                                cTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGGTATCAATCACCGTGGTTTCCATGTCGTGCATCTCGGCGAACTGGCGCGTATCGTTGTGGTTCAGCGCATGGCTGAAGACGGTATTGTgcg                                    >  1:24781/1‑251 (MQ=37)
                                                                                                                                                                                                                       tATCCGGCCTACAAAATCTCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGcca                               <  1:77961/249‑1 (MQ=255)
                                                                                                                                                                                                                        aTCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGGGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGccacca                            <  2:810264/251‑1 (MQ=255)
                                                                                                                                                                                                                            ggCCTAAAAAATCGTAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTGACACACCTGATGCGACTACTCAGCGACGAATTCCGTAATACACTTCGTTCCAGCGCATCGCGTCTTTAAACGCTGGCAGGCGGGTAGCGTTATCAATCACCGGGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGa                        <  2:285197/251‑1 (MQ=255)
                                                                                                                                                                                                                             gCCTACAAAATCGCAGAGTGTAGGCCTGATAAGACGCGTAAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGCATCGTTTTCAATCACCGTGATTTCAAATTCGTGCATCTCGGCGAACTGGGGCGTAACGTTGAGGTTCCACGCACGGCTTGAGACGGTATGGTGCGCGccccccgaggg                       >  2:595543/1‑242 (MQ=255)
                                                                                                                                                                                                                              ccTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCCCCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCTCATATCGTTGAGGGTCCGAGCATGGCTGCAGACGGTATGGTGCGCGCCCCCCGCGAgg                      >  2:347618/1‑251 (MQ=255)
                                                                                                                                                                                                                                   aaaaACGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCa                 <  2:135556/251‑1 (MQ=255)
                                                                                                                                                                                                                                    aaaTCGCCGCGTGTAGGCCTGATAAGACGCGACAGCGTCGCAACAGGCGTTACGCACCTGATGCGACTACGTAGCGACGAAATCCGTAATACACTTCGTTCCATCGCAGCGCGTCTTTAAACGCTTGCAGGCGGGTATCGGTGTCAATCACCGTGATTGCAAGGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGATGATCCAc                <  2:205031/251‑1 (MQ=255)
                                                                                                                                                                                                                                    aaaTCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGCTGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTTTCAATCACCGAGATTTATATGTCGGGCATCTCTGCGAACTGGCGCCTATCGTTGGGGGTCAGCGCATGGCTTGAGACGGGATGGCGCGCGccccccgcgcgagtacac                >  2:758493/1‑235 (MQ=255)
                                                                                                                                                                                                                                     aaTCGCAGCGTGTAGGCCTGATAAGCCGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATTGTTATCAATCACCGTGATTTTAATGTCGTGACTCCCGTGGAACTGGCGCATATCGATGGGGGTCAGCGCACGGCTGAAGACGGTATGGGgcgccacccccacgaggagccccc               >  2:512931/1‑231 (MQ=255)
                                                                                                                                                                                                                                        cGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACCCTTCGTTCCAGCGCATCGCGTCTATAAACGCTGGCAGGCTGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGATGTTCAGCGCATGGCTGAAGACGGTTTGGTGCGCGCCACCAGCGAGGATCCACGCtt            <  1:584342/251‑1 (MQ=255)
                                                                                                                                                                                                                                         gCAGCGTGTAGGCCTGATAATACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCCCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTc           <  1:911505/251‑1 (MQ=255)
                                                                                                                                                                                                                                         gCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTGTCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTc           <  1:496030/251‑1 (MQ=255)
                                                                                                                                                                                                                                          cAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCg          <  1:6759/251‑1 (MQ=255)
                                                                                                                                                                                                                                           aGCGTGTATGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGACATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTATAAATAACCGTGATTTCAATGTCGGGCAACCCGGCGAACATGCGCCTATTGTTGAGGTTCAGCGACTGGCCGAAGAAGGACTGGTTCGCGCCACCAGCCAGGacacacgccttcg         >  2:167525/1‑239 (MQ=255)
                                                                                                                                                                                                                                           aGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGCTGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGACGGGTCTCGTTATCAATCACCGTGATTTACATGTCGTGCAACTCGGCGAACTGGAGCATGACGTTGATGGTCAGGCCTAGGCTGAAGAATGTATGGGTCGAGCCACACGCGAGGATCCACGccttgg         >  2:213804/1‑246 (MQ=255)
                                                                                                                                                                                                                                           aGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAGCGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGAGCCTCGCTTCgg         <  2:708176/251‑1 (MQ=255)
                                                                                                                                                                                                                                            gCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAAACACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCgg         >  1:449839/1‑250 (MQ=255)
                                                                                                                                                                                                                                                  aGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACGTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGt  <  2:632408/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
GATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTGTCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGT  >  NC_000913/66538‑67029

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: