Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077793,404,0730TG27.8% 152.3 / 16.8 72G245G (GGT→GGGY75_RS16905oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (9/43);  new base G (20/0);  total (29/43)
Fisher's exact test for biased strand distribution p-value = 3.21e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGCGCTGCTGCATAAGTTGGGTTATCAGGTCGTTGCCGTTTCCGGTCGCGAAAGTACCCATGAATATCTGAAAAGTTTAGGTGCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATTCTATACCCAGGCGGCAAAAGAGATATCTCTGTCAGAGGCACCGAACTTTGCCGAGGCCATCATTAATAACCAGATCCAGGGTCGCACGCTGGTGAAGGTTAACTAACCATTTAGC  >  NC_007779/3403830‑3404322
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                            
ggCGCTGCTGCATAAGTTGGGTTATCAGGTCGTTGCCGTTTCCGGTCGCGAAAGTACCCATGAATATCTGAAAAGTTTAGGTGCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGGTTTAct                                                                                                                                                                                                                                                        >  1:653131/1‑250 (MQ=255)
  cgcTGCTGCATAAGTTGGGTTATCAGGTCGTTGCCGTTTCCGGTCGCGAAAGTACCCATGAATATCTGAAAAGTTTAGGTGCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTActct                                                                                                                                                                                                                                                      <  2:906476/250‑1 (MQ=255)
   gctgctGCATAAGTTGGGTTATCAGGTCGTTGCCGTTTCCGGTCGCGAAAGTACCCATGAATATCTGAAAAGTTTAGGTGCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGTTTTAATCTg                                                                                                                                                                                                                                                     >  1:917298/1‑250 (MQ=255)
    ctgctgCATAAGTTGGGTTATCAGGTCGTTGCCGTTTCCGGTCGCGAAAGTACCCATGAATATCTGAAAAGTTTAGGTGCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGc                                                                                                                                                                                                                                                    >  1:496055/1‑250 (MQ=255)
    ctgctgCATAAGTTGGGTTATCAGGTCGTTGCCGTTTCCGGTCGCGAAAGTACCCATGAATATCTGAAAAGTTTAGGTGCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGGGGCGGCCTGTGGTCTGGCGGGGGGTTTTACTCTGc                                                                                                                                                                                                                                                    >  2:722819/1‑250 (MQ=255)
      gctgcATAAGTTGGGTTATCAGGTCGTTGCCGTTTCCGGTCGCGAAAGTACCCATGAATATCTGAAAAGTTTAGGTGCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTGCCa                                                                                                                                                                                                                                                  >  2:497065/1‑250 (MQ=255)
       ctgcATAAGTTGGGTTATCAGGTCGTTGCCGTTTCCGGTCGCGAAAGTACCCATGAATATCTGAAAAGTTTAGGTGCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTGCCaa                                                                                                                                                                                                                                                 >  1:1001215/1‑250 (MQ=255)
        tgcATAAGTTGGGTTATCAGGTCGTTGCCGTTTCCGGTCGCGAAAGTACCCATGAATATCTGAAAAGTTTAGGTGCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAAc                                                                                                                                                                                                                                                <  2:216782/250‑1 (MQ=255)
              aGTTGGGTTATCAGGTCGTTGCCGTTTCCGGTCGCGAAAGTACCCATGAATATCTGAAAAGTTTAGGTGCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGt                                                                                                                                                                                                                                          <  2:987754/250‑1 (MQ=255)
               gTTGGGTTATCAGGTCGTTGCCGTTTCCGGTCGCGAAAGTACCCATGAATATCTGAAAAGTTTAGGTGCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTc                                                                                                                                                                                                                                         <  2:1065512/250‑1 (MQ=255)
                 tGGGTTATCAGGTCGTTGCCGTTTCCGGTCGCGAAAGTACCCATGAATATCTGAAAAGTTTAGGTGCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCAt                                                                                                                                                                                                                                       <  2:854943/250‑1 (MQ=255)
                 tGGGTTATCAGGTCGTTGCCGTTTCCGGTCGCGAAAGTACCCATGAATATCTGAAAAGTTTAGGTGCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTACTCTGCCAACCACGGTCAt                                                                                                                                                                                                                                       >  1:44679/1‑250 (MQ=255)
                    gTTATCAGGTCGTTGCCGTTTCCGGTCGCGAAAGTACCCATGAATATCTGAAAAGTTTAGGTGCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGcc                                                                                                                                                                                                                                    <  2:526207/250‑1 (MQ=255)
                       aTCAGGTCGTTGCCGTTTCCGGTCGCGAAAGTACCCATGAATATCTGAAAAGTTTAGGTGCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGGTTTACTCTGCCAACCACGGTCATGCCAtt                                                                                                                                                                                                                                 >  1:1118388/1‑250 (MQ=255)
                             tCGTTGCCGTTTCCGGTCGCGAAAGTACCCATGAATATCTGAAAAGTTTAGGTGCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCt                                                                                                                                                                                                                           <  2:702098/250‑1 (MQ=255)
                              cGTTGCCGTTTCCGGTCGCGAAAGTACCCATGAATATCTGAAAAGTTTAGGTGCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTg                                                                                                                                                                                                                          <  1:258869/250‑1 (MQ=255)
                                    cGTTTCCGGTCGCGAAAGTACCCATGAATATCTGAAAAGTTTAGGTGCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAAt                                                                                                                                                                                                                    <  1:950890/250‑1 (MQ=255)
                                     gTTTCCGGTCGCTAAAGTACCCATGAATATCTGAAAAGTTTAGGTGCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATg                                                                                                                                                                                                                   <  2:873354/250‑1 (MQ=255)
                                       ttCCGGTCGCGAAAGTACCCATGAATATCTGAAAAGTTTAGGTGCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTc                                                                                                                                                                                                                 <  2:655315/250‑1 (MQ=255)
                                         ccGGTCGCGAAAGTACCCATGAATATCTGAAAAGTTTAGGTGCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCg                                                                                                                                                                                                               <  2:299600/250‑1 (MQ=255)
                                          cGGTCGCGAAAGTACCCATGAATATCTGAAAAGTTTAGGTGCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGt                                                                                                                                                                                                              <  1:1000171/250‑1 (MQ=255)
                                            gTCGCGAAAGTACCCATGAATATCTGAAAAGTTTAGGTGCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGttt                                                                                                                                                                                                            <  1:10281/250‑1 (MQ=255)
                                             tCGCGAAAGTACCCATGAATATCTGAAAAGTTTAGGTGCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGGGGCGGCCTGTGGTCTGGCGGGGGGGGTTTCTCTGCCCACCACGGTCAAGCCCATTATTTTGCGGTATGTGCGtttt                                                                                                                                                                                                           >  2:891593/1‑249 (MQ=255)
                                                cgAAAGTACCCATGAATATCTGAAAAGTTTAGGTGCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCaa                                                                                                                                                                                                        <  1:418971/250‑1 (MQ=255)
                                                   aaGTACCCATGAATATCTGAAAAGTTTAGGTGCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAggg                                                                                                                                                                                                     <  1:931102/250‑1 (MQ=255)
                                                    aGTACCCATGAATATCTGAAAAGTTTAGGTGCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTGCCAAACACCGTGCTGCCCTTTTTTCTGCGGGATGTGCGTTTGCAAgggg                                                                                                                                                                                                    >  2:440734/1‑250 (MQ=255)
                                                         ccATGAATATCTGAAAAGTTTAGGTGCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGGCACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGGTTTCTCTGCCCAACACGGGCCTGCCATTTTTTCTGCGGCAAGTCCGTTTTCAAGgggggggt                                                                                                                                                                                               >  2:120799/1‑248 (MQ=255)
                                                               atatCTGAAAAGTTTAGGTGCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTa                                                                                                                                                                                        <  2:667947/250‑1 (MQ=255)
                                                                      aaaaGTTTAGGTGCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTACTTTGCCAACCACCGTCACGTCATTTATTTTGCGTAATGTCCGTTTGCAAGGGGTGGATTTAGCAATGACgt                                                                                                                                                                                 >  1:233333/1‑249 (MQ=255)
                                                                             tAGGTGCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGa                                                                                                                                                                          <  1:297200/250‑1 (MQ=255)
                                                                              aGGTGCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTACTTTTCCCACCACGGTCATGCCCTTTTTTCTGCGTAATGTCCTGTTGCGAGGGGGGGATTCAGTAACGGTGCCCCCAGaa                                                                                                                                                                         >  2:830615/1‑250 (MQ=255)
                                                                                gTGCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAAcg                                                                                                                                                                       <  1:918167/250‑1 (MQ=255)
                                                                                 tGCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAAcgc                                                                                                                                                                      <  2:636668/250‑1 (MQ=255)
                                                                                 tGCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGGTTTTCTCTGTCAACCCCGGGCATGCCCTTTTTTCTGGGTAATGTCCGTTTTCAAGGGGGGGAGTCTGATAAGACGCCCCCCGAAagc                                                                                                                                                                      >  2:104661/1‑250 (MQ=255)
                                                                                  gCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGGGGCGGCCTGTGGGCTGGCGGGGGGGTTTTTTCTGTCAACCACGGTCATGCCCTTTTTTTTGCGTTATGTTCGTTTTCAAGGGGGGGGTTTAGTAAAGAACCCACCAGaaaccc                                                                                                                                                                     >  2:258142/1‑246 (MQ=255)
                                                                                    tAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGGGGTCTGGCGGGGGGGTTTTTTCTGCCCACCCCCGCCACGCCACTTTTTTTGCGCTATGCTCGGTTGCAAGGGGGGGTGTTCATGAATACGCCGCCACAAcgcgcc                                                                                                                                                                   >  2:221948/1‑247 (MQ=255)
                                                                                     aGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCcgcg                                                                                                                                                                  <  2:955591/250‑1 (MQ=255)
                                                                                     aGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGCGGCGGCCCGTGGTCTGGCGGGTGGGTTTTCTCTGCCCACCAAGGGGATGCCCTTTATTTTGTGGAATGTTCGTGTGGAAAGGGGGGGGTCAGTAATGACGCCACCACAAAGCcccg                                                                                                                                                                  >  2:845047/1‑250 (MQ=255)
                                                                                           gtTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGGGGCGGCCTGTGGGCTGGCGGGGGGTTTTTACTCTCCCACCCCCGGCCTGCCCTTTTTTTTGCGGTATGTTCCGTTGGAAGGGGGGGGTTCACGAATGACCCCCCCCTAACGCCGCGCACaaa                                                                                                                                                            >  2:530639/1‑249 (MQ=255)
                                                                                                 ccTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGc                                                                                                                                                      <  2:362852/250‑1 (MQ=255)
                                                                                                   tCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTACTCTGCCAACCACGGGCATGCCATTTATTCTGCGGAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAg                                                                                                                                                    >  1:489934/1‑250 (MQ=255)
                                                                                                    cGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGc                                                                                                                                                   <  1:54052/250‑1 (MQ=255)
                                                                                                     gTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACATAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGGTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGGAATGACGCCACCAGAACGCCGCGCACCAGCCTGGCAGCg                                                                                                                                                  >  1:643144/1‑250 (MQ=255)
                                                                                                     gTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGGTTTTCTCTGCCCACCACGGGCATGCCCTTTTTTTTGCGGAATTTTCGTTTTGCAGGGGGGGAGTTCGTAATTACGCCCCCCAAAACCCGCGCCCAACCCCGGCAGCg                                                                                                                                                  >  2:330593/1‑250 (MQ=255)
                                                                                                                 ccGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGa                                                                                                                                      <  2:316613/250‑1 (MQ=255)
                                                                                                                      tCCCTTCCTCTCGAAAAACAGTTCTGGGCTGGGGTAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGTGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAATCCTGGCAGCGACTGGTCTCCGAtttgc                                                                                                                                 <  2:592641/250‑3 (MQ=255)
                                                                                                                            ccTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTACTCTGCCCACCACGGGCCTGCCATTTTTTTTGCGGAATGTCCGTTTGCAAGGGGGGGGATTAGGAATGGAACCCCCAAAAACCCGCGCACAAGCCTGGCAGCGACCGGTGGCCCGTTTAcccggat                                                                                                                           >  2:1110406/1‑245 (MQ=255)
                                                                                                                               ctGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGAGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTGCCAAACACGGTCATGCCATTTATTCTGCGGAATGTCCGTTTGCAAGGGGGGGGTTCAGTAATGACGCCACCAGAACGCCGCGCACAAACCTGGCAGCGACTGGTCGCCGATTTTCCCGAATCAt                                                                                                                        >  2:344118/1‑250 (MQ=255)
                                                                                                                                tGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTGCCAACCACGGGCCTGCCCTTTTTTCTGCGGTATGTCCGCTTGGAAGGGGGGGGTTCAGTAATGACGCCACCAGAACGCCGCGCACAAGGCCGGGCGGGCCTGGTCGCCGGTTTACCGGGATCAtt                                                                                                                       >  2:938620/1‑250 (MQ=255)
                                                                                                                                  gAAAAACAAGTCTGGGCGGGGTCAAGTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGGCATGCCATTTATGCTGCGTAATGGCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCAACAGAACGCCGCGCAAAAGCCTGGCAGCGACTGGTCGCCGAATTACCGGAATCATGCt                                                                                                                     <  2:815873/250‑1 (MQ=255)
                                                                                                                                      aaCAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTACTCTTCCAACCACGGGCATTGCATTTATTCTGCGGAATGTCCCGTTTCAAGGGGGGGGTTCAGTAATGACGCCACCAGAAACCCGCGCACAAGACCGGGAGCGACTGGTCGCCGGCTTTACAGAATCATTTtataa                                                                                                                 >  2:649567/1‑249 (MQ=255)
                                                                                                                                        cAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTTGCAAACACGGGCATGCCCTTTTTTTTGCGGAATGTCCGGTTTCAAGGGGGGGGTTCAGGAATGACGCCCCCAGAACACCGCGCCCAAGCCCGGGAGCGACTGGTGGCCGCTTTACCCGGATCATTCTATAacc                                                                                                               >  2:611363/1‑250 (MQ=255)
                                                                                                                                            tCTGGGCTGGGGCAATTGACCCCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATTCTATACCCAggc                                                                                                          <  1:517196/250‑1 (MQ=255)
                                                                                                                                               gggctgggGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACGGGGCGCTGATTTACCGGAATCATTCTATACCCAggcggc                                                                                                      <  2:648632/250‑1 (MQ=255)
                                                                                                                                                  ctgggGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATTCTATACCCAGGCGGCaaa                                                                                                   <  1:103346/250‑1 (MQ=255)
                                                                                                                                                  ctgggGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTACTCTGCCAACCACGGGCATGCCCTTTATTTTGCGGAATGTCCGTTTGCAAGGGGGGGGTTCAGTAATGACGCCCCCAGAAAGCCGCGCACAAGCCTGGCCGCGGCTGGTGGCCGGTTTTCCGGAAACCTTCTATAACCAGGGGGCcaa                                                                                                   >  2:176690/1‑250 (MQ=255)
                                                                                                                                                    ggggCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATTCTATACCCAGGCGGCAAAag                                                                                                 <  1:497065/250‑1 (MQ=255)
                                                                                                                                                    ggggCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATTCTATACCCAGGCGGCAAAag                                                                                                 <  2:781558/250‑1 (MQ=255)
                                                                                                                                                    ggggCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGGGGCGGCCTGTGGTCTGGCGGGGGGGTTTACTCTGCCAAACACGGGCATGCCATTTATTTTGCGTAATGTCCGTTTGCAAGGGGGGGGTTCAGGAATGACGCCACCAGAAAGCCGCGCACAAGCCCGGGAGCGGAGGGTCGCCGATTTACCGGGATCATTCTATACCCAGGCGGCGAAag                                                                                                 >  2:467097/1‑250 (MQ=255)
                                                                                                                                                       gCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTGCCAACCACGGGCATGCCATTTATTCTGCGGAATGTCCGTTTGCAAGGGGGGGGTTCAGTAATGACGCCACCAGAAAGCCGCGCAAAAGCCTGGCAGCGACTGGTCGCCGATTTAACGGAATCCTTCTAAACCCCAGCGGCCAAaggag                                                                                              >  2:1085930/1‑247 (MQ=255)
                                                                                                                                                              acacCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATTCTATACCCAGGCGGCAAAAGAGATATctct                                                                                       <  2:192352/250‑1 (MQ=255)
                                                                                                                                                                  cGTTGGCGACAAAGTTGTGGCAAAAGTGGTGGGGCAAATGGATTACGGCGGCCGCGTGGCGGCCTGTTGTCTGGCGGGTGGTTTTTCTCTGCCAAACACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCCCCCGAACGGCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATTCTATACCCAGGCGGCAAAAGAGATATCTCTGTCa                                                                                   <  2:44679/250‑1 (MQ=255)
                                                                                                                                                                        cGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATTCTATACCCAGGCGGCAAAAGAGATATCTCTGTCAGAGGCa                                                                             <  1:36582/250‑1 (MQ=255)
                                                                                                                                                                             aaGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATTCTATACCCAGGCGGCAAAAGAGATATCTCTGTCAGAGGCACCGaa                                                                        <  2:908340/250‑1 (MQ=255)
                                                                                                                                                                             aaGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTACTCTGCCCACCACGGGCATTCCCTTTTTTCTGCGTAATGTCCGTTTTCCAGGGGGGGGGTTAGGAATGAAGCCACCCGAAAGCCGCGGCCAAGCCTGGCAGCGACTGGTCGCCGGTTTACCCGGAACATTCTATAACCACGCGGCCAAAAGCAATATCCTCTCTGAGGGAcccga                                                                        >  2:332842/1‑247 (MQ=255)
                                                                                                                                                                              aGTGCTGGCAAAAGTGCTGGCGCAAATGAATTTCGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATTCTATACCCAGGCGGCAAAAGAGATATCTCTGTCAGAGGCACCGAAc                                                                       <  2:37814/250‑1 (MQ=255)
                                                                                                                                                                                        aaaGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGGCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGGTTCAGGAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGGCGCCGATTTACCGGAATCATTCTATACCCAGGCGGCAAAAGAGATATCTCTGTCAGAGGCACCGCACTTTGCCGAgg                                                             >  1:981090/1‑250 (MQ=255)
                                                                                                                                                                                           gTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGTGACTGGTCGCCGATTTACCGGAATCATTCTATACCCAGGCGGCAAAAGAGATATCTCTGTCAGAGGCACCGAACTTTGCCGAGGCca                                                          <  2:524503/250‑1 (MQ=255)
                                                                                                                                                                                                ggCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTACTCTGCCAACCACGGGCATGCCATTTATTCTGGGGAATGTCCGTTTGCAAGGGGGGGGTTCAGTAATGACGCCACCAGAAAGCCGCGCACAAGCCCGGCAGCGACTGGTCGCCGATTTACCGGGATCATTCTATACCCACGCGGGAAAAGAGAAATCTCTGTGAGAGGCACCGAACATTTCCGAGGCCATCAtt                                                     >  2:806201/1‑250 (MQ=255)
                                                                                                                                                                                                  cgcAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGGTTTTCTCTGCCAAACACGGGCATGGCATTTTTTTTGCGGAATGTCCGTTTGCAAGGGGGGGGTTCAGAAATGACGCCACCAGAACGCCGCGCCCAAGACTGGCAGCGAGTGGTCGCCGCGTTACCGGAATCATTATATACCCAAGAGGGCAAAAAGAGATCTCTCTCAGGGGCACAAAACATTGCCGGGGCCATGATcaa                                                   >  2:396272/1‑250 (MQ=255)
                                                                                                                                                                                                     aaaTGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTACTCTGCCAACCACGGTCATGCCATTTTTTTTGCGTAATGTCCGTTTTCAAGGGGGGGGTTCAGGAATTTAGCCACCAGAAAGCCGCGCACCAACCTGGGAGGGACTGGGCGCCGGTTTACCGGAAACAATCTACACCCAGGGGGGCAAAAGAGAATTTCCGGTAAAGGGGACGAACTTTGCCGAGGCCAGCAAtaaaaa                                                >  2:964744/1‑250 (MQ=255)
                                                                                                                                                                                                        tGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATTCTATACCCAGGCGGCAAAAGAGATATCTCTGTCAGAGGCACCGAACTTTGCCGAGGCCATCATTAATAACCa                                             <  2:1118388/250‑1 (MQ=255)
                                                                                                                                                                                                         gAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAGTCATTCTATACCCAGGCGGCAAAAGAGATATCTCTGTCAGAGGCACCGAACTTTGCCGAGGCCATCATTAATAACCAg                                            <  1:325760/250‑1 (MQ=255)
                                                                                                                                                                                                               cggcggCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATTCTATACCCAGGCGGCAAAAGAGATATCTCTGTCAGAGGCACCGAACTTTGCCGAGGCCATCATTAATAACCAGATCCAg                                      <  2:19066/250‑1 (MQ=255)
                                                                                                                                                                                                               cggcggCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATTCTATACCCAGGCGGCAAAAGAGATATCTCTGTCAGAGGCACCGAACTTTGCCGAGGCCATCATTAATAACCAGATCCAg                                      >  2:875853/1‑250 (MQ=255)
                                                                                                                                                                                                                  cggcTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATTCTATACCCAGGCGGCAAAAGAGATATCTCTGTCAGAGGCACCGAACTTTGCCGAGGCCATCATTAATAACCAGATCCAGGGt                                   >  1:451000/1‑250 (MQ=255)
                                                                                                                                                                                                                     cTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATTCTATACCCAGGCGGCAAAAGAGAAATCTCTGTCAGAGGCACCGAACTTTGCCGAGGCCATCATTAATAACCAGATCCAGGGTCGc                                >  1:466234/1‑250 (MQ=255)
                                                                                                                                                                                                                              ggcCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATTCTATACCCAGGCGGCAAAAGAGATATCTCTGTCAGAGGCACCGAACTTTGCCGAGGCCATCATTAATAACCAGATCCAGGGTCGCACGCTGGTg                       <  2:496055/250‑1 (MQ=255)
                                                                                                                                                                                                                              ggcCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATTCTATACCCAGGCGGCAAAAGAGATATCTCTGTCAGAGGCACCGAACTTTGCCGAGGCCATCATTAATAACCAGATCCAGGGTCGCACGCTGGTg                       <  1:938620/250‑1 (MQ=255)
                                                                                                                                                                                                                                     ggTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATTCTATACCCAGGCGGCAAAAGAGATATCTCTGTCAGAGGCACCGAACTTTGCCGAGGCCATCATTAATAACCAGATCCAGGGTCGCACGCTGGTGAAGGTta                <  1:732321/250‑1 (MQ=255)
                                                                                                                                                                                                                                      gTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATTCTATACCCAGGCGGCAAAAGAGATATCTCTGTCAGAGGCACCGAACTTTGCCGAGGCCATCATTAATAACCAGATCCAGGGTCGCACGCTGGTGAAGGTtaa               >  1:831306/1‑250 (MQ=255)
                                                                                                                                                                                                                                          ggCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATTCTATACCCAGGCGGCAAAAGAGATATCTCTGTCAGAGGCACCGAACTTTGCCGAGGCCATCATTAATAACCAGATCCAGGGTCGCACGCTGGTGAAGGTtaactaa           >  1:422015/1‑250 (MQ=255)
                                                                                                                                                                                                                                            cGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATTCTATACCCAGGCGGCAAAAGAGATATCTCTGTCAGAGGCACCGAACTTTGCCGAGGCCATCATTAATAACCAGATCCAGGGTCGCACGCTGGTGAAGGTTAACTAAcc         >  1:452490/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    ttttACTCTTCCAACCACGGTCATGCCATTTATTCTCCGTAATTTCCGTTTGCACTGGGTGGATTCCGTAATGACGCCACCAGAACGCCGCGCACAATCCTTGCCGCGACTGGTCTCTGCTTTACCGTCATCATTCTATACCCATGCGGCAACATTGCTATCTCTGCCAGCTGCCCCGCACTTTCCCGCGGCCATCATTAATAACCAGCTCCAGGGTCGCACGCTCGTGTCGTTTAACTCACCATCTacc  >  2:570516/1‑248 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                            
GGCGCTGCTGCATAAGTTGGGTTATCAGGTCGTTGCCGTTTCCGGTCGCGAAAGTACCCATGAATATCTGAAAAGTTTAGGTGCTAGCCGTGTTCTCCCTCGTGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCATTCTATACCCAGGCGGCAAAAGAGATATCTCTGTCAGAGGCACCGAACTTTGCCGAGGCCATCATTAATAACCAGATCCAGGGTCGCACGCTGGTGAAGGTTAACTAACCATTTAGC  >  NC_007779/3403830‑3404322

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