Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,709,1790AC39.1% 89.6 / 12.7 46L54V (TTG→GTG) unknownN‑acetylmannosamine‑6‑phosphate 2‑epimerase (EC 5.1.3.9)
Reads supporting (aligned to +/- strand):  ref base A (27/1);  new base C (0/18);  total (27/19)
Fisher's exact test for biased strand distribution p-value = 6.74e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCGCCATGGCAAGACAGCCCATTTCATGAATTTTTTTAACCGCACTTTCAATCTCTACCGGGCGCGGGCGACGAGTGCCATCCACGGCGATAATGTCCGCTCCGGCGTGCGCCAAGTCTTCAATGTCCTGTAAAAA‑C‑GGGGTGATA‑CGTACCGGGCTGTCGGGAAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGGCGCACGGCTTTTAAATTGTCCACACCTTCAATGCGCAACCCCGCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTGCGGCTGCATTTGGCAGTCATTCTAGTCAATGGAGACATTCTACTAAAATAAAA  >  1/1708933‑1709409
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                      
ccGCCATGGCAAGACAGCCCATTTCATGAATTTTTTTAACCGCACTTTCAATCTCTACCGGGCGCGGGCGACGAGTGCCATCCACGGCGATAATGTCCGCTCCGGCGTGCGCCAAGTCTTCAATGTCCTGTAAAAA‑C‑GGGGTGATA‑CGTACCGGGCTGTCGGGAAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGGCGCACGGCTTTTAAATTGTCCACACCTTCAATGCGCAAccc                                                                                                                                                                                                                                     >  1:307802/1‑250 (MQ=255)
      tGGCAAGACAGCCCATTTCATGAATTTTTTTAACCGCACTTTCAATCTCTACCGGGCGCGGGCGACGAGTGCCATCCACGGCGATAATGTCCGCTCCGGCGTGCGCCAAGTCTTCAATGTCCTGTAAAAA‑C‑GGGGTGATA‑CGTACCGGGCTGTCGGGAAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGGCGCACGGCTTTTAAATTGTCCACACCTTCAATGCGCAAccccgccgc                                                                                                                                                                                                                               >  1:309815/1‑250 (MQ=255)
      tGGCAAGACAGCCCATTTCATGAATTTTTTTAACCGCACTTTCAATCTCTACCGGGCGCGGGCGACGAGTGCCATCCACGGCGATAATGTCCGCTCCGGCGTGCGCCAAGTCTTCAATGTCCTGTAAAAA‑C‑GGGGTGATA‑CGTACCGGGCTGTCGGGAAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGGCGCACGGCTTTTAAATTGTCCACACCTTCAATGCGCAAccccgccgc                                                                                                                                                                                                                               >  1:297658/1‑250 (MQ=255)
           aGACAGCCCATTTCATGAATTTTTTTAACCGCACTTTCAATCTCTACCGGGCGCGGGCGACGAGTGCCATCCACGGCGATAATGTCCGCTCCGGCGTGCGCCAAGTCTTCAATGTCCTGTAAAAA‑C‑GGGGTGATA‑CGTACCGGGCTGTCGGGAAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGGCGCACGGCTTTTAAATTGTCCACACCTTCAATGCGCAAccccgccgccccgc                                                                                                                                                                                                                          >  1:531574/1‑250 (MQ=255)
                gCCCATTTCATGAATTTTTTTAACCGCACTTTCAATCTCTACCGGGCGCGGGCGACGAGTGCCATCCACGGCGATAATGTCCGCTCCGGCGTGCGCCAAGTCTTCAATGTCCTGTAAAAA‑C‑GGGGTGATA‑CGTACCGGGCTGTCGGGAAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGGCGCACGGCTTTTAAATTGTCCACACCTTCAATGCGCAACCCCGCCGCCCCGCCGATaa                                                                                                                                                                                                                    >  1:406630/1‑251 (MQ=255)
                                   ttAACCTCATTTTCAACCTCTCCCGGGCGCGGGCGACGAGTCCCTCCCACGGCGATAATGTCCGCCCGGGCGTGCGCAAAGTCTTCAATGTCCTGTAAAAA‑C‑GGGGTGATC‑CGTCCCGGCCTGTGGGGAAGCTCCCGTTTGACGATCCCGATGACCGGCACATTTCCGGTGGGGCGCAGGGTTTTTAAATTGTCCACCCCTTCAATGCCCACCCCCGCCGCCCCGCCGATAACGGAAGCCTGCGCCATagc                                                                                                                                                                                                 <  1:338110/251‑1 (MQ=255)
                                   taacccGCACTTTCACTCTCTCCCGGGCGCGGGCGACGAGTCCCATCCACGGCGATAATGCCCGCTCCGGCGTGCGCAAAGTCTTCAATGTCCTGAAAAAA‑C‑GGGGTGATA‑CGTACCGGGCTGTCGGGAAGACCCCGTTTGACGATCCCGATGACCGGCACTTTTCCGGTGGGGCGCACGGTTTTTAAATTGCCCACCCCTTCAATGCGCACCCCCGCCGCCCCGCCGATAACGGAAGCCTGCGCCATagc                                                                                                                                                                                                 <  1:127139/247‑1 (MQ=255)
                                              ttCAATCTCTACCGGGCGCGGGCGACGAGTGCCATCCACGGCGATAATGTCCGCTCCGGCGTGCGCCAAGTCTTCAATGTCCTGTAAAAA‑C‑GGGGTGATA‑CGTACCGGGCTGTCGGGAAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGGCGCACGGCTTTTAAATTGTCCACACCTTCAATGCGCAACCCCGCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAAttt                                                                                                                                                                                      >  1:163960/1‑251 (MQ=255)
                                              ttCAATCTCTACCGGGCGCGGGCGACGAGTGCCATCCACGGCGATAATGTCCGCTCCGGCGTGCGCCAAGTCTTCAATGTCCTGTAAAAA‑C‑GGGGTGATA‑CGTACCGGGCTGTCGGGAAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGGCGCACGGCTTTTAAATTGTCCACACCTTCAATGCGCAACCCCGCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAAttt                                                                                                                                                                                      >  1:114210/1‑251 (MQ=255)
                                                   tctcTACCGGGCGCGGGCGACGAGTGCCATCCACGGCGATAATGTCCGCTCCGGCGTGCGCCAAGTCTTCAATGTCCTGTAAAAA‑C‑GGGGTGATA‑CGTACCGGGCTGTCGGGAAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGGCGCACGGCTTTTAAATTGTCCACACCTTCAATGCGCAACCCCGCCGCCCCGCCGACAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGt                                                                                                                                                                                 >  1:132594/1‑251 (MQ=255)
                                                   cccctcccGGGCGGGGGCGACGAGTCCCACCCACGGCGATAATGTCCGCCCCGGCGGGCGCAAATTCTCCAATGTCCTG‑AAAAA‑CGGGGGTGATA‑CGTCCCGGGCTGGGGGGAAGACCCCGTTTGACGATCCCGATGACCGCCACATTTCCGGTCGGGCGCACGGTTTTAAAATTGTCCCCCCCTTAAATGCGCACCCCCGCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGt                                                                                                                                                                                 <  1:313372/245‑1 (MQ=255)
                                                                 gggcgACGAGTGCCATCCACGGCGATAATGTCCGCTCCGGCGTGCGCCAAGTCTTCAATGTCCTGTAAAAA‑C‑GGGGTGATA‑CGTACCGGGCTGTCGGGAAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGGCGCACGGCTTTTAAATTGTCCACACCTTCAATGCGCAACCCCGCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGAc                                                                                                                                                                    >  1:79716/1‑250 (MQ=255)
                                                                      aggaGTCCCACCCAGGGCGATAATGCCCGCCCGGGCGTGCGCCAAGTTTTCAATGTCCT‑TAAAAA‑CGGGGGTGATA‑CGTACGGGGCTGTCGGGAAGATCCCTTTTGACGATCCCGATGACCGGCACATTTCCGGGCGGGCGCACGGTTTTTAAATTGCCCCCCCCTTAAATGCGCACCCCCGCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACcatca                                                                                                                                                               <  1:511414/248‑1 (MQ=255)
                                                                          gTGCCATCCACGGCGATAATGTCCGCTCCGGCGTGCGCCAAGTCTTCAATGTCCTGTAAAAA‑C‑GGGGTGATA‑CGTACCGGGCTGTCGGGAAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGGCGCACGGCTTTTAAATTGTCCACACCTTCAATGCGCAACCCCGCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCa                                                                                                                                                           >  1:18533/1‑250 (MQ=255)
                                                                           cgccacccgcgggaaTAATGCCCCCTCCGGCGGGCCCAAAGTCTTCATTTTCCTGTAAAAC‑C‑GGGGTGATA‑CGTCCGGGGCTTGCGGGAAGCTCCCTTTTGACAATCCCGATGCCCGGCACTTTTACGGTGGGGCGCACGGTTTTTAATTTGTCCCCCCCTTCATTGCCCACCCCCGCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCAc                                                                                                                                                          <  1:594410/236‑1 (MQ=255)
                                                                                      gCGATAATGTCCGCTCCGGCGTGCGCCAAGTCTTCAATGTCCTGTAAAAA‑C‑GGGGTGATA‑CGTACCGGGCTGTCGGGAAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGGCGCACGGCTTTTAAATTGTCCACACCTTCAATGCGCAACCCCGCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAaga                                                                                                                                              >  1:339153/1‑251 (MQ=255)
                                                                                      gCGATAATGTCCGCTCCGGCGTGCGCCAAGTCTTCAATGTCCTGAAAAAA‑C‑GGGGTGATA‑CGTCCCGGGCTGTGGGGAAGATCCCGTTTGCCGATACCGATGATCGGCACATTTCCGGTCGGGCGCCCGGTTTTTAAATTGCCCACCCCTTCAATGCGCACCCCCGCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAag                                                                                                                                               <  1:133700/250‑1 (MQ=255)
                                                                                         aTAATGTCCGCTCCGGCGTGCGCCAAGTCTTCAATGTCCTGTAAAAA‑C‑GGGGTGATA‑CGTACCGGGCTGTCGGGAAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGGCGCACGGCTTTTAAATTGTCCACACCTTCAATGCGCAACCCCGCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAgg                                                                                                                                            >  1:466906/1‑250 (MQ=255)
                                                                                                 cGCTCGGGCGTGCGACAATTCTTCATTGTCCTATAAAAG‑C‑GGGGTGATA‑CGTACCGGGCTTTGGGGAATACCCCGTTTGACGATACCGATGCTCGGCACTTTTCCGGTCGGGCGCGCGGTTTTTAAATTGTCCCCCCCTTCAATGCGCACCCCCGCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATc                                                                                                                                   <  1:178169/251‑1 (MQ=255)
                                                                                                            cccccAAGTCTTAAATGTCCT‑TAAAAACC‑GGGGTG‑TACCGTCCGGGGCTGTCGGGAAGATCCCGTTTGCCGATCCCGAGGACCGGCACTTTTACGGTCGGGCGCCCGGTTTTAAAATTGCCCCCCCCTCCAAGGCCCACCCCCGCCCCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAg                                                                                                                         <  1:79821/247‑1 (MQ=255)
                                                                                                                                            gggTGATA‑CGTACCGGGCTGTCGGGAAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGGCGCACGGCTTTTAAATTGTCCACACCTTCAATGCGCAACCCCGCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAAtttt                                                                                          >  1:501619/1‑251 (MQ=255)
                                                                                                                                                gATA‑CGTACCGGGCTGTCGGGAAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGGCGCACGGCTTTTAAATTGTCCACACCTTCAATGCGCAACCCCGCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTgag                                                                                      >  1:370186/1‑251 (MQ=255)
                                                                                                                                                      gTACCGGGCTGTCGGGAAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGGCGCACGGCTTTTAAATTGTCCACACCTTCAATGCGCAACCCCGCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAg                                                                                  >  1:242431/1‑250 (MQ=255)
                                                                                                                                                                gccGGGAAGACCCCTTTTGACGATCCCGATGATCGGCACATTTCCGGTCGGGCGCCCGGCTTTAAATTTGCCCACCCCTTCAATGCGCACCCCCCCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTtgatg                                                                                                    <  1:413444/220‑1 (MQ=255)
                                                                                                                                                                 tCGGGAAGATCCCGTTTGACAATACCGATGATCGGCTCATTTACGGGGGGGCGCAGGCCTTTTAAATTGTCCACCCCTTAAATGCGCACCCCCGCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAg                                                                       <  1:406742/250‑1 (MQ=255)
                                                                                                                                                                        gATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGGCGCACGGCTTTTAAATTGTCCACACCTTCAATGCGCAACCCCGCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAActc                                                                 >  1:140030/1‑249 (MQ=255)
                                                                                                                                                                                    cGATACCGATGATCGGCACATTTACGGTCGGGCGCACGGCTTTTAAATTGTCCACACCTTCAATGCGCAACCCCGCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAActc                                                                 >  1:77860/1‑237 (MQ=255)
                                                                                                                                                                                    cGATACCGATGACCGGCACATTTGCGGTGGGGCGCCCGGTTTTTAAATTGCCCCCCCCTTCAATGCCCACCCCCGCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTGCGGc                                                   <  1:199503/251‑1 (MQ=255)
                                                                                                                                                                                                      cttttgcggtggggcgccgggtttttAAATTGTCCACACCTTCAATGCCCACCCCCGCCGTCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTGCGGCTGCATTTGGCAGTCATTc                                 <  1:549506/229‑1 (MQ=255)
                                                                                                                                                                                                      caTTTACGGTCGGGCGCACGGCTTTTAAATTGTCCACACCTTCAATGCGCAACCCCGCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCttttt                                                         >  1:518667/1‑227 (MQ=255)
                                                                                                                                                                                                      caTTTACGGTCGGGCGCACGGCTTTTAAATTGTCCACACCTTCAATGCGCAACCCCGCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTGCGGCTGCATTTGGCAGTCAtt                                  >  1:520836/1‑250 (MQ=255)
                                                                                                                                                                                                       aTTTACGGTCGGGCGCACGGCTTTTAAATTGTCCACACCTTCAATGCGCAACCCCGCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTGCGGCTGCATTTGGCAGTCAtt                                  >  1:153676/1‑249 (MQ=255)
                                                                                                                                                                                                           cgggTCGGGCGCACGGTTTTAAAATTGCCCCCACCTCCAATGCGCACCCCCGCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTGCGGCTGCATTTGGCAGTCATTCTAGTc                            <  1:449847/249‑1 (MQ=255)
                                                                                                                                                                                                           ccgggcGGGCGCCCGGTTTTTAATTTGCCCCCCCCTTCAATGCCCACCCCCCCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTGCGGCTGCATTTGGCAGTCATTCTAGTc                            <  1:108741/246‑1 (MQ=255)
                                                                                                                                                                                                           ccGGTGGGGCGCACGGTTTTTAAATGACCCCCCCCTAAAATGCGCAACCCCGCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTGCGGCTGCATTTGGCAGTCATTCTAGTc                            <  1:64573/250‑1 (MQ=255)
                                                                                                                                                                                                           ccGGTCGGGCGCCCGGCTTTAAAATTGCCCCCCCCTTCAAGGCGCACCCCCGCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTGCGGCTGCATTTGGCAGTCATTCTAGTc                            <  1:495702/250‑1 (MQ=255)
                                                                                                                                                                                                            cGGTCGGGCGCACGGCTTTTAAATTGTCCACACCTTCAATGCGCAACCCCGCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTGCGGCTGCATTTGGCAGTCATTCTAGTCa                           >  1:8781/1‑251 (MQ=255)
                                                                                                                                                                                                                cGGGCGCACGGCTTTTAAATTGTCCACACCTTCAATGCGCAACCCCGCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTGCGGCTGCATTTGGCAGTCATTCTAGTCAATgg                       >  1:406998/1‑251 (MQ=255)
                                                                                                                                                                                                                          gCTTTTAAATTGTCCACACCTTCAATGCGCAACCCCGCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTGCGGCTGCATTTGGCAGTCATTCTAGTc                            >  1:154240/1‑236 (MQ=255)
                                                                                                                                                                                                                          gCTTTTAAATTGTCCACACCTTCAATGCGCAACCCCGCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTGCGGCTGCATTTGGCAGTCATTCTAGTCAATGGAGACATTcta             >  1:253367/1‑251 (MQ=255)
                                                                                                                                                                                                                           tttttAATTTGTCCACCCCTTAAAGGCCCACCCCCGCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTGCGGCTGCATTTGGCAGTCATTCTAGTCAATGGAGACATTctac            <  1:354706/250‑1 (MQ=255)
                                                                                                                                                                                                                            ttttAAATTGCCCCCACCTTCAATGCGCACCCCCCCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTGCGGCTGCATTTGGCAGTCATTCTAGTCAATGGAGACATTctact           <  1:74599/251‑1 (MQ=255)
                                                                                                                                                                                                                            tttaaaaTTGTCCCCCCCTTCAATGCGCACCCCCCCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTGCGGCTGCATTTGGCAGTCATTCTAGTCAATGGAGACATTctact           <  1:537226/247‑1 (MQ=255)
                                                                                                                                                                                                                             tttAAATTGTCCACACCTTCAATGCGCAACCCCGCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTGCg                                                     >  1:576206/1‑208 (MQ=255)
                                                                                                                                                                                                                                     gTCCACACCTTCAATGCGCAACCCCGCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGTACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTGCGGCTGCATTTGGCAGTCATTCTAGTCAATGGAGACATTCTACtaaaataaaa  >  1:350386/1‑251 (MQ=255)
                                                                                                                                                                                                                                            ccTTCAATGCGCAACCCCGCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTGCGGCTGCATTTGGCAGTCATTc                                 >  1:340614/1‑213 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                      
CCGCCATGGCAAGACAGCCCATTTCATGAATTTTTTTAACCGCACTTTCAATCTCTACCGGGCGCGGGCGACGAGTGCCATCCACGGCGATAATGTCCGCTCCGGCGTGCGCCAAGTCTTCAATGTCCTGTAAAAA‑C‑GGGGTGATA‑CGTACCGGGCTGTCGGGAAGATCCCGTTTGACGATACCGATGATCGGCACATTTACGGTCGGGCGCACGGCTTTTAAATTGTCCACACCTTCAATGCGCAACCCCGCCGCCCCGCCGATAACGGAAGCCTGCGCCATAGCAGCGACAATTTCCGGTTTATCCATCGGACCATCATCCACCGGCTGACAAGAGGCAATTAATCCGTTCTGAAGTTGGGCGAAAACGGCTTGATGAGATAATTTTTGAGATAGCTGTGACATAGAAACTCCTCTTTTTTGCGGCTGCATTTGGCAGTCATTCTAGTCAATGGAGACATTCTACTAAAATAAAA  >  1/1708933‑1709409

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