Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,514,336 A→T 7.6% intergenic (‑28/+17) fecB ← / ← fecA ferric citrate ABC transporter periplasmic binding protein/TonB‑dependent outer membrane ferric citrate transporter and signal transducer; ferric citrate extracelluar receptor; FecR‑interacting protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,514,3360AT 7.6% 187.2 / 5.3 66intergenic (‑28/+17)fecB/fecAferric citrate ABC transporter periplasmic binding protein/TonB‑dependent outer membrane ferric citrate transporter and signal transducer; ferric citrate extracelluar receptor; FecR‑interacting protein
Reads supporting (aligned to +/- strand):  ref base A (29/32);  new base T (3/2);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.07e-01

GCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTGAACG  >  NC_000913/4514048‑4514636
                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                                            
gcgcgTTCCGACGGACTGCCACTGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTaa                                                                                                                                                                                                                                                                                                  <  2:138775/301‑1 (MQ=255)
cctcgTTCCGTCGGCCTGCCCCGGTTTCAGGTGTGCACGCACTTCTGGCAGGATGCGGTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCTTGCCGGGCtttta                                                                                                                                                                                                                                                                                                  <  1:487018/298‑3 (MQ=255)
  gcgTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCTTGCCGGGCttttagc                                                                                                                                                                                                                                                                                                <  2:532166/301‑5 (MQ=255)
  gcgTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAaca                                                                                                                                                                                                                                                                                                <  2:200596/301‑1 (MQ=255)
   cgTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCTTGCCGGGCttttaaga                                                                                                                                                                                                                                                                                               >  2:487664/1‑296 (MQ=255)
        ggacggacgGCCACGGTTTGAGGTGCGGACGAACTGTGGTCAGGATGCGTGTTGCATCGTGATCGTCGGCAATACCGATCGGGATGACGTCCACGTCGGCCAGCGCATCGGCGTACGGGAGGTCCAGCCCCTCAATCCGTTGTGGCTTTTTTTCGTGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACAGTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAg                                                                                                                                                                                                                                                                                          <  1:500080/300‑1 (MQ=255)
        cGACGGACTGCCACGTTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAg                                                                                                                                                                                                                                                                                          <  2:554496/301‑1 (MQ=255)
        cGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCCAAAAGACAACGGATAAATGCAAACAAAATCAACATCCAGCTAAAAACCCGCACAGCCGGGCCTTaaaccaacaa                                                                                                                                                                                                                                                                                          >  1:212178/1‑293 (MQ=255)
        cGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACAtaaa                                                                                                                                                                                                                                                                                          >  2:262902/1‑298 (MQ=255)
                  ccACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAAca                                                                                                                                                                                                                                                                                >  2:380204/1‑300 (MQ=255)
                           ctggTGCGCAGGTATTTTGGGGAGGTTGCGTTTGGCATCGTTATCGTCTGCAATACCGATCGGGATGACGGCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCa                                                                                                                                                                                                                                                                       <  1:338632/299‑1 (MQ=255)
                           cAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCa                                                                                                                                                                                                                                                                       <  2:512393/301‑1 (MQ=255)
                             ggTGGGCACGCCGTTCGGGCAGGATGGGTTTTGCTTCGTTATCGTCGGCCATACTGTTCGGTATGCCGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCata                                                                                                                                                                                                                                                                     <  1:112739/301‑1 (MQ=255)
                                              ataacgatgcaaaacgCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTg                                                                                                                                                                                                                                                    >  2:377056/16‑301 (MQ=255)
                                               gCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGAGGCTgg                                                                                                                                                                                                                                                   >  2:631973/1‑301 (MQ=255)
                                                   ggtGCGTTTTGACTCGTTAGCGGTGGCAATACCGATTGGGATGACGGCCAGGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGAACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTg                                                                                                                                                                                                                                               <  1:316166/299‑1 (MQ=255)
                                                    aTGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGc                                                                                                                                                                                                                                              >  2:504202/1‑301 (MQ=255)
                                                                        gtcgGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTtgttgttgt                                                                                                                                                                                                                          >  2:492858/1‑301 (MQ=255)
                                                                          cgGCAATGCCGGTGGGGGTGACGTCCACGGCGGCCAGCTCAGCGGCGAACGAGAGTTCCAGCACCACAGTCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTgtcg                                                                                                                                                                                                                        <  1:557800/301‑1 (MQ=255)
                                                                               atcccGGTCGGGATGACGGCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATa                                                                                                                                                                                                                   <  1:379708/298‑1 (MQ=255)
                                                                                               cGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCCACGACCCCTGCCTATACAGCGTGCGCGTCTGGCCTGCATAGAACCCTTTGTTGTTGTCGTCATAAGAGCCGATTAaataa                                                                                                                                                                                                   >  2:301023/1‑300 (MQ=255)
                                                                                                  ccACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTcc                                                                                                                                                                                                <  2:212421/301‑1 (MQ=255)
                                                                                                   cACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCAAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTTCATAGATGCCTTTGTTGTTGTCGTCATAAGCGCGGATAAagcacacct                                                                                                                                                                                               >  1:610576/1‑294 (MQ=255)
                                                                                                         ggcCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCCTATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATACGCGCGGATGAAGTAGTCATGGACGa                                                                                                                                                                                         >  2:207710/1‑301 (MQ=255)
                                                                                                             aGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTGTTTCGAGTGTAAACGTGCCTTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGAt                                                                                                                                                                                     <  2:2665/301‑1 (MQ=255)
                                                                                                                  gctcttccgatctagTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGtttt                                                                                                                                                                                <  2:25611/288‑1 (MQ=255)
                                                                                                                     ggCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGCAGTAGTACTGGTCGAAGATGTTTTTca                                                                                                                                                                             >  2:485533/1‑301 (MQ=255)
                                                                                                                     ggCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGAAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAAAATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTGGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCAGGTGGAAGATGTTTTTca                                                                                                                                                                             >  2:448347/1‑301 (MQ=255)
                                                                                                                      gcggacgtgctttCCAGTACCTCATTCCGTTTGTGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGGGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAACACGTCTAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGACCTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTcac                                                                                                                                                                            <  1:447783/290‑1 (MQ=255)
                                                                                                                       cGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAACCGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCCTATACAGCGTGCGCGGCTGGCCTGCATAGATCCCTTTGCTGTTGTCGTCATAAGAGCGGATGAAGTAGTCATGGGCGAAACAGTTTTTcacc                                                                                                                                                                           >  1:116052/1‑300 (MQ=255)
                                                                                                                               agTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGttttagatcggaagagc                                                                                                                                                                   >  1:25560/1‑288 (MQ=255)
                                                                                                                                         ccacAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTTCGCGGCTGGCCTGCATCGATGCCTTTGTTGTTGTCGTCATAAGAGCTGAACAAGTAATCCTGGTCGAAAATGTTTTTCACCACGAATGCAAAGTTCAGa                                                                                                                                                         >  1:560212/1‑301 (MQ=255)
                                                                                                                                           acaATTCGTTTGTGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATc                                                                                                                                                       <  1:631096/301‑1 (MQ=255)
                                                                                                                                             attCCGGTGTGTCGTTTTTTCGGGTGTAAACGTGCCGTGTTCGTCCTGAACCGGGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACAGAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTg                                                                                                                                                     <  1:122755/299‑1 (MQ=255)
                                                                                                                                                ccGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACGTCAGAACGTCAGCGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCa                                                                                                                                                  <  1:262634/301‑1 (MQ=255)
                                                                                                                                                         ggTTTTTTCGATTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAGTGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGc                                                                                                                                         <  1:279335/300‑1 (MQ=255)
                                                                                                                                                             tttgCGGGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGAAAAGGGGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCCTATACAGCGTGCGCGGCTGGCCTGCATAGAGGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCaaa                                                                                                                                     <  1:207467/301‑1 (MQ=255)
                                                                                                                                                              ttgggatggtaacctgggcgggtTCGTCCGGAACAGTGGCGGCAAAGGCGTGGCTGTTCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACAGTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACAGCAGAACTTCAACGACCCCGGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAg                                                                                                                                    <  1:626929/280‑1 (MQ=255)
                                                                                                                                                                tCGGGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATGCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTc                                                                                                                                  <  1:503530/301‑1 (MQ=255)
                                                                                                                                                                                                    gcgtcAAAGGCGTTGCTGTTTACCGGTAGCTGGCCGGCAAAAAGAAATCGGATAAATGCCAACATACTCACCTTCCAGTTACAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCAGATACAGCGTGCGCTGCTGGCCTGCATAGATGCCGTTGGTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCggg                                                                                              <  1:484932/301‑1 (MQ=255)
                                                                                                                                                                                                    gaggcAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCggg                                                                                              <  2:560926/299‑1 (MQ=255)
                                                                                                                                                                                                      ggcAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGTCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACACCACGCGTGCGACCCAGAGCATAAAGACAGGaa                                                                                            >  2:386224/1‑301 (MQ=255)
                                                                                                                                                                                                        cAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTACTAAGAGCGGATGAAGTAGTCCTGGTCGCAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCACCCCGGTCGCCCCAGAGCATGAAGcacgcaccg                                                                                          >  1:419955/1‑293 (MQ=255)
                                                                                                                                                                                                               gTGGCTGATCCACAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATGCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCgg                                                                                   <  1:376560/301‑1 (MQ=255)
                                                                                                                                                                                                                     gATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTCGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCAGGCCAAAGTCATAAGCCACGCGTGCGCCACAGAGACAGAAGCCCGGAAAGCGGCCGGGAcacc                                                                             >  1:486821/1‑298 (MQ=255)
                                                                                                                                                                                                                              cagcagGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGgcgc                                                                    <  2:443855/301‑1 (MQ=255)
                                                                                                                                                                                                                                 cagGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGACCGGTACTGCCGTCGGCCctcc                                                                 >  1:541670/1‑300 (MQ=255)
                                                                                                                                                                                                                                  agGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACGTCAACGACCCCTGCATATCAAGCGTGCGCGGCTGGCCTGCATAGATGCCTTGGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCGGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGctctc                                                                <  1:492212/301‑1 (MQ=255)
                                                                                                                                                                                                                                        gCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTAATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCAGGTACTGCCGTCGGCGCTCTCTTTCAc                                                          >  2:631141/1‑301 (MQ=255)
                                                                                                                                                                                                                                                  aaCGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGACGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAACGTCATACGCCACGCGTGCGCCCCCGAGCATGCAGCCGGTAATGCGGACGATACTGCCGTCgcccaccccttcccacgtattgccc                                                >  1:508984/1‑277 (MQ=255)
                                                                                                                                                                                                                                                               cAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGCCGTTATCCGCAAAct                                   >  1:247388/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                    taataaCATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTAAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCCCTTTCACCGTATTGGCGTTATACGCAAActggctg                              >  2:545884/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                    tAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGCTGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAAATCTGCCATCTGCGGGCCAAAGTCATACCCCACGCGTGCGCCCCACATCATGACGCCGAGAATGAGTCCGGTACTACCGTCGGCGCTCTCTTTCACCATATTTGCATTATCCGCAAAatagctt                              >  2:545886/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                       tCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGctggac                           >  2:41708/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                        ccctgtccgcttAAAAGCCCGGCAAGTCGGGCGTTAACACAGCAGAGCTTCAACGACCCCTGCGTATACAGGGTGCGCGGCTGGCTTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAATTAGTCCTGGTCGAAGATGGTTTTCGCACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCAACGCGTGCGCGCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTATCACCGTATTGGCGTTATCCGCAAACTGGctggact                          <  1:545211/290‑1 (MQ=255)
                                                                                                                                                                                                                                                                        aagcttCCAGCGAAAGGCTCGGCAAGGCGGTGGTTAACAATTCAGAACTTGAACGACATCTTCAGATACAGTGGGGGCGGCTGGACTGCATATATGCTGTTGTTGGTGTCGTCATAATGTAGGATGTAGTCGTCCGTGTCGAAGATGTTTTTCATGCCGAATGCGAGGTTCAGATCTGCCATCTGCGGGCCAAAGTGATACGCGACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTAGGCGCTCTCTTTCAACGTATTGGCGTTATCCGCAAACTGGctggact                          <  1:545213/297‑1 (MQ=255)
                                                                                                                                                                                                                                                                         acaTTCCAGCTAAAATGTCGGCAAGTCGGGTGTTAACATATCATGACTGCAACGACCCGTGAATATACAGCGTGCGCGGCTGGCCTGCATAGAGGTCGTTGTTGTTGTCGTCAGAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCAGCTGCGGGCCAAAGTCAGAGGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGctggactg                         <  1:45228/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                         acaTTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGctggactg                         >  2:212516/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                          aagtCCAGCGTAAAGCCTGTCAAGCCGGGCGTTAACACATCAGAGCTGCAACGACCTCTGCATTTACCGCGTGGGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTAGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTATCAGACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGTCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGctggactgg                        <  1:630264/298‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  ttAAAGGCCCCGCCAGCCGGGCGTTAACACATGAGTACTTACACGACCCCGGCATTTACAGCGTGCGCGGCTGGCCTGGATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTTGTCGAAGATGTTTTTCACACCGAATGCCAGGTGCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCt                <  1:401046/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                   tAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCta               >  2:401560/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                     aacgCCCGGCTTGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTACACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAACGCCATACGCCACGTGTGCGCCCCAGAGCATGAAGCAGGGAATGCGGCCGGTACTGACGTCGGCGCTCACCTCACCCGTATTAGCGTAATCCGCAAACTCGctggaccgacaatccctcat             >  1:624181/4‑287 (MQ=255)
                                                                                                                                                                                                                                                                                      acgCCCGGCTTGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGttt            >  2:360510/3‑300 (MQ=255)
                                                                                                                                                                                                                                                                                      aaGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTc            >  2:126735/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                        gCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACACGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCaa          >  2:443693/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                             gCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCTTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTGa     <  2:247655/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                aGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCCTTATCCGCAAACTGGCTGGACTTGAAAACGCTGACCAACTTGAAcc  >  1:269980/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                                            
GCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTACTGCCGTCGGCGCTCTCTTTCACCGTATTGGCGTTATCCGCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTGAACG  >  NC_000913/4514048‑4514636

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑