Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,516,632 T→C T16A (ACC→GCC)  fecA ← ferric citrate outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,516,6320TC100.0% 289.4 / NA 90T16A (ACC→GCC) fecAferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (36/54);  total (36/54)

CCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTC  >  NC_000913/4516385‑4516878
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
cccAGCGGTTTTTCCTGCAGTCCGCTGCCGGCCAGCAGTTGTTGCAGGCCGTTCTCGACGTCGTAATCGCCGTGCAGGCCGTTTCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAAGCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGGGCTATATTAACCTGTGCAGCAAACGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGtt                                                                                                                                                                                                                                                     <  1:556816/251‑1 (MQ=255)
 ccAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGtt                                                                                                                                                                                                                                                     >  1:41902/1‑250 (MQ=255)
      ggTTTTGCCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCaa                                                                                                                                                                                                                                               <  1:343771/251‑1 (MQ=255)
       gTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCaaa                                                                                                                                                                                                                                              <  1:130724/251‑1 (MQ=255)
        ttttACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAACCGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCaa                                                                                                                                                                                                                                               >  2:501837/1‑249 (MQ=255)
            cccTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGgtgt                                                                                                                                                                                                                                         <  1:258939/250‑1 (MQ=255)
            aCCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGgtgt                                                                                                                                                                                                                                         >  2:352098/1‑251 (MQ=255)
             ccTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGagtt                                                                                                                                                                                                                                        >  1:510253/1‑251 (MQ=255)
             ccTGCAGGCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGAGCGGCATACTGATTGAGCGCTTTATCGAGCGATACCGGTGCGAGATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGagtt                                                                                                                                                                                                                                        >  1:204493/1‑251 (MQ=255)
                  aGTCCGCTGCCGTCCAGCAGTTGTTGCAGGACGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGtttt                                                                                                                                                                                                                                   <  2:626723/251‑1 (MQ=255)
                   ggccGCTGCCGCCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTa                                                                                                                                                                                                                                  <  2:520047/249‑1 (MQ=255)
                       gCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGCGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGa                                                                                                                                                                                                                               >  1:369819/1‑250 (MQ=255)
                              tCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGCAAAAcgcg                                                                                                                                                                                                                       >  1:424629/1‑251 (MQ=255)
                              tCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGAGTTGTTTTACGAAAAAcgcg                                                                                                                                                                                                                       >  1:97697/1‑251 (MQ=255)
                               ccAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTGTTACGAAAAAcgcg                                                                                                                                                                                                                       <  2:576028/250‑1 (MQ=255)
                                 agcagTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGGTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCTAATGGCGGGAACCAAAGGTGTTGTTTTACGAAAAACGCGTaa                                                                                                                                                                                                                    <  2:577365/251‑1 (MQ=255)
                                  gcagTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAAc                                                                                                                                                                                                                   >  2:143767/1‑251 (MQ=255)
                                     gttgttGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGc                                                                                                                                                                                                                <  2:567870/251‑1 (MQ=255)
                                        gttgCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTc                                                                                                                                                                                                             <  1:47943/251‑1 (MQ=255)
                                             gggCCGCTCGCGACGTCGGAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATAcc                                                                                                                                                                                                        <  2:310856/250‑1 (MQ=255)
                                              ggCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTATGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGAAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATCTTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCCAAGGTGTTGTTTTACGAAAAACGCGTAAAGGCGACATACCt                                                                                                                                                                                                       >  2:606600/1‑251 (MQ=255)
                                               gCCGCTCTCGACGTCTTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCtt                                                                                                                                                                                                      <  2:492967/251‑1 (MQ=255)
                                                           gtcgtAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCcatcatcat                                                                                                                                                                                          >  1:516392/1‑251 (MQ=255)
                                                             cgtAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGATGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCAtt                                                                                                                                                                                         <  1:570491/250‑1 (MQ=255)
                                                             cgtAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCAttt                                                                                                                                                                                        <  1:10550/251‑1 (MQ=255)
                                                                   cGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTttgttg                                                                                                                                                                                  <  2:515616/251‑1 (MQ=255)
                                                                    gCCGTGCATGCCGTTGCTCTTCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTtgttgt                                                                                                                                                                                 <  1:502815/251‑1 (MQ=255)
                                                                     ccGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTtgtg                                                                                                                                                                                >  1:104003/1‑251 (MQ=255)
                                                                            ggCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCt                                                                                                                                                                         <  2:190684/251‑1 (MQ=255)
                                                                             gCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCt                                                                                                                                                                         >  2:269596/1‑250 (MQ=255)
                                                                             gCCGTTGCCCTGCTTGCCGCGCCTCAGGCTGGCGTCAACCGAGAGTGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGTTGCGATAGTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTACCCAAAGGGGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTa                                                                                                                                                                        <  1:448879/251‑1 (MQ=255)
                                                                              ccGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTAt                                                                                                                                                                       <  2:616936/251‑1 (MQ=255)
                                                                                      ctGCTTGCCGCGCGTCAGGCTGGCGGCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCATCCCCATCATCACTTTTGTTGTGTTCAGCTATGAGTCg                                                                                                                                                                 <  2:103586/249‑1 (MQ=255)
                                                                                            gCCGGGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGGGTGAGGCATACTGATGGAGCGCTTTATCGAGCGATCCCGGTGCGATAGTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCCGGCGAATTGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGGCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGaa                                                                                                                                                         <  2:502800/251‑1 (MQ=255)
                                                                                            gCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACgaga                                                                                                                                                          <  1:525971/250‑1 (MQ=255)
                                                                                            gCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGaa                                                                                                                                                         <  2:475581/251‑1 (MQ=255)
                                                                                             ccGCGCGTCAGGGTTGCGTCAACCGTGATGGCAAAGCCGTTGTGTGCGGCATACTGATTGAGCGCTTTAGCGAGCGATCGCGGTGCGATATGAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCGCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTGCTCCATCATCATTGTTGTTGTGTTCAGCTATGAGTCGAACGGGAAg                                                                                                                                                        <  2:70078/251‑1 (MQ=255)
                                                                                                   gTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTaaa                                                                                                                                                  <  1:209111/251‑1 (MQ=255)
                                                                                                      cggCTGGCGTCACCCTAGAGGGTAAATCCGCTGTGGGCGGCATACTGATTGAGCGCGTGATCGAGCGATCCCGGTGCGATAGTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGGTAACCAAAGGTGTTTTTTTACGAAAATCGCGTAACGGCGTCAGACCTTCCCCATCATCATTGTTGTTGTGTTCAGCTATGATTCTAACGAGAAGGGTaaaca                                                                                                                                                <  2:274674/249‑3 (MQ=255)
                                                                                                              gTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCg                                                                                                                                       <  1:306583/251‑1 (MQ=255)
                                                                                                               ccAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATGGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTTATACCTTCCCCATCAGCAATTTTGTTGGGGTCAGCTATGAGTCGAACGAGAAGGGTAAAATGGACAATCGa                                                                                                                                      <  2:273988/250‑1 (MQ=255)
                                                                                                                 aaCCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGaa                                                                                                                                     >  1:585251/1‑250 (MQ=255)
                                                                                                                 aaCCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCACGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGaaa                                                                                                                                    <  2:104003/251‑1 (MQ=255)
                                                                                                                  aCCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAATCGAGAGACCGGCCAGCTGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAGCTGCGGCAGACCTTCCCCATCATCACTTTGGTCGTGTACAGCTATGAGTGGAACGTGAAGGGTAAAAATAAGAATCGAAAt                                                                                                                                   <  2:29176/251‑1 (MQ=255)
                                                                                                                  aCCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAAt                                                                                                                                   >  1:393056/1‑251 (MQ=255)
                                                                                                                     gcgagGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGCAAATCGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTCAACAGGACAATCGAAATaa                                                                                                                                 >  2:332959/3‑250 (MQ=255)
                                                                                                                      agagGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATa                                                                                                                                  >  1:177002/1‑248 (MQ=255)
                                                                                                                       gagGGTAAATCCGCTGTGTGCGGCATACTGAGTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATTCCTTCCGCATCATCATTTTTGTTGTGGTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGaa                                                                                                                              <  2:365059/251‑1 (MQ=255)
                                                                                                                        agGGTAAATCCGCTGGGTGCGGCATACTGATTGAACGCTTTAGCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAat                                                                                                                             <  2:411881/251‑1 (MQ=255)
                                                                                                                             aaaTCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAattatt                                                                                                                         <  1:154788/250‑1 (MQ=255)
                                                                                                                              aaTCCGCTGTGTTCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGGTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTAtttt                                                                                                                       <  2:122393/251‑1 (MQ=255)
                                                                                                                              aaTCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTAtttt                                                                                                                       <  1:355891/251‑1 (MQ=255)
                                                                                                                               agccGCTGTGTGCGGCATACTGATTGAGCGTTATATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCGAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGTCTCAGGCGAATGGCGTTAACCAAATGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTATAAAGGACAATCGAAATAAGAATTATTTTc                                                                                                                      <  2:350386/249‑1 (MQ=255)
                                                                                                                                               aTACTGATTGAGCGCTTTATCGAGCGCTCACGGTGCGCTATTAACTTGTGCAGCAAAAGCGGAAACCGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCTTACCTTCCCCATCATCATTTTTGTTGTGTTCCGCTATGAGTCGAACGAGCAGGGTAAAAAGGACAATCGAAATAAGCATTATTTTCCTTACAGTGGTGAAAt                                                                                                      >  2:293665/1‑251 (MQ=255)
                                                                                                                                                   ggATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTATCCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTGACAGTGGTGAAATGttg                                                                                                  <  2:545981/250‑2 (MQ=255)
                                                                                                                                                       tGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATc                                                                                               <  1:585231/250‑1 (MQ=255)
                                                                                                                                                        gAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAAGCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCa                                                                                             <  1:585936/251‑1 (MQ=255)
                                                                                                                                                          gcgcTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGt                                                                                           <  1:39725/251‑1 (MQ=255)
                                                                                                                                                          gcgcTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTGCCCCAGCATCATTGTTGTTGTGTGCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGt                                                                                           <  1:570623/251‑1 (MQ=255)
                                                                                                                                                               ttATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTAcc                                                                                        <  2:496753/249‑1 (MQ=255)
                                                                                                                                                                    gAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAAt                                                                                 >  1:256967/1‑251 (MQ=255)
                                                                                                                                                                     aGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAAt                                                                                 >  1:434808/1‑250 (MQ=255)
                                                                                                                                                                        gATCCCGGTGCGATATTCACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGGGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGAc                                                                             >  2:434488/1‑251 (MQ=255)
                                                                                                                                                                         aTCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACt                                                                            >  2:172550/1‑251 (MQ=255)
                                                                                                                                                                          tCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGCAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTg                                                                           >  2:25243/1‑251 (MQ=255)
                                                                                                                                                                           cccGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTg                                                                           <  1:19598/250‑1 (MQ=255)
                                                                                                                                                                           cccGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTg                                                                           <  1:606600/250‑1 (MQ=255)
                                                                                                                                                                            ccGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGaa                                                                         <  2:144724/251‑1 (MQ=255)
                                                                                                                                                                               gTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAAttt                                                                      >  2:117553/1‑251 (MQ=255)
                                                                                                                                                                                   gATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAAc                                                                  <  2:60753/251‑1 (MQ=255)
                                                                                                                                                                                       ttAACCTGTGAAGCAAAAGCGGAAAACGAGAGACCGGCCAACGGCAGCAGGCTCAGACGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGCAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACggg                                                               >  2:116438/1‑250 (MQ=255)
                                                                                                                                                                                              gtgCAGCAAAAGCGGAAAACGAGAGGCCGGCCAGCGGCAGCAGGCTCAGGCGAATTGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCGTCCCCATCATCATTGTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTGCCTTACAGTGGTGAAATGTTTATCCAGTACCTCGTAATAGACTGAATTTTGACGGGAAGCGttt                                                       <  2:29637/251‑1 (MQ=255)
                                                                                                                                                                                                gcagcaAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTg                                                      <  1:161274/250‑1 (MQ=255)
                                                                                                                                                                                                     aaaaGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGAt                                                <  2:434808/251‑1 (MQ=255)
                                                                                                                                                                                                     aaaaGCGGAAAACGAGAGACCGGCCAGCGGCAGAAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGACTTTTAACGGGAAGCGTTTGCGCGAt                                                >  1:127782/1‑250 (MQ=255)
                                                                                                                                                                                                       aaGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAg                                                                                                                                <  1:469113/169‑1 (MQ=255)
                                                                                                                                                                                                       aaGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAg                                                                                                                                >  2:469113/1‑169 (MQ=255)
                                                                                                                                                                                                       aaGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATa                                               >  1:155109/1‑250 (MQ=255)
                                                                                                                                                                                                       aaGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCAAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGaga                                               >  2:471336/1‑248 (MQ=255)
                                                                                                                                                                                                            gAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTc                                         >  1:353941/1‑251 (MQ=255)
                                                                                                                                                                                                             aaaaCGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGACTGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGTCTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTc                                         >  2:596494/1‑250 (MQ=255)
                                                                                                                                                                                                                 cGAGAGACCGGCCAGCGGCAGCAAGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCTATAACGTTCAGGAt                                    <  2:72254/251‑1 (MQ=255)
                                                                                                                                                                                                                      gaCCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAACGTGTTGTTTTACGAAGAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAATCGTTTGCGCAATAACGTTCAGGATCGCAt                               <  1:484617/251‑1 (MQ=255)
                                                                                                                                                                                                                       aCCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATc                              >  1:139365/1‑251 (MQ=255)
                                                                                                                                                                                                                        ccGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCg                             >  1:84042/1‑251 (MQ=255)
                                                                                                                                                                                                                          ggCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGt                           <  1:14290/251‑1 (MQ=255)
                                                                                                                                                                                                                                 ggCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTc                    <  1:352098/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gcagcaGGCTCAGGCGTATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGACCGAGAAGGGTAAAAAGGACAATCGAACTAAGAATTATTTTCCTTACAGTGGTGACATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGTTCGCTTCGGTATTTTTc                    >  2:105218/1‑250 (MQ=255)
                                                                                                                                                                                                                                                 gAATGGCGTTAACCATAGGGGTTGTTTTACGCAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGTGAACGTccggc    >  2:248494/1‑248 (MQ=255)
                                                                                                                                                                                                                                                   aTGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGCAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTc  >  1:4240/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
CCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTC  >  NC_000913/4516385‑4516878

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