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,6320TC96.7% 267.6 / ‑0.6 92T16A (ACC→GCC) fecAferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (39/48);  total (43/49)
Fisher's exact test for biased strand distribution p-value = 5.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.58e-01

CCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTC  >  NC_000913/4516385‑4516878
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
cccAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGtt                                                                                                                                                                                                                                                      <  1:183592/251‑1 (MQ=255)
cccAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGAGtt                                                                                                                                                                                                                                                      >  1:636325/1‑251 (MQ=255)
 ccAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTa                                                                                                                                                                                                                                                     <  1:215049/251‑1 (MQ=255)
     cGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAAcc                                                                                                                                                                                                                                                  >  1:588174/1‑250 (MQ=255)
       gtattACCTACAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGCCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGCGAGGGTAAATACGCTGTGTGCGGCATACTGATTGAGCGCTTTATAGAGCGCTCCCGGTGCGATATTAACCTGTTAAGCCAAAGAGGAAAACCAAAGAACGGCCAACGGCCGGCGGCGCAGTAGAATGGCGTTAACCac                                                                                                                                                                                                                                                >  2:467011/4‑249 (MQ=255)
             ccTGCAGGCCGCTGCCGTCCAGCAGTGGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAAGCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGtgtt                                                                                                                                                                                                                                         <  2:132593/251‑1 (MQ=255)
                gccgTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCACCCGCGAGGGTAAATCCGCTGTGTGCGGCCTACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCCGCAACAGCGGACATCGAGCTACCGGCCCGCGGCAGCAGGCTCAGGCGAATGGCGTTaatccacggagttttt                                                                                                                                                                                                                                      >  2:593709/4‑237 (MQ=255)
                   gTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTa                                                                                                                                                                                                                                   <  1:311642/251‑1 (MQ=255)
                            cGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGAGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGATAACCAAAGGTGTTGTTTTACGAAAAAcg                                                                                                                                                                                                                          >  2:446564/1‑251 (MQ=255)
                               ccAGCAGTTGTTGCAGGCCGCCCTCGCCGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTAGCGAGCGAGCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAAcgcg                                                                                                                                                                                                                        <  2:279266/250‑1 (MQ=255)
                                 agcagTTGATGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGGAAATCAGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTTCGGAAAActcgga                                                                                                                                                                                                                      >  2:376160/1‑245 (MQ=255)
                                 agcagTGGTTGCAGGCCGCTCTCGCCGTCGTAATCGCCGTTCAGGCCGTTGCTCTGCTTGCCGCGCGGCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTaa                                                                                                                                                                                                                     <  2:18541/251‑1 (MQ=255)
                                        gttgCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTc                                                                                                                                                                                                              <  1:446564/251‑1 (MQ=255)
                                          tgCAGGCCGCTCTCGACGTCGTAATCGACGAGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGGAAATCCGCTGCGAGCGGAACACTGATGGAGCGCTATAGCGAGCTATACCGGTGCGGTATCAAACTGTGAAGCAAGAGCGGACCCCGCGAGACCGGCAAGCGGCAGAAGGCTCAGGCGAATGGCGCTGAACAAAGATGTGGTTTTCCGaaaacagcaatgcggtgtact                                                                                                                                                                                                            >  2:18559/1‑234 (MQ=255)
                                            cAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATAc                                                                                                                                                                                                          >  2:363394/1‑251 (MQ=255)
                                            cAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATAc                                                                                                                                                                                                          >  1:351799/1‑251 (MQ=255)
                                               gCCGCGCTCGACGTCGTAATCGACGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACGGATTGAGCGATTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCtt                                                                                                                                                                                                       <  1:680429/251‑1 (MQ=255)
                                                 cGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGCGGGTACATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCAtcacttcc                                                                                                                                                                                                     >  1:385596/1‑244 (MQ=255)
                                                  gCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTAAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTccc                                                                                                                                                                                                    >  1:462605/1‑251 (MQ=255)
                                                       cGACGTCGTAATCGCCGTGCAGGCCGGTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCcatca                                                                                                                                                                                               <  1:15085/251‑1 (MQ=255)
                                                        gACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCcatcat                                                                                                                                                                                              <  1:432532/251‑1 (MQ=255)
                                                                              ccGTTTCTCTGTTTGCGGCGCGTCAGGCTGTCGTCAACCGTGAGGGTAAATCCGCTGGGGGCGGCATACTGATTGAGCGCTTTATCGAGCGAGCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCTAATGGCGTTAACCAAAGGTGTTGTTTTTCGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCATCTAt                                                                                                                                                                        <  2:162599/251‑1 (MQ=255)
                                                                               cGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGAGCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATg                                                                                                                                                                       <  2:266109/251‑1 (MQ=255)
                                                                                        gCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCTATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAAc                                                                                                                                                               <  1:438002/250‑1 (MQ=255)
                                                                                        gCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAAc                                                                                                                                                               <  1:36220/250‑1 (MQ=255)
                                                                                               gcgcgTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAggg                                                                                                                                                       <  1:656625/251‑1 (MQ=255)
                                                                                               gcgcgTCAGGCTGGCGGCAACCGTGAGGGTAAATCCGCGGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGGCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAggg                                                                                                                                                       <  2:679730/251‑1 (MQ=255)
                                                                                                 gcgTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTa                                                                                                                                                     <  1:289388/251‑1 (MQ=255)
                                                                                                  cgTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTaa                                                                                                                                                    <  2:476941/251‑1 (MQ=255)
                                                                                                  cgTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGGGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTaa                                                                                                                                                    <  1:92503/251‑1 (MQ=255)
                                                                                                     cAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTTTTGTGTTCAGCTATGATTCGAACGAGAAGGGAaaaa                                                                                                                                                  >  2:623328/1‑250 (MQ=255)
                                                                                                     cAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAAcagatcggaagagc                                                                                                                                                  <  1:593709/250‑14 (MQ=255)
                                                                                                     cAGGCTGGCGTCAACCGAGAGGGTAAAACCGCTGTGTTCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTaaaaa                                                                                                                                                 <  1:308028/251‑1 (MQ=255)
                                                                                                        gCTGTCGTCAACCGAGAGGGTAAATCCGATGTGTGCGTCATACTGATTGTGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCATCTGGCTCAGGCGAATGGCGTTAACCAAAGGGGTTGTTTTACGAAAAACTCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTTTGAGTCGAACGAGAAGGGTAAAAAgta                                                                                                                                              <  1:627419/251‑3 (MQ=255)
                                                                                                        gCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATTCCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGa                                                                                                                                              <  1:397450/251‑1 (MQ=255)
                                                                                                        gCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGa                                                                                                                                              <  1:227217/251‑1 (MQ=255)
                                                                                                        gCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGa                                                                                                                                              <  1:627394/251‑1 (MQ=255)
                                                                                                         cTGGCGGCAACCGCGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGAGCGAGATTAACATGTGAAGCAAAAGCGGAAAACGAGAGACCGGCCATCGGCAGCAGGCTCAGGCGAAGGGCGGTAACCAAAGGGGTTGGTTTACGAAAAACGCGTAAAGGCGTCATACCTTGACCTTCATGCTTTTTGTTGTGTTCAGCAATGAGTAGAAAGAGAAGGGAAAAAAGGAc                                                                                                                                             >  2:179737/1‑251 (MQ=255)
                                                                                                           ggCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACaa                                                                                                                                           <  1:349064/251‑1 (MQ=255)
                                                                                                           ggCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGAGTGAGCGCTTTATCGAGCGATCCCGTTGCGATATTAACCGGTGCAGCAAAAGCGGAAATCGAGAGACCGGCCAGCGGCAGCAGGCTCAGTCGAATGGCGTTAACCAAATGTGTTGTTTTTCGAAAAACGCGTAACGGCGTCATACCTGCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGGCGAACGAGAAGGTTAAAAAGGACa                                                                                                                                            <  2:620071/250‑1 (MQ=255)
                                                                                                            gcggcAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCGCAGGCGAATGGCGTTAACCAAATGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAAt                                                                                                                                          <  2:613719/247‑1 (MQ=255)
                                                                                                            gCGTCAACCGAGAGGGTAAATCCGCTGGGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCTAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTTTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAAt                                                                                                                                          <  1:26030/251‑1 (MQ=255)
                                                                                                              gTCAACCGAGAGGGTAAATCCGATGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGAGTTAACCAAAGGGGTTGTTTTACGAAAACCGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTacacaagataatcg                                                                                                                                        >  2:108412/1‑238 (MQ=255)
                                                                                                                    cGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATa                                                                                                                                   >  1:598189/1‑250 (MQ=255)
                                                                                                                     gagaGGGGAAATCCGCTGTGGGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAg                                                                                                                                 <  2:296369/251‑1 (MQ=255)
                                                                                                                        agGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTAt                                                                                                                                                                        <  2:176203/209‑1 (MQ=255)
                                                                                                                        agGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTAt                                                                                                                                                                        >  1:176203/1‑209 (MQ=255)
                                                                                                                         gggTAAATCCGCTGTGTGCGGCATTCTGATTGAGCGCTGTATCGAGCGATCCCGGTGCGTTAGTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGAGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAatt                                                                                                                             <  2:530163/251‑1 (MQ=255)
                                                                                                                         gggTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGCGCGATCCCGGTGCGATATTAACCTGGGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGAAGCAGGATCAGGAGAATGGCGGTAACCAAAGGTGTTGTTTTACGAAAAACGCGCAAGGGCGTCATACCTTCCCCATCATCATTTTTGATGAGTTCAGCTATATGTAGACAGAGAATGGGAAAAAGGTCAATCGAGATAAGAatt                                                                                                                             >  2:404157/1‑251 (MQ=255)
                                                                                                                           gTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCg                                                                                                                                        >  1:241214/1‑238 (MQ=255)
                                                                                                                           gTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCg                                                                                                                                        <  2:241214/238‑1 (MQ=255)
                                                                                                                             aaaTCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTAttt                                                                                                                         >  1:449960/1‑251 (MQ=255)
                                                                                                                                 ccGCTGTGTGCGGTATTCTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTcc                                                                                                                      <  1:456839/250‑1 (MQ=255)
                                                                                                                                                   tGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACTAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGAGTTGTTTTACGAAAAACGCGTCCCGGCGTCATa                                                                                                                                                                                                           >  2:484802/1‑147 (MQ=255)
                                                                                                                                                   tGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATa                                                                                                                                                                                                           <  1:484802/147‑1 (MQ=255)
                                                                                                                                                     aTTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTAt                                                                                                 >  1:15954/1‑251 (MQ=255)
                                                                                                                                                      ttGAGCGCTTTATCGAGCGATCCCGGAGCGATATTAACCTGTGCAGCACAAGCGGAAAACGAGAGACCGGCCAGCAGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAACACGCGTAACGGCGTCATACCTTCCCCCTCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGGGAAGGGTAAAAAGGACAATAGAAATAAGAATTATTTTCCTTACAGTGGGGAAATGATTaac                                                                                                >  2:45141/1‑249 (MQ=255)
                                                                                                                                                           cgcTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACaa                                                                                                                                           >  1:252249/1‑203 (MQ=255)
                                                                                                                                                           cgcTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACaa                                                                                                                                           <  2:252249/203‑1 (MQ=255)
                                                                                                                                                            gcTTTATCGCGCGATCCCGGCGCGATATTACCCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGATGTGTTCAGCCATGAGTCGAAAGTGAAGGGTAAAACGGACAATCGAAATAAGCATTATTTTCCTTACAGTGGTGCATAGTTTATCCATTAc                                                                                          >  2:260928/1‑251 (MQ=255)
                                                                                                                                                              tttATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGACCGAGAAGGGTAAAAAGGACAATCGAAATAAGACTTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACcg                                                                                       >  2:88817/1‑251 (MQ=255)
                                                                                                                                                                    gAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGCGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAACGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTTTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTCACGCGTAAt                                                                                 >  1:277653/1‑251 (MQ=255)
                                                                                                                                                                       ggATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTGAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATGTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGa                                                                              <  1:591904/250‑1 (MQ=255)
                                                                                                                                                                        gATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGCCCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACgaga                                                                                                                                                           >  2:191804/1‑174 (MQ=255)
                                                                                                                                                                        gATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACgaga                                                                                                                                                           <  1:191804/174‑1 (MQ=255)
                                                                                                                                                                             cggttcgattatacccggtgCAGCAAAAGCGGAAAACGAGAGATCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCTTTAACTAACGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTTTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGTGTAATAGTCTGAAt                                                                        <  2:534329/234‑1 (MQ=255)
                                                                                                                                                                                tGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGACTGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATAATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAAAATTATTTTCCTTACAGGGGTGAAATGTTTATCAAGTACCGCGTAATAGATTGAAtttt                                                                     >  2:562725/1‑251 (MQ=255)
                                                                                                                                                                                  cGATATTAACCTGTGCAGCAAAAGCGGAAACCGAGCGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAACGGAGTTGTTTTACGAAACACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGCGCAGGGAAAACAGGACAATCGAACAAATAATTCTTTTCCTTACAGAGGGGAAATGTTTGTCCATACCCGCGTATTGGACTGAAttata                                                                    >  2:201667/1‑247 (MQ=255)
                                                                                                                                                                                     tatTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAAc                                                                                                                                                               <  2:568219/157‑1 (MQ=255)
                                                                                                                                                                                     tatTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAAc                                                                                                                                                               >  1:568219/1‑157 (MQ=255)
                                                                                                                                                                                        tAACCTGTGCAGCAAAAGCGGAAAATGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGa                                                              <  1:503906/250‑1 (MQ=255)
                                                                                                                                                                                                 cagcaAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAATTAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTgcg                                                    >  2:362099/1‑251 (MQ=255)
                                                                                                                                                                                                   gcaAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCg                                                                                                                                                                  >  1:235602/1‑140 (MQ=255)
                                                                                                                                                                                                   gcaAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCg                                                                                                                                                                  <  2:235602/140‑1 (MQ=255)
                                                                                                                                                                                                           ggAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTTCGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGtt                                          <  2:42684/251‑1 (MQ=255)
                                                                                                                                                                                                              aaaCGAGCGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTc                                         >  2:623318/1‑249 (MQ=255)
                                                                                                                                                                                                               aaCGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAACACGCGTAACGGCGTCATAACTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGATTATCCAGTACCGAGTAATAAACTGACTTTTAACGGGCAGCGTTTGCGCGATacagttcatg                                      >  2:583683/1‑242 (MQ=255)
                                                                                                                                                                                                               aaCGAGAGACCGGCCAGAGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGCCTGAATTTTAACGGGAAGCGTTTGAGCGATAACGTTCAgg                                      >  2:170390/1‑251 (MQ=255)
                                                                                                                                                                                                                aCGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATTGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTgc                                                     <  2:362916/235‑1 (MQ=255)
                                                                                                                                                                                                                aCGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTgc                                                     >  1:362916/1‑235 (MQ=255)
                                                                                                                                                                                                                       aCCGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATc                              <  1:618975/251‑1 (MQ=255)
                                                                                                                                                                                                                        ccGGCCAGCGGCAGCAGGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTACAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGAACAAGACTGAATTTTAACGGAAAACGTTTGCGCGATAACGTTCAGGATCGCATCg                             >  2:174508/1‑251 (MQ=255)
                                                                                                                                                                                                                           gCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGt                           >  1:333333/1‑250 (MQ=255)
                                                                                                                                                                                                                                 gtcagcagGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACACGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTc                    <  1:427050/249‑1 (MQ=255)
                                                                                                                                                                                                                                       agGCTCAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGCGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGAGTTTGCGCGATAACGTTAAGGATCGCATCGGTATTTTTCAGCggg              >  2:1303/1‑251 (MQ=255)
                                                                                                                                                                                                                                            cAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGt         <  1:363394/251‑1 (MQ=255)
                                                                                                                                                                                                                                            cAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGt         >  2:617258/1‑251 (MQ=255)
                                                                                                                                                                                                                                            cAGGCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACAGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCATCGGGagagg         >  2:553296/1‑247 (MQ=255)
                                                                                                                                                                                                                                               gCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTGTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTccc      <  1:623318/251‑1 (MQ=255)
                                                                                                                                                                                                                                               gCGAATGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAAAGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTAGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGAATCGGTATTTTTCAGCGGGAACGTccc      >  1:410395/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  aaTGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCt   <  1:530879/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   aTGGCGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGCAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTc  >  1:444207/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 < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: