Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1154,2040GT41.4% 116.8 / 21.1 60C609F (TGC→TTC) unknownRibonucleotide reductase of class III (anaerobic), large subunit (EC 1.17.4.2)
Reads supporting (aligned to +/- strand):  ref base G (3/31);  new base T (24/0);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 2.28e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCACCAAACAATGGCGCGCCGAAACGGGTTATGCCTTCAGTTTGTATTCCACTCCAAGCGAAAACCTGTGCGATCGCTTCTGCCGTTTAGATAC‑CAAACAATTCGGCGTCATCGACGGCGTGACCGACAAAGGCTACTACACCAAC‑AGCTACCATTT‑GGACGTGGAGAAAAAAGTCAAT‑CCATACG‑ACAAACTGGACTTTGAAATGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTCATTTGCTACGGC‑GAATACCCAAACATCCAGCACAACCTCAAAGCACTGGAAGACGTGTGGAATTACAGCTACGACCGCGTGCCTTACTACGGCACCAACACGCCTATCGACGAATGCTACGAATGC‑GGTTTCACCGGC‑G‑AATTTGAATGCACCAGCAAAGGCTTC‑ACCTG‑CCCGAAATGCGGCAACCACGACAGC‑GAAAAAGTCTCCGTCACCCGACGGGTTTGCGGCTACCT  >  1/153964‑154432
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                           
ccaccaAACAATGGCGCGCCGAAACGGGTTATGCCTTCAGTTTGTATTCCACTCCAAGCGAAAACCTGTGCGATCGCTTCTGCCGTTTAGATAC‑CAAACAATTCGGCGTCATCGACGGCGTGACCGACAAAGGCTACTACACCAAC‑AGCTACCATTT‑GGACGTGGAGAAAAAAGTCAAT‑CCATACG‑ACAAACTGGACTTTGAAATGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTCATTTTCTACGGC‑Gaa                                                                                                                                                                                                                                  >  1:25331/1‑251 (MQ=255)
     aaaCAATGGCGCGCCGAAACGGGTTATGCCTTCAGTTTGTATTCCACTCCAAGCGAAAACCTGTGCGATCGCTTCTGCCGTTTAGATAC‑CAAACAATTCGGCGTCATCGACGGCGTGACCGACAAAGGCTACTACACCAAC‑AGCTACCATTT‑GGACGTGGAGAAAAAAGTCAAT‑CCATACG‑ACAAACTGGACTTTGAAATGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTCATTTTCTACGGC‑GGATAcc                                                                                                                                                                                                                              >  1:550916/1‑250 (MQ=255)
     aaaCAATGGCGCGCCGAAACGGGTTATGCCTTCAGTTTGTATTCCACTCCAAGCGAAAACCTGTGCGATCGCTTCTGCCGTTTAGATAC‑CAAACAATTCGGCGTCATCGACGGCGTGACCGACAAAGGCTACTACACCAAC‑AGCTACCATTT‑GGACGTGGAGAAAAAAGTCAAT‑CCATACG‑ACAAACTGGACTTTGAAATGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTCATTTTCTACGGC‑GAATAcc                                                                                                                                                                                                                              >  1:539831/1‑250 (MQ=255)
                  ccaaaccggtttATGCCTTCAGTTTTTATCCCCCTCCAAGCAAAACCCTGTGCGATCGCTCCGGCCGTTTAAATAC‑CAAACATTTCGGCGTCATCGACGGCGTGACCGACAAAGGCTACTCCACCAAC‑ACCTACCATTT‑GGACGTGAAAAAAAAAGTCAAT‑CCATACG‑ACAAACTGGACTTTAAAAGGG‑CTTACCCGCCACTTGCCAGCGGCGGTTTCATTTGCTACGGC‑GAATACCCAAACATCCAGcac                                                                                                                                                                                                                <  1:457566/241‑1 (MQ=255)
                         gggTTATGCCTTCAGTTTGTATTCCACTCCAAGCGAAAACCTGTGCGATCGCTTCTGCCGTTTAGATAC‑CAAACAATTCGGCGTCATCGACGGCGTGACCGACAAAGGCTACTACACCAAC‑AGCTACCATTT‑GGACGTGGAGAAAAAAGTCAAT‑CCATACG‑ACAAACTGGACTTTGAAATGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTTATTTTCTACCGC‑GGATACCCCAACATcccccccaacctc                                                                                                                                                                                                          >  1:195596/1‑239 (MQ=255)
                          ggTTATGCCTTCAGTTTGTATTCCACTCCAAGCGAAAACCTGTGCGATCGCTTCTGCCGTTTAGATAC‑CAAACAATTCGGCGTCATCGACGGCGTGACCGACAAAGGCTACTACACCAAC‑AGCTACCATTT‑GGACGTGGAGAAAAAAGTCAAT‑CCATACG‑ACAAACTGGACTTTGAAATGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTCATTTGCTACGGC‑GAAAAACCCAACAACCAGCACaaacccaa                                                                                                                                                                                                        >  1:271371/1‑245 (MQ=255)
                             tATGCCTTCATTTTGTTTTCCCCCCCAAGCAAAACCCTGTGCGATCGCTTCTGCCGTTTAGATAC‑CAAACAATTCGGCGTCATCGACGGCGTGACCGACAAAGGCTACTACACCAAC‑AGCTACCATTT‑GGACGTGGAAAAAAAAGTCAAT‑CCATACG‑ACAAACTGGATTTTAAAATGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTCATTTGCTACGGC‑GAATACCCAAACATCCAGCACAACCTCAAAGc                                                                                                                                                                                                     <  1:174006/251‑1 (MQ=255)
                              ttGCCTTCAGTTTGTATCCCCCCCCAAGCAAAACCCGGTGCGATCGCTTCTGCCTTTTAGATAC‑AAAACAATTCGGCGTCATCGCCGGCGTGACCGAAAAAGGCTACTACCCCAAC‑AGCTCCCATTT‑GGACGGGGAAAAAAAAGCCAAC‑CCATACG‑AAAACCGGGATTTTGAAATGG‑CGTCCCCGCCACTTGCCAGCGGCGGTTTCATTTGCTACGGC‑GAATACCCAAACATCCAGCACAACCTCAAAGCa                                                                                                                                                                                                    <  1:120745/250‑1 (MQ=255)
                              aTGCCTTCAGTTTGTATTCCACTCCAAGCGAAAACCTGTGCGATCGCTTCTGCCGTTTAGATAC‑CAAACAATTCGGCGTCATCGACGGCGTGACCGACAAAGGCTACTACACCAAC‑AGCTACCATTT‑GGACGTGGAGAAAAAAGTCAAT‑CCATACG‑ACAAACTGGACTTTGAAATGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTCATTTTGTACGGC‑GGAAACACAAAAATCCAACACAACCTCAAAgaa                                                                                                                                                                                                    >  1:379317/1‑249 (MQ=255)
                               tGCCTTCAGTTTGTATTCCACCCAAAGCAAAACCCTGTGCGATCGCTTCTGCCTTTTAGATAC‑AAAACAATCGGGCGTCATCGAGGGCGTGACCGAAAAGGGCTACTACACCAAC‑AGCTACCATTT‑GGACGGGGAGAAAAAAGTAAAC‑CCATACG‑AAAAACGGGACTTTGAAATGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTCATTTGCTACGGC‑GAATACCCAAACATCCAGCACAACCTCAAAGCAc                                                                                                                                                                                                   <  1:572822/251‑1 (MQ=255)
                                gCCTTCAGTTTGTATTCCACTCCAAGCGAAAACCTGTGCGATCGCTTCTGCCGTTTAGATAC‑CAAACAATTCGGCGTCATCGACGGCGTGACCGACAAAGGCTACTACACCAAC‑AGCTACCATTT‑GGACGTGGAGAAAAAAGTCAAT‑CCATACG‑ACAAACTGGACTTTGAAATGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTTCTTTTCTTCGGC‑GGATACCCAAACATCCAGCCCAAcccccaagcca                                                                                                                                                                                                   >  1:52413/1‑241 (MQ=255)
                                gCCTTCAGTTTGTATTCCACTCCAAGCGAAAACCTGTGCGATCGCTTCTGCCGTTTAGATAC‑CAAACAATTCGGCGTCATCGACGGCGTGACCGACAAAGGCTACTACACCAAC‑AGCTACCATTT‑GGACGTGGAGAAAAAAGTCAAT‑CCATACG‑ACAAACTGGACTTTGAAATGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTCATTTTCTTCGGCGGAAAACCC‑AACAACCCGCACAAACCCAAAGCAc                                                                                                                                                                                                   >  1:319187/1‑250 (MQ=255)
                                  cTTCAGTTTGTATCCCACCCAAAGCAAAACCCTGTGCGATCGCTTCTGCCTTTAAGATAC‑AAAACAATTCGGCGTCATCGCCGGCGGGACCGAAAAAGGCTACTACACCAAC‑AGCTACCATTG‑GGACGTGGAAAAAAAAGTCAAT‑CCAAACG‑ACAACCGGGATTTTGAAATGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTCATTTGCTACGGC‑GAATACCCAAACATCCAGCACAACCTCAAAGc                                                                                                                                                                                                     <  1:290041/246‑1 (MQ=255)
                                    tCAGTTTGTTTCCCACTCCAAGCGAAACCCTGTGCGATCTCTTCTGCCGTTTAGATAC‑CAAACAATTCGGCGTCATCGCCGGGGTGACCAAAAAAGGCTACTACACCAAC‑AGCTCCCATTG‑GGACGGGGAGAAAAAAGCCAAT‑CCATACA‑ACAACCTGGACTTTGAAATGG‑CGTCCCCGCCACTTGCCAGCGGCGGTTTCATTTGCTACGGC‑GAATACCCAAACATCCAGCACAACCTCAAAGCACTGGaa                                                                                                                                                                                              <  1:360047/251‑1 (MQ=255)
                                       ttttttATTCCCCTCCAAGCAAAAACCTGTGCGATCGCTTCTGCCGTTTAGATAC‑AAAACAATCCGGCGTCATCGACGGCGTGACCGACAAAGGCTCCTACACAAAC‑AGCTCCAATTT‑GGACGGGGAAAAAAAAGCCAAC‑CCATACG‑ACAAACTGGATTTTAAAAGGG‑CGTCCCCGCCACTTGCCAGCGGCGGTTTCATTTGCTACGGC‑GAATACCCAAACATCCAGCACAACCTCAAAGCACTGGAAGAc                                                                                                                                                                                           <  1:78743/246‑1 (MQ=255)
                                         ttGTATTCCACTCCAAGCGAAAACCTGTGCGATCGCTTCTGCCGTTTAGATAC‑CAAACAATTCGGCGTCATCGACGGCGTGACCGACAAAGGCTACTACACCAAC‑AGCTACCATTT‑GGACGTGGAGAAAAAAGTCAAT‑CCATACG‑ACAAACTGGACTTTGAAATGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTCATTTTCTACCGC‑GGATACCCCAACATCCAACACAACCCCAAAGCACCGGAAGACgg                                                                                                                                                                                         >  1:251004/1‑250 (MQ=255)
                                                     caaaggaaaaccctgtgcgtccccTTCGGCCGTTTAGATAC‑AAACCAATTCGGGGCCACCGGCGGGGGCCCCGAAAAAGGCTACTACCCCAAC‑AGCACCCTTTT‑GGAGGGGGAAAAAAAATAAAAC‑CCAAACG‑AAAACCTGGA‑TTTAAAAAGGCCTTACCCGCCACTTGCCAGCGGCGGTTTCATTTGCTACGGC‑GAATACCCAAACATCCAGCACAACCTCAAAGCACTGGAAGACGTGTGGAATTACAg                                                                                                                                                                             <  1:641494/228‑1 (MQ=255)
                                                        aGCGAAAACCTGTGCGATCGCTTCTGCCGTTTAGATAC‑CAAACAATTCGGCGTCATCGACGGCGTGACCGACAAAGGCTACTACACCAAC‑AGCTACCATTT‑GGACGTGGAGAAAAAAGTCAAT‑CCATACG‑ACAAACTGGACTTTGAAATGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTCATTTTCTACGGC‑GGAAACCCAAAAATCCAACACAACCTCCAAAAACAGGCAGACGGGTGGAATTACaacca                                                                                                                                                                          >  1:485855/1‑247 (MQ=255)
                                                         gCGAAAACCTGTGCGATCGCTTCTGCCGTTTAGATAC‑CAAACAATTCGGCGTCATCGACGGCGTGACCGACAAAGGCTACTACACCAAC‑AGCTACCATTT‑GGACGTGGAGAAAAAAGTCAAT‑CCATACG‑ACAAACTGGACTTTGAAATGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTTATTTTCTACGGC‑GGAAACCCCAAAAACCAGCACAACCTCAAAACACTGGaaaaagggggggaataaaacca                                                                                                                                                                          >  1:219842/1‑230 (MQ=255)
                                                         gCGAAAACCTGTGCGATCGCTTCTGCCGTTTAGATAC‑CAAACAATTCGGCGTCATCGACGGCGTGACCGACAAAGGCTACTACACCAAC‑AGCTACCATTT‑GGACGTGGAGAAAAAAGTCAAT‑CCATACG‑ACAAACTGGACTTTGAAATGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTCATTTTCTACGGC‑GAATACCCCAACAACCCGCACAAACTCAAAACACTGGAAGACGGGTgggaataaagcta                                                                                                                                                                          >  1:199476/1‑239 (MQ=255)
                                                            aaacccTGGGCGATCGCTTCTGCCGTTTA‑AAACACAAAAAATTCGGGGTCACCGACGGCGTGCCCCAAAAAGCCTCCCACACCAAC‑AGCTACCATTG‑GGACGTGGAAAAAAAAGTAAAC‑CCAAACA‑AAAAACGGGACTTTGAAATGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTCATTTGCTACGGC‑GAATACCCAAACATCCAGCACAACCTCAAAGCACTGGAAGACGTGTGGAATTACAGCTACGAc                                                                                                                                                                      <  1:562011/247‑1 (MQ=255)
                                                              aaCCTGTGCGATCGCTTCTGCCGTTTAGATAC‑CAAACAATTCGGCGTCATCGACGGCGTGACCGACAAAGGCTACTACACCAAC‑AGCTACCATTT‑GGACGTGGAGAAAAAAGTCAAT‑CCATACG‑ACAAACTGGACTTTGAAATGG‑CGTACCCGCCACTTGCCAGCGGCGGGTTCATTTTCTACGGC‑GAATAACCAAAAATCCCGCACAACCCCAAAGAACTGGAAGACGTGGGGGATTACAGCTACGAcc                                                                                                                                                                     >  1:391923/1‑250 (MQ=255)
                                                                ccTGTGCGATCGCTTCTGCCGTTTAGATAC‑CAAACAATTCGGCGTCATCGACGGCGTGACCGACAAAGGCTACTACACCAAC‑AGCTACCATTT‑GGACGTGGAGAAAAAAGTCAAT‑CCATACG‑ACAAACTGGACTTTGAAATGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTCATTTTCTACGGC‑GGAAACCCAAACAACCACCACAACCTCCAAACACTGGAAGACGGGTGGAAATACAGCTACGACcgc                                                                                                                                                                   >  1:370636/1‑250 (MQ=255)
                                                                           tcTTCTGCCGTTTAGATAC‑CAAACAAGTCGGGGTCACCGACGGCGTAACCAACAAAGGCTACTCCACCAAC‑ACCACCCATTG‑GGACGGGGAGAAAAAAGTCAAT‑CCAAACG‑ACAAGCTGGATTTTAAAATGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTCATTTGCTACGGC‑GAATACCCAAACATCCAGCACAACCTCAAAGCACTGGAAGACGTGTGGAATTACAGCTACGACCGCGTGCCTtactac                                                                                                                                                       <  1:37848/250‑1 (MQ=255)
                                                                           ccTTCTGCCTTTTAAATAC‑CAAACAATGCGGCGTCATCGGGGGCGTGACCAAAAAAGCCTCCTACCCCAAC‑ACCTCCCATTT‑GGAGGTGGAAAAAAAAGTAAAC‑CCATACG‑ACAAACGGGATTTTGAAATGG‑CGTCCCCGCCACTTGCCAGCGGCGGTTTCATTTGCTACGGC‑GAATACCCAAACATCCAGCACAACCTCAAAGCACTGGAAGACGTGTGGAATTACAGCTACGACCGCGTGCCTtactac                                                                                                                                                       <  1:410641/250‑1 (MQ=255)
                                                                           ccTTCTGCCTTTAAAATAC‑AAAAAAATCCGGCGTCATCGACGGCGTGACCGAAAAAGCCTACACCACCAAC‑AACTACAATTG‑GGACGTGGAAAAAAAAGTCAAC‑CCATCCG‑ACAACCTGGACTTTAAAAGGG‑CGTCCCCGCCACTTGCCAGCGGCGGTTTCATTTGCTACGGC‑GAATACCCAAACATCCAGCACAACCTCAAAGCACTGGAAGACGTGTGGAATTACAGCTACGACCGCGTGCCTtactac                                                                                                                                                       <  1:533164/250‑1 (MQ=255)
                                                                           ccTTCTGCCGTTTAGATCC‑AAAACAATTCGGCGTCCTCGACGGCGTGCCCAACAAAGGCTACTACCCCAAC‑AGCTCCCATTG‑GGAGGGGGAGAAAAAAGTCAAT‑CCATACG‑ACAACCTGGACTTTGAAATGG‑CGTCCCCGCCACTTGCCAGCGGCGGTTTCATTTGCTACGGC‑GAATACCCAAACATCCAGCACAACCTCAAAGCACTGGAAGACGTGTGGAATTACAGCTACGACCGCGTGCCTtactac                                                                                                                                                       <  1:241764/250‑1 (MQ=255)
                                                                           ccTTCTGCCGTTTAGATAC‑CAAACAATTCGGCGTCATCGACGGCGTGACCGAAAAAGGCTACTCCACCAAC‑AGCTCCCATTT‑GGACGTGGAAAAAAAAGTCAAT‑CCATACG‑ACAAACTGGACTTTGAAATGG‑CGTCCCCGCCACTTGCCAGCGGCGGTTTCATTTGCTACGGC‑GAATACCCAAACATCCAGCACAACCTCAAAGCACTGGAAGACGTGTGGAATTACAGCTACGACCGCGTGCCTtactac                                                                                                                                                       <  1:401313/250‑1 (MQ=255)
                                                                            cTTCTGCCGTTTAGATAC‑CAAACAATTCGGCGTCATCGACGGCGTGACCGACAAAGGCTACTACACCAAC‑AGCTACCATTT‑GGACGTGGAGAAAAAAGTCAAT‑CCATACG‑ACAAACTGGACTTTGAAATGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTTATTTATCACGGG‑CAATACCCCAAAATCCCGCACAACCCCACAAAACTGGAAGAAGTGAGGGAATACACCTACCACCGCGTGCCCtactac                                                                                                                                                       >  1:613080/1‑250 (MQ=255)
                                                                            cTTCTGCCGTTTAAATCC‑AAAACAATTCGGCGTCATCGGCGGCGTGCCCGACAAAGGCTACTCCA‑CAACAAGCT‑CCATTTGGGACGGGGAAAAAAAATTCAAT‑CCATACG‑ACAAACTGGATTTTGAAAGGG‑CGTCCCCGCCACTTGCCAGCGGCGGTTTCATTTGCTACGGC‑GAATACCCAAACATCCAGCACAACCTCAAAGCACTGGAAGACGTGTGGAATTACAGCTACGACCGCGTGCCTtactac                                                                                                                                                       <  1:579608/250‑1 (MQ=255)
                                                                                    gTTTAGATAC‑CAAACAATTCGGCGTCATCGACGGCGTGACCGACAAAGGCTACTACACCAAC‑AGCTACCATTT‑GGACGTGGAGAAAAAAGTCAAT‑CCATACG‑ACAAACTGGACTTTGAAATGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTCATTTTCTACGGC‑GAATACCCAAAAATCCAGCACAACCTCAAAACACTGGAAGACGTGTGGAATTACAACTACGACCGCGGGCCTTACTACGGCACCaa                                                                                                                                               >  1:133485/1‑250 (MQ=255)
                                                                                     tttAGATAC‑CAAACAATTCGGCGTCATCGACGGCGTGACCGACAAAGGCTACTACACCAAC‑AGCTACCATTT‑GGACGTGGAGAAAAAAGTCAAT‑CCATACG‑ACAAACTGGACTTTGAAATGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTCATTTTCTACGGC‑GAATACCCAAAAATCCAGCACAACCTCAAAGCACTGGAAGACGGGTGGAATTACAGCTACGACCGCGCGCCCTACTACGGCACCAaca                                                                                                                                             >  1:29606/1‑251 (MQ=255)
                                                                                             acaaaaCAATTCGGCGTCACCGACGGCGTGACCGAAAAAGGCTCCTACACCAAC‑AGCTCCCTTTG‑GGACGGGGAAAAAAAAGTCAAT‑CCATACG‑ACAAACTGGCCTTTGAAAGGG‑CGTCCCCGCCACTTGCCAGCGGCGGTTTCATTTGCTACGGC‑GAATACCCAAACATCCAGCACAACCTCAAAGCACTGGAAGACGTGTGGAATTACAGCTACGACCGCGTGCCTTACTACGGCACCAACACGCCTAt                                                                                                                                      <  1:104945/248‑1 (MQ=255)
                                                                                               cAAACAATTCGGCGTCATCGACGGCGTGACCGAAAAAGGCTACTACACCAAC‑AGCTACCATTT‑GGCCGTGGAGAAAAAAGTCAAC‑CCATACG‑AAAAACTGGATTTTGAAAGGG‑CGTCCCCGCCACTTGCCAGCGGCGGTTTCATTTGCTACGGC‑GAATACCCAAACATCCAGCACAACCTCAAAGCACTGGAAGACGTGTGGAATTACAGCTACGACCGCGTGCCTTACTACGGCACCAACACGCCTATcg                                                                                                                                    <  1:184704/251‑1 (MQ=255)
                                                                                                aaaCAATTCGGCGTCATCGAGGGCGTGACCGAAAAAGGCTACTCCACCAAC‑AGCTACCATTT‑GGACGTGGAAAAAAAAGTCATT‑CCATCCG‑ACAAACTGGATTTTAAAAGGG‑CGTCCCCGCCACTTGCCAGCGGCGGTTTCATTTGCTACGGC‑GAATACCCAAACATCCAGCACAACCTc                                                                                                                                                                                                          <  1:23330/180‑1 (MQ=255)
                                                                                                 aaCAATTCGGCGTCATCGAGGGCGTGACCGACAAAGGCTACTACACCAAC‑AGCTCCCATTG‑GGAGGTGGAGAAAAAAGTCAAT‑CCATACG‑AAAAACTGGACTTTGAAAGGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTCATTTGCTACGGC‑GAATACCCAAACATCCAGCACAACCTCAAAGCACTGGAAGACGTGTGGAATTACAGCTACGACCGCGTGCCTTACTACGGCACCAACACGCCTATcgac                                                                                                                                  <  1:585932/251‑1 (MQ=255)
                                                                                                                    cGGCGTGCCCGAAAAAGGCTACTACACCAAC‑AGCTACCATGT‑GGCGGGGGAAAAAAAACT‑AATGCCA‑GCGCAAAAACGGGTTTTAGAAGTGG‑TGTACCCGCCACTTGCCAGCGGCGGTTTCATTTGCTACGGC‑GAATACCCAAACATCCAGCACAACCTCAAAGCACTGGAAGACGTGTGGAATTACAGCTACGACCGCGTGCCTTACTACGGCACCAACACGCCTATCGACGAATGCTACGAATGC‑GGtt                                                                                                              <  1:309662/251‑1 (MQ=255)
                                                                                                                          gACCGACAAAGGCTACTACACCAAC‑AGCTACCATTT‑GGACGTGGAGAAAAAAGTCAAT‑CCATACG‑ACAAACTGGACTTTGAAATGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTCATTTTCTACGGC‑GAATACCCAAACATCCAGCACAACCTCAAAGCACTGGAAGACGGGTGGGAATACAACCACGACCCCGCGCCCTACTACGGCACCAACACGCCTATCGACGAATGCTACGAATGC‑GGGTTCACCg                                                                                                        >  1:227452/1‑251 (MQ=255)
                                                                                                                              gACAAAGGCTACTACACCAAC‑AGCTACCATTT‑GGACGTGGAGAAAAAAGTCAAT‑CCATACG‑ACAAACTGGACTTTGAAATGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTCATTTTTTACGGC‑GAATAACCCAACCTCCCCCACAACCTCAAAGCAATGGAAGACGTGGGGAATTATAGCTACGACCCCCTGCCCTACTACGGCACCCACACCCCCATTGACGAATGCTACGAATGG‑GGGTTCACCGGC‑gg                                                                                                   >  1:47937/1‑250 (MQ=255)
                                                                                                                                   aGGCTACTACACCAAC‑AGCTACCATTT‑GGACGTGGAGAAAAAAGTCAAT‑CCATACG‑ACAAACTGGACTTTGAAATGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTTATTTTCTACGGC‑GAAAACCCCAACAACCCGCACAACCCCAAAACACTGGAAGACGTGTGGGAATACAGCCACGACCCCGGGCCTTACTACCGCACCAACACCCCCATCGACGAAAGCTACGAATGC‑GGGTTTACCGGC‑G‑GAtttt                                                                                             >  1:53780/1‑250 (MQ=255)
                                                                                                                                    ggCTACTACACCAAC‑AGCTACCATTT‑GGACGTGGAGAAAAAAGTCAAT‑CCATACG‑ACAAACTGGACTTTGAAATGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTCATTTGCTACGGC‑GAATACCCAAACATCCAGCACAACCTCAAAACACTGGAAGACGTGTGGAATTACAGCTACGACCCCGTGCCTTACTACGGCACCAACACGCCCATCGACGAATGCTACGAATGC‑GGTTTCACCGGC‑G‑AATTTGa                                                                                            >  1:175322/1‑251 (MQ=255)
                                                                                                                                                  aacacctccctttgggACGGGGAAAAAAAAGTCAAT‑CCATACG‑ACAAACTGGTTTTTGAAATGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTCATTTGCTACGGC‑GAATACCCAAACATCCAGCACAACCTCAAAGCACTGGAAGACGTGTGGAATTACAGCTACGACCGCGTGCCTTACTACGGCACCAACACGCCTATCGACGAATGCTACGAATGC‑GGTTTCACCGGC‑G‑AATTTGAATGCACCAGCaa                                                                                <  1:324889/237‑1 (MQ=255)
                                                                                                                                                     cctcccATTT‑GGGGGTGGAAAAAAAAGTCAAC‑CCACACA‑ACAACCTGTATTTTAAAGTGG‑CGCCCCCGCCACTTGCCAGCGGCGGTTTCATTTGCTACGGC‑GAATACCCAAACATCCAGCACAACCTCAAAGCACTGGAAGACGTGTGGAATTACAGCTACGACCGCGTGCCTTACTACGGCACCAACACGCCTATCGACGAATGCTACGAATGC‑GGTTTCACCGGC‑G‑AATTTGAATGCACCAGCAAAGGc                                                                            <  1:148113/247‑1 (MQ=255)
                                                                                                                                                      cTACCATGT‑GGACGTGGAGAAAAAAGTCAAC‑CCAGACG‑ACAACCTGGCCTTGGAAAGGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTCATTTGCTACTGC‑GAATACCCAAACATCCAGCACAACCTCAAAGCACTGGAAGACGTGTGGAATTACAGCTACGACCGCGTGCCTTGCTACGGCACCAACACGCCTATCGACGAATGCTACGAATGC‑GGTTTCACCGGC‑G‑AATTTGAATGCACCAGCAAAGGc                                                                            <  1:311523/250‑1 (MQ=255)
                                                                                                                                                         ccATTT‑GGACGTGGAGAAAAAAGTCAAT‑CCATACG‑ACAAACTGGACTTTGAAATGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTCATTTTCTACGGC‑GAATACCCAAACATCCCGCACAACCTCAAAGCACTGGAAGACGTGTGGGATTACAGCTACGACCCCGTGCCTTACTACGGCACCAACACGCCTATCGACGAATGCTACGAAAGC‑GGGTTTACCGGC‑G‑GATTTTAATGCACCAGCAAAGGCCTC‑a                                                                       >  1:12781/1‑251 (MQ=255)
                                                                                                                                                          cATTT‑GGACGTGGAGAAAAAAGTCAAT‑CCATACG‑ACAAACTGGACTTTGAAATGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTTATTTTCTACGGC‑GAATAACCCAACATCCAACACAACCCCAAAACACTGGAAGAAGGGTGGGATTACAGCTAAGACCCCGTGCCTTACTACGGCCCCAACACCCCCATcgac                                                                                                                                  >  1:51446/1‑195 (MQ=255)
                                                                                                                                                           aTTT‑GGACGTGGAGAAAAAAGTCAAC‑CCATACG‑ACAAACTGGTTTTTAAAATGG‑CGTCCCCGCCACTTGCCAGCGGCGGTTTCATTTGCTACGGC‑GAATACCCAAACATCCAGCACAACCTCAAAGCACTGGAAGACGTGTGGAATTACAGCTACGACCGCGTGCCTTACTACGGCACCAACACGCCTATCGACGAATGCTACGAATGC‑GGTTTCACCGGC‑G‑AATTTGAATGCACCAGCAAAGGCTTC‑Acc                                                                     <  1:143688/251‑1 (MQ=255)
                                                                                                                                                                       gaaaaaaaaGTCAAC‑CCATCCG‑ACAACCGGGATTTTGAAATGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTCATTTGCTACGGC‑GAATACCCAAACATCCAGCACAACCTCAAAGCACTGGAAGACGTGTGGAATTACAGCTACGACCGCGTGCCTTACTACGGCACCAACACGCCTATCGACGAATGCTACGAATGC‑GGTTTCACCGGC‑G‑AATTTGAATGCACCAGCAAAGGCTTC‑ACCTG‑CCCGAAAt                                                          <  1:33811/247‑1 (MQ=255)
                                                                                                                                                                         gaAAAAAGTCAAT‑CCATACG‑ACAAACTGGACTTTGAAATGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTCACTTCCTACCGC‑GAATACCCAAAAAACCCGCACAACCCCCAAACACACGAAGACCTGTGGAACTACAGCTACGCCCCCCCGCCCTACTACCGCACCCACACCCCCACCGACGCACGCTACGAAAGA‑GGTTTCACCCGC‑C‑AATTTGTATGCACCCAAAAACACCTC‑ACCTG‑CCCCAAATGCg                                                       >  1:194412/1‑251 (MQ=255)
                                                                                                                                                                                    aT‑CCATACG‑ACAAACTGGAGTTTAAAAGGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTCATTTGCTACGGC‑GAATACCCAAACATCCAGCACAACCTCAAAGCACTGGAAGACGTGTGGAATTACAGCTACGACCGCGTGCCTTACTACGGCACCAACACGCCTATCGACGAATGCTACGAATGC‑GGTTTCACCGGC‑G‑AATTTGAATGCACCAGCAAAGGCTTC‑ACCTG‑CCCGAAATGCGGCAACCACGAc                                            <  1:352651/251‑1 (MQ=255)
                                                                                                                                                                                            cG‑ACAAACTGGACTTTGAAATGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTCATTTTCTACGGC‑GGATACCCAAACATCCAGCACAAACTCAAAACACTGGAAGACGTGTGGAATTACAGCTACGACCCCGGGCCTTACTACGGCACCAACACGCCTATCGACGAATGCTACGAATGC‑GGGTTTACCGGC‑GAAATTTG‑ATGCAACAGCAAAAGCTTC‑ACCCT‑CCCCAAAAGCGGGAACCACGACAGC‑Gaaa                                    >  1:341539/1‑251 (MQ=255)
                                                                                                                                                                                              aaCAACCTGGTCTTTGAAAGGG‑CGTCCCCGCCACTTGCCAGCGGCGGTTTCATTTGCTACGGC‑GAATACCCAAACATCCAGCACAACCTCAAAGCACTGGAAGACGTGTGGAATTACAGCTACGACCGCGTGCCTTACTACGGCACCAACACGCCTATCGACGAATGCTACGAATGC‑GGTTTCACCGGC‑G‑AATTTGAATGCACCAGCAAAGGCTTC‑ACCTG‑CCCGAAATGCGGCAACCACGACAGC‑Gaaaa                                   <  1:461044/250‑1 (MQ=255)
                                                                                                                                                                                                         cTTTGAAATGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTTATTTTCTACGGG‑GAATACCCCAAAAACCAGCACAACCCCAAAACACCGGAAGACGGGTGGGAATACAACTACGACCCCGTGCCTTACTACGGCCCCAACACCCCCATCGACGAATGCCACGAAAGC‑GGTTTCACCCGC‑GCAATTT‑AATGCACCACCAAAAGCTTC‑Accccccccaaa                                                            >  1:261818/1‑208 (MQ=255)
                                                                                                                                                                                                           ttGAAATGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTTATTTTCTAAGGG‑GAAAAACCAAACAACCCACACAACCTCCAAACACTTGAAGACGGGTGGGAATACAACTACCACCCCCCGCCTTTCTACCGCACCAACACGCCCATCGACGAAAGCTACCAAAGC‑GGGTTCACCGGCGG‑AATTT‑AATGCACCCGCAAAAGCCTC‑CCCCGCCCCAAAATGC‑GCAACCCCGACAAC‑Gaaaaaaaccccctcac                       >  1:512170/1‑239 (MQ=255)
                                                                                                                                                                                                                   g‑CGTCCCCGCCACTTGCCAGCGGCGGTTTCATTTGCTACGGC‑GAATACCCAAACATCCAGCACAACCTCAAAGCACTGGAAGACGTGTGGAATTACAGCTACGACCGCGTGCCTTACTACGGCACCAACACGCCTATCGACGAATGCTACGAATGC‑GGTTTCACCGGC‑G‑AATTTGAATGCACCAGCAAAGGCTTC‑ACCTG‑CCCGAAATGCGGCAACCACGACAGC‑GAAAAAGTCTCCGTCACCCGACggg               <  1:551750/250‑1 (MQ=255)
                                                                                                                                                                                                                        aCCCGCCACTTGCCAGCGGCGGTTTCATTTTGTACGGC‑GGATACCCCAACATCCCACACAACCTCCAAACACTGGAAGACGGGTGGGAATACAGCTACGACCCCCGGCCCTACTACGGCACCAACACCCCCATCGACGAAAGCTACGAATGCGGGTTTC‑CCGGC‑G‑GAATTGAATGCACCAGCAAAAGCTTCTCCCTC‑CCCG‑AATGCGGCAACCAAGACAGCGGAAAAAG‑CTCCGTCccccccccgggttg           >  1:52557/1‑240 (MQ=255)
                                                                                                                                                                                                                             ccACTTGCCAGCGGCGGTTTCATTTGCTACGGC‑GAATACCCAAACATCCAGCACAACCTCAAAGCACTGGAAGACGTGTGGAATTACAGCTACGACCGCGTGCCTTACTACGGCACCAACACGCCTATCGACGAATGCTACGAATGC‑GGTTTCACCGGC‑G‑AATTTGAATGCACCAGCaa                                                                                <  1:626149/179‑1 (MQ=255)
                                                                                                                                                                                                                                cTTGCCAGCGGCGGTTTCATTTGCTACGGC‑GAATACCCAAACATCCAGCACAACCTCAAAGCACTGGAAGACGTGTGGAATTACAGCTACGACCGCGTGCCTTACTACGGCACCAACACGCCTATCGACGAATGCTACGAATGC‑GGTTTCACCGGC‑G‑AATTTGAATGCACCAGCAAAGGCTTC‑ACCTG‑CCCGAAATGCGGCAACCACGACAGC‑GAAAAAGTCTCCGTCACCCGACGGGTTTGCGGCTACCt  <  1:243494/251‑1 (MQ=255)
                                                                                                                                                                                                                                  tGCCAGCGGCGGTTTCATTTGCTACGGC‑GAATACCCAAACATCCAGCACAACCTCAAAGCACTGGAAGACGTGTGGAATTACAGCTACGACCGCGTGCCTTACTACGGCACCAACACGCCTATCGACGAATGCTACGAATGC‑GGTTTCACCGGC‑G‑AATTTGAATGCACCAGCAAAGGCTTC‑ACCTG‑ccc                                                               <  1:221374/189‑1 (MQ=255)
                                                                                                                                                                                                                                               ttCATTTGCTACGGC‑GAATACCCAAACATCCAGCACAACCTCAAAGCACTGGAAGACGTGTGGAATTACAGCTACGACCGCGTGCCTTACTACGGCACCAACACGCCTATCGACGAATGCTACGAATGC‑GGTTTCACCGGC‑G‑AATTTGAATGCACCAGCAAAGGCTTC‑Acc                                                                     >  1:468806/1‑171 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                           
CCACCAAACAATGGCGCGCCGAAACGGGTTATGCCTTCAGTTTGTATTCCACTCCAAGCGAAAACCTGTGCGATCGCTTCTGCCGTTTAGATAC‑CAAACAATTCGGCGTCATCGACGGCGTGACCGACAAAGGCTACTACACCAAC‑AGCTACCATTT‑GGACGTGGAGAAAAAAGTCAAT‑CCATACG‑ACAAACTGGACTTTGAAATGG‑CGTACCCGCCACTTGCCAGCGGCGGTTTCATTTGCTACGGC‑GAATACCCAAACATCCAGCACAACCTCAAAGCACTGGAAGACGTGTGGAATTACAGCTACGACCGCGTGCCTTACTACGGCACCAACACGCCTATCGACGAATGCTACGAATGC‑GGTTTCACCGGC‑G‑AATTTGAATGCACCAGCAAAGGCTTC‑ACCTG‑CCCGAAATGCGGCAACCACGACAGC‑GAAAAAGTCTCCGTCACCCGACGGGTTTGCGGCTACCT  >  1/153964‑154432

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