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.

GATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTGTCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGT  >  NC_000913/66538‑67029
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
gaTCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGCCGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGCATGCCGGATGCGCTTTGCTTATCCGGCCTCCAAAATCGCAGCGTGTCGGCCTGAt                                                                                                                                                                                                                                                     >  2:404430/1‑250 (MQ=18)
    ttctAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCACTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAggc                                                                                                                                                                                                                                                >  2:414524/1‑249 (MQ=255)
      ctAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGTCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCTCATCAGGCGTTGAATGACGGATGCGCTTTGCTTATCCGGCCTACCAAATCGGAACGCGTAGGCCTGATAAGACTCGCCAGCGTCGCATAC‑GACGTTGCATGCCGGATGCGCTTTGCTTATCCGGCCTACGAAATCGCAGCGTGTAGGCCTGATATGAcgc                                                                                                                                                                                                                                              >  2:173395/1‑251 (MQ=255)
        aataTGTTGACTCCTTCGTGCCAGATGCGCTTTGCTTATCCGGCCTACAAAAGCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGAGGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGAcgc                                                                                                                                                                                                                                              <  1:438982/246‑1 (MQ=21)
         aCATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTACGCCTGCTAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGCTAAGACGCGTc                                                                                                                                                                                                                                           >  2:241907/1‑251 (MQ=255)
                ggcgccTTCGTGCCGGATGCGCGTTGCTCAACCCGGCTACAAAATCGCAGCGTGTAGGCCGGATAAGAAGCGCCAGCGTCGAAGCAGGCGTTGAAGGCCGGATGCGCTTTGCGTATCCGGACTACAAAATCGGAACGCGTAGGCCTGAGAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCg                                                                                                                                                                                                                                    <  2:196713/247‑1 (MQ=18)
                     cTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAGTCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACTCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCC‑GGCGTTAAATGCCGGATGCGCTTTGCTTATCCGGCCTACCAAATCGCAGCGCGTAGGCCTGATAAGACGCGGCAGCGgcgccac                                                                                                                                                                                                                                >  2:708766/1‑247 (MQ=18)
                       tCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCCGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACCAAATCGGAACGAGTAGGCCTGATAAGCAGCGCCAGCGTCGCA‑CACGGCGTTGAATGCCGGATGCGCTTTGCTTATCCTGCCGACAACCTCGAAGCGCGTAGGCCTGATGAGACGCGTCAGCGTCGTATCAgg                                                                                                                                                                                                                             >  2:908744/1‑251 (MQ=17)
                       tCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGTCCTGATAAGACGCGTCAGCGTCGCATCa                                                                                                                                                                                                                               <  1:883853/249‑1 (MQ=255)
                        cGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGATTATCCGGCCTACAAAATCGCCGCGTGTGGGCCTGTTCAGACGCGTCAGCGTCGCATCAggg                                                                                                                                                                                                                            >  2:929328/1‑250 (MQ=18)
                           gCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGtt                                                                                                                                                                                                                         >  1:306699/1‑251 (MQ=255)
                           gCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGtt                                                                                                                                                                                                                         >  1:135556/1‑251 (MQ=255)
                             cGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTAc                                                                                                                                                                                                                       <  1:8746/251‑1 (MQ=255)
                               gATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACac                                                                                                                                                                                                                     <  1:556778/251‑2 (MQ=255)
                               gATGCGCTGTGCTCGTCCGGCCGACAAATGCGCAGTGTTGAGGCCTGATAAGACGCGCCCGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCGTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGGCGCGCCAGTGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACac                                                                                                                                                                                                                     <  2:524771/251‑2 (MQ=17)
                                    gcTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACAcctg                                                                                                                                                                                                                <  2:413114/251‑3 (MQ=255)
                                       ttGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGAt                                                                                                                                                                                                              <  1:176282/250‑1 (MQ=255)
                                       ttGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGAt                                                                                                                                                                                                              <  1:791142/250‑1 (MQ=255)
                                        tGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGc                                                                                                                                                                                                            >  1:313327/1‑251 (MQ=255)
                                             aTCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACTCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCAGACCTACAAAATCGCAGCGTGTAGGCCTGATCAGAAGCGTCCGCGTCGCATCAGGCGGTACACAcctgcttcgaccc                                                                                                                                                                                                       >  2:734619/1‑240 (MQ=18)
                                             aTCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCTGCCTACAAAATCGCAGCGTGTAGGCCTGCTAAGACGCGTCGGGGTCGCATCAGGCGTTACACACCTGATGCgact                                                                                                                                                                                                        >  2:47930/1‑247 (MQ=21)
                                             aTCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAACGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACTCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCgact                                                                                                                                                                                                        <  2:442209/250‑4 (MQ=21)
                                              tCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTAc                                                                                                                                                                                                      >  1:95149/1‑251 (MQ=255)
                                                 ggCCTACAAATGCGCAGCGTGTAGGCCTAAGAAGACGCGACAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCGTAGCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTa                                                                                                                                                                                                   <  2:512221/251‑1 (MQ=18)
                                                          aaTCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGaaa                                                                                                                                                                                          <  1:655068/251‑1 (MQ=255)
                                                           aTCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGaaat                                                                                                                                                                                         >  1:511647/1‑250 (MQ=255)
                                                                 gCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTa                                                                                                                                                                                    >  1:271955/1‑250 (MQ=18)
                                                                     gtAGGCCGGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGGCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATac                                                                                                                                                                                <  2:420781/250‑1 (MQ=18)
                                                                              aTAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGCCGAAATCCGTAATACACTTCGTTcc                                                                                                                                                                      >  1:612003/1‑251 (MQ=18)
                                                                               tAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGCTGCGCTTTGCTTATCCGGCCTACAAAATCGCACCGTGTATGCCTGCTAAGACGCGTCAGCGTCGCACCAGGCGTTACACACCTGATGCGCCTACTTAGCGACGAAATCCGTAATCCCCTTCTTTccc                                                                                                                                                                     >  2:693596/1‑250 (MQ=18)
                                                                                aaGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCca                                                                                                                                                                     <  1:970517/250‑1 (MQ=18)
                                                                                aaGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCca                                                                                                                                                                     <  1:524073/250‑1 (MQ=18)
                                                                                       gcccgCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCGTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCcagcgcagcg                                                                                                                                                             <  2:341194/247‑1 (MQ=18)
                                                                                            cGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCt                                                                                                                                                        >  1:987664/1‑251 (MQ=18)
                                                                                                gCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTa                                                                                                                                                     >  1:475848/1‑250 (MQ=18)
                                                                                                    cAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCTTGTAGGCCTGCTCAGACTCTTCAGCGTCGCATCAGGCGTTACACACCTGCTGCGACTACATATCGACGACATCCGTAATACCCTTCGTTCCAGCGCAGCGCGGCTTTACACGc                                                                                                                                                >  2:6759/1‑251 (MQ=17)
                                                                                                     aGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGCTAAGACGAGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATACGTAATACACTTCGTTCCCGCGCAGCGCGTCTTTAAACgat                                                                                                                                               >  2:575496/1‑249 (MQ=18)
                                                                                                     aGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCt                                                                                                                                               >  1:439771/1‑251 (MQ=255)
                                                                                                      ggCGTTGAATGTCGGCTGCGCTTTGCTTACCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCTCATCA‑GGCGTTGAAGGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTg                                                                                                                                              <  2:135323/251‑1 (MQ=18)
                                                                                                            gAATGCCGGTAGCTCTTTGCTTATCCGGCCAACCAAAACGGAACGCGGAGGCCTGAAAAGACGCGCCGGCGTCGCATCA‑GGTGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAAGCCGTAATACACTGCGTTCCAGCGCAGCGCGTTTTTAAACGCTGGCAGGc                                                                                                                                        <  2:571900/251‑1 (MQ=18)
                                                                                                            gAATGCCGGATTCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGGAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAgg                                                                                                                                         <  2:799294/250‑1 (MQ=18)
                                                                                                                   ggATGCGCTTTGCTTATCCGGCCTACAAGATCGGACCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGCATGCCGGATGCGCTTTGCTTATCCGGCCTACAAACTCGCAGCGTGTAGGCCTGATAAGACGCGTCAGGGTCGCATCAGGCATTACACCCCTTATGCGACTCCTTAGCGACACAATCCGAAATACGCTTGGTTCCAGCGCCTGGCGTCTTTAAACGCTAGCAGGCgtggtgc                                                                                                                                 >  2:475929/1‑247 (MQ=9)
                                                                                                                    gATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACCCCCCTGATGCGACTACTTAGCGTCGAAATCCGTAATAAACTTCGTTCCAGCGCACCGCGTCTTTAACCGCTGGCAGGCggttatcg                                                                                                                                >  2:223981/1‑244 (MQ=17)
                                                                                                                          cTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGCCGAAATCCGTAATACACTTCGTTCCAGAGCAGCGCGTCTTTAAACGCTGGCAGGCgggtattcttatc                                                                                                                           >  1:537598/1‑238 (MQ=18)
                                                                                                                          cTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTatc                                                                                                                           >  1:372730/1‑250 (MQ=18)
                                                                                                                            ttGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGTCGCGTCAGCGTCGCATCAGACGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGCTCCAGCGCAGCGCGTCTTTAAACGCTGACAGGCggatatcgtttacaat                                                                                                                        >  2:13493/1‑236 (MQ=18)
                                                                                                                              gCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGGCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTatcaatca                                                                                                                      <  1:177175/251‑1 (MQ=21)
                                                                                                                              gCTTAGCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTatcaatca                                                                                                                      <  1:600919/251‑1 (MQ=21)
                                                                                                                               cTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTCCACACCTGATGCGCCTACTTAGCGACGACATCCGTAATACCCTTCGTTCCAGCGCAGCGCGACTTTAAACGCTGGCCGGCggggatcgttgtcactcac                                                                                                                     >  2:691235/1‑233 (MQ=17)
                                                                                                                                ttATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGAACGTTATCAATCAcc                                                                                                                    >  1:158268/1‑251 (MQ=18)
                                                                                                                                 tcaccgtccagcAAAATCGGAACGCGTAGGCCTGCTTAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATTCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTGAGCGACGAAATCCGTAATTCACTTCGTTCCAGCGCCGCGCGTCTTTAAACGCTGTCAGGCGGGTATCGTTATCAATCAc                                                                                                                     <  2:67473/238‑1 (MQ=18)
                                                                                                                                      ggCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGCCTAATTAGCGACGAAATCCGTAATAAACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCggggatcggttttaatcaacctggtt                                                                                                              >  2:878655/1‑226 (MQ=18)
                                                                                                                                       cccTACAAAATCGTTACGCGTAGGCCTGATAATGCGCGGCGGCGTCGCAGCA‑GGCGTTGAATGCCGGATGCGATTTGCTTACACGGGCTACAAAATCGCATCGTGTAGGCCTGATAAGACGCGTCAGCGTAGCATCAGACGTTACACACCTTAGGCGTCTATTTAGCGACGTAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTTGCAGGCGGGTATCGTTATCAATCTCCGTTAttt                                                                                                             <  2:835880/250‑1 (MQ=255)
                                                                                                                                           aCAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCaa                                                                                                          <  1:575701/250‑1 (MQ=21)
                                                                                                                                                      acgcgTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCAtct                                                                                              <  1:738691/250‑1 (MQ=25)
                                                                                                                                                        gcgTAGGCCTGATAAGACGCGCCAGCTTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCCCCGTGATTTCAATGTCGTGCATCTCg                                                                                            <  1:215389/251‑1 (MQ=25)
                                                                                                                                                        gcgTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCg                                                                                            <  1:172253/251‑1 (MQ=25)
                                                                                                                                                                  gATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTACTACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCATTCACCGTGATTTCAATGTCGTGCATCTCGGTGAACTGgc                                                                                  >  1:515350/1‑251 (MQ=25)
                                                                                                                                                                   aTAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGgcg                                                                                 <  1:459363/251‑1 (MQ=25)
                                                                                                                                                                   aTAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGgcg                                                                                 <  1:621931/251‑1 (MQ=25)
                                                                                                                                                                   aTAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGgcg                                                                                 <  1:47930/251‑1 (MQ=25)
                                                                                                                                                                            gcCAGCGTCGCATCA‑GGCGTTGAATGCCGGCTGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGAGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCCGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTATCAATCACCGTACTTTTCATGGCGTGCATATCGGCGGAATGGCGCATCTCGtt                                                                        >  1:756626/1‑251 (MQ=18)
                                                                                                                                                                            gcCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGtt                                                                        <  1:693596/251‑1 (MQ=37)
                                                                                                                                                                               aGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGGTCATCTCGGCGAACTGGCGCATATccttt                                                                       >  1:300319/1‑245 (MQ=25)
                                                                                                                                                                                gCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAgg                                                                    <  1:414524/251‑1 (MQ=37)
                                                                                                                                                                                ctcttcgcctcctgCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGCCGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTTTCCATCCCCGTGATTTCCATGTCGTGCATCTCGGCGAACTGGCGCATATCCTTggg                                                                     >  1:168354/14‑249 (MQ=25)
                                                                                                                                                                                  gTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTTCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTATCAATCAACGTGATTTCAATGTCGTGCATCTCGGCGAAATGGCGCCAATCGTTGAGGtt                                                                  >  1:217463/1‑251 (MQ=25)
                                                                                                                                                                                      cATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGGTATAAATCACAGTGGTTTCAATGTAGTGCATCTCGGCGAAATGGGGCATATCGGTGAGGttaagg                                                              >  1:844684/1‑247 (MQ=21)
                                                                                                                                                                                      cATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGAACGTTATCAATCACCGGGATTTTTATATCGTGCATCTCGGCCGACTGGCGCATATTGTTGTGGttaggc                                                              >  2:175604/1‑247 (MQ=21)
                                                                                                                                                                                         cA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCa                                                           <  1:929328/251‑1 (MQ=37)
                                                                                                                                                                                         cA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTATCAATCCCCGGGATTTCAATGTCGTGGATCTCCGCGAACTGGCGCATATCGTTGAGGTTCAgcgc                                                            >  1:429794/1‑247 (MQ=25)
                                                                                                                                                                                               gTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTATCCATCACCGTGATTTCAATGTCGTACATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAgcgccttgct                                                      >  1:599797/1‑242 (MQ=37)
                                                                                                                                                                                                 ttgaggCCGGATGCGCTTTGCGTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCATGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTGATCAATCACCGTGATTTCAAGGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTg                                                    <  2:440771/245‑1 (MQ=37)
                                                                                                                                                                                                   aaTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGCCGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGTTAGTTATAAAACACCGGGGTTTCAATGTTGTGTATCTCGGCGCACTGGCGCCTATCCTTTAGGTGACGTGCAT‑GGccgcc                                                  >  2:601392/1‑246 (MQ=16)
                                                                                                                                                                                                    aTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGGTCAGCGCAT‑GGCTGa                                                   >  1:461821/1‑248 (MQ=37)
                                                                                                                                                                                                      gCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACg                                              <  2:313327/251‑1 (MQ=38)
                                                                                                                                                                                                        cggatggtcTTTGATTATCCGGCCTACAAAAACGCAGCGGGTAGGCCCGATAAGACGCGGCAGCGTCGCATCAGTCGTTACACATCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTAGCGTTATCAAGCACCGTGATTTTAATGTCGTGCATCTCGGCGAACTGTCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAATACGGt                                            <  2:687803/243‑1 (MQ=255)
                                                                                                                                                                                                           attagcTTTGGTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTCCTTAGCGACGAAATTCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATg                                         <  2:327825/247‑1 (MQ=37)
                                                                                                                                                                                                             gcgcTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATgg                                        <  1:137351/250‑1 (MQ=38)
                                                                                                                                                                                                              cgcTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTGCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTg                                      <  1:218024/251‑1 (MQ=38)
                                                                                                                                                                                                               gcTTGGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTgc                                     <  2:372730/251‑1 (MQ=38)
                                                                                                                                                                                                                cTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGGTATCAATCACCGTGGTTTCCATGTCGTGCATCTCGGCGAACTGGCGCGTATCGTTGTGGTTCAGCGCAT‑GGCTGAAGACGGTATTGTgcg                                    >  1:24781/1‑251 (MQ=37)
                                                                                                                                                                                                                       tATCCGGCCTACAAAATCTCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGcca                               <  1:77961/249‑1 (MQ=255)
                                                                                                                                                                                                                        aTCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGGGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGccacca                            <  2:810264/251‑1 (MQ=255)
                                                                                                                                                                                                                            ggCCTAAAAAATCGTAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTGACACACCTGATGCGACTACTCAGCGACGAATTCCGTAATACACTTCGTTCCAGCGCATCGCGTCTTTAAACGCTGGCAGGCGGGTAGCGTTATCAATCACCGGGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGa                        <  2:285197/251‑1 (MQ=255)
                                                                                                                                                                                                                             gCCTACAAAATCGCAGAGTGTAGGCCTGATAAGACGCGTAAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGCATCGTTTTCAATCACCGTGATTTCAAATTCGTGCATCTCGGCGAACTGGGGCGTAACGTTGAGGTTCCACGCAC‑GGCTTGAGACGGTATGGTGCGCGccccccgaggg                       >  2:595543/1‑242 (MQ=255)
                                                                                                                                                                                                                              ccTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCCCCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCTCATATCGTTGAGGGTCCGAGCAT‑GGCTGCAGACGGTATGGTGCGCGCCCCCCGCGAgg                      >  2:347618/1‑251 (MQ=255)
                                                                                                                                                                                                                                   aaaaACGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCa                 <  2:135556/251‑1 (MQ=255)
                                                                                                                                                                                                                                    aaaTCGCCGCGTGTAGGCCTGATAAGACGCGACAGCGTCGCAACAGGCGTTACGCACCTGATGCGACTACGTAGCGACGAAATCCGTAATACACTTCGTTCCATCGCAGCGCGTCTTTAAACGCTTGCAGGCGGGTATCGGTGTCAATCACCGTGATTGCAAGGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGATGATCCAc                <  2:205031/251‑1 (MQ=255)
                                                                                                                                                                                                                                    aaaTCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGCTGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTTTCAATCACCGAGATTTATATGTCGGGCATCTCTGCGAACTGGCGCCTATCGTTGGGGGTCAGCGCAT‑GGCTTGAGACGGGATGGCGCGCGccccccgcgcgagtacac                >  2:758493/1‑235 (MQ=255)
                                                                                                                                                                                                                                     aaTCGCAGCGTGTAGGCCTGATAAGCCGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATTGTTATCAATCACCGTGATTTTAATGTCGTGACTCCCGTGGAACTGGCGCATATCGATGGGGGTCAGCGCAC‑GGCTGAAGACGGTATGGGgcgccacccccacgaggagccccc               >  2:512931/1‑231 (MQ=255)
                                                                                                                                                                                                                                        cGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACCCTTCGTTCCAGCGCATCGCGTCTATAAACGCTGGCAGGCTGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGATGTTCAGCGCAT‑GGCTGAAGACGGTTTGGTGCGCGCCACCAGCGAGGATCCACGCtt            <  1:584342/251‑1 (MQ=255)
                                                                                                                                                                                                                                         gCAGCGTGTAGGCCTGATAATACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCCCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTc           <  1:911505/251‑1 (MQ=255)
                                                                                                                                                                                                                                         gCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTGTCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTc           <  1:496030/251‑1 (MQ=255)
                                                                                                                                                                                                                                          cAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCg          <  1:6759/251‑1 (MQ=255)
                                                                                                                                                                                                                                           aGCGTGTATGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGACATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTATAAATAACCGTGATTTCAATGTCGGGCAACCCGGCGAACATGCGCCTATTGTTGAGGTTCAGCGACT‑GGCCGAAGAAGGACTGGTTCGCGCCACCAGCCAGGacacacgccttcg         >  2:167525/1‑239 (MQ=255)
                                                                                                                                                                                                                                           aGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGCTGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGACGGGTCTCGTTATCAATCACCGTGATTTACATGTCGTGCAACTCGGCGAACTGGAGCATGACGTTGATGGTCAG‑GCCTAGGCTGAAGAATGTATGGGTCGAGCCACACGCGAGGATCCACGccttgg         >  2:213804/1‑246 (MQ=255)
                                                                                                                                                                                                                                           aGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAGCGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGAGCCTCGCTTCgg         <  2:708176/251‑1 (MQ=255)
                                                                                                                                                                                                                                            gCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAAACACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCgg         >  1:449839/1‑250 (MQ=255)
                                                                                                                                                                                                                                                  aGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACGTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGt  <  2:632408/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
GATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTGTCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGT  >  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: ‑