Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 72,892 G→T H12Q (CAC→CAA thiQ ← thiamine ABC transporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091372,8920GT97.5% 342.2 / ‑3.4 121H12Q (CAC→CAAthiQthiamine ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base T (59/58);  total (61/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

Only 100 of 121 total aligned reads displayed.

CCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAG  >  NC_000913/72649‑73138
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      
cccGATGTTCTGCGCGACCGTCAGGTGACTGCACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACTGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGAAAATGTTGGTaa                                                                                                                                                                                                                                                  >  1:332812/1‑250 (MQ=255)
   gATGTTCTGCGCGACCGTCAGGTGACTGAACAGGATGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCACGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAAAGCATCGGCAAATGTTGGTAAAGcc                                                                                                                                                                                                                                              >  2:819980/1‑251 (MQ=255)
      gttccgcgCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCTACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAg                                                                                                                                                                                                                                            <  2:646560/243‑1 (MQ=255)
      gTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGGGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAgg                                                                                                                                                                                                                                           <  2:368970/251‑1 (MQ=255)
       ttCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGCTCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAggg                                                                                                                                                                                                                                          >  1:915478/1‑250 (MQ=255)
          tgcgCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGat                                                                                                                                                                                                                                       >  1:382230/4‑251 (MQ=255)
            cgCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGACGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGatat                                                                                                                                                                                                                                     <  2:746210/250‑1 (MQ=255)
               gACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAATGCATCGGCAAATGTTGGTAAAGCCAGGTGATAtc                                                                                                                                                                                                                                    >  1:32640/1‑249 (MQ=255)
               gACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATAtcag                                                                                                                                                                                                                                  <  1:735009/251‑1 (MQ=255)
                       ggggacTGAACAGGTAGTTTACCCGCAACAGCATCGTCACCAGACGGCGTGACGGCTGCGACGTGGTGCGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCTCCACCTGCTCGCCGCGTGCTACCTTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCCGtt                                                                                                                                                                                                                            <  2:281128/244‑1 (MQ=255)
                       ggTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGCCGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGtttt                                                                                                                                                                                                                          >  1:507172/1‑251 (MQ=255)
                        gtgacTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGtttt                                                                                                                                                                                                                          >  1:76012/5‑250 (MQ=255)
                           acTGACCAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACa                                                                                                                                                                                                                      <  2:170745/250‑1 (MQ=255)
                           acTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGCCGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCTCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAAAGTTGGTAAAGACAGGTGATATCAGACAGTTTTaaa                                                                                                                                                                                                                       >  2:27655/2‑249 (MQ=255)
                              gAACAGGTTGTTTTCACGAACAAGCATCCACACCGGACGGCGGGACGGCGGCGTCGTGGGGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTGCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCTCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGAGATCAGTCAGTTTTAACAtc                                                                                                                                                                                                                    <  2:45213/250‑2 (MQ=255)
                               aaCAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACAtctc                                                                                                                                                                                                                  <  1:981778/251‑4 (MQ=255)
                                    gttgttTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGGCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGcccc                                                                                                                                                                                                             <  1:942638/251‑1 (MQ=255)
                                     ttgttTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCg                                                                                                                                                                                                            >  1:824724/1‑251 (MQ=255)
                                         tttcCTGAAACAGCATCGACACCGGAAGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAACAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTGAAGCCAGGACATATCAGTCAGTTTTAACATCTCCCCCCGGt                                                                                                                                                                                                          >  2:629233/3‑249 (MQ=255)
                                           tcCTGAAACAGCATCGACGCCGGACGGCGTTACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGt                                                                                                                                                                                                       <  2:331123/250‑1 (MQ=255)
                                               gAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGttttt                                                                                                                                                                                                   <  1:154573/250‑1 (MQ=255)
                                               gAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTc                                                                                                                                                                                                  >  1:328612/1‑251 (MQ=255)
                                               gAAACAGCATCGACACCGGACGGCGTGACGGCGGCGGCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTc                                                                                                                                                                                                  <  1:511756/251‑1 (MQ=255)
                                                    aGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCa                                                                                                                                                                                             <  1:19534/251‑1 (MQ=255)
                                                      cATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACg                                                                                                                                                                                           <  1:857216/251‑1 (MQ=255)
                                                        ccGACACCGGACGGCGTGACGGCGTCGTCGTGGTGTGATCAACGCCGGCGATAGTCAGCGAACTGTTGGCCGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTGTTACCTGCGCCGTTTGGCCCGAGGATCTCCACCTGCTCGCCGCGTTCTGCCGTTAAGCTAAAACGCATCTGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTGTCAATCACGGt                                                                                                                                                                                         <  2:983528/250‑1 (MQ=255)
                                                          gACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCCCCTGCTCGCCGCGTTCTACCGTTACGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGAAAACAGTCAGATTTAACATCTCGCCCCGGTATGTTTTCattacaggaa                                                                                                                                                                                        >  2:245499/1‑241 (MQ=255)
                                                             acCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAAc                                                                                                                                                                                     >  1:257442/1‑250 (MQ=255)
                                                               cGGACGGCGTGACTGCGGCGTCGTGGTGTGATCAACGCCGTCGTTAGTCAGCGAGCCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGTTCTCCACCTGCTCGCCTCGTTCTACCGTTAAGCTAAAACGCATCGGCATATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACAGCTCGTCCCGGTAGTTTTTCAATCACGGTAAAcagc                                                                                                                                                                                  <  2:596178/251‑1 (MQ=255)
                                                                ggACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAAcagca                                                                                                                                                                                 <  1:90919/251‑1 (MQ=255)
                                                                ggACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAAcagc                                                                                                                                                                                  >  1:419630/1‑250 (MQ=255)
                                                                   cGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGAAAACAGCAGa                                                                                                                                                                               >  1:845781/1‑250 (MQ=255)
                                                                    ggCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGaa                                                                                                                                                                              >  1:321959/1‑250 (MQ=255)
                                                                    ggCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGAGATATCAGTCAGTTTTAAACTCTCGCCCCGGTAGTTTTTCAATCACggataaaagcagca                                                                                                                                                                              >  2:181801/1‑238 (MQ=255)
                                                                                       ggtgtgATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAAt                                                                                                                                                           <  1:114224/249‑1 (MQ=255)
                                                                                          tggatcaacgctgtggATAGTCAGCGAACCGCTGGCTGGCGGCAGAAAACCGACGATCAAGTTCAGCAGGGGACTTGTCCCTGCGCCGATTGGCCCTAGGATCGCCTCCTGCTCGCCGCGTGCTACCGTTAAGCTAAAACGCATCGTCAAATGTTTGTAAAGCCAGGTGATATCAGTCAGGTTTAACATCTCTCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTaa                                                                                                                                                        <  2:889146/235‑1 (MQ=255)
                                                                                              tcaaCGCCGTCGATAGTCAGAGAACCGCTGGCTGGCGTCAGAAAACCGGCGATAAAGTTCAGCAGGGAACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGACAGGGGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCACTCACGGAATACAGCAGAAAACAGTGTAGCAGAAGAATTAAAGCgg                                                                                                                                                   >  2:463193/4‑251 (MQ=255)
                                                                                                aaCGCCGTCGATAGTCGGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTg                                                                                                                                                 >  1:214994/2‑251 (MQ=255)
                                                                                                aaCGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGGTCAGCAGGGTACTTTTACCTGCGCCGCGTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTg                                                                                                                                                 <  2:741318/250‑1 (MQ=255)
                                                                                                   gccgtCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGCTATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGAcc                                                                                                                                              >  2:71854/5‑251 (MQ=255)
                                                                                                        cgcaaGTCACCGCACCGCTGGCTTGCGACAGAAACCCTGCGATCAAGTTCAGCAGGGTTCTTTTACCTGCCCCGCTTGGCCCGAGGATCGCCCCCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCAGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTTGCAGCAGAATTAACGCGGTGCCCGcgcc                                                                                                                                         <  2:865158/247‑4 (MQ=255)
                                                                                                        cGATAGCCAGCGAACCGCGGGCTGGCGTCAGAAAACCTGCGATCAAGTTCAGCAGGGTACTTGTCCCTGCGCCGCTTGGCCCGAGTATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCTCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGcgcc                                                                                                                                         <  1:545862/251‑4 (MQ=255)
                                                                                                         gATAGACAGCGCACCGATGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGGACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATAACGGTAAACAGCCGAAACCAGAATAGCAGCAGAATTAACGCGGTTACCGcgccc                                                                                                                                        >  2:491606/1‑247 (MQ=255)
                                                                                                            aGTCAGCGAACGGCTGGCTGGCGTCAGAAAACCGGCGAACAAGTTCAGCAGGGTACTTTTCCCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTc                                                                                                                                      <  2:660944/250‑1 (MQ=255)
                                                                                                                 gCGAACCGCTGGCTGGAGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTTCTTTTACCTGCGCCGCTTGGCCAGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCGCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTTCTGGCt                                                                                                                                <  1:384796/251‑1 (MQ=255)
                                                                                                                  cGAACCGCAGGCGGGCGGCAGAAAACCGGCGATCAAGGTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCAGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTg                                                                                                                               <  2:139372/251‑1 (MQ=255)
                                                                                                                       cggtggccggcGTCAGAAAACCGGCGAGCAAGTGCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTGCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATa                                                                                                                          <  2:672346/248‑1 (MQ=255)
                                                                                                                        gctggctgGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAg                                                                                                                         <  1:712461/251‑1 (MQ=255)
                                                                                                                         ctggctggCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAgg                                                                                                                        >  1:122646/1‑251 (MQ=255)
                                                                                                                           tgctggcGTCAGAAAACCGGCGATCAAGTTCAGCATGGTACTTTTCCCTGCTCCGCTTTGCCCGAGGATCGCCACCTGCTCGCTGCGTTCTACCGTTAAGCTAAAACTCATCGGCAAATGTTGGTAAAGCCAGGGGATATCAGTCAGTTTTAACATCGCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAAGTAACGCGGTGACCGCGCCGTCCTGGCTGCGATatga                                                                                                                       <  2:599934/249‑4 (MQ=255)
                                                                                                                           ggctggcGTCAGAAAACCGGCGATCAAGTTCAGCAGGTTACTTTTACCTTCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTTAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAg                                                                                                                      <  2:896566/251‑1 (MQ=255)
                                                                                                                               ggcGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCa                                                                                                                   <  1:910172/250‑1 (MQ=255)
                                                                                                                                gcGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTCCTTTTACCTGCGCCGCTTGGCCGGAGGATCGCCGCCGGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAAt                                                                                                                 <  1:703931/251‑1 (MQ=255)
                                                                                                                                      gAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGGATTTTTTCAATCACGTTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTg                                                                                                           <  1:432302/251‑1 (MQ=255)
                                                                                                                                         aaCCGGCGAGCAAGTTCAGCAGGGTACTGTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTa                                                                                                        <  1:449732/251‑1 (MQ=255)
                                                                                                                                           ccggagtgcaatttCAGCAGGTTACTTTTCCCTGCGTCGCTTGGACCGAGGATCGCCACCTGCTCGCCGCGTTCTAACGTTAAGCTAAACCGCATCGGCAAATGTTGGTAAAGCCATGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAACCAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTaga                                                                                                      <  2:277324/239‑1 (MQ=255)
                                                                                                                                             ggCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAgaga                                                                                                    <  2:382230/251‑1 (MQ=255)
                                                                                                                                                    aagttCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGATGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACg                                                                                             >  1:358735/4‑251 (MQ=255)
                                                                                                                                                         cagcagGGTACTTTTACCTGTGCCGCTTGGCCCGAGGATCGCCCCCTGCTCGCCGTGTTCGACCGTGAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGTTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTATCGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAg                                                                                         <  2:751007/250‑1 (MQ=255)
                                                                                                                                                         cagcagGGTACTTTTACCTGCGCCGCTTGGCACGAGGATCGCCATCTGCTCGCCGCGATCTTCCGTTAAGCTAAAACGCATCGGCAAATTTTGGTAAAGCCAGGGGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAgg                                                                                        <  1:686216/251‑1 (MQ=255)
                                                                                                                                                             ggggtccttctccctgcgcCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAACCAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGttc                                                                                    <  2:242983/236‑3 (MQ=255)
                                                                                                                                                             agGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACCGAGTAGCAGCAGAATTAACGAGGTGACCGCGCCGTCCTGGCTGCGGTAGGAGCCAATTTGCTGGGAGAGATAAAACGGCAGTGTGc                                                                                    >  2:569243/1‑251 (MQ=255)
                                                                                                                                                                  aCTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGaaa                                                                               <  1:93437/251‑1 (MQ=255)
                                                                                                                                                                   cTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACAGTTAAGCTAATACGCATCGGCAAATGTTGGTAAAGCCAGGCGATAACAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGGGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGCTAGGAGCCAATTAGCTGGTAGAGATAAAACGGCAGGGGGCGGaca                                                                               >  2:379146/1‑248 (MQ=255)
                                                                                                                                                                   cTTTTACCTGCGCCGCTCGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAACTGTTGGTAAAGCAAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGAACGCGCGGCCCTGGCTGCGAAAGGACCCAATTTGCTGGATGAGATAAAACGGCAGGGTGCgggaa                                                                               >  2:504572/1‑247 (MQ=255)
                                                                                                                                                                     tttACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAatc                                                                             >  1:782198/1‑250 (MQ=255)
                                                                                                                                                                     tttACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGAAAACAGCAGAAACCAGAGTAGCCGCAGAATTCACGCGGTGCCCGCGCCGTCCTAGCTGCGAGAGGAGCCAATCTGCTGGTAGAGAGAAAATGGCAGGGCGCGCAAata                                                                             >  2:776457/1‑249 (MQ=255)
                                                                                                                                                                      ttACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAatca                                                                            <  2:289102/250‑1 (MQ=255)
                                                                                                                                                                        acctgcgcCGCTTGGCCCGAGGATCGCCACCTGCACGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCACGTGATATCAGTCAGTTTTAACATCTAGCCCCGGTAGTTTTTCAAACACGGTAAACAGCAGTAAACAGCGTAGCAGCAGAATCAACGAGGTGACCGCACCGACTTTGCTGCGATAGGAGCCAATTGGATGGTAGAGATAAAACGGCAGAGGGCGGAAatctt                                                                           >  2:614265/5‑247 (MQ=255)
                                                                                                                                                                         cctgcgcCGCTTGGCCCGAGGATCGCCACCTGCTAGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGt                                                                        >  1:109604/4‑251 (MQ=255)
                                                                                                                                                                                 gcTTGGACCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCAAATTTGCTGGAAGAGATAAAACGGCAGGGGGCGGGAATCATCGTTACAGaa                                                                 >  1:522693/1‑250 (MQ=255)
                                                                                                                                                                                    gggCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACACCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGTAGCCAATCTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCAGCGTTACCGaacaac                                                             <  1:133372/250‑1 (MQ=255)
                                                                                                                                                                                     ggCCCGAGGATCGCCACCGGCGCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGTAAATCAGCGTTACCGAACAACg                                                            <  2:147182/251‑1 (MQ=255)
                                                                                                                                                                                         cGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCCGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGGTAGGAGCCAATTTGCTGGTAGAGATAACACGGCAGGGGGCGGAAATCATCGTTACCGAACAACGcccc                                                        >  2:784634/1‑249 (MQ=255)
                                                                                                                                                                                          gctgatcggtcacTGCTCACCGCCTTCGACCGTTAAGCTAAAACGCACCGGCAAACGTTGGTAAAGCCAGTTGATATCAGGCAGTTTTACAATCTCCGCCCGGTAGTTTTTTAATCACTGTAAACAGCTGAAAACAGAGGGGCAGCAGAATTAACGCGGTGACCGCGCCGTTCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATACAACGGCAGTGTGCGGATATCATCGTTACCGACCAACGcaatt                                                       <  2:671123/239‑5 (MQ=255)
                                                                                                                                                                                          gAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGAGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTACACAGCAGAAAACAGAGTAGAAGCAGACTTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGccacc                                                       >  1:887249/1‑251 (MQ=255)
                                                                                                                                                                                           aGGATCGCCACCAGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGccacca                                                      <  2:591522/251‑1 (MQ=255)
                                                                                                                                                                                            ggatggccccctgcatgccgcGTTCGACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCcaccac                                                     <  2:19153/235‑1 (MQ=255)
                                                                                                                                                                                            ggATCGCCCCCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCcaccac                                                     <  2:702408/251‑1 (MQ=255)
                                                                                                                                                                                            ggATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAACACGCATCGGCCAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCcaccac                                                     >  1:374965/1‑251 (MQ=255)
                                                                                                                                                                                                  ccACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCaaa                                               <  1:697980/251‑1 (MQ=255)
                                                                                                                                                                                                    aCCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAAt                                             <  1:48590/251‑1 (MQ=255)
                                                                                                                                                                                                        gCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGGGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCAc                                          >  2:703281/1‑250 (MQ=255)
                                                                                                                                                                                                               ccgttctcacGTTAAGCTAAAACTCATCCGCAAACGTTGGTAAATCCCGGTGATATCAGCTAGTTTTAACATCTCGCCCCGGTAGGTTTTCAATCACGGTAAACAGCAGAACACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCTCGCCAAAATCACCAATCGAc                                  <  2:754285/242‑1 (MQ=255)
                                                                                                                                                                                                                 gTTCTACCGTAAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGAc                                  >  1:187090/1‑249 (MQ=255)
                                                                                                                                                                                                                   tctaCCGTTAAGCTAAAACGCCTCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTAGCCCCGGTAGTTTTTCAATACCGGGAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGGGCGGGAATCATCGTTACCGAACAACGCAACCACGCCAAAATCACCAATCGACAGcc                              >  2:822751/4‑250 (MQ=255)
                                                                                                                                                                                                                        cGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGagtagca                                                                                                                                                                  >  2:107079/1‑114 (MQ=255)
                                                                                                                                                                                                                        cGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGagtagca                                                                                                                                                                  <  1:107079/114‑1 (MQ=255)
                                                                                                                                                                                                                         gTTAAGCTAAAACGCATCGGCAAATATTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGAAGTTTTTCAATCACGGTAAACAGCAGAAAACAGCGTAGCAGCAGAATTAACGCGGTGACCTCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCTGAGTGCGAGAATACTCGTTACCGAACAACGCCACCACGACAAAAACAACAATCGaatccacgcagg                        >  2:663097/1‑240 (MQ=255)
                                                                                                                                                                                                                             aGCTAAAACGCATCGGCAAATGTGGGTAAAGCCAGGTGATATCAGTCAGGTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCaaa                    <  2:311053/251‑1 (MQ=255)
                                                                                                                                                                                                                               cTAAAACGCATCGGCAAATGTTGGTAAAGCCCGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGccggcaaa                    >  1:756129/1‑242 (MQ=255)
                                                                                                                                                                                                                                tAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATATGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCaaaagc                 >  1:887357/1‑248 (MQ=255)
                                                                                                                                                                                                                                tAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGAAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCACTCGACAGCACGCAGGCaaaagc                 >  1:366806/1‑248 (MQ=255)
                                                                                                                                                                                                                                 aaaaCGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGcc                >  1:752537/1‑251 (MQ=255)
                                                                                                                                                                                                                                  aaaCGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCCGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCaaaaaccc               >  1:527090/1‑246 (MQ=255)
                                                                                                                                                                                                                                  aaaCGCATCGGCAAATGTTGGGAAAGCCAGGTGATATCAGTCAGGTTTAACATCTCGCCCCGGTAGTTTGTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCa               <  2:195512/251‑1 (MQ=255)
                                                                                                                                                                                                                                   aaCGCATCGGCAAATGTTGGTAAAGACAGGTGATATCAGTCAGTTTTAACATCTCGCCCAGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCATCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGTTAAAACGGCAGGGTGAGGAAATCATCGTTACCGAACAAAGCCACCACGACAAAATCACCAAACGAAAACACGccggaaaaagcca               >  2:342550/1‑238 (MQ=255)
                                                                                                                                                                                                                                       cATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCg                                                                                                                                             >  1:308551/1‑120 (MQ=255)
                                                                                                                                                                                                                                       cATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCg                                                                                                                                             <  2:308551/120‑1 (MQ=255)
                                                                                                                                                                                                                                               aaTGTTGGTAAAGACAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGAAGCAGCAGAATTAACGCGGTGACCGCGCCGTCATGGCTGCGATCGGAGCCAATTTGCTGGTAGCGATAAACCGGCAGGATGCGGAAATCATCGTTACCGCACAAAGCACCCACGCCAAAATCACAACACGACAGCACGAAGGCAAAAGCCacacccagcccccg  >  2:928198/1‑238 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      
CCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAG  >  NC_000913/72649‑73138

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: