Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,664,8960TG36.7% 155.6 / 15.5 80G169G (GGT→GGG) yfhRputative peptidase
Reads supporting (aligned to +/- strand):  ref base T (9/41);  new base G (28/1);  total (38/42)
Fisher's exact test for biased strand distribution p-value = 1.58e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCCGGAAATATGTCCGCCCACTGGCCGCTGGTCAGTTGGTTACCCGAGCGTAATTTCAACGTTTTTATGTTTGATTATCGCGGGTTTGGTAAATCAAAAGGCACGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTAT‑ATCGCCAGCGTCAGCCCGATCCCGCTTTTACTCATTCACGGTAAAGCTGATCACGTTATCCCATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCG  >  NC_000913/2664651‑2665145
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                           
gCCGGAAATATGTCCGCCCACTGGCCGCTGGTCAGTTGGTTACCCGAGCGTAATTTCAACGTTTTTATGTTTGATTATCGCGGGTTTGGTAAATCAAAAGGCACGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATc                                                                                                                                                                                                                                                         <  2:611977/250‑1 (MQ=255)
               gCCCACTGGCCGCTGGTCAGTTGGTTACCCGAGCGTAATTTCAACGTTTTTATGTTTGATTATCGCGGGTTTGGTAAATCAAAAGGCACGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGc                                                                                                                                                                                                                                         <  2:130122/251‑1 (MQ=255)
                  cACTGGCCGCTGGTCAGTTGGTTACCCGAGCGTAATTTCAACGTTTTTATGTTTGATTATCGCGGGTTTGGTAAATCAAAAGGCACGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGGTGGTTTTTGTTGGGGTGATGGGGAAGGCATACGGGCggg                                                                                                                                                                                                                                      >  1:247693/1‑250 (MQ=255)
                    cTGGCCGCTGGTCAGTTGGTTGCCCGAGCGTAATTTCAACGTTTTTATGTTTGATTATCGCGGGTTTGGTAAATCAAAAGGCACGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGGGGATATTCT‑GGGTGTTATTTGTCGGGGGGGTTGTGGAGGCATACGTGCGGTgg                                                                                                                                                                                                                                    >  2:474420/1‑250 (MQ=255)
                      ggctgCTGGTCAGTTGGTTACCCGAGCGTAATTTCACCGTTTTTATGTTTGATTATCGCGGGTTTGGGCAATGAAAAGGCACGCCGTCCCCGGCCGAATTGCTGGACGATACGCAAAGAGCCATCACTGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTATTCGGGCAGCGTATTGGCGGAGCGGCTATTCTTGGAGGTT‑TTTGTCGGGGGGAGAGTGAAGGGATTCGTGCGGGGATc                                                                                                                                                                                                                                  >  2:338142/5‑251 (MQ=255)
                      ggCCGCTGGTCAGTTGGTTACCCGAGCGTAATTTCAACGTTTTTATGTTTGATTATCGCGGGTTTGGTAAATCAAAAGGCACGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTTTTTGGCGGGGGGGTTGTGAAGGCATAAGTGCGGGGATc                                                                                                                                                                                                                                  >  1:98113/1‑251 (MQ=255)
                      ggCCGCTGGTCAGTTGGTTACCCGAGCGTAATTTCAACGTTTTTATGTTTGATTATCGCGGGTTTGGTAAATCAAAAGGCACGCCGTCCCACGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTTCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑TGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATc                                                                                                                                                                                                                                  <  1:41614/251‑1 (MQ=255)
                       gCCGCTGGTCAGTTGGTTACCCGAGCGTAATTTCAACGTTTTTATGTTTGATTATCGCGGGTTTGGTAAATCAAAAGGCACGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATcc                                                                                                                                                                                                                                 <  1:224605/251‑1 (MQ=255)
                       gCCGCTGGTCAGTTGGTTACCCGAGCGTAATTTCAACGTTTTTATGTTTGATTATCGCGGGTTTGGTAAATCAAAAGGCACGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATcc                                                                                                                                                                                                                                 <  1:26502/251‑1 (MQ=255)
                         cGCTGGTCAGTTGGTTACCCGAGCGTAATTTCAACGTTTTTATGTTTGATTATCGCGGGTTTGGTAAATCAAAAGGCACGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTc                                                                                                                                                                                                                               <  2:345872/251‑1 (MQ=255)
                                          cccGAGCGTAATTTCAACGTTTTTATGTTTGATTATCGCGGGTTTGGTAAATCAAAAGGCACGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGGGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTTCctc                                                                                                                                                                                                              >  1:234014/1‑251 (MQ=255)
                                                    aTTTCAACGTTTTTATGTTTGATTATCGCGGGTTTGGTAAATCAAAAGGCACGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGGGGTCGTGAAGGCATACGTGCGGGGGTTCCCGACTCCCCATTTGCCTCTTAttaaacc                                                                                                                                                                                                    >  2:21274/1‑245 (MQ=255)
                                                    aTTTCAACGTTTTTATGTTTGATTATCGCGGGTTTGGTAAATCAAAAGGCACGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGGTTTTGGTCGGGGGGGTCGTGCAGGCATACGGGCGGGGATCCTCGACTCCACATTTGCCTCTTTTGCAAcc                                                                                                                                                                                                    >  1:61186/1‑251 (MQ=255)
                                                        cAACGTTTTTATGTTTGGTTATCGCGGGTTTGGTAAATCAAAAGGCACGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTTT‑GGGTGTTATTTGTCGGGGGGGTCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTTTGCAACCATcc                                                                                                                                                                                                >  1:410081/1‑250 (MQ=255)
                                                        cAACGTTTTTATGTTTGATTATCGCGGGTTTGGTAAATCAAAAGGCACGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATc                                                                                                                                                                                                 <  2:595167/250‑1 (MQ=255)
                                                        cAACGTTTTTATGTTTGATTATCGCGGGTTTGGTAAATCAAAAGGCACGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑g                                                                                                                                                                                               <  1:474420/251‑1 (MQ=255)
                                                          aCGTTTTTATGTTTGATTATCGCGGGTTTGGTAAATCAAAAGGCACGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGGGATCGTGAAAGCCTAAGCGGGGGGATCCTCGACTCCCCATTTTCCTCTTATTCGCCCATC‑Gcc                                                                                                                                                                                             >  1:288848/1‑251 (MQ=255)
                                                                 tATGTTTGATTATCGCGGGTTTGGTAAATCAAAAGGCACGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑Gccaaccaa                                                                                                                                                                                       <  2:98113/250‑1 (MQ=255)
                                                                     tttGATTATCGCGGGTTTGGTAAATCAAAAGGCACGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGGAAATTCT‑GGATGTTATTTGTCGGGGGGGTCGTGGAGGCATACGGGCGGGGATCCTCGACTCCACATTTTCCTCTTATTTAACCATC‑Gcccaccacacaa                                                                                                                                                                                   >  2:173077/1‑245 (MQ=255)
                                                                        gATTATCGCGGGTTTGGTAAATCAAAAGGCACGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGATcc                                                                                                                                                                                <  1:261877/250‑1 (MQ=255)
                                                                           tATCGCGGGTTTGGTAAATCAAAAGGCACGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTTCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATAATTT‑GGGTGTTATTGGTCGGGGGGATCGTGAAGGGATACGAGGGGTGGTCCTCCAACCCACCTTTTCCTCTTTTTCAACCATC‑GCCACCCAAATGaacacctg                                                                                                                                                                            >  2:667660/1‑244 (MQ=255)
                                                                                cgGGTTTGGTAAATCAAAAGGCACGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGGATATTCT‑GGATGTTTTTTGTCGGGGGGGTCGTGAAGGCATACGTGCGGTGATCCTCGAATCCACATTTTCCTCTTTTGCCACCATC‑CCCCAACAAAAGATCCCCCGCAggg                                                                                                                                                                       >  2:541088/1‑249 (MQ=255)
                                                                                cgGGTTTGGTAAATCAAAAGGCACGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGATCCCCGGCAGTg                                                                                                                                                                       <  1:234591/251‑1 (MQ=255)
                                                                                  ggTTTGGTAAATCAAAAGGCACGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGATCCCCGGCAGTGGc                                                                                                                                                                     <  2:81449/251‑1 (MQ=255)
                                                                                   gTTTGGTAAATCAAAAGGCACGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGACCCCCGGCAGTGGCt                                                                                                                                                                    <  1:460412/251‑1 (MQ=255)
                                                                                    tttGGTAAATCAAAAGGCACGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGGCGGGGGGATCGTGAAGGCATACGTGCGGGGATCCTCGACTCCACATTTTCCTCTTATGCAACCATC‑GCCAACCAAATGATCCCCGCCAGTGGCtc                                                                                                                                                                   >  2:564159/1‑250 (MQ=255)
                                                                                      tGGTAAATCAAAAGGCACGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGGGGATATTTT‑GGGTGTTATTTTTCGGGGGGGTTGTGATCGCATACGTGTGGGGATTCTCGACTCCACATTTTCCTCTTTTTCCACCATCTCCCAACCAA‑TGATccccccccggggccaat                                                                                                                                                                 >  2:46027/1‑238 (MQ=255)
                                                                                          aaaTCAAAAGGCACGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGATCCCCGGCAGTGGCtacttac                                                                                                                                                              <  2:410081/250‑1 (MQ=255)
                                                                                            aTCAAAAGGCACGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTa                                                                                                                                                                                                           <  1:527439/204‑1 (MQ=255)
                                                                                                 aaGGCACGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGGGAATATTCT‑GGATGTTATTTGGCGGGGGGGTCGGGAAGGGATACGTGCGGGGATCCTCGACTCCACATTTTGCTCTTTTTCAACCATC‑GCCCACCAAACGATTCCCGGCAGGGGCCACTTTCTtgatga                                                                                                                                                       >  2:361236/1‑250 (MQ=255)
                                                                                                  aGGCACGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTCTTGGTCGGGGGGAGCGTGAAGGCATATGGGCGGTGATCCTCGACTCCACCTTTGCCTCTTATGCCAACATC‑GCCCACCCAATGATCCCCCGCAGCGGCTCCTTAATTGATgagg                                                                                                                                                     >  1:595648/1‑250 (MQ=255)
                                                                                                   ggCACGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGGATATTCT‑TGACGTTATTGGTCGGGGGGATCGGGAAGGCATACGTGCGGGGATCCTCGACTCCACATTTGCCTCTTATGCAACCACC‑GCCAACCCAACCATCCCCCGCCGTGGCtcatttcctgcttgatg                                                                                                                                                    >  2:162290/1‑235 (MQ=255)
                                                                                                                aGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGGATATTCT‑GGGTGTTTTTGGTCGGGGGGGTCGTGGAGGCATACGGGCGGGGGTTCTCGACTCCACATTTGTCTCTTTTGCAAACATC‑GCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGagaggttaaacggggg                                                                                                                                        >  2:140453/1‑238 (MQ=255)
                                                                                                                aGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTTTTGGTCGGGGGGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCCCATTTGCCTCTTATGCAACCATC‑GCCCACCAAATGATc                                                                                                                                                                                 >  1:565852/1‑209 (MQ=255)
                                                                                                                aGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGATc                                                                                                                                                                                 <  2:565852/209‑1 (MQ=255)
                                                                                                                aGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGGGGTCGTGAAGGCATACGTGCGGGGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCCACCAAATGATCCCCGGCAGGGGCTACTTTCTTGATGAGAGTTACAgcggcgg                                                                                                                                        >  2:49059/1‑249 (MQ=255)
                                                                                                                          gCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGGATATTTTTGGA‑GGTATTTGTCGGGGGGATCGTGGAGGCATAAGGGCGGTGATCCTTGACTCCCCATTTGCCCCTTATGCAACCAAC‑GCCCACCAAATGATCCCCGGCAGTGGCTCCTTACTTTATGAAAGTTACAACGGCGAAAATTAT‑ATCg                                                                                                                            >  2:597548/1‑251 (MQ=255)
                                                                                                                               aCGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGACCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTAT‑ATCGCCAGc                                                                                                                       <  1:580009/251‑1 (MQ=255)
                                                                                                                                       cAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAACCCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTAT‑ATCGCCAGCGTCAGccc                                                                                                               <  2:465055/251‑1 (MQ=255)
                                                                                                                                        aaaGTGCCATCAATGTGGAGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGGATTTTCT‑TGATGTTATTGGTCGGGGGGATCCTGAAGGGATACGTGCGGTGGTCCCCCCCGCCACAATTGCCTCTTTTACAACCCTC‑GCCACCCCAATTACCCCCCGCAGGGGCTACTTTCTTGCTGTGGGTTACACCGGCGAAAACTAT‑ATCGCCaaagggagcccc                                                                                                              >  1:3057/1‑240 (MQ=255)
                                                                                                                                             gCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTAT‑ATCGCCAGCGTCAGCCCGATccc                                                                                                         <  1:140453/251‑1 (MQ=255)
                                                                                                                                             gCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTAT‑ATCGCCAGCGTCAGCCCGATcc                                                                                                          <  1:648044/250‑1 (MQ=255)
                                                                                                                                               cATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTAT‑ATCGCCAGCGTCAGCCCGATCCCGc                                                                                                       <  1:361236/251‑1 (MQ=255)
                                                                                                                                                aTCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCATCCATC‑GCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTAT‑ATCGCCAGCGTCAGCCTGATCCCGCt                                                                                                      <  1:479095/251‑1 (MQ=255)
                                                                                                                                                aTCAATGTGGAGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGGGCTGTTCGGGCAGAGTATTGGCGGGGGGGCTATTCTTGGA‑GTTTTTGGTACGGGAGATCGTGATGGCATACGTGGGGTGAGCCTCCACTCCCCATTTGCCTCTTATTCAACCCTC‑TCCAACCAAATGAACCCCGGCCGTGGCTACCTACCTTATGAGAGTTACATCCGCGAAAAATATTAT‑TCCAGCGTCACCCCGATccccct                                                                                                      >  1:235586/1‑248 (MQ=255)
                                                                                                                                                  cAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGGTGTTATTGGGCGGGGGGGTCGGGGAGGAATACGTGCGGTGATCCCCGACTCCACATTTTCCTCTTTTGCAACCAAC‑GCCAACCAAATGATCCCCGGCAGTGGGTACTTACTTGATTAGAGGTACAACCGCGACAATTAT‑ATCGCCAACGGCAACCCGATCCCGcctt                                                                                                    >  2:624645/1‑248 (MQ=255)
                                                                                                                                                   aaTGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGACCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTAT‑ATCGCCAGCGTCAGCCCGATCCCGCtttt                                                                                                   <  2:288848/251‑1 (MQ=255)
                                                                                                                                                        ggtgCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGGATATTTT‑GGGTGTTATTGGGCGGGGGGGTCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACCTTTGCCTCTTATGCAACCATC‑GCCAACCAAATGATCCCCCGCAGTGGGTACTTACTTGATGAGAGTTACAGCGGCGAAAATTAT‑ATCGCCAGCGGCCGCCCGATCCCGCCTTTtctcc                                                                                              >  1:577923/1‑250 (MQ=255)
                                                                                                                                                           gcgcCATCGCAGTGATGTAACCCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTAT‑ATCGCCAGCGTCAGCCCGATCCCGCTTTTACTCATTc                                                                                           <  1:404968/251‑1 (MQ=255)
                                                                                                                                                               cATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAAAATTCT‑GGGTGTTATTTGTCGGGGGGATCGTGGAGGCCTACGTGCGGGGATCCTCGACTCCCCATTTGCCTCTTTTGCAACCCTC‑CCCAACCAAAAGACCCCCCGCAGTGGCTACTTACTTGATGAGAGTTACAACCGCGAAAATTAT‑ATCGCCCGCGTCCACCCGAACCCCCTTTTACTCATTCCCgg                                                                                       >  2:279606/1‑251 (MQ=255)
                                                                                                                                                                agcGCAGTGATGTAAACCCACAACGTCTGGGGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGTCGAAAATTAT‑ATCGCCAGCGTCAGCCCGATCCCGCTTTTACTCATTCACGGt                                                                                      <  2:235586/249‑1 (MQ=255)
                                                                                                                                                                  cGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGATCCCTGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTAT‑ATCGCCAGCGTCAGCCCGATCCCGCTTTTACTCATTCACGGTa                                                                                     <  2:234014/250‑1 (MQ=255)
                                                                                                                                                                      gTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGCATATTCT‑GGGTGTTTTTGGCCCGGGGGATCGTGAAGGCAAACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTAT‑ATCCCCAGCGTCCGCCCCATCCCCCTTTTACTCATTCACGGGaaaacc                                                                                >  2:193215/1‑248 (MQ=255)
                                                                                                                                                                           gTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTAT‑ATCGCCAGCGTCAGCCCGATCCCGCTTTTACTCATTCACGGTAAAGCTGATc                                                                            <  1:46027/250‑1 (MQ=255)
                                                                                                                                                                               aCCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGTGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTAT‑ATCGCCAGCGTCAGCCCGATCCCGCTTTTACTCATTCACGGTAAAGCTGATCACGtn                                                                       <  2:305114/251‑2 (MQ=255)
                                                                                                                                                                               aCCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGCATATTCT‑GGATGTTATTGGTCGGGGGGATCGGGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTTTGCaa                                                                                                                                                                                                      >  1:465951/1‑126 (MQ=255)
                                                                                                                                                                               aCCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCaa                                                                                                                                                                                                      <  2:465951/126‑1 (MQ=255)
                                                                                                                                                                               aCCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTAT‑ATCGCCAGCGTCAGCCCGATCCCGCTTTTACTCATTCACGGTAAAGCTGATCACGtt                                                                       <  1:21274/251‑1 (MQ=255)
                                                                                                                                                                                 ccACAACGGCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTAT‑ATCGCCAGCGTCAGCCCGATCCCGCTTTTACTCATTCACGGTAAAGCTGATCACGTTAt                                                                     <  2:61186/251‑1 (MQ=255)
                                                                                                                                                                                    caACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGTCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCATCCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAACGGCTAAAATTAT‑ATCGCCAGCGTCATCCCGATCCCCCTTTTACTCATTCCCGGTAAAGCTGATCACGTTATccc                                                                  <  2:679377/251‑1 (MQ=255)
                                                                                                                                                                                        gTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTAT‑ATCGCCAGCGTCAGCCCGATCCCGCTTTTACTCATTCACGGTAAAGCTGATCACGTTATCCCATgg                                                              <  1:564159/251‑1 (MQ=255)
                                                                                                                                                                                            ggTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTAT‑ATCGCCAGCGTCAGCCCGATCCCGCTTTTACTCATTCACGGTAAAGCTGATCACGTTATCCCATGgcagc                                                          <  1:624645/251‑1 (MQ=255)
                                                                                                                                                                                               gCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑Gccaacca                                                                                                                                                                                        <  2:284820/123‑1 (MQ=255)
                                                                                                                                                                                               gCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGGGATCGTGAAGGCATACGTGCGGGGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑Gccaacca                                                                                                                                                                                        >  1:284820/1‑123 (MQ=255)
                                                                                                                                                                                                    gcGGGCAGAGTATTGGCGGGGCGAATAGGCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTAT‑ATCGCCAGCGTCAGCCCGATCCCGCTTTTACTCATACACGGTAAAGCTGATCCCGTTATCCCATAGCAGCACAGCGaa                                                  <  2:49992/250‑1 (MQ=255)
                                                                                                                                                                                                      gggCAGAGTATTGGCGGGGCGAATATTTT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGGTCCTCGACTCCACATTTTCCTCTTATGCAACCATC‑GCCCACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTAT‑ATCGCCAGCGGCAGCCCCATCCCGCTTTTTCTCATTTACGGTAAAACTGGTCACGGTATCCCATTGCAGCACAGCggaa                                                 >  1:375786/1‑247 (MQ=255)
                                                                                                                                                                                                      gggCAGAGTATTGGCGGGGCGAATATTCT‑GGGTGTTATTGGTCGGGGGGGTCGGGAAGGCATACGTGCGGGGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTTATGAGAGTTACAACGGCGGAAATTAT‑ATCGCCAGCGTCAGCCCGATCCCCCTTTTTCTCATTCACGGTAAAGCTGATCACGTTATACCATGGCAGcaccacgaaa                                                 >  1:463127/1‑243 (MQ=255)
                                                                                                                                                                                                         cAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTAT‑ATCGCCAGCGTCAGCCCGATCCCGCTTTTACTCATTCACGGTAAAGCTGATCACGTTATCCCATGGCAGcaca                                                       <  1:529005/241‑1 (MQ=255)
                                                                                                                                                                                                               aTTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTAT‑ATCGCCAGCGTCAGCCCGATCCCGCTTTTACTCATTCACGGTAAAGCTGATCACGTTATCCCATGGCAGCACAGCGAAAAGTTGTATAg                                       <  1:49059/251‑1 (MQ=255)
                                                                                                                                                                                                                 tGGCGGGGCGAATATTCT‑GGATGTTATTGGGCGTGGGGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGACCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTAT‑ATCGCCAGCGTCAGCCCGATCCCGCTTTTACTCATTCACGGTAAAGCTGATCACGTTATCCCATGGCAGCACAGCGAAAAGTTGTATAGc                                      <  2:3057/250‑1 (MQ=255)
                                                                                                                                                                                                                        gCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGCAGGCATACGTGCGGTGATCCTCGACTCCACGTTTGCCTCTTATGCAACAATC‑GCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTAT‑ATCGCCAGCGTCCGCCCGATCCCGCTTTTACTCATTCCCGGTAAAGCTGATCACGTTATCCCATGGCAGCACAGCGCAAAGTTGTATAGCCTTGCaa                               >  2:701217/1‑250 (MQ=255)
                                                                                                                                                                                                                          gAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTAT‑ATCGCCAGCGTCAGCCCGATCCCGCTTTTACTCATTCACGGTAAAGCTGATCACGTTATCCCATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCaaa                              >  2:234216/1‑249 (MQ=255)
                                                                                                                                                                                                                                    ggaagTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTAT‑ATCGCCAGCGTCAGCCCGATCCCGCTTTTACTCATTCACGGTAAAGCTGATCACGTTATCCCATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCaa                      <  1:162290/244‑1 (MQ=255)
                                                                                                                                                                                                                                    ggATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCATATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTAT‑ATCGCCAGCGTCAGCCCGATCCCGCTTTTACTCATTCACGGTAAAGCTGATCACGTTATCCCATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAAc                   <  1:597548/251‑1 (MQ=255)
                                                                                                                                                                                                                                        gTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTAT‑ATCGCCAGCGTCAGCCCGATCCCGCTTTTACTCATTCACGGTAAAGCTGATCACGTTATCCCATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGtt               <  1:173077/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 cGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTAT‑ATCGCCAGCGTCAGCCCGATCCCGCTTTTACTCATTCACGGTAAAGCTGATCACGTTATCCCATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAAt      >  2:60397/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  ggggTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTAT‑ATCGCCAGCGTCAGCCCGATCCCGCTTTTACTCATTCACGGTAAAGCTGATCACGTTATCCCATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATc     >  1:31901/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    ggTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTAT‑ATCGCCAGCGTCAGCCCGATCCCGCTTTTACTCATTCACGGTAAAGCTGATCACGTTATCCCATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATccc   >  2:213531/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    ggTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTAT‑ATCGCCAGCGTCAGCCCGATCCCGCTTTTACTCATTCACGGTAAAGCTGATCACGTTATCCCATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATccc   >  2:296442/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     gTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTAT‑ATCGCCAGCGTCAGCCCGATCCCGCTTTTACTCATTCACGGTAAAGCTGATCACGTTATCCCATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCg  >  1:288105/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                           
GCCGGAAATATGTCCGCCCACTGGCCGCTGGTCAGTTGGTTACCCGAGCGTAATTTCAACGTTTTTATGTTTGATTATCGCGGGTTTGGTAAATCAAAAGGCACGCCGTCCCAGGCCGGATTGCTGGACGATACGCAAAGTGCCATCAATGTGGTGCGCCATCGCAGTGATGTAAACCCACAACGTCTGGTGCTGTTCGGGCAGAGTATTGGCGGGGCGAATATTCT‑GGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATC‑GCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAATTAT‑ATCGCCAGCGTCAGCCCGATCCCGCTTTTACTCATTCACGGTAAAGCTGATCACGTTATCCCATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCG  >  NC_000913/2664651‑2665145

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