Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,867,454 2 bp→AA coding (97‑98/993 nt) rpoS ← RNA polymerase sigma factor RpoS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,867,4540TA100.0% 351.1 / NA 91Q33M (CAG→TTG) rpoSRNA polymerase sigma factor RpoS
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (45/46);  total (45/46)
*NC_0009132,867,4550GA98.9% 333.6 / NA 91Q33M (CAG→TTG) rpoSRNA polymerase sigma factor RpoS
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (44/45);  total (45/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACT  >  NC_000913/2867216‑2867691
                                                                                                                                                                                                                                               ||                                                                                                                                                                                                                                            
aCGATTGCCATAACGGCGGGCAATTTTTATCACCAGACGCAAGGTTCTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCCCGCAGTGCGCGACGCGCCAAAAAAACTTCTTCTTCGGCCGGTAACAGGGGGGAAAAACCAATCTCACCCAGGGAAAGATGAGTCGCGGCCCACACACCCTGTGTGGCTCCCTGCGGTAACAGTTCCTCTTCGGCCCAATCGGTTTCACTGGG‑GTCCAATTCTACTAAgg                                                                                                                                                                                                                                   >  2:546629/1‑251 (MQ=255)
     tGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCtt                                                                                                                                                                                                                               <  2:301268/250‑1 (MQ=255)
        cATAACGGCGGGCAATTTTTTCCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAAATTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACCCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑GTCCAAATCTACTAAGGCCTTTTct                                                                                                                                                                                                                           >  1:815082/1‑250 (MQ=255)
             cggcggGCAATTTTTACCACCAGACGCAACGTACTCTCGATCATCCGGCGGCGAGATGCGACATCGCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCCCTGGG‑TGCCAAGTCTATTAAGGCCTTTTCGTCaat                                                                                                                                                                                                                      <  2:44313/251‑2 (MQ=255)
               gcggGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCaaaa                                                                                                                                                                                                                     >  2:92605/1‑250 (MQ=255)
                    cAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTc                                                                                                                                                                                                                >  2:82093/1‑250 (MQ=255)
                     aaTTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCaa                                                                                                                                                                                                              <  2:778354/251‑1 (MQ=255)
                        ttttACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTc                                                                                                                                                                                                           >  2:396277/1‑251 (MQ=255)
                            accaccAGACGAAAGTTACTCGCGAGCATCCGGCGGCGAGAGGCGACATCTCCACGCAGGGCGCGAAGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGtt                                                                                                                                                                                                       <  1:478445/251‑1 (MQ=255)
                                 caGATGCAATTTACTCCCGTTCATCCGGCGGCGATAGGCGACATCTCCACGCATTGCGCGACGCGCAAAATAAACGGCTGCGTCGTCCGTTAACAGTGGGGAATAACCTCTCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTTCGATAACAGTTCCTCTTCGGCCAAATCGTTATCGCTGGG‑TTCCAATTCTACTAAGGCCTGTTCGTCAAAAACCTCAACTCCGGTCtcat                                                                                                                                                                                                  <  2:438665/251‑1 (MQ=255)
                                 caGACGCAAGTTACTCTCGATCAGCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtcat                                                                                                                                                                                                  <  2:815082/251‑1 (MQ=255)
                                  aGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtcatc                                                                                                                                                                                                 <  2:310843/251‑1 (MQ=255)
                                   gACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtcatc                                                                                                                                                                                                 >  2:183091/1‑250 (MQ=255)
                                          gTTACTCTCGATAATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCACAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTAACCAAGGCAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGAGAACAGTTCCTCTTCGTCCAACACGTTATCACTGGG‑TTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCATCCTCATCAAATTCCg                                                                                                                                                                                         >  2:385215/1‑251 (MQ=255)
                                           ttACTCTCGATCATCCAGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTCCTAAGGCCTTTTCGTCACAAACCTCAACTCCGTTCTCATCAAATTCCGc                                                                                                                                                                                        >  2:569323/1‑251 (MQ=255)
                                              ctctCGATCATCCGGCGGCAAGAGGCGACATATCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTACAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTACTCAAATTCCGCATc                                                                                                                                                                                     >  2:724033/1‑251 (MQ=255)
                                                ctcGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGGTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCagc                                                                                                                                                                                     <  2:765443/249‑3 (MQ=255)
                                                  cGATCATCCGGCGGCGAGAGGCGACAGCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCTTTCTCATCAAATTCCGCATCtt                                                                                                                                                                                   <  2:820355/249‑1 (MQ=255)
                                                   gATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCAt                                                                                                                                                                                <  1:372565/251‑1 (MQ=255)
                                                       atcCGGCGGCGAGAGGCGACATCTCCCCGCAGTGCGCGACGCGCAAAATAAACGTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATATCAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTaa                                                                                                                                                                            <  1:555767/251‑1 (MQ=255)
                                                          cggcggCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATc                                                                                                                                                                         <  1:773481/251‑1 (MQ=255)
                                                            gcggcgAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCAt                                                                                                                                                                       <  1:601816/251‑1 (MQ=255)
                                                              ggcgAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATg                                                                                                                                                                      >  1:345463/1‑250 (MQ=255)
                                                                 gagaGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAg                                                                                                                                                                                                                                    <  2:626224/185‑1 (MQ=255)
                                                                 gagaGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAg                                                                                                                                                                                                                                    >  1:626224/1‑185 (MQ=255)
                                                                 gagaGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACt                                                                                                                                                                  <  1:517046/251‑1 (MQ=255)
                                                                 gagaGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACt                                                                                                                                                                  >  2:110794/1‑251 (MQ=255)
                                                                     ggCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCGAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTc                                                                                                                                                               >  1:702791/1‑250 (MQ=255)
                                                                        gACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCCTCAAATTCCGCATCTTCATTTAATTCATGAACTTTCAGCg                                                                                                                                                           >  2:146978/1‑251 (MQ=255)
                                                                         aCATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGt                                                                                                                                                          <  1:279720/251‑1 (MQ=255)
                                                                                  cGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACt                                                                                                                                                 >  1:282053/1‑251 (MQ=255)
                                                                                          gcgACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGt                                                                                                                                         <  2:124446/251‑1 (MQ=255)
                                                                                           cgACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTg                                                                                                                                        <  1:190300/251‑1 (MQ=255)
                                                                                              cgcgCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCt                                                                                                                                     >  2:20995/1‑251 (MQ=255)
                                                                                               gcgcAAAATAAACTTCTTCTGCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTc                                                                                                                                    <  1:50904/251‑1 (MQ=255)
                                                                                                  cAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCCGACTCATAAGGTGGCTCCt                                                                                                                                  >  2:496021/1‑250 (MQ=255)
                                                                                                    aaaTAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGAGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAAGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTACATCATGAACTTTCAGCGCATTCTGACTCATACGGGGGCTCCTAcc                                                                                                                               >  2:33705/1‑251 (MQ=255)
                                                                                                     aaTAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTAccc                                                                                                                              >  2:248348/1‑251 (MQ=255)
                                                                                                      aTAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCg                                                                                                                             >  1:665708/1‑251 (MQ=255)
                                                                                                             tcttcttcGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCACCACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCCAAAACCTCAACTCCGTTCTCATCAAATTCCGCCTCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGCTccc                                                                                                                      >  2:186985/1‑251 (MQ=255)
                                                                                                                tcttcGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATccc                                                                                                                      <  1:555989/248‑1 (MQ=255)
                                                                                                                tcttcGGACGGTAACAGGTGTGAATAACCAAGCTCGCCAAGGTAAAGGTGAATCGCGTCCAACACACGCTGTGGGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCATCGTATTCTGACTCATAAGGTGGCCCCTGCCCGTGATCCCtt                                                                                                                    <  2:843099/250‑1 (MQ=255)
                                                                                                                   tcGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTg                                                                                                                   <  1:228859/248‑1 (MQ=255)
                                                                                                                      gCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGaa                                                                                                             <  1:20995/251‑1 (MQ=255)
                                                                                                                        cGTTAACAGTGGTGAATAACCAATCGCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACa                                                                                                           <  1:587559/251‑1 (MQ=255)
                                                                                                                         gTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACAt                                                                                                          >  1:177374/1‑251 (MQ=255)
                                                                                                                         gTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACAt                                                                                                          <  2:183884/251‑1 (MQ=255)
                                                                                                                         gTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACAt                                                                                                          >  1:137240/1‑251 (MQ=255)
                                                                                                                           tAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGGCCAACACACGCTGTGTGGCTGCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACAtt                                                                                                         <  2:32907/250‑1 (MQ=255)
                                                                                                                              cAGTGGTGAATAACCCATCTCACCAAGGTAACGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCACATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCACCTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAGGGTGGCTCCTACCCGTGCTCCCTTGACGGAACATTCa                                                                                                       >  2:173299/1‑249 (MQ=255)
                                                                                                                               aGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGc                                                                                                    >  2:419281/1‑251 (MQ=255)
                                                                                                                               aGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGc                                                                                                    >  1:286644/1‑251 (MQ=255)
                                                                                                                                    tgAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCaaa                                                                                                 >  1:255984/1‑249 (MQ=255)
                                                                                                                                     gAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATCAGGTGGCTCCTACCAGTGATCCCTTGACGGACCATTCAAGCAAAAGc                                                                                              >  2:71907/1‑251 (MQ=255)
                                                                                                                                        tAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTg                                                                                           >  1:130598/1‑251 (MQ=255)
                                                                                                                                          aCCAATCTCACTAAGATACAGCTGAGCCCCGTCTGACACACGCTGTCAGGCTCCGAGAGATAACAGTTCCACTTCAACCAACTCGTTATCACTGGT‑CTCCAATTCTACTAAGGCCTTTTCGTAACAAACCTCACCTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGCACTTCCAACGTATTCTGACTCATAAGGAGGCTCCTACCCGCGATCCCTTGCCGGAACATTCAAGCAAGAGCCTGGt                                                                                         >  2:321576/1‑251 (MQ=255)
                                                                                                                                            cAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTc                                                                                       >  1:627711/1‑251 (MQ=255)
                                                                                                                                             aaTCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTcc                                                                                      <  1:18942/251‑1 (MQ=255)
                                                                                                                                                ctcACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCt                                                                                                                                                                                    >  1:171223/1‑154 (MQ=255)
                                                                                                                                                ctcACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCt                                                                                                                                                                                    <  2:171223/154‑1 (MQ=255)
                                                                                                                                                     ccaGGTAAAGCGGAGTCGCATCCCACACCCGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTACGGCCAAATAGTTATCACTGGG‑ATCCAATTCTACTAAGTCCTTTTCGTCACAAACCTCAACTCCGTTCTCATAAAAATCAGCATCTTCATTTACATCATGAACTCTCAGCGTCGTCTGACTAATAAGTGGGCTCCTACCCGCGATCCCTTGACGGATCATTCACGCAAAAGCCTGGTTCCGCCGAttt                                                                              >  2:29990/3‑251 (MQ=255)
                                                                                                                                                      aaGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGAttt                                                                              >  2:373636/1‑250 (MQ=255)
                                                                                                                                                          tAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGc                                                                         >  1:141364/1‑251 (MQ=255)
                                                                                                                                                            aaGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATACGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTAtagccg                                                                       >  2:11518/1‑246 (MQ=255)
                                                                                                                                                             aGCTGAGGCGCGCCCAAGACACGCTGTGTGGCCCCCAGCGATAACAGTTCCTCTTCGGCCAAATCGTGAGC‑CTGGGGTTCCAATTTTACAAAGGTCTTTGCGTCAAAAAACTCATCTCCGTTCTCATCAAATTCCGCATCTTCAGTAAAATCATGACCTTTCAGCGTATTCTGACTCATAAGGTGCCTCTTCCCCGGGATCCCTTGACGGAACATAAAAACAAAAGCCTGGTTCCGCCGATTTATCGCTGc                                                                      <  1:607008/251‑1 (MQ=255)
                                                                                                                                                                       cggccAACACACGCTGTGTGGCTCCCTGCGATAACATTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAAc                                                            <  2:494602/248‑1 (MQ=255)
                                                                                                                                                                            aaCACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGc                                                       <  2:70309/251‑1 (MQ=255)
                                                                                                                                                                               ccacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCggg                                                    <  2:171586/250‑1 (MQ=255)
                                                                                                                                                                                cacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGt                                                   >  2:345089/1‑251 (MQ=255)
                                                                                                                                                                                cacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGt                                                   <  1:724033/251‑1 (MQ=255)
                                                                                                                                                                                 acGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCCACCCGTGATCCCTTGACGGAACATTCAAGCAAAGGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGtt                                                  <  2:157808/251‑1 (MQ=255)
                                                                                                                                                                                  cgctttttgtctccctccTATACCAGTTCCTCTTCGGCCAATTCGTTCTCCCTGGT‑TTCCAATTCTTCTAACGCCTTTTCTTCAAAAACCTCACCTCCGTTCTCCTCAAATTCCGCATCTTCTTTTAAATCCTGACCTTTCAGCGTATTCTGACTCATAAGGTGCCTCCTCCCCGTGATCCCTTGACTGACCCTTCAAGCAAAAGCCTGTTTCCTCCGATTTATCGCTGCGGCAAATAACGCAGCGGGttt                                                 <  2:548618/241‑1 (MQ=255)
                                                                                                                                                                                  cGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGttt                                                 <  2:498651/251‑1 (MQ=255)
                                                                                                                                                                                  cGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGttt                                                 <  1:322159/251‑1 (MQ=255)
                                                                                                                                                                                      gtgtgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTAc                                               <  1:247554/249‑1 (MQ=255)
                                                                                                                                                                                      gtgtgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACgg                                             <  1:551160/251‑1 (MQ=255)
                                                                                                                                                                                        gtgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCTCTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCTTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACGTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTGGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGa                                            <  2:766488/250‑1 (MQ=255)
                                                                                                                                                                                               ccTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGTGATTACGGATTTCCCCt                                    <  1:835260/251‑1 (MQ=255)
                                                                                                                                                                                                  gCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGt                                 <  2:282053/251‑1 (MQ=255)
                                                                                                                                                                                                     aTAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAAc                              >  1:319026/1‑251 (MQ=255)
                                                                                                                                                                                                       aaCAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTCTCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACg                             >  2:41081/1‑250 (MQ=255)
                                                                                                                                                                                                       aaCAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGa                            <  1:323438/251‑1 (MQ=255)
                                                                                                                                                                                                          aGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAAtt                         >  1:830346/1‑251 (MQ=255)
                                                                                                                                                                                                                 ctTCGGCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCaaa                    <  1:345089/249‑1 (MQ=255)
                                                                                                                                                                                                                     ggCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAACAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCa               >  2:69308/1‑250 (MQ=255)
                                                                                                                                                                                                                      gCCAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCa               >  2:554042/1‑249 (MQ=255)
                                                                                                                                                                                                                        cAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTcc                                                                                      <  1:126329/176‑1 (MQ=255)
                                                                                                                                                                                                                        cAAATCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTcc                                                                                      >  2:126329/1‑176 (MQ=255)
                                                                                                                                                                                                                         aaaTCGTTATCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCgt          >  1:443736/1‑251 (MQ=255)
                                                                                                                                                                                                                                 aTCACTGGG‑TTCCATTTCTTCTAAGGCCTTCTCGTCCAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTGCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCTGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCATGCGTGTTGAAc   <  2:505742/250‑1 (MQ=255)
                                                                                                                                                                                                                                 aTCACTGGG‑TTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACt  >  2:249056/1‑251 (MQ=255)
                                                                                                                                                                                                                                               ||                                                                                                                                                                                                                                            
ACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGG‑TTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACT  >  NC_000913/2867216‑2867691

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