Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,514,287 G→C L8V (CTT→GTT)  fecB ← ferric citrate ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,514,2870GC97.9% 315.2 / NA 97L8V (CTT→GTT) fecBferric citrate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base C (40/54);  minor base A (2/0);  total (42/55)
Fisher's exact test for biased strand distribution p-value = 1.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAG  >  NC_000913/4514041‑4514532
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
gCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACaaaa                                                                                                                                                                                                                                                   <  1:254642/251‑1 (MQ=255)
 cggcgcgcgcGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTGCGGGCAGGATGCGTGTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGGCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCATCACAATCCGTTGTGGCGTTTTTTCGAGTGGAAACGTGCCGTGTTCGTCCTGAACCGGGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAaaaacaaaac                                                                                                                                                                                                                                                  <  1:267624/249‑1 (MQ=255)
   gcgcgcgcGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTGGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCTCGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCTTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACg                                                                                                                                                                                                                                                 <  2:378714/250‑1 (MQ=255)
      cgcgcgTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGCTCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGAt                                                                                                                                                                                                                                              >  2:250336/1‑250 (MQ=255)
      cgcgcgTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGGTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGGTTTTTCGAGTGGAAACGTGCCGGGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAATCGGATa                                                                                                                                                                                                                                             <  1:78200/251‑1 (MQ=255)
          ggtcccGACGGACTGCCACGGTTTCATGTGCGCTCGCACTTCGGGCAGGATGCGTTTTGCATCGTTACCGTCGGCAATACCGAGCGGGCTGACGTCCACGTCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTAGTCCTGAAACGTGGCGGCAAAGGCGTGGCGGATCGCCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATg                                                                                                                                                                                                                                         <  2:520983/247‑1 (MQ=255)
          cgatCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGAGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGAACACCAGCAGCAGGCCTGCCAAACCAAAACGGATAAATg                                                                                                                                                                                                                                         >  2:464713/4‑251 (MQ=255)
          cgTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCTGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATg                                                                                                                                                                                                                                         <  2:508489/251‑1 (MQ=255)
                gACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACa                                                                                                                                                                                                                                   <  1:127654/251‑1 (MQ=255)
                  gggACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATa                                                                                                                                                                                                                                 <  1:86915/250‑1 (MQ=255)
                   ggACTGCCACGGTTTCAGGTGCGCACGGACTGCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGGTCCAGCACTACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTTATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATaa                                                                                                                                                                                                                                <  1:607996/251‑1 (MQ=255)
                     aCTGCCTCGGTTTCAGGTGCGCACGCACTGCGGGCAGGATGCGTTTGGCATCGTGAGCGTCGGCAATACCGATCTGGCTGACGGCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCTCCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGTCGGCAAATGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCTACATAATc                                                                                                                                                                                                                              <  2:140073/251‑1 (MQ=255)
                     aCTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACAGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGAGTAAACGTGCCGTGTTCGTCCTGAAACGTGGCGGCAAAGGCGTGGCTGATCACCAGAAGCAGTCCTGCAAAAACAAAACGGATAAATGACAACATAATc                                                                                                                                                                                                                              >  1:361904/1‑251 (MQ=255)
                     aCTGCCACGATTTCAGGTGCGCACGCGCTTCGGGCAGGATGCGTTTTGCATCGTTCTCGTCGGCAATACCGAACGGGCTGACGTCCACGGCGGCGAACGCCTCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCCAAGGCGAGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATCAATGCCAACAAAATc                                                                                                                                                                                                                              >  2:80050/1‑251 (MQ=255)
                          cACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGTCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGTCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACAtt                                                                                                                                                                                                                         <  1:438877/251‑1 (MQ=255)
                                   aGGTGCGCACGCACTTCGGACAGGATGCGTCTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGTGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGCGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTaa                                                                                                                                                                                                                 >  2:610651/1‑250 (MQ=255)
                                    ggTGCGCACTCACTTCGGGCAGGATGCGTTGTGCATCGTTATCGTCGGCAATACCGATCGGGCTAACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTaaaa                                                                                                                                                                                                               <  2:482532/251‑1 (MQ=255)
                                    ggTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTaaaa                                                                                                                                                                                                               <  1:171749/251‑1 (MQ=255)
                                    ggTGCGCACGCACTTCGGGCAGGATGCGTTTGGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCAGCGGCGAACGAGAGTTCCAGCGCCACAATCCGTTGTGGCGTTTTTGCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTTCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTaaaa                                                                                                                                                                                                               <  1:610291/251‑1 (MQ=255)
                                          cacgcacTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCgg                                                                                                                                                                                                         >  1:244249/1‑251 (MQ=255)
                                                  tCGGTCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGcc                                                                                                                                                                                                   <  2:338194/249‑1 (MQ=255)
                                                  tCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCg                                                                                                                                                                                                  <  2:130343/250‑1 (MQ=255)
                                                   cGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCgg                                                                                                                                                                                                 <  2:318783/250‑1 (MQ=255)
                                                      gCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGGATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGt                                                                                                                                                                                             <  1:414753/251‑1 (MQ=255)
                                                        aGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGtt                                                                                                                                                                                            >  1:127014/1‑250 (MQ=255)
                                                           aTGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCCCATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAac                                                                                                                                                                                         >  2:35611/1‑250 (MQ=255)
                                                               ggtttGCATCGTTATCGTCGGCAATACCGATCGGGCTTACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGGTCCAGCACCACAATCCGTTGTTGCGTTTTTGCGAGGGTAAACGTGCCGTGGTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGAGCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATTCCAACAGAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATc                                                                                                                                                                                    <  2:257881/249‑1 (MQ=255)
                                                                ttttGCATGGTTATCGTCGGAAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCAGATGTGGCGGTTTTTCGAGTGTAAACGTGCCCTGTTCGTCCTGAACAGGGGCGGCAAAGGCGTGGCTGATAACCGGCAGCAGGCCTGCAAAAACAAAACGCGTAAATGCCAACATATTCACATTCCAGCTAACAGACCGGCAAGCCCGGCGTTAACACATcc                                                                                                                                                                                   >  2:513969/1‑250 (MQ=255)
                                                                     cATCGGTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGCACCGTGGCGGCAAAGGCGTGGCTGATCCCCAGCAGCCGGCCTGCAAAACCAAAACGGATAAATGCCACAATACTCACATTCCATCTAAAATCACGGCAAGGCgggtgtaaaccccgcagaact                                                                                                                                                                              >  1:514319/1‑233 (MQ=255)
                                                                         gTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTGTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACAGTCCAGCTaaaa                                                                                                                                                                                                               <  1:515204/214‑1 (MQ=255)
                                                                         gTTATCGTCAGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCCGCACCACAATCCGTTGTGGCGTTTTTTCGCGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTaaaa                                                                                                                                                                                                               >  2:515204/1‑214 (MQ=255)
                                                                           tATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAAc                                                                                                                                                                         <  2:150891/250‑1 (MQ=255)
                                                                                   gcgaTACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGc                                                                                                                                                                <  1:524643/248‑1 (MQ=255)
                                                                                    cAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCCGCAGCCGGCCTGCAAAAAAAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATAAGAACttacacgaacactgc                                                                                                                                                                >  2:210756/1‑237 (MQ=255)
                                                                                     aaTACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCat                                                                                                                                                              <  1:419154/251‑1 (MQ=255)
                                                                                          cGATCGGGCTGACGTCCACGGCGGCCAGCGCATCTGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGTCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGTCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACAGTCCAGCTAAAAGCCCTGCAAGCCGGGCGTTAACACATCATAACGTCAACTACCTCTGCATATACa                                                                                                                                                         <  2:309866/251‑1 (MQ=255)
                                                                                            aTCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGc                                                                                                                                                       >  1:589727/1‑251 (MQ=255)
                                                                                            aTCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGc                                                                                                                                                       >  2:53184/1‑251 (MQ=255)
                                                                                                 gctGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGATCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCCGGCCTGCAAAAACAAAACGGAAAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGAACACTGCATATACAGCGTgcc                                                                                                                                                  >  2:163887/3‑250 (MQ=255)
                                                                                                    gACGTCCACTGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCATAAGCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCATGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTTCATATACCGCGTgcgt                                                                                                                                                 <  2:377592/249‑2 (MQ=255)
                                                                                                               cggcCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCACTGCATATACAGCGAGCGCGGCTGGCCTGCAt                                                                                                                                    >  1:557277/1‑251 (MQ=255)
                                                                                                                ggcCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCAt                                                                                                                                    >  2:542950/1‑250 (MQ=255)
                                                                                                                   catcccctCGGCGAACGAGAGTTCCAGCCCCACAATCCGTTGGGGCGTTTTTTCCAGTGTAAACGTGTCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGAGCACCAGCATCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCTCATTCCCGCAAAAAGCTCGGCAAGCCGTGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGAt                                                                                                                                <  2:228354/244‑1 (MQ=255)
                                                                                                                        cATCGGCGAACGAGAGTGCCAGCCCCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACTTGCCGTGTTCGTCCTTAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCTGGCCTGCAAAAACAAAACTGATAAATGCCAACATAATCACAGTCCAGCTAAAAGCCCGTTAAGCCGCGCGGTAACGCGTCAGAGCGTCAACTACCCCTGTATATACAACGTGCGCGGCTTGCCCGCATATAGGCCtt                                                                                                                           <  2:88123/251‑1 (MQ=255)
                                                                                                                         aTCGGCGAACGAGAGTTCCAGCACCACAATCCGTGGTGGCGTTTTTTCGAGTGTAAACGTGGCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGGATATACAGCGTGCGCGGCGGGCCTGCAGAGATGCCttt                                                                                                                          <  2:227444/251‑1 (MQ=255)
                                                                                                                           cGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCTTGCGCGGCTGGCCTGCATAGATGCCTttgt                                                                                                                        <  1:180979/251‑1 (MQ=255)
                                                                                                                                aaCGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTctgttgttg                                                                                                                   >  2:355755/1‑251 (MQ=255)
                                                                                                                                 aCGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTtgttgttgt                                                                                                                  <  2:196441/251‑1 (MQ=255)
                                                                                                                                               accacAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCgg                                                                                                    <  1:222675/251‑1 (MQ=255)
                                                                                                                                               aacacaATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAACAACCCAACGGCTAAATGCCAACATACTCACATTCCAGCTAAACGCCCGGCAAGCCGGGCGTTAACACCTCAGAACTTCAACGACCTCAGCATATACAGCGTGCGCGGCTGTCCTGCATAGATGCCTTCGTTGTTGTCGTCATa                                                                                                           >  1:609455/3‑244 (MQ=255)
                                                                                                                                                         gTTGTGGCGTTTTTGCGAGTGTAAACGTGCCGTGTTCGTCCTGATCCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCAGAAGAGCGGATGAagtagt                                                                                          <  2:22826/251‑1 (MQ=255)
                                                                                                                                                          ttGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGCTGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTc                                                                                         >  2:566196/1‑251 (MQ=255)
                                                                                                                                                           tgtgGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGACGTAGTcc                                                                                        >  2:225308/1‑251 (MQ=255)
                                                                                                                                                              ggCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTTCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTg                                                                                      <  2:319455/250‑1 (MQ=255)
                                                                                                                                                              ggCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTg                                                                                      <  1:194085/250‑1 (MQ=255)
                                                                                                                                                               gCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAACACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGAGTTAACCCATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTggg                                                                                    >  1:614614/1‑250 (MQ=255)
                                                                                                                                                                cGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTc                                                                                   <  1:542950/251‑1 (MQ=255)
                                                                                                                                                                    ttttCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGaa                                                                                <  1:250336/250‑1 (MQ=255)
                                                                                                                                                                    ttttCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTGCAACGACCCCTGCATATACAGCGTGCTCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGaa                                                                                <  1:282776/250‑1 (MQ=255)
                                                                                                                                                                     tttcgcgtcaaaaCGTGACGTGTTCGTCCCGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGAAGGCCTGCAAAAAAAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCACGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGCCCTGCATCGATGCCTTTGTTGTTGTCGTCATAAGCACGGATGAAGAAGTCCTGGTCggcg                                                                               >  2:521762/11‑247 (MQ=255)
                                                                                                                                                                            tgtAAACGTGCAGTGTTCGACCTGAACCGTGGCGGCAAAGTCGCGGATGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATAACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTCACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGCGCGCATGACATAGTCCTGGTCGAAGATGttttt                                                                       >  2:384610/1‑251 (MQ=255)
                                                                                                                                                                             gtAAACGTGCCGTGTTCTTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTc                                                                      <  1:2199/251‑1 (MQ=255)
                                                                                                                                                                               aaaCGTGCCGTGTTCGTCCTGAACCTTGGCGGCAAAGGCGTGGCTGATCACAAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAAGCACATTCCAGCTAAAAGCCCGGCAAGCCGTGCGTTAACACATCAGAACTTCAACTACCCCTTCATATACAGCGTGCGCGGCTGGCCTGCAt                                                                                                                                    <  2:377602/187‑1 (MQ=255)
                                                                                                                                                                               aaaCGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCCGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCCCCTTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCCTATAAAGCGTGCGCGGCTGgcatgcat                                                                                                                                    >  1:377602/1‑187 (MQ=255)
                                                                                                                                                                                 aCGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCacac                                                                  >  1:107057/1‑251 (MQ=255)
                                                                                                                                                                                  cGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACAGTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCAGATACAGCGTGCGCGGCTGGCCTGCATATATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACAcc                                                                 <  2:217135/251‑1 (MQ=255)
                                                                                                                                                                                    ggCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGa                                                               <  2:244249/250‑1 (MQ=255)
                                                                                                                                                                                      ccgggtTCGGCCTGAACCGTGGCGGCAAAGGCGTGGCTGAGCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCAAATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCGTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGGCGAAGATGTTTTTCACGCCGAAt                                                             <  2:425783/247‑1 (MQ=255)
                                                                                                                                                                                      ccGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGGGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAAt                                                             <  2:519979/251‑1 (MQ=255)
                                                                                                                                                                                          gtTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCa                                                         <  1:473336/251‑1 (MQ=255)
                                                                                                                                                                                            tCGTCCTGAACCGTGGCGCCCAAGGCGTGGCTGACCACCAGCAGCACGACTGAAACAACAAAACGGATAAATGCCAACAAAATCACATTCCAGCTAAACGCCAGGCACGCCGGGAGTTAACACATCAGAAATTCAAAGACCCCACCACATCCCGCGTGCGCGGCTGGCCTGAATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGACGAAGTAGTCCTGGTCGTAGATGTTTTTCACACCGGCTGCCAgg                                                       >  2:403175/1‑251 (MQ=255)
                                                                                                                                                                                              gTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATACGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCCGGtt                                                     >  2:166365/1‑251 (MQ=255)
                                                                                                                                                                                                ccTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAacac                                                                                                                                                                                       <  1:131757/119‑1 (MQ=255)
                                                                                                                                                                                                ccTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAacac                                                                                                                                                                                       >  2:131757/1‑119 (MQ=255)
                                                                                                                                                                                                ccTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCa                                                   <  1:53184/251‑1 (MQ=255)
                                                                                                                                                                                                 cTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTc                                                    >  2:162218/1‑249 (MQ=255)
                                                                                                                                                                                                   gAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGAt                                                <  1:468114/251‑1 (MQ=255)
                                                                                                                                                                                                     aCCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCt                                              <  1:116620/251‑1 (MQ=255)
                                                                                                                                                                                                       cGTGGCGGCAAAGGCGTTGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGACGTTTTTCCCACCGAATGCCAGGTTCAGATCTGc                                            >  2:497601/1‑251 (MQ=255)
                                                                                                                                                                                                       cGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGc                                            >  1:530962/1‑251 (MQ=255)
                                                                                                                                                                                                        ggggctgcAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGTGCGTTAACACATCAGAACTCCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCGTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTGTCACACCGAATGCCAGGTTCAGATCTGcc                                           <  1:513969/245‑1 (MQ=255)
                                                                                                                                                                                                               cAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCATCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCg                                    <  1:579115/251‑1 (MQ=255)
                                                                                                                                                                                                                aaaGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGCTAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCgg                                   >  1:296643/1‑251 (MQ=255)
                                                                                                                                                                                                                 aaGGCGTTGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCggg                                  <  2:4704/251‑1 (MQ=255)
                                                                                                                                                                                                                            gATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCAt                        >  1:258687/1‑250 (MQ=255)
                                                                                                                                                                                                                               cACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGc                    >  1:393630/1‑251 (MQ=255)
                                                                                                                                                                                                                                  cagcagcagGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCat                 <  2:127014/251‑2 (MQ=255)
                                                                                                                                                                                                                                   agcagcagGCCTGCAAAAACAAAACGGATCAATGCCAACATAATCACATTCAAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTttgt                                                                                                                        >  2:224182/1‑147 (MQ=255)
                                                                                                                                                                                                                                   agcagcagGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTttgt                                                                                                                        <  1:224182/147‑1 (MQ=255)
                                                                                                                                                                                                                                     cagcagGCCTGCAACAACCAACCGGATAAATGCCAACATAATCACATTCCCGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGCCCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCATTTGCTGTTGTCGTCATAAGAGCGGATGAGGTAGTCCTGGTCGAAGCTGTTTTTCACACCGTATGCCAGGTTCAGATCTGACATCTGCGGGCCAAAGTCATACGCCAcgcg              >  1:104946/1‑251 (MQ=255)
                                                                                                                                                                                                                                        ccggCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCCTATACAGCGTGCGCGGCTGGCCTGCATCGATGCCTTTGTTGTTGTCGTCATAAGAGCGGGTGACGTAGTCCTGGTCGAAGATGTTTTTCACACCGCATGCCAGGGTCAGTTCTGCCATCTGCGGGCCAAAGTCATCCGCAACGCGTgc           >  2:260792/3‑251 (MQ=255)
                                                                                                                                                                                                                                        cagGCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTgc           >  1:408177/1‑251 (MQ=255)
                                                                                                                                                                                                                                           gCCTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACTCCACGCGTGCGcc        <  1:497998/251‑1 (MQ=255)
                                                                                                                                                                                                                                             cTGCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCCTAGATGCCTTTGTTGTTGTCGTCATAAGCGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGACTGCCAGGTTCAGATCTGCCATCTGCGCGCCACAGTCATCCGCCACGCGTGCGcccc      >  2:213092/1‑251 (MQ=255)
                                                                                                                                                                                                                                               gCAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCa     <  2:107057/250‑1 (MQ=255)
                                                                                                                                                                                                                                                cAAAAACAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGcccc      >  1:137454/1‑248 (MQ=255)
                                                                                                                                                                                                                                                 aaaaaCAAAACGGATAAATGCCAACATAATCACATTACAGCTAAAAGCCCGGCAAGCCGGGCGTTAACAAATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGATTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCACCag    >  1:380673/1‑249 (MQ=255)
                                                                                                                                                                                                                                                 aaaaaCAAAACGGATAAATGCCAACATAATCACATGCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACGCATCAGAACTTCAACGACCCCGGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCagag  <  2:134774/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
GCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAG  >  NC_000913/4514041‑4514532

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: ‑