Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,761,620 T→G 5.3% K641Q (AAG→CAG)  yfjK ← radiation resistance protein; DEAD/H helicase‑like protein; CP4‑57 putative defective prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,761,6200TG 5.3% 790.8 / 31.3 227K641Q (AAG→CAG) yfjKradiation resistance protein; DEAD/H helicase‑like protein; CP4‑57 putative defective prophage
Reads supporting (aligned to +/- strand):  ref base T (108/107);  new base G (6/6);  total (114/113)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.20e-03

Only 100 of 237 total aligned reads displayed.

ATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAGATATTCCGAATAACTTTCAGCACTTAAATAGCCGGCAAGTTTTGCGATCATGACATCCGTAGGCTTAACAGGGACCCCATTGCGACTGAACGTATGGGGAACCAGTC  >  NC_000913/2761350‑2761892
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                
     tGTCGGCTAATACCGAGCGCGAGCCTGATAAAAAGCTGATCAACGTACCCAATGGATCTACTCCAAATGGCCTGAGTATACCGAAGATACTGGCTCAACTCATCCACATCAGGCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCGGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCGTTACAGCgtg                                                                                                                                                                                                                                                                         <  4:366588/275‑1 (MQ=255)
                    aGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTAc                                                                                                                                                                                                                                                                               >  4:333218/1‑254 (MQ=255)
                    aGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCCCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACAtc                                                                                                                                                                                                                                                          >  2:460258/1‑275 (MQ=255)
                    aGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCACCAAGTCTACCAGTAGGTGTAACNTCTGAACGGGAATGCCAACCTATTCTAAGGTATTGACACCTGGTGGTAGAGGTAAACTGTGTAACCCCAATTTGACAGGTNGGAAAACGACTTGTCTACGAATACTNCCCGTTACAGCGTGTTTTTCAAcccaacc                                                                                                                                                                                                                                                          >  3:94373/1‑269 (MQ=255)
                      cgcgcgcCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGCATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATTGAAACTTTCGAACGCAAATTTGACATGTGTGTAATCGACTTTGCTGCGTTTACCCCGCTTTAAAGCGTGGTGTTTAACTACAtcct                                                                                                                                                                                                                                                        >  4:487460/1‑275 (MQ=255)
                      cgcgcgcCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTGTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACAtcct                                                                                                                                                                                                                                                        >  4:503306/1‑275 (MQ=255)
                        cgcgcCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCg                                                                                                                                                                                                                                                                           <  1:51960/254‑1 (MQ=255)
                        cgcgcCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCg                                                                                                                                                                                                                                                                           <  2:333219/254‑1 (MQ=255)
                                    aCATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAg                                                                                                                                                                                                                                           <  2:487461/274‑1 (MQ=255)
                                                  tACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGt                                                                                                                                                                                                                            >  2:313581/1‑275 (MQ=255)
                                                   gcccAACCGATCTATCCCAAAGGTCCTTAGTAACCGGAAGATGCCTGCTACACTGATCCACATAAGTCTCGTCGGAAAATTCCAGCAAGTCGACCAGTCGGTGTAACGTTTGAAGGGGAATGCCAATCTCTTATAAGGCAGTGACACCTGGTGGTAGATGGAAACTGTCGAACGCCAATTTGACATGTGTGTAATCTACTTTGCTACGTATACCACGTTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGta                                                                                                                                                                                                                           <  4:358506/274‑1 (MQ=255)
                                                     ccAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTa                                                                                                                                                                                                                                                              >  2:160203/1‑238 (MQ=255)
                                                         tCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACAtc                                                                                                                                                                                                                                                          <  4:160202/238‑1 (MQ=255)
                                                         actagtcaCTCCAAATGTCCTTGGTATCCCGAAGATACTGGCTCAACTCATCCACATCCGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCGCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAAAGACATACTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCg                                                                                                                                                                                                                     <  4:105021/268‑1 (MQ=255)
                                                          cGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGTCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGGCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATCTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCg                                                                                                                                                                                                                     <  2:274123/274‑1 (MQ=255)
                                                          agATCTACTCCGCATGGCCCGGGCATCCGGGAGAGAATGGGTCAACTCGGTCGCATAAACTTCGACAGGAAAATCAGGCAAGTCTACCAGTCGGTTTAACGCATGAATGGGAAGGCCAATATCTTATAAGGGATTGACACCTGGTGGTAGATGGAAACTTTCGAAGGCCAAGTTGACATGTGGGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGGGGTGTTTAACTACATCCTCCCTCAAGCTGAGCACCTTAGGTAGGGTATAGGCg                                                                                                                                                                                                                     <  4:386088/273‑1 (MQ=255)
                                                                       aTGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTAt                                                                                                                                                                                                              <  1:71959/268‑1 (MQ=255)
                                                                                    cccGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATcc                                                                                                                                                                                          >  4:419883/1‑274 (MQ=255)
                                                                                        aaGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGa                                                                                                                                                                                                                                            <  4:212092/221‑1 (MQ=255)
                                                                                         aGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGCGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCACCATCAACCTGAGCACCTTAGGCAGGGTATAGCCGAAGGTATTGGTTATGAGCTTGAAATCAttgtt                                                                                                                                                                                     >  2:9909/1‑275 (MQ=255)
                                                                                             aCTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCTCGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTAAAATCAttgttggag                                                                                                                                                                                 >  4:284052/1‑272 (MQ=255)
                                                                                                 gCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCaa                                                                                                                                                                             >  4:318638/1‑275 (MQ=255)
                                                                                                      aCTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGa                                                                                                                                                                         >  2:536003/1‑274 (MQ=255)
                                                                                                                                       gCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCttt                                                                                                                                                                                                                                                                                 <  6:157291/137‑1 (MQ=255)
                                                                                                                                       gCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTGTAAAGCgt                                                                                                                                                                                                                                                                          >  8:543671/1‑144 (MQ=255)
                                                                                                                                       gCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTGTAAAGCgt                                                                                                                                                                                                                                                                          <  7:543671/144‑1 (MQ=255)
                                                                                                                                            tCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGc                                                                                                                                  >  1:109214/1‑275 (MQ=255)
                                                                                                                                                ccAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTGTAAAGCgt                                                                                                                                                                                                                                                                          <  7:1046095/135‑1 (MQ=255)
                                                                                                                                                ccAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTGTAAAGCgt                                                                                                                                                                                                                                                                          >  8:1046095/1‑135 (MQ=255)
                                                                                                                                                    tCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAg                                                                                                                                                                                                                                                                             <  8:169877/128‑1 (MQ=255)
                                                                                                                                                      ggTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAggg                                                                                                                                                                                                                             >  3:18778/1‑174 (MQ=255)
                                                                                                                                                            aCGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCaa                                                                                                                                                                                                                                                 >  7:453716/1‑148 (MQ=255)
                                                                                                                                                                  gAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCCGa                                                                                                            >  4:156896/1‑275 (MQ=255)
                                                                                                                                                                   aaTTGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGGATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAg                                                                                                             <  1:99504/273‑1 (MQ=255)
                                                                                                                                                                      gggAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAg                                                                                                             <  4:409658/270‑1 (MQ=255)
                                                                                                                                                                         aaTGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCt                                                                                                                                                                                                                                      <  8:612734/146‑1 (MQ=255)
                                                                                                                                                                          aTGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTGTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAGATATTCCGa                                                                                                    <  2:503307/275‑1 (MQ=255)
                                                                                                                                                                             ccAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAGATATTCCGAATa                                                                                                 >  2:339633/1‑275 (MQ=255)
                                                                                                                                                                               aaTCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtaggtagg                                                                                                                                                                                                                              >  7:964532/1‑148 (MQ=255)
                                                                                                                                                                                aTCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgt                                                                                                                                                                                                                                                                          >  8:199910/1‑103 (MQ=255)
                                                                                                                                                                                aTCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgt                                                                                                                                                                                                                                                                          <  7:199910/103‑1 (MQ=255)
                                                                                                                                                                                aTCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtaggtagg                                                                                                                                                                                                                              <  7:2715/147‑1 (MQ=255)
                                                                                                                                                                                 tctcTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGgtgt                                                                                                                                                                                                                                                                     <  6:1030786/107‑1 (MQ=255)
                                                                                                                                                                                      tctAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTa                                                                                                                                                                                                                                                                  <  5:1059337/105‑1 (MQ=255)
                                                                                                                                                                                      tctAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTa                                                                                                                                                                                                                                                                  >  6:1059337/1‑105 (MQ=255)
                                                                                                                                                                                          aGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATc                                                                                                                        <  1:151160/239‑1 (MQ=255)
                                                                                                                                                                                             cATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTta                                                                                                                                                                                                                                    >  8:33325/1‑128 (MQ=255)
                                                                                                                                                                                               ttGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTa                                                                                                                                                                                                                                                                                <  7:745700/82‑1 (MQ=255)
                                                                                                                                                                                               ttGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTAc                                                                                                                                                                                                                                                                               <  8:629343/83‑1 (MQ=255)
                                                                                                                                                                                               ttGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTTCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGACGGCAATGGAAAGAGTACGTTTCTACCCTTCAGAAATCCACTTGAGTTAATAATCGATCTGATTTTTCAGATTTTCCGAATACCTTCAAGCACTTAAAt                                                                                >  4:475200/1‑274 (MQ=255)
                                                                                                                                                                                               ttGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGaaa                                                                                                                                               >  2:557902/1‑211 (MQ=255)
                                                                                                                                                                                               ttGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGaaa                                                                                                                                               >  2:358769/1‑211 (MQ=255)
                                                                                                                                                                                                 gACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCTACTTTGCTGCGTATACCACGCTTTACAGCGTGCTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTg                                                                                                                                                                                                            <  8:839117/148‑1 (MQ=255)
                                                                                                                                                                                                  acacCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGa                                                                                                                                                                                                                    >  7:165289/1‑139 (MQ=255)
                                                                                                                                                                                                   canCTGGNGGTAGATGGANACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCa                                                                                                                                           <  4:557901/211‑1 (MQ=255)
                                                                                                                                                                                                   cacCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGta                                                                                                                                                                                                                           <  2:239066/131‑1 (MQ=255)
                                                                                                                                                                                                    acCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAGATATTCCGAATAACTTTCAGCACTTAAATAGCCgg                                                                          <  2:284053/275‑1 (MQ=255)
                                                                                                                                                                                                     ccTGTTGTTAGATGGAAACTTTCTAAATCNNTTTTGACAATGGTGTAATCGACTTTGCTGTNTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGANGGTATTGGTTATGATCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAGATATTCCGAATAACTTTCAGCACTTAAATAGCCGGc                                                                         <  3:69850/275‑1 (MQ=255)
                                                                                                                                                                                                     ccTGGTGGTAGATGGAAATTTTCGAACGCCAATTTGAGATGTGTGTTATGGACTTTTCTGCGTATTCCCCTCTTTTTAGCGGGGTGTTTAACTACATCCTCCATCAAACGTAGCACCTTAGGTAGGGTTTAGCCGAAGGTATTGGTCATGAGCTTCAAATCAATGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATGTTTCAGATATTCCGAATAACTATCAGCACGTAAATAGCCGGc                                                                         <  4:9909/275‑1 (MQ=255)
                                                                                                                                                                                                     ccTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTa                                                                                                                                                                                                        >  6:245960/1‑148 (MQ=255)
                                                                                                                                                                                                            gtAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCGTAGGTAGGGTATAGCCGAAGGAATGGGTTATGAGCTTCAAATCAttgtt                                                                                                                                                                                     >  2:29385/1‑160 (MQ=255)
                                                                                                                                                                                                               gATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGaa                                                                                                                                                                                                                   >  7:615018/1‑127 (MQ=255)
                                                                                                                                                                                                               gATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCa                                                                                                                                                                                               >  7:570418/1‑147 (MQ=255)
                                                                                                                                                                                                                ttGGAAACTTTGGAACGCCAATTTAGCATGTGTGTAATCGACTTTGCTTCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAgg                                                                                                                                                                                 <  4:29384/159‑1 (MQ=255)
                                                                                                                                                                                                                   gAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTaa                                                                                                                                                                                                                                                                 >  5:633019/1‑77 (MQ=255)
                                                                                                                                                                                                                   gAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTaa                                                                                                                                                                                                                                                                 <  6:633019/77‑1 (MQ=255)
                                                                                                                                                                                                                   gAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTta                                                                                                                                                                                                                                    >  5:87805/1‑106 (MQ=255)
                                                                                                                                                                                                                   gAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTta                                                                                                                                                                                                                                    <  6:87805/106‑1 (MQ=255)
                                                                                                                                                                                                                   gAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTAtt                                                                                                                                                                                                             >  6:484473/1‑129 (MQ=255)
                                                                                                                                                                                                                     aaCTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtaggtag                                                                                                                                                                                                                               <  8:456341/109‑1 (MQ=255)
                                                                                                                                                                                                                     aaCTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtaggt                                                                                                                                                                                                                                 <  8:947727/107‑1 (MQ=255)
                                                                                                                                                                                                                                     tttGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTAc                                                                                                                                                                                                                                                                               >  5:115929/1‑45 (MQ=255)
                                                                                                                                                                                                                                          cATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgtgg                                                                                                                                                                                                                                                                        >  8:219630/1‑47 (MQ=255)
                                                                                                                                                                                                                                          cATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgtgg                                                                                                                                                                                                                                                                        <  7:219630/47‑1 (MQ=255)
                                                                                                                                                                                                                                            tgtgtgTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAGATATTCCGAATAACTTTCAGCACTTAAATAGCCGGCAAGTTTTGCGATCATGACATCCGTAGGCTTAACAGGGAc                                  >  2:202412/1‑275 (MQ=255)
                                                                                                                                                                                                                                               gtgtAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtaggta                                                                                                                                                                                                                                >  8:442785/1‑82 (MQ=255)
                                                                                                                                                                                                                                               gtgtAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAGATATTCCGAATAACTTTCAGCACTTAAATAGCCGGCAAGTTTTTCGATCATGACa                                                      >  1:23051/1‑252 (MQ=255)
                                                                                                                                                                                                                                                tgtAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAGATATTCCGAATAACTTTCAGCACTTAAATCGCCGGCAAGTTTTGCGATCATGACTTCCGTAGCCTTAACAGGGAccccc                              >  4:257916/1‑274 (MQ=255)
                                                                                                                                                                                                                                                  tAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGttt                                                                                                                                                                                                                                                                   >  5:439097/1‑44 (MQ=255)
                                                                                                                                                                                                                                                  tAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGttt                                                                                                                                                                                                                                                                   <  6:439097/44‑1 (MQ=255)
                                                                                                                                                                                                                                                  tAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATATAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAGATATTCCGAATAACTTTCAGCACTTAAATAGCCGGCAAGTTTTGCGATCATGACATCCGTAGGCTTAACAGGGACCCCAtt                            <  2:374187/275‑1 (MQ=255)
                                                                                                                                                                                                                                                   aaTCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGtat                                                                                                                                                                                                                          >  5:357254/1‑84 (MQ=255)
                                                                                                                                                                                                                                                   aaTCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGtat                                                                                                                                                                                                                          <  6:357254/84‑1 (MQ=255)
                                                                                                                                                                                                                                                   aaTCGACTTTGCTGCGTATACCACGCTTTACAACGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtaggta                                                                                                                                                                                                                                <  6:955263/78‑1 (MQ=255)
                                                                                                                                                                                                                                                      cGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAGATATTCCGAATAACTTTCAGCACTTAAATAGCCGGCAAGTTTTGCGATCATGACATCCGTAGGCTTAACAGGGACCCCATTGCGa                        <  2:181294/275‑1 (MQ=255)
                                                                                                                                                                                                                                                          tttGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAGATATTCCGAATAACTTTCAGCACTTAAATAGCCGGCAAGTTTTGCGATCATGACATCCGTAGGCTTAACAGGGACCCCATTGCGACTGa                    >  2:336725/1‑275 (MQ=255)
                                                                                                                                                                                                                                                           ttGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGa                                                                                                                                                                                                                                            <  7:357326/58‑1 (MQ=255)
                                                                                                                                                                                                                                                           ttGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGa                                                                                                                                                                                                                                            >  8:357326/1‑58 (MQ=255)
                                                                                                                                                                                                                                                           ttGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAg                                                                                                                                                                                                                  <  8:675317/84‑1 (MQ=255)
                                                                                                                                                                                                                                                           ttGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAg                                                                                                                                                                                                                  >  7:675317/1‑84 (MQ=255)
                                                                                                                                                                                                                                                                 cGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAGATATTCCGAATAACTTTCAGCACTTAAATAGCCGGCAAGTTTTGCGATCATGACATCCGTAGGCTTAACAGGGACCCCATTGCGACTGAACGTATg             <  4:462524/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                       ccACGCTGTAAAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCt                                                                                                                                                      >  7:900703/1‑132 (MQ=255)
                                                                                                                                                                                                                                                                         aCGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGa                                                                                                                                                                                   <  5:17648/101‑1 (MQ=255)
                                                                                                                                                                                                                                                                         aCGCTGTAAAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGtt                                                                                                                                                              >  7:387356/1‑122 (MQ=255)
                                                                                                                                                                                                                                                                          cGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCAt                                                                                                                                                                                         <  5:96423/94‑1 (MQ=255)
                                                                                                                                                                                                                                                                           gCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAGATATTCCGAATAACTTTCAGCACTTAAATAGCCGGCAAGTTTTGCGATCATGACATCCGTTGGCTTAACAGGGACCCCATTCCGACTGAACGTATGGGGAAcc      >  4:163643/1‑272 (MQ=255)
                                                                                                                                                                                                                                                                            cTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCg                                                                                                                                     >  5:54781/1‑144 (MQ=255)
                                                                                                                                                                                                                                                                             tttACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAGATATTCCGAATAACTTTCAGCACTTAAATAGCCGGCAAGTTTTGCGATCATGACATCCGTAGGCTTAACAGGGACCCCATTGCGACTGAACGTATGGGGAACCAGTc  >  4:338170/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                             tttACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAGATATTCCGAATAACTTTCAGCACTTAAATAGCCGGCAAGTTTTGCGATCATGACATCCGTAGGCTTAACAGGGACCCCATTGCGACTGAACGTATGGGGAACCAGTc  <  2:475201/274‑1 (MQ=255)
                                                                                                                                                                                                                                                                             tttACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAGATATTCCGAATAACTTTCAGCACTTAAATAGCCGGCAAGTTTTGCGATCATGACATCCGTAGGCTTAACAGGGACCCCATTGCGACTGAACGTATGGGGAACCAGTc  <  4:202411/274‑1 (MQ=255)
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                
ATCAATGTCGTCTAATACCGAGCGCGCGCCTGATAAACATCTGATCAACGTACCCAATCGATCTACTCCAAATGTCCTGAGTATCCCGAAGATACTGGCTCAACTCATCCACATCAGCCTCGTCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATAATCGATCTGATTTTTCAGATATTCCGAATAACTTTCAGCACTTAAATAGCCGGCAAGTTTTGCGATCATGACATCCGTAGGCTTAACAGGGACCCCATTGCGACTGAACGTATGGGGAACCAGTC  >  NC_000913/2761350‑2761892

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: