Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1188,0800AC43.3% 92.5 / 14.9 68V80G (GTC→GGC) unknown4‑hydroxy‑2‑oxoglutarate aldolase (EC 4.1.3.16) @ 2‑dehydro‑3‑deoxyphosphogluconate aldolase (EC 4.1.2.14)
Reads supporting (aligned to +/- strand):  ref base A (26/12);  new base C (0/29);  total (26/42)
Fisher's exact test for biased strand distribution p-value = 1.03e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGCCGGTCGGCATAATTTGCAACTGTGCGTAAGGGCCTAACAAGGCTTTAATCATTTTCACGCCACCAGAGGCTTCCGCCGGGAAGAATTTCACGGTTTCAATACCGACGTCCAACGCGGCTTCAATTGCCATCGGATTGTTTACCCCCGGTGTGATGGGTAACGCTAAATCCTGACATAATTGCACGATTTTA‑GGATTAAATCCCGGTGTGACGACCACATCTGCACCGCTGTTTTTA‑GCGGCGACGACCTGTTCGGCGGTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAATATCTTCCGCGTTGTCTAAAGCAATCACCGGCACGACTTTTAATTGACGAAATTTTTCAATGAGTTGGTTAGTTGTGTATGACATAAATC  >  1/187832‑188323
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    
cgccgGTCGGCATAATTTGCAACTGTGCGTAAGGGCCTAACAAGGCTTTAATCATTTTCACGCCACCAGAGGCTTCCGCCGGGAAGAATTTCACGGTTTCAATACCGACGTCCAACGCGGCTTCAATTGCCATCGGATTGTTTACCCCCGGTGTGATGGGTAACGCTAAATCCTGACATAATTGCACGATTTTA‑GGATTAAATCCCGGTGTGACGACCACATCTGCACCGCTGTTTTTA‑GCGGCGACGAcc                                                                                                                                                                                                                                                    <  1:720187/251‑1 (MQ=255)
    ggTCGGCATAATTTGCAACTGTGCGTAAGGGCCTAACAAGGCTTTAATCATTTTCACGCCACCAGAGGCTTCCGCCGGGAAGAATTTCACGGTTTCAATACCGACGTCCAACGCGGCTTCAATTGCCATCGGATTGTTTACCCCCGGTGTGATGGGTAACGCTAAATCCTGACATAATTGCACGATTTTA‑GGATTAAATCCCGGTGTGACGACCACATCTGCACCGCTGTTTTTA‑GCGGCGACGACCTGtt                                                                                                                                                                                                                                                >  1:684993/1‑251 (MQ=255)
      tCGGCATAATTTGCAACTGTGCGTAAGGGCCTAACAAGGCTTTAATCATTTTCACGCCACCAGAGGCTTCCGCCGGGAAGAATTTCACGGTTTCAATACCGACGTCCAACGCGGCTTCAATTGCCATCGGATTGTTTACCCCCGGTGTGATGGGTAACGCTAAATCCTGACATAATTGCACGATTTTA‑GGATTAAATCCCGGTGTGACGACCACATCTGCACCGCTGTTTTTA‑GCGGCGACGACCTGTTc                                                                                                                                                                                                                                               <  1:780119/250‑1 (MQ=255)
              aTTTGCAACTGTGCGTAAGGGCCTAACAAGGCTTTAATCATTTTCACGCCACCAGAGGCTTCCGCCGGGAAGAATTTCACGGTTTCAATACCGACGTCCAACGCGGCTTCAATTGCCATCGGATTGTTTACCCCCGGTGTGATGGGTAACGCTAAATCCTGACATAATTGCACGATTTTA‑GGATTAAATCCCGGTGTGACGACCACATCTGCACCGCTGTTTTTA‑GCGGCGACGACCTGTTCGGCGGTCaa                                                                                                                                                                                                                                      >  1:110613/1‑251 (MQ=255)
                   cAACTGTGCGTAAGGGCCTAACAAGGCTTTAATCATTTTCACGCCACCAGAGGCTTCCGCCGGGAAGAATTTCACGGTTTCAATACCGACGTCCAACGCGGCTTCAATTGCCATCGGATTGTTTACCCCCGGTGTGATGGGTAACGCTAAATCCTGACATAATTGCACGATTTTA‑GGATTAAATCCCGGTGTGACGACCACATCTGCACCGCTGTTTTTA‑GCGGCGACGACCTGTTCGGCGGTCAATACgg                                                                                                                                                                                                                                 <  1:860565/251‑1 (MQ=255)
                   cAACTGTGCGTAAGGGCCTAACAAGGCTTTAATCATTTTCACGCCACCAGAGGCTTCCGCCGGGAAGAATTTCACGGTTTCAATACCGACGTCCAACGCGGCTTCAATTGCCATCGGATTGTTTACCCCCGGTGTGATGGGTAACGCTAAATCCTGACATAATTGCACGATTTTA‑GGATTAAATCCCGGTGTGACGACCACATCTGCACCGCTGTTTTTA‑GCGGCGACGACCTGTTCGGCGGTCAATACgg                                                                                                                                                                                                                                 <  1:851324/251‑1 (MQ=255)
                      cTGTGCGTAAGGGCCTAACAAGGCTTTAATCATTTTCACGCCACCAGAGGCTTCCGCCGGGAAGAATTTCACGGTTTCAATACCGACGTCCAACGCGGCTTCAATTGCCATCGGATTGTTTACCCCCGGTGTGATGGGTAACGCTAAATCCTGACATAATTGCACGATTTTA‑GGATTAAATCCCGGTGTGACGACCACATCTGCACCGCTGTTTTTA‑GCGGCGACGACCTGTTCGGCGGTCAATACGG‑TTc                                                                                                                                                                                                                             <  1:797171/251‑1 (MQ=255)
                                   ccTAACAAGGCTTTAATCATTTTCACGCCACCAGAGGCTTCCGCCGGGAAGAATTTCACGGTTTCAATACCGACGTCCAACGCGGCTTCAATTGCCATCGGATTGTTTACCCCCGGTGTGATGGGTAACGCTAAATCCTGACATAATTGCACGATTTTA‑GGATTAAATCCCGGTGTGACGACCACATCTGCACCGCTGTTTTTA‑GCGGCGACGACCTGTTc                                                                                                                                                                                                                                               >  1:300054/1‑221 (MQ=255)
                                                        ttcccGCCACCAGAGGCTCCCGCCGGAAAGAATTCCACGTTTTCAATCCCGACGTCAACCGCGGCTTCAATTGCCATCGGTTTGTTTCCCCCCGGGGTGATGGGTAACGCTAAATCCTGACATAATTGCACGTTTTTA‑GGATTAAACCCCGGTGTGACGCCCACATCTGCCCCGCTGTTTTTA‑GCGGCGACGCCCTGTTCGGCGGTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACg                                                                                                                                                                                           <  1:762579/247‑1 (MQ=255)
                                                                accaGAGGCTTCCGCCGGGAAGAATTTCACGGTTTCAATACCGACGTCCAACGCGGCTTCAATTGCCATCGGATTGTTTACCCCCGGTGTGATGGGTAACGCTAAATCCTGACATAATTGCACGATTTTA‑GGATTAAATCCCGGTGTGACGACCACATCTGCACCGCTGTTTTTA‑GCGGCGACGACCTGTTCGGCGGTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATaaa                                                                                                                                                                                    >  1:37011/1‑250 (MQ=255)
                                                                 ccaGAGGCTTCCGCCGGGAAGAATTTCACGGTTTCAATACCGACGTCCAACGCGGCTTCAATTGCCATCGGATTGTTTACCCCCGGTGTGATGGGTAACGCTAAATCCTGACATAATTGCACGATTTTA‑GGATTAAATCCCGGTGTGACGACCACATCTGCACCGCTGTTTTTA‑GCGGCGACGACCTGTTCGGCGGTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACg                                                                                                                                                                                  >  1:669332/1‑251 (MQ=255)
                                                                 ccaGAGGCTCCCGCCGGGAAGAATTTCACGGTTTCAATCCCGACGTCCAACGCGGCTTCAATTGCCATCGGATTGTTTCCCCCGGGTGTGATGGGTAACGCTAAATCCTGACATAATTGCACGATTTTA‑GGATTAAACCCCGGTGTGACGACCACATCTGCACCGCTGTTTTTA‑GCGGCGACGCCCTGTTCGGCGGTAAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACg                                                                                                                                                                                  <  1:67357/251‑1 (MQ=255)
                                                                       gCTTCCGCCGGGAAGAATTTCACGGTTTCAATACCGACGTCCAACGCGGCTTCAATTGCCATCGGATTGTTTCCCCCCGGTGTGATGGGTAACGCTAAACCCTGACATAATTGCACGATTTTA‑GGATTAAATCCCGGTGTGACGACCACATCTGCCCCGCTGTTTTTA‑GCGGCGACGCCCTGTTCGGCGGTCAATCCGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATgg                                                                                                                                                                             <  1:596714/250‑1 (MQ=255)
                                                                                                tttCAATACCGACGTCCAACGCGGCTTCAATTGCCATCGGATTGTTTACCCCCGGTGTGATGGGTAACGCTAAATCCTGACATAATTGCACGATTTTA‑GGATTAAATCCCGGTGTGACGACCACATCTGCACCGCTGTTTTTA‑GCGGCGACGACCTGTTcgg                                                                                                                                                                                                                                             <  1:420344/162‑1 (MQ=255)
                                                                                                                           cAATTGCCATCGGATTGTTTACCCCCGGTGTGATGGGTAACGCTAAATCCTGACATAATTGCACGATTTTA‑GGATTAAATCCCGGTGTGACGACCACATCTGCACCGCTGTTTTTA‑GCGGCGACGACCTGTTCGGCGGTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCAttt                                                                                                                        >  1:635130/1‑251 (MQ=255)
                                                                                                                                           gTTTCCCCCCGGTGTGATGGGTAACGCTAAATCCTGACATAATTGCACGTTTTTA‑GGATTAAATCCCGGTGTGACGCCCACATCTGCCCCGCTGTTTTTA‑GCGGCGACGCCCTGTTCGGCGGTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACg                                                                                                                                                                                           <  1:123059/168‑1 (MQ=255)
                                                                                                                                           gTTTACCCCCGGTGTGATGGGTAACGCTAAATCCTGACATAATTGCACGATTTTA‑GGATTAAATCCCGGTGTGACGACCACATCTGCACCGCTGTTTTTA‑GCGGCGACGACCTGTTCGGCGGTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGgtgt                                                                                                           >  1:105833/1‑248 (MQ=255)
                                                                                                                                               ccccccGGTGTGGTGGGTAACGCTAAATCCTGACATAATTGCTGGTTTTTA‑GGATAAAATCCGGGTGCGACGCCCACATCTGCACGGTTGTTTTGC‑GCGGGGACGCCCTGTTCGGCGGTAAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCag                                                                                                    <  1:766050/250‑1 (MQ=255)
                                                                                                                                                   cgggTGTGATGGGTACCGATAAATCCTGACATATTTGCACGATTTTA‑GAATTAAATCCGGGTGTGACGACCACATCTGCCCCGCTGTTTTTG‑GCGGCGACGTCCTGTTGGGCGTTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGt                                                                                                <  1:234711/249‑1 (MQ=255)
                                                                                                                                                   ccGGTGTGATGGGTAACGCTAAATCCTGACATAATTGCACGATTTTA‑GGATTAAATCCCGGTGTGACGACCACATCTGCACCGCTGTTTTTA‑GCGGCGACGACCTGTTCGGCGGTCAATACGG‑TTcc                                                                                                                                                                                                                            <  1:135347/127‑1 (MQ=255)
                                                                                                                                                     ggTGTGATGGGTAACGCTAAATCCTGACATAATTGCACGATTTTA‑GGATTAAATCCCGGTGTGACGACCACATCTGCACCGCTGTTTTTA‑GCGGCGACGACCTGTTCGGCGGTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTaa                                                                                              >  1:221183/1‑251 (MQ=255)
                                                                                                                                                      gtgtGATGGGTAACGCTAAATCCTGACATAATTGCACGATTTTA‑GGATTAAATCCCGGTGTGACGACCACATCTGCACCGCTGTTTTTA‑GCGGCGACGACCTGTTCGGCGGTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTaaa                                                                                             >  1:821097/1‑251 (MQ=255)
                                                                                                                                                       tgtgATGGGTAACGCTAAATCCTGACATAATTGCACGTTTTTA‑GGATTAAATCCGGGTGTGACGCCCACATCTGCCCCGTTTTTTTTG‑GCGGCGACGCCCTGTTCGGCGGTCAATTCGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTaaaa                                                                                            <  1:89117/251‑1 (MQ=255)
                                                                                                                                                        gtgATGGGTAACGCTAAATCCTGACATAATTGCACGTTTTTA‑GGATTAAACCCGGGTGTGACGCCCACATCTGCACCGTTTTTTTTA‑GCGGCGACGCCCTGTTCGGCGGTCAATCCGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACg                                                                                                                                                                                  <  1:624749/164‑1 (MQ=255)
                                                                                                                                                        gtgATGGGTAACGCTAAATCCTGACATAATTGCACGATTTTA‑GGATTAAATCCCGGTGTGACGACCACATCTGCACCGCTGTTTTTA‑GCGGCGACGACCTGTTCGGCGGTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAat                                                                                           >  1:492177/1‑251 (MQ=255)
                                                                                                                                                         tgATGGGTAACGCTAAACCCTGACATAATTGCACGTTTTTA‑GGATAAAACCCCGGTGTGACGACCACATCTGCACCGTTGTTTTGA‑GCGGCGACGCCCTGTTCGGCGGTAAATCCGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAata                                                                                          <  1:603235/251‑1 (MQ=255)
                                                                                                                                                           aTGGGTAACGCTAAATCCTGACATAATTGCACGATTTTA‑GGATTAAATCCCGGTGTGACGACCACATCTGCCCCGCTTTTTTTA‑GCGGCGACGACCTGTTCGGCGGTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAatat                                                                                         <  1:51953/250‑1 (MQ=255)
                                                                                                                                                               gaaactcaaaaTCCTGACATAATTGCGCG‑TTTTAGGGATTAAATCCGGGTGGGGCGACCACGGCTGCCCCGGT‑TTTTTAGGGGGGGGCGCCGTGTTCGGCGGTAAATGGGG‑TTGCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAATATCTTc                                                                                     <  1:323961/242‑1 (MQ=255)
                                                                                                                                                                   cGCTAAATCCTGACATAATTGCACGATTTTA‑GGATTAAATCCCGGTGTGACGACCACATCTGCACCGCTGTTTTTA‑GCGGCGACGACCTGTTCGGCGGTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTgc                                                                                                                                                                      >  1:96945/1‑165 (MQ=255)
                                                                                                                                                                             tgaaaTAATTGCCCGTTTTTA‑GAATAAAACCCGGGTGTGACGCCCACATCTCCACCGCTGTTTTTA‑GCGGCGACGCCCTGTCCGGCGGTAAAT‑CGGTTTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAATATCTTCCGCGTTGTCTAAAGc                                                                      <  1:3422/247‑1 (MQ=255)
                                                                                                                                                                             tgaaaTAATTGCACGTTTTTA‑GGTTTAAACCCGGGTGTGACGCCCACATCTCCCCCGCTGTTTTTA‑GCGGCGACGCCCTGTTCGGCGGTAAATCCGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAATATCTTCCGCGTTGTCTAAAGc                                                                      <  1:871357/247‑1 (MQ=255)
                                                                                                                                                                             taaCATAATTGCCCGTTTTTA‑GGATTAACCCCGGGTGGGACGACCACTTCTGCCCCGTTGTTTTTA‑GCGGCGACGCCCTGTTCGGCGGTCAATACGT‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAATATCTTCCGCGTTGTCTAAAGc                                                                      <  1:370272/249‑1 (MQ=255)
                                                                                                                                                                             tGACATAATTGCACGTTTTTA‑GGATTAAATCCCGGTGTGACGCCCACATCTGCACCGCTTTTTTGA‑GCGGCGACGCCCTGTTCGGCGGTCAATACGG‑TCCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAATATCTTCCGCGTTGTCTAAAGc                                                                      <  1:521170/251‑1 (MQ=255)
                                                                                                                                                                             tGACATAATTGCACGTTTTTA‑GGATAAACCCCCGGTGTGACGACCCCATCTGCACCGTTTTTTTTA‑GTGGCGACGCCCTGTTCGGCGGTAAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAATATCTTCCGCGTTGTCTAAAGc                                                                      <  1:692794/251‑1 (MQ=255)
                                                                                                                                                                             tGACATAATTGCACGATTTTA‑GGATAAAATCCCGGTGTGACGACCACATCTGCCCCGCTGTTTTGA‑GCGGCGACGCCCTGTGCGGCGGTCAATCCGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAATATCTTCCGCGTTGTCTAAAGc                                                                      <  1:448568/251‑1 (MQ=255)
                                                                                                                                                                               aCATAATGGCCCGTTTTTA‑GGATAAAACCCGGGTGTGACGCCCACATCTGCCCCGTTTTTTTTA‑GCGGCGACGCCCTGTTGGGCGGTCAATGCGT‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAATATCTTCCGCGTTGTCTAAAGcc                                                                     <  1:813227/250‑2 (MQ=255)
                                                                                                                                                                                    tttccACGTTTTTA‑GAATTAACTCCCGGTGTGACGCCCACATCTGCCCCGTTTTTTTTG‑GCGGCGACGCCCTGTTCGGCGGTAAATACGG‑TCCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAATATCTTCCGCGTTGTCTAAAGCAATCAcc                                                               <  1:38822/247‑1 (MQ=255)
                                                                                                                                                                                          cGATTTTA‑GGATTAAATCCCGGTGTGACGACCACATCTGCACCGCTGTTTTTA‑GCGGCGACGACCTGTTCGGCGGTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGCGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAatat                                                                                         >  1:887079/1‑219 (MQ=255)
                                                                                                                                                                                            aTTTTA‑GGATTAAATCCCGGTGTGACGACCACATCTGCACCGCTGTTTTTA‑GCGGCGACGACCTGTTCGGCGGTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGcggcgggaaccgcatttgac                                                                                                                                                   >  1:626921/1‑145 (MQ=255)
                                                                                                                                                                                            aTTTTA‑GGATTAAATCCCGGTGTGACGACCACATCTGCACCGCTGTTTTTA‑GCGGCGACGACCTGTTCGGCGGTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGcc                                                                                                      >  1:702180/1‑204 (MQ=255)
                                                                                                                                                                                            aTTTTA‑GGATTAAATCCCGGTGTGACGACCACATCTGCACCGCTGTTTTTA‑GCGGCGACGACCTGTTCGGCGGTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAATATCTTc                                                                                     >  1:393981/1‑221 (MQ=255)
                                                                                                                                                                                            aTTTTA‑GGATTAAATCCCGGTGTGACGACCACATCTGCACCGCTGTTTTTA‑GCGGCGACGACCTGTTCGGCGGTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAATATCTTCCGCGTTGTCTAAAGCAATCACCGGCACGa                                                        >  1:218791/1‑250 (MQ=255)
                                                                                                                                                                                            aTTTTA‑GGATTAAACCCCGGTGTGACGACCACATCTGCCCCGCTGTTTTTA‑GCGGCGACGACCTGTTCGGCGGTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAg                                                                                                                                                         <  1:822718/153‑1 (MQ=255)
                                                                                                                                                                                                     ataaaacccgggTGTGACGCCCACATCTCCCCCGCTTTTTTTA‑GCGGCGACGCCTTGTTCGGCGGTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGAAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAATATCTTCCGCGTTGTCTAAAGCAATCACCGGCACGACTTTTAAt                                                <  1:120869/240‑1 (MQ=255)
                                                                                                                                                                                                     aTTAAACCCGGGTGTGACGCCCACATCTGCCCCGCTTTTTTGA‑GCGGCGACGACCTGTTCGGCGGTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAATATCTTCCGCGTTGTCTAAAGCAATCACCGGCACGACTTTTAAt                                                <  1:356018/250‑1 (MQ=255)
                                                                                                                                                                                                          accccGGTGTGACGACCACATCTGCCCCGCTGTTTTTA‑GCGGCGACGCCCTGTTCGGCGGTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGtaat                                                                                                                                                                                       <  1:437763/108‑1 (MQ=255)
                                                                                                                                                                                                                gtgtGACGACCACATCTGCACCGCTGTTTTTA‑GCGGCGACGACCTGTTCGGCGGTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAATATCTTCCGCGTTGTCTAAAGCAATCACCGGCACGACTTTTAATTGACGAAAtttt                                    >  1:228753/1‑251 (MQ=255)
                                                                                                                                                                                                                   tgACGCCCACATCTGCCCCGTTTTTTTTA‑GGGGCGACGCCCTGTTCGGCGGTAAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAATATCTTCCGCGTTGTCTAAAGCAATCACCGGCACGACTTTTAATTGACGAAATTTTTc                                  <  1:633058/250‑1 (MQ=255)
                                                                                                                                                                                                                       gacCACATCTGCACCGCTGTTTTTA‑GCGGCGACGACCTGTTCGGCGGTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAATATCTTCCGCGTTGTCTAAAGCAATCACCGGCACGACTTTTAATTGACGAAATTTTTCAATGa                             >  1:699139/1‑251 (MQ=255)
                                                                                                                                                                                                                            caTCTGCACCGCTGTTTTTA‑GCGGCGACGCCCTGTTCGGCGGTAAATGCGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAATATCTTCCGCGTTGTCTAAAGCAATCACCGGCACGACTTTTAAt                                                <  1:660075/227‑1 (MQ=255)
                                                                                                                                                                                                                            caTCTGCACCGCTGTTTTTA‑GCGGCGACGACCTGTTCGGCGGTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAATATCTTCCGCGTTGTCTAAAGc                                                                      >  1:41604/1‑205 (MQ=255)
                                                                                                                                                                                                                            caTCTGCACCGCTGTTTTTA‑GCGGCGACGACCTGTTCGGCGGTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAATATCTTCCGCGTTGTCTAAAGc                                                                      >  1:246339/1‑205 (MQ=255)
                                                                                                                                                                                                                            caTCTGCACCGCTGTTTTTA‑GCGGCGACGACCTGTTCGGCGGTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAATATCTTCCGCGTTGTCTAAAGCAATCACCGGCACGACTTTTAATTGACGAAATTTTTCAATGAGTTgg                        >  1:94823/1‑251 (MQ=255)
                                                                                                                                                                                                                            caTCTGCACCGCTGTTTTTA‑GCGGCGACGACCTGTTCGGCGGTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGGGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAATATCTTCCGCGTTGTCTAAAGCAATCACCGGCACGACTTTTAATTGACGAAATTTTTCAATGAGTTgg                        >  1:536669/1‑251 (MQ=255)
                                                                                                                                                                                                                            caTCTGCACCGCTGTTTTTA‑GCGGCGACGACCCGTTCGGCGGTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGTCAAGGTGTCTGCCAGAGGTAAAATATcgtc                                                                                     >  1:377938/1‑187 (MQ=255)
                                                                                                                                                                                                                              tctgccccgttttttttA‑GCGGCGACGCCCTGTTCGGCGGTAAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAATATCTTCCGCGTTGTCTAAAGCAATCACCGGCACGACTTTTAATTGACGAAATTTTTCAATGAGTTGGtt                      <  1:537892/239‑1 (MQ=255)
                                                                                                                                                                                                                              tCTGCCCCGCTTTTTTTA‑GCGGCGACGCCCTGTTCGGCGGTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAATATCTTCCGCGTTGTCTAAAGCAATCACCGGCACGACTTTTAATTGACGAAATTTTTCAATGAGTTGGtt                      <  1:100644/251‑1 (MQ=255)
                                                                                                                                                                                                                               cTGCACCGCTGTTTTTA‑GCGGCGACGACCTGTTCGGCGGTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGCGTCTGCCAGAGGTAAAATATCTTCCGCGTTGTCTAAAGCAATCACCGGCACGACTTTTAATTGACGAAATTTTTCAATGAGTTGGTTa                     >  1:601679/1‑251 (MQ=255)
                                                                                                                                                                                                                               cTGCACCGCTGTTTTTA‑GCGGCGACGACCTGTTCGGCGGTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGCGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAATATCTTCCGCGTTGTCTAAAGCAATCACCGGCACGACTTTTAATTGACGAAATTTTTCAATGAGTTGGTTa                     >  1:503181/1‑251 (MQ=255)
                                                                                                                                                                                                                                tgcgccgttttttttA‑GCGGCGACGCCCTGTTCGGCGGTAAATAGGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAATATCTTCCGCGTTGTCTAAAGCAATCACCGGCACGACTTTTAATTGACGAAATTTTTCAATGAGTTGGTTAg                    <  1:537589/241‑1 (MQ=255)
                                                                                                                                                                                                                                tccacgtcTGTTTTTA‑GGGGCGACGCCCTTTTGGGCGTTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAATATCTTCCGCGTTGTCTAAAGCAATCACCGGCACGACTTTTAATTGACGAAATTTTTCAATGAGTTGGTTAg                    <  1:613349/244‑1 (MQ=255)
                                                                                                                                                                                                                                   cccGCTGTTTTTA‑GCGGCGACGACCTGTTCGGGGGTAAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAATATCTTCCGCGTTGTCTAAAGCAATCACCGGCACGACTTTTAATTGACGAAATTTTTCAATGAGTTGGTTAGTtg                 <  1:846764/250‑1 (MQ=255)
                                                                                                                                                                                                                                      gCTGTTTTTA‑GCGGCGACGACCTGTTCGGCGGTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAATATCTTCCGCGTTGTCTAAAGCAATCACCGGCACGACTTTTAATTGACGAAATTTTTCAAtcatttctttcggcgcgt              >  1:220499/1‑234 (MQ=255)
                                                                                                                                                                                                                                        tttttttA‑GGGGCGACGCCCTGTTGGGCGGTAAATACGT‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAATATCTTCCGCGTTGTCTAAAGCAATCACCGGCACGACTTTTAATTGACGAAATTTTTCAATGAGTTGGTTAGTTGTGTAt            <  1:293184/249‑1 (MQ=255)
                                                                                                                                                                                                                                        tttttttA‑GCGGCGACGCCGTTTGCGGGGGTCAATCCGT‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAATATCTTCCGCGTTGTCTAAAGCAATCACCGGCACGACTTTTAATTGACGAAATTTTTCAATGAGTTGGTTAGTTGTGTAt            <  1:150432/249‑1 (MQ=255)
                                                                                                                                                                                                                                        tttttttA‑GCGGCGACGACCTGTTCGGCGGTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAATATCTTCCGCGTTGTCTAAAGCAATCACCGGCACGACTTTTAATTGACGAAATTTTTCAATGAGTTGGTTAGTTGTGTAt            <  1:274998/249‑1 (MQ=255)
                                                                                                                                                                                                                                              ttggcggcgACGCCCTGTTCGGCGGTAAATAGGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAATATCTTCCGCGTTGTCTAAAGCAATCACCGGCACGACTTTTAATTGACGAAATTTTTCAATGAGTTGGTTAGTTGTGTATGACAt       <  1:687798/248‑1 (MQ=255)
                                                                                                                                                                                                                                                    gcgACGCCTTGTTCGGCGGTCAATCCGT‑TCCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAATATCTTCCGCGTTGTCTAAAGc                                                                      <  1:69075/182‑1 (MQ=255)
                                                                                                                                                                                                                                                         gacCTGTGGGGCGGTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAATATCTTCCGCGTTGTCTAAAGCAATCACCGGCACGACTTTTAATTGACGAAATTTTTCAATGAGTTGGTTAGTTGTGTATGACATAAATc  <  1:40073/245‑1 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    
CGCCGGTCGGCATAATTTGCAACTGTGCGTAAGGGCCTAACAAGGCTTTAATCATTTTCACGCCACCAGAGGCTTCCGCCGGGAAGAATTTCACGGTTTCAATACCGACGTCCAACGCGGCTTCAATTGCCATCGGATTGTTTACCCCCGGTGTGATGGGTAACGCTAAATCCTGACATAATTGCACGATTTTA‑GGATTAAATCCCGGTGTGACGACCACATCTGCACCGCTGTTTTTA‑GCGGCGACGACCTGTTCGGCGGTCAATACGG‑TTCCCGCCGCAATCAGAAATTCCGGGCGCTGTTGACGTAATAAACGAATGGCTTCTGCGGCGGCAGGGGAGCGAAAGGTGATTTCCGCCACCGGCAAACCATTTTGCGCCAAGGTGTCTGCCAGAGGTAAAATATCTTCCGCGTTGTCTAAAGCAATCACCGGCACGACTTTTAATTGACGAAATTTTTCAATGAGTTGGTTAGTTGTGTATGACATAAATC  >  1/187832‑188323

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: ‑