Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,748,6340AC42.7% 121.2 / 11.9 75V366G (GTC→GGC) unknownPhosphoglucosamine mutase (EC 5.4.2.10)
Reads supporting (aligned to +/- strand):  ref base A (38/5);  new base C (0/32);  total (38/37)
Fisher's exact test for biased strand distribution p-value = 2.79e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATTAAATCCTAGTTGGCTTCTTTTACTGCTGCTGCAATTTCTTCTGCACATTGCTGTGCCAAAGTGCCATCTTCACATTCCACCATGATGCGGATAAGCGGTTCGGTACCTGATTTACGCAACAAAATTCGCCC‑TTTGCCCGCCAA‑CCGTTTTT‑CCATCTCGGCAGCAA‑CCATTTTCACCGCACCACTTTCCAACGGATTTT‑CACCACCGGCAAAACGTACATTGATTAGCACTT‑GCGGG‑AATAATTTGACCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGACCGGAGTTTTCTCCGCCTAAGCCCCAACCGCGTTCCTGCATTTTTTCCAATACATAACGATCGCCCACATTGGCACGCACAAATGGAATCGCTAATTGTTTTAACGCCAATTCAAGGCTCATGTTGCTC  >  1/1748385‑1748876
                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  
aTTAAATCCTAGTTGGCTTCTTTTACTGCTGCTGCAATTTCTTCTGCACATTGCTGTGCCAAAGTGCCATCTTCACATTCCACCATGATGCGGATAAGCGGTTCGGTACCTGATTTACGCAACAAAATTCGCCC‑TTTGCCCGCCAA‑CCGTTTTT‑CCATCTCGGCAGCAA‑CCATTTTCACCGCACCACTTTCCAACGGATTTT‑CACCACCGGCAAAACGTACATTGATTAGCACTT‑GCGGG‑AATAATTTGa                                                                                                                                                                                                                                                    >  1:461582/1‑250 (MQ=255)
      tCCTAGTTGGCTTCTTTTACTGCTGCTGCAATTTCTTCTGCACATTGCTGTGCCAAAGTGCCATCTTCACATTCCACCATGATGCGGATAAGCGGTTCGGTACCTGATTTACGCAACAAAATTCGCCC‑TTTGCCCGCCAA‑CCGTTTTT‑CCATCTCGGCAGCAA‑CCATTTTCACCGCACCACTTTCCAACGGATTTT‑CACCACCGGCAAAACGTACATTGATTAGCACTT‑GCGGG‑AATAATTTGACCGCAc                                                                                                                                                                                                                                              >  1:692563/1‑250 (MQ=255)
                  tcttTTACTGCTGCTGCAATTTCTTCTGCACATTGCTGTGCCAAAGTGCCATCTTCACATTCCACCATGATGCGGATAAGCGGTTCGGTACCTGATTTACGCAACAAAATTCGCCC‑TTTGCCCGCCAA‑CCGTTTTT‑CCATCTCGGCAGCAA‑CCATTTTCACCGCACCACTTTCCAACGGATTTT‑CACCACCGGCAAAACGTACATTGATTAGCACTT‑GCGGG‑AATAATTTGACCGCACCCGCCAATTCAtt                                                                                                                                                                                                                                 >  1:297857/1‑251 (MQ=255)
                             tgctgcATTTTCTTCTGCACATTGCTGTGCCAAAGTGCCATCTTCACATCCCCCCATGAGGCGGATAAGCGGTTGGGTCCCTGTTTTACGCAACAAAATTCGCCT‑TTTGCCCGCCAC‑CCGTTTTC‑CCATCTCGGCAGCAC‑CCATTTTCCCCGCCCCTTTTTCCAACGGTTTTT‑CCCCACCGGCAAAACGTACATTGATTAGCACTT‑GCGGA‑AATATTTTGCCCGCACCCGCCAATTCATTAAGGGACAAAc                                                                                                                                                                                                                      <  1:323698/251‑1 (MQ=255)
                             tgctgcAATTTCTTCTGCACATTGCTGTGCCAAAGTGCCATCTTCACATTCCACCATGATGCGGATAAGCGGTTCGGTACCTGATTTACGCAACAAAATTCGCCC‑TTTGCCCGCCAA‑CCGTTTTT‑CCATCTCGGCAGCAA‑CCATTTTCACCGCACCACTTTCCAACGGATTTT‑CACCACCGGCAAAACGTACATTGATTAGCACTT‑GCGGG‑AATAATTTGACCGCACCCGCCAATTCATTAAGGGACaaa                                                                                                                                                                                                                       >  1:255271/1‑250 (MQ=255)
                               ctgcATTTTCTTCTGCACATTGCTGTGCCAAAGTGCCATCTTCACTTCCCCCCATGAGGCGGATAAGCGGTTCGTTCCTTGTTTTACGAAAAAAAATTCCCCC‑TTTCCCCGCCAA‑CCTTTTTC‑CCATCTCGGCAGCAC‑CCATTTCCCCCGCCCCCTTTTCCAACGGTTTTC‑CCCCCCCGGCAAAACGTACATTGATTAGCCCTT‑GCGGG‑AATAATTTGCCCGCCCCCGCCAATTCATTAAGGGACAAACGa                                                                                                                                                                                                                    <  1:601241/251‑1 (MQ=255)
                                             gCACATTGCTGTGCCAAAGTGCCATCTTCACATTCCACCATGATGCGGATAAGCGGTTCGGTCCCTGATTTACGAAAAAAAATTCGCCC‑TTTGCCCCCCAC‑CCGTTTTT‑CCATCTCGGCAGCAA‑CCATTTTCCCCGCCCCCCTTTCCAACGGTTTTT‑CCCCCCCGGCAAAACGTACATTGATTAGCACTG‑GCGGG‑AATAATTTGCCCGCCCCCGCCAATTCATTAAGGGACAAACGATGTTTTGc                                                                                                                                                                                                            <  1:502087/245‑1 (MQ=255)
                                             gCACATTGCTGTGCCAAAGTGCCATCTTCACATTCCACCATGATGCGGATAAGCGGTTCGGTACCTGATTTACGCAACAAAATTCGCCC‑TTTGCCCGCCAA‑CCGTTTTT‑CCATCTCGGCAGCAA‑CCATTTTCACCGCACCACTTTCCAACGGATTTT‑CACCACCGGCAAAACGTACATTGATTAGCACTT‑GCGGG‑AATAATTTGACCGCACCCGCCAATTCATTAAGGGACTAACg                                                                                                                                                                                                                     >  1:495404/1‑236 (MQ=255)
                                             gCACATTGCTGTGCCAAAGTGCCATCTTCACATTCCACCATGATGCGGATAAGCGGTTCGGTACCTGATTTACGCAACAAAATTCGCCC‑TTTGCCCGCCAA‑CCGTTTTT‑CCATCTCGGCAGCAA‑CCATTTTCACCGCACCACTTTCCAACGGATTTT‑CACCACCGGCAAAACGTACATTGATTAGCACTT‑GCGGG‑AATAATTTGACCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCg                                                                                                                                                                                                       >  1:528423/1‑250 (MQ=255)
                                                       gtgCCAAAGTGCCATCTTCACATTCCACCATGATGCGGATAAGCGGTTCGGTACCTGATTTACGCAACAAAATTCGCCC‑TTTGCCCGCCAA‑CCGTTTTT‑CCATCTCGGCAGCAA‑CCATTTTCACCGCACCACTTTCCAACGGATTTT‑CACCACCGGCAAAACGTACATTGATTAGCACTT‑GCGGG‑AATAATTTGACCGCACCCGCCAATTCATTaa                                                                                                                                                                                                                               >  1:680560/1‑216 (MQ=255)
                                                       gtgCCAAAGTGCCATCTTCACATTCCACCATGATGCGGATAAGCGGTTCGGTACCTGATTTACGCAACAAAATTCGCCC‑TTTGCCCGCCAA‑CCGTTTTT‑CCATCTCGGCAGCAA‑CCATTTTCACCGCACCACTTTCCAACGGATTTT‑CACCACCGGCAAAACGTACATTGATTAGCACTT‑GCGGG‑AATAATTTGACCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACg                                                                                                                                                                                            >  1:265613/1‑251 (MQ=255)
                                                        tgcaaaaGTGCATTCTTCCCATCCCACCATGATGCGGATAGGGGGTTCGGTACCTGATTTACGCAAAAAAATTCCCCC‑TTTGCCCGCCAC‑CCTTTTTC‑CCATCTCGGCAGCAC‑CCATTTCCCCCGCCCCCCTTTCCAACGGATTTT‑CCCCCCCGGAAAAACGTACATTGATTAGCACTG‑GCGGA‑AATATTTTGCCCCCCCCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACgg                                                                                                                                                                                           <  1:695254/247‑1 (MQ=255)
                                                              aGTGCCATCTTCACATTCCACCATGATGCGGATAAGCGGTTCGGTACCTGATTTACGCAACAAAATTCGCCC‑TTTGCCCGCCAA‑CCGTTTTT‑CCATCTCGGCAGCAA‑CCATTTTCACCGCACCACTTTCCAACGGATTTT‑CACCACCGGCAAAACGTACATTGATTAGCACTT‑GCGGG‑AATAATTTGACCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAagag                                                                                                                                                                                     >  1:144901/1‑251 (MQ=255)
                                                               gTGCCATCTTCACATTCCACCATGATGCGGATAAGCGGTTCGGTACCTGATTTACGCAACAAAATTCGCCC‑TTTGCCCGCCAA‑CCGTTTTT‑CCATCTCGGCAGCAA‑CCATTTTCACCGCACCACTTTCCAACGGATTTT‑CACCACCGGCAAAACGTACATTGATTAGCACTT‑GCGGG‑AATAATTTGACCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGcc                                                                                                                                                                                                  >  1:394737/1‑237 (MQ=255)
                                                                     tCTTCACATTCCACCATGATGCGGATAAGCGGTTCGGTACCTGATTTACGCAACAAAATTCGCCC‑TTTGCCCGCCAA‑CCGTTTTT‑CCATCTCGGCAGCAA‑CCATTTTCACCGCACCACTTTCCAACGGATTTT‑CACCACCGGCAAAACGTACATTGATTAGCACTT‑GCGGG‑AATAATTTGACCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCaat                                                                                                                                                                              >  1:822806/1‑251 (MQ=255)
                                                                      cTTCACATTCCACCATGATGCGGATAAGCGGTTCGGTACCTGATTTACGCAACAAAATTCGCCC‑TTTGCCCGCCAA‑CCGTTTTT‑CCATCTCGGCAGCAA‑CCATTTTCACCGCACCACTTTCCAACGGATTTT‑CACCACCGGCAAAACGTACATTGATTAGCACTT‑GCGGG‑AATAATTTGACCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCaata                                                                                                                                                                             >  1:153358/1‑251 (MQ=255)
                                                                         cacaTCCCCCCAGGATGCGGAAAAGCGGTTCGTTACCTGATCCCCGCAACAAAATTCGCCC‑TTTCCCCCCCAC‑CCGTTTTC‑CCATCCCGGCAGCAA‑CCATTTCCCCCCCCCCACTTTCCAACGGATTTC‑CCCCCCCGGAAAAACGTACATTGATTAGCACTT‑GGGGG‑AATAATTTGACCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATg                                                                                                                                                                          <  1:848077/251‑1 (MQ=255)
                                                                          ccatccccccaTGATGCGAATAGGCGTTCCGGTCCCTGATTTACCCAACAAAATCCCCCC‑TTTCCCCCCCAC‑CCTTTTTC‑CCATCCCGGCACCAA‑CCATTTCCCCCGCCCCCCTTCCCAACGAATTTC‑CCCCCCCGCAAAAACGTACATTGATTAGCACTT‑GGGGA‑AATATTTTGACCGCCCCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGc                                                                                                                                                                         <  1:386430/243‑1 (MQ=255)
                                                                          acaTTCCACCATGATGCGGATAAGCGGTTCGGTACCTGATTTACGCAACAAAATTCGCCC‑TTTGCCCGCCAA‑CCGTTTTT‑CCATCTCGGCAGCAA‑CCATTTTCACCGCACCACTTTCCAACGGATTTT‑CACCACCGGCAAAACGTACATTGATTAGCACTT‑GCGGG‑AATAATTTGACCGCACCCGCCAATTc                                                                                                                                                                                                                                    >  1:67853/1‑192 (MQ=255)
                                                                          acaTTCCACCATGATGCGGATAAGCGGTTCGGTACCTGATTTACGCAACAAAATTCGCCC‑TTTGCCCGCCAA‑CCGTTTTT‑CCATCTCGGCAGCAA‑CCATTTTCACCGCACCACTTTCCAACGGATTTT‑CACCACCGGCAAAACGTACATTGATTAGCACTT‑GCGGG‑AATAATTTGACCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATg                                                                                                                                                                          >  1:362872/1‑250 (MQ=255)
                                                                           cacccccccaTGAGGCGGATAAGCGGTTCGGTCCCTGATTTCCGCAAAAAAATCCCCCC‑TCTCCCCCCCAC‑CCGTTTTC‑CCTTCCCGGCAGCAC‑CCTTTTTCCCCGCCCCTTTTTCCAACGGTTTTC‑CCCCCCCGGAAAAACGTACATTGATTAGCCTTT‑GCGGG‑AATATTTGGCCCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATAc                                                                                                                                                                                             <  1:444897/223‑1 (MQ=255)
                                                                             ttCCACCATGATGCGGATAAGCGGTTCGGTACCTGATTTACGCAACAAAATTCGCCC‑TTTGCCCGCCAA‑CCGTTTTT‑CCATCTCGGCAGCAA‑CCATTTTCACCGCACCACTTTCCAACGGATTTT‑CACCACCGGCAAAACGTACATTGATTAGCACTT‑GCGGG‑AATAATTTGACCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGc                                                                                                                                                                                          >  1:727898/1‑231 (MQ=255)
                                                                             ttCCACCATGATGCGGATAAGCGGTTCGGTACCTGATTTACGCAACAAAATTCGCCC‑TTTGCCCGCCAA‑CCGTTTTT‑CCATCTCGGCAGCAA‑CCATTTTCACCGCACCACTTTCCAACGGATTTT‑CACCACCGGCAAAACGTACATTGATTAGCACTT‑GCGGG‑AATAATTTGACCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAagag                                                                                                                                                                                     >  1:620410/1‑236 (MQ=255)
                                                                               cccccaTGATGCGGATAAGGGGTTGGGTCCCTGATTTACGCAACAAAATTCCCCC‑TTT‑CCCCCCAA‑CCCTTTTTCCCATCGCGGCAGCAA‑CCATTTCCCCCGCACCATTTCCCAACGGATTTC‑CCCCCCCGACAAACCGTACATTGATTAGCACTT‑GCGGA‑AATAATTTGCCCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATcc                                                                                                                                                                    <  1:881159/248‑1 (MQ=255)
                                                                                    atgatgCGGATAAGCGGTTCGTTCCCTGTTTTCCCCAACAAAATCCGCCC‑TTTGCCCCCCAC‑CCTTTTTT‑CCATCTCGGCAGCAA‑CCATTTTCACCGCCCCACTTTCCAACGGTTTTT‑CCCCCCCGGCAAAACGTACATTGATTAGCACTT‑GGGGG‑AATAATTGGCCCGCCCCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGt                                                                                                                                                               <  1:696005/251‑1 (MQ=255)
                                                                                          gggATAAGGGGTTCGGTCCCGGATTTACGAAAAAAAATCCGCCC‑TTTCCCCGCCAC‑CCTTTTTC‑CCATCTCGGCAGCAC‑CCATTTTCACCGCCCCACTTTCCAACGGATTTT‑CCCCCCCGGAAAAACGTACATTGATTAGCACTT‑GCGGA‑AATAATTTGCCCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAgtgt                                                                                                                                                          <  1:263714/249‑1 (MQ=255)
                                                                                            gATAAGCGGTCCGGTACCTGTTTTACGCAAC‑AAATTCCCCCTTTTCCCCGCAAA‑CCTTTTTT‑CCATCTCGGCAGCAA‑CCATTTTCACCGCACCTCTTCCCAACGGATTTT‑CCCCCCCGGCAAAACGTACATTGATTAGCACTT‑GCGGA‑AATAATTTGACCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTg                                                                                                                                                        <  1:378110/250‑1 (MQ=255)
                                                                                                 gCGGTTCGGTCCCTGATTTACGCAAAAAATTTCCCCC‑TTTCCCCGCCAA‑CCGTTTTT‑CCTTCTCGGAAGCAC‑CCATTTTCCCCGCCCCCCTTTCCACCGATTTTT‑CCCCCCCGGCAAAACGTACTTTGATTAGCTTTT‑GCGGG‑AAAATTTTGCCCCCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATc                                                                                                                                                   <  1:590654/250‑1 (MQ=255)
                                                                                                    gTTCGGTCCCTGATTTACGCAACAAATTCCGCCC‑TTTGCCCGCCAC‑CCGTTTTT‑CCATCTCGGCAGCAA‑CCATTTCCCCCGCCCCTCTTCCCAACGGTTTTT‑CCCCCCCGGCAAAACGTACATTGATTAGCCCTT‑GCGGA‑AATAATTTGCCCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGa                                                                                                                                                     <  1:276302/245‑1 (MQ=255)
                                                                                                    gTTCGGTACCTGATTTACGCAACAAAATTCGCCC‑TTTGCCCGCCAA‑CCGTTTTT‑CCATCTCGGCAGTAA‑CCATTTTCACCGCACCACTTTCCAACGGATTTT‑CACCACCGGCAAAACGTACATTGATTAGCACTT‑GCGGG‑AATAATTTGACCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGt                                                                                                                                                               >  1:503501/1‑235 (MQ=255)
                                                                                                    gTTCGGTACCTGATTTACGCAACAAAATTCGCCC‑TTTGCCCGCCAA‑CCGTTTTT‑CCATCTCGGCAGCAA‑CCATTTTCACCGCACCACTTTCCAACGGATTTT‑CACCACCGGCAAAACGTACATTGATTAGCACTT‑GCGGG‑AATAATTTGACCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCa                                                                                                                                               >  1:806186/1‑251 (MQ=255)
                                                                                                     ttCGGTACCTGATTTACGCAAAAAACTTCGCCC‑TTTCCCCGCCAA‑CCTTTTTT‑CCTTCTCGGCAGCAC‑CCTTTTCCCCCGCCCCACTTCCCAACGGTTTTT‑CCCCCCCGGAAAAACGTACATTGATTAGCACTT‑GCGGG‑AATATTTTGCCCGCCCCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCaa                                                                                                                                              <  1:408505/251‑1 (MQ=255)
                                                                                                         gTACCTGATTTACGCAACAAAATTCGCCC‑TTTGCCCGCCAA‑CCGTTTTT‑CCATCTCGGCAGCAA‑CCATTTTCACCGCACCACTTTCCAACGGATTTT‑CACCACCGGCAAAACGTACATTGATTAGCACTT‑GCGGG‑AATAATTTGACCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGat                                                                                                                                          >  1:173063/1‑251 (MQ=255)
                                                                                                               gATTTCCGCAACAAATTTCCCCC‑TTTCCCCGCCAC‑CCTTTTTC‑CCATCTCGGCAGCAA‑CCTTTTTCACCGCCCCACTTTCCAACGGATTTT‑CCCCCCCGGCAAAACGTACATTGATTAGCCCTT‑GCGGA‑AATATTTTGCCCGCCCCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATAc                                                                                                                                                                                             <  1:334183/194‑1 (MQ=255)
                                                                                                                 ttcaCGCAAAAAATTTCGCCC‑TTTCCCCGCCAC‑CCTTTTTT‑CCATCTCGCCA‑CAACCCATTTCCACCCCCCCACTTCCCAGGGGATTTG‑CCCCCCCGGAAAACCGTATATGAATTAGCCTTT‑GCGGT‑AATAATTTGCCCGCCCCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGAc                                                                                                                                  <  1:580680/248‑1 (MQ=255)
                                                                                                                              aatccccc‑TTTCCCCCCCAACCCTTTTTT‑CCATCCCGCCAGCAA‑CCATTTCCCCCCCCCCTCTTTCCAA‑GGTTTTTCCCCCCCCGCCAAAACGTAAATTGTTTAGCACTT‑GGGGA‑AATAATTTGCCCGCCCCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGACCGGAGTTTtctc                                                                                                                      <  1:887610/246‑1 (MQ=255)
                                                                                                                                   ccc‑TTTGCCCGCCAA‑CCGTTTTT‑CCATCTCGGCAGCAA‑CCATTTTCACCGCACCACTTTCCAACGGATTTT‑CACCACCGGCAAAACGTACATTGATTAGCACTT‑GCGGG‑AATAATTTGACCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGACCGGAGTTTTCTCCGCCTa                                                                                                                >  1:207815/1‑251 (MQ=255)
                                                                                                                                                    ccGTTTTT‑CCATCTCGGCAGCAA‑CCATTTTCACCGCACCACTTTCCAACGGATTTT‑CACCACCGGCAAAACGTACATTGATTAGCACTT‑GCGGG‑AATAATTTGACCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGcc                                                                                                                                                                                                  >  1:382335/1‑154 (MQ=255)
                                                                                                                                                    ccGTTTTT‑CCATCTCGGCAGCAA‑CCATTTTCACCGCACCACTTTCCAACGGATTTT‑CACCACCGGCAAAACGTACATTGATTAGCACTT‑GCGGG‑AATAATTTGACCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGACCGGAGTTTTCTCCGCCTAAGCCCCAACCGCGtt                                                                                                 >  1:789962/1‑251 (MQ=255)
                                                                                                                                                               aTCTCGGCAGCAA‑CCATTTTCACCGCACCACTTTCCAACGGATTTT‑CACCACCGGCAAAACGTACATTGATTAGCACTT‑GCGGG‑AATAATTTGACCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGACCGGAGTTTTCTCCGCCTAAGCCCCAACCGCGTTCCTGCAtttt                                                                                       >  1:552205/1‑251 (MQ=255)
                                                                                                                                                                           aaccttttccccccccccccTTTCCACCGGATTTC‑CCCCCCCGACAAACCTTACATTAATTAGCACTT‑GCGGT‑AATATTTTGCCCCCCCCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGACCGGAGTTTTCTCCGCCTAAGCCCCAACCGCGTTCCTGCATTTTTTCCAATACAt                                                                            <  1:748007/232‑1 (MQ=255)
                                                                                                                                                                           a‑CCATTTTCACCGCACCACTTTCCAACGGATTTT‑CACCACCGGCAAAACGTACATTGATTAGCACTT‑GCGGG‑AATAATTTGACCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGACCGGAGTTTTCTCCGCCTAAGCCCCAACCGCGTTCCTGCATTTTTTCCAATACATa                                                                           >  1:83442/1‑251 (MQ=255)
                                                                                                                                                                             ccATTTTCACCGCACCACTTTCCAACGGATTTT‑CACCACCGGCAAAACGTACATTGATTAGCACTT‑GCGGG‑AATAATTTGACCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGc                                                                                                                                                                         >  1:346297/1‑156 (MQ=255)
                                                                                                                                                                             ccATTTTCACCGCACCACTTTCCAACGGATTTT‑CACCACCGGCAAAACGTACATTGATTAGCACTT‑GCGGG‑AATAATTTGACCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGACCGGAGTTTTCTCCGCCTAAGCCCCAACCGCGTTCCTGCATTTTTTCCAATACATaa                                                                          >  1:25412/1‑251 (MQ=255)
                                                                                                                                                                              ctttttCCCCCCACCCCTTTCCAACGGATTTC‑CCCCACCGGAAAACCGTACATTGATTAGCACTT‑GCGGA‑AATATTTTGCCCGCCCCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGACCGGAGTTTTCTCCGCCTAAGCCCCAACCGCGTTCCTGCATTTTTTCCAATACATAAc                                                                         <  1:764161/249‑1 (MQ=255)
                                                                                                                                                                              catttcccccGCCCCACTTCCCACCGTATTTC‑CCCCCCCGGCAAACCGTACATTGATTAGCA‑TT‑GGGGGAAATAATTTGCCCGCCCCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGACCGGAGTTTTCTCCGCCTAAGCCCCAACCGCGTTCCTGCATTTTTTCCAATACATAAc                                                                         <  1:593860/243‑1 (MQ=255)
                                                                                                                                                                                ttttCACCGCCCCACTTCCCAAC‑GATTTTCCCCCCCCGGAAAACCGTACATGGATAAGCACTT‑GGGGA‑AATATTTGGCCCGCCCCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGAAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGACCGGAGTTTTCTCCGCCTAAGCCCCAACCGCGTTCCTGCATTTTTTCCAATACATAACGa                                                                       <  1:545841/251‑1 (MQ=255)
                                                                                                                                                                                tttcccccccaccccTTCCCAACGTATTTT‑CACCACCGGCAAAACGTACATTGATTAGCCTTT‑GCGGA‑AATAATTTACCCGCCCCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGACCGGAGTTTTCTCCGCCTAAGCCCCAACCGCGTTCCTGCATTTTTTCCAATACATAACGa                                                                       <  1:802472/242‑1 (MQ=255)
                                                                                                                                                                                aagagacacccccccttcccAACGGTTTTC‑CCCCACCGGCAAAACGTACATTGATTAGCCCTG‑GGGGA‑AATAATTTGCCCCTCCCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGACCGGAGTTTTCTCCGCCTAAGCCCCAACCGCGTTCCTGCATTTTTTCCAATACATAACGa                                                                       <  1:426980/233‑1 (MQ=255)
                                                                                                                                                                                 tttcccccccccacTTCCCAACGGTTTTT‑CACCCCCGGCAAAACGTAAATTGATTAGCACTT‑GCGGG‑AATATTTTGCCCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGACCGGAGTTTTCTCCGCCTAAGCCCCAACCGCGTTCCTGCATTTTTTCCAATACATAACGAt                                                                      <  1:236489/241‑1 (MQ=255)
                                                                                                                                                                                      cccccccacTTCCCAAGGGATTTC‑CCCCCCCGGAAAAACGTAAATTGATTAGCCCTT‑GGGGG‑AATATTTTGCCCGCCCCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGACCGGAGTTTTCTCCGCCTAAGCCCCAACCGCGTTCCTGCATTTTTTCCAATACATAACGATCGccc                                                                 <  1:817936/246‑1 (MQ=255)
                                                                                                                                                                                      ccGCACCACTTTCCAACGGATTTT‑CACCACCGGCAAAACGTACATTGATTAGCACTT‑GCGGG‑AATAATTTGACCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGcc                                                                                                                                                                                                           >  1:82079/1‑113 (MQ=255)
                                                                                                                                                                                      ccGCACCACTTTCCAACGGATTTT‑CACCACCGGCAAAACGTACATTGATTAGCACTT‑GCGGG‑AATAATTTGACCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGt                                                                                                                                                               >  1:755501/1‑157 (MQ=255)
                                                                                                                                                                                           ccacTTTCCAACGGATTTT‑CACCACCGGCAAAACGTACATTGATTAGCACTT‑GCGGG‑AATAATTTGACCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGACCGGAGTTTTCTCCGCCTAAGCCCCAACcg                                                                                                     >  1:469941/1‑210 (MQ=255)
                                                                                                                                                                                            ccatttcccagcggttttccccccccGGCAAAACGTACATTGATTACCCTTT‑GGGGG‑AATATTTGGCCCCCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGACCGGAGTTTTCTCCGCCTAAGCCCCAACCGCGTTCCTGCATTTTTTCCAATACATAACGATCGCCCACAtt                                                            <  1:693826/227‑1 (MQ=255)
                                                                                                                                                                                                             ttcccccccgaaaaaaCGTAAATTGATTAGCCCTT‑GCGGA‑AATAATTGGCCCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGACCGGAGTTTTCTCCGCCTAAGCCCCAACCGCGTTCCTGCATTTTTTCCAATACATAACGATCGCCCACATTGGCACGCACAAATGGaa                                           <  1:571122/239‑1 (MQ=255)
                                                                                                                                                                                                             tccccccccGGCAAAACGTAAATTGATTAGCTCTG‑GGGGA‑AAAAATTGGACCACCCCCGCCAATTCATTAAGGGACAAACGATGTTTGGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGACCGGAGTTTTCTCCGCCTAAGCCCCAACCGCGTTCCTGCATTTTTTCCAATACATAACGATCGCCCACATTGGCACGCACAAATGGaa                                           <  1:232486/244‑1 (MQ=255)
                                                                                                                                                                                                             tcccccaccGACAAACCGAACATTGATTA‑CCCTTGGCGGA‑AATATTTTGCCCCCCCCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGACCGGAGTTTTCTCCGCCTAAGCCCCAACCGCGTTCCTGCATTTTTTCCAATACATAACGATCGCCCACATTGGCACGCACAAATGGaa                                           <  1:221571/247‑1 (MQ=255)
                                                                                                                                                                                                                      ggCAAAACGTACATTGATTAGCACTT‑GCGGG‑AATAATTTGACCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGACCGGAGTTTTCTCCGCCTAAGCCCCAACCGCGTTCCTGCATTTTTTCCAATACATAACGATCGCCCACATTGGCACGCACAAATGGAATCGCTAAtt                                  >  1:378101/1‑251 (MQ=255)
                                                                                                                                                                                                                           aaCGTACTTTGATTAGCGCTT‑GCGGG‑AATATTTTGCCCCCCCCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGataata                                                                                                                                      <  1:683729/146‑1 (MQ=255)
                                                                                                                                                                                                                           aaCGTACATTGATTAGCACTT‑GCGGG‑AATAATTTGACCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGAc                                                                                                                                  >  1:415059/1‑150 (MQ=255)
                                                                                                                                                                                                                           aaCGTACATTGATTAGCACTT‑GCGGG‑AATAATTTGACCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGACCGGAGTTTTCTCCGCCTAAGCCCCAACcg                                                                                                     >  1:754394/1‑179 (MQ=255)
                                                                                                                                                                                                                           aaCGTACATTGATTAGCACTT‑GCGGG‑AATAATTTGACCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGACCGGAGTTTTCTCCGCCTAAGCCCCAACCGCGTTCCTGCATTTTTTCCAATACATAACGATCGCCCACATTGGCACGCACAAATGGAATCGCTAATTGtttt                             >  1:663253/1‑251 (MQ=255)
                                                                                                                                                                                                                               taaatggATTAGCCCTT‑GGGGA‑AATAATTGGACCGCCCCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGACCGGAGTTTTCTCCGCCTAAGCCCCAACCGCGTTCCTGCATTTTTTCCAATACATAACGATCGCCCACATTGGCACGCACAAATGGAATCGCTAATTGTTTTAACg                         <  1:871854/245‑1 (MQ=255)
                                                                                                                                                                                                                                     gATTAGCACTT‑GCGGG‑AATAATTTGACCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGACCGGAGTTTTCTCCGCCTAAGCCCCAAc                                                                                                       >  1:673875/1‑167 (MQ=255)
                                                                                                                                                                                                                                     gATTAGCACTT‑GCGGA‑AATAATTTGCCCCCCCCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGACCGGAGTTTTCTCCGCCTAAGCCCCAACCGCGTTCCTGCATTTTTTCCAATACATAACGATCGCCcac                                                               <  1:309667/207‑1 (MQ=255)
                                                                                                                                                                                                                                        tAGCACTT‑GCGGG‑AATAATTTGACCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGACCGGAGTTTTCTCCGCCTAAGCCCCAACCGCGTTCCTGCATTTTTTCCAATACATAACGATCGCCCACATTGGCACGCACAAATGGAATCGCTAATTGTTTTAACGCCAATTCaa                >  1:259924/1‑251 (MQ=255)
                                                                                                                                                                                                                                          gCACTT‑GCGGG‑AATAATTTGACCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGACCGGAGTTTTCTCCGCCTAAGCCCCAACCGCGTTCCTGCATTTTTTCCAATACATAACGATCGCCCACATTGGCACGCACAAATGGAATCGCTAATTGTTTTAACGCCAATTCAAgg              >  1:816010/1‑251 (MQ=255)
                                                                                                                                                                                                                                           ctcTT‑GCGGA‑AATAATTCGCCCGCCCCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGACCGGAGTTTTCTCCGCCTAAGCCCCAACCGCGTTCCTGCATTTTTTCCAATACATAACGATCGCCCACATTGGCACGCACAAATGGAATCGCTAATTGTTTTAACGCCAATTCAAgg              <  1:865449/248‑1 (MQ=255)
                                                                                                                                                                                                                                             ccctgggggaaatattTTGCCCGCCCCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGACCGGAGTTTTCTCCGCCTAAGCCCCAACCGCGTTCCTGCATTTTTTCCAATACATAACGATCGCCCACATTGGCACGCACAAATGGAATCGCTAATTGTTTTAACGCCAATTCAAGGc             <  1:463925/241‑1 (MQ=255)
                                                                                                                                                                                                                                               t‑GCGGG‑AATAATTTGACCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGata                                                                                                                                         >  1:128671/1‑123 (MQ=255)
                                                                                                                                                                                                                                                   cggaaataatTTGCCCGCCCCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGACCGGAGTTTTCTCCGCCTAAGCCCCAACCGCGTTCCTGCATTTTTTCCAATACATAACGATCGCCCACATTGGCACGCACAAATGGAATCGCTAATTGTTTTAACGCCAATTCAAGGCTCATGt       <  1:193543/247‑1 (MQ=255)
                                                                                                                                                                                                                                                    ggaaatattTTGCCCCCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGACCGGAGTTTTCTCCGCCTAAGCCCCAACCGCGTTCCTGCATTTTTTCCAATACATAACGATCGCCCACATTGGCACGCACAAATGGAATCGCTAATTGTTTTAACGCCAATTCAAGGCTCATGt       <  1:353515/247‑1 (MQ=255)
                                                                                                                                                                                                                                                        ataatTTGCCCCCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGACCGGAGTTTTCTCCGCCTAAGCCCCAACCGCGTTCCTGCATTTTTTCCAATACATAACGATCGCCCACATTGGCACGCACAAATGGAATCGCTAATTGTTTTAACGCCAATTCAAGGCTCATGTTGCTc  <  1:377965/251‑1 (MQ=255)
                                                                                                                                                                                                                                                         taatTTGACCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGACCGGAGTTTTCTCCGCCTAAGCCCCAACCGCGTTCCTGCATTTTTTCCAATACATAACGa                                                                       >  1:407163/1‑181 (MQ=255)
                                                                                                                                                                                                                                                          aatTTGACCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGACCGGAGTTTTCTCCGCCTAAGCCCCAACCGCGTTCCTGCATTTTTTCCAATACATAACGATCGCCCACATTGGCACGCACAAATGGAATCGCTAATTGTTTTAACGCCAATTCAAGGc             >  1:614497/1‑238 (MQ=255)
                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  
ATTAAATCCTAGTTGGCTTCTTTTACTGCTGCTGCAATTTCTTCTGCACATTGCTGTGCCAAAGTGCCATCTTCACATTCCACCATGATGCGGATAAGCGGTTCGGTACCTGATTTACGCAACAAAATTCGCCC‑TTTGCCCGCCAA‑CCGTTTTT‑CCATCTCGGCAGCAA‑CCATTTTCACCGCACCACTTTCCAACGGATTTT‑CACCACCGGCAAAACGTACATTGATTAGCACTT‑GCGGG‑AATAATTTGACCGCACCCGCCAATTCATTAAGGGACAAACGATGTTTTGCCATCGCAGCCAATACGGCAAGAGAAGCAATAATGCCATCCCCGGTAGTGTTGCGATCGGCAATGATAATATGACCGGAGTTTTCTCCGCCTAAGCCCCAACCGCGTTCCTGCATTTTTTCCAATACATAACGATCGCCCACATTGGCACGCACAAATGGAATCGCTAATTGTTTTAACGCCAATTCAAGGCTCATGTTGCTC  >  1/1748385‑1748876

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