Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,760,3840AC31.1% 41.2 / 13.4 45V289G (GTG→GGG) selAselenocysteine synthase
Reads supporting (aligned to +/- strand):  ref base A (25/6);  new base C (4/10);  total (29/16)
Fisher's exact test for biased strand distribution p-value = 1.70e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGAT‑CATCTCTTTTTTACCAAC‑AATAATTCCTGCC‑TGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTC‑ACCAGACTGACGCCCGCCGC‑AATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGCCTGACGATAATCATTCGCGT  >  NC_000913/3760151‑3760630
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                       
ggATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGGT‑CATCTCTTTTTTACCCAA‑AAACATTCCTGCC‑TGCGGCCCGCCTAAAAAATTTTCGCCGGCAACACTC‑ACCCGACTGAcgcccccc                                                                                                                                                                                                                                        >  2:315912/1‑248 (MQ=255)
     tGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGAT‑CATCGCTTTTTTACCAAC‑AAAAATTCCTGCCCTG‑GGCCCCCCTCACAAATTGTCGCCGGGGAAACTC‑ACCAGACTAGCGCCCGCCCC‑aa                                                                                                                                                                                                                                   >  1:445833/1‑250 (MQ=255)
                  cgctgcgGGTAAGCAGGCGCAGGGTCGGTAATTTTTAACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGAT‑CATCTCTTTTTTACCAAC‑AATAATTTATGCC‑TGCGGCCCGCCTAACAACTTTTCGACGGAGAAACTC‑CCCCCAATGACGCCCCCCCC‑AAACAACTCATGCggg                                                                                                                                                                                                                     >  2:354476/1‑250 (MQ=255)
                   gctgcgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGCGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGAT‑CATCTCTTTTTTACCAAC‑AATAATTCCTGCC‑TGCGGCCCCGCTAACAACTTGTCGCCGGGGGAACTC‑CCCCGACCGACCCCCGCCCC‑aaacacaccctgcggc                                                                                                                                                                                                                     >  2:333342/1‑236 (MQ=255)
                             ggcgggCGCAGGGTCGGTATTTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCCCGCTTCAGCGGGTGGTTTTGCAGGCGGGCGAT‑CATCTCTTTTTTCCCAAC‑AATAATTCCTGCC‑TGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACCC‑ACCAGCCTGACGCCCGCCGC‑AATCAACTCCTGCGGCATTGGCTCttt                                                                                                                                                                                                          <  1:317187/247‑1 (MQ=255)
                                aGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGAT‑CATCTCTTTTTTACCCAC‑AATAATTCCTGCC‑TGCGGCCCGCCTTACAACTTGTCGCCGGCGAAACTC‑CACAGACTGACGCCCGCCGC‑AATCAACTCCTGCGGCCTTGGCTCTTTCgg                                                                                                                                                                                                       >  1:19507/1‑251 (MQ=255)
                                aGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGAT‑CATCTCTTTTTTACCAAC‑AATAATTCCTGCC‑TGCGGCCCGCCCAACAACTTGTCGCCGGAGAAACTC‑ACCAGACTGACGCCCGCCGC‑AATCAAATCCTGCGGCATTGGCTCTTTCgg                                                                                                                                                                                                       >  1:280472/1‑251 (MQ=255)
                                     cAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGAT‑CATCTCTTTTTTACCAAC‑AATAATTCCTGCC‑TGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTC‑ACCAGACTGACGCCCGCCGC‑AATCAACTCCTGCGGCATTGGCTCTTTTGGCCAAc                                                                                                                                                                                                  >  1:204234/1‑251 (MQ=255)
                                     cAGGGCCGGTAATTTTCCCCTAAGAGCTTCGGGGTGAAAATAAAGACGAAACGTGGCTCCCAGCGCCGCGAGGGTCTTTTTACCCGCGCGTAATGCCCGCTTCAGCGGGTGGCTTTGCAGGCGGGCGAT‑CATCTCTTTTTTCCCAAC‑AATAATCCCTGCC‑TGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTC‑ACCAGACTGACGCCCGCCGC‑AATCAACTCCTGCGGCATTGGCTCTTTCGGCAAAc                                                                                                                                                                                                  <  1:178137/251‑1 (MQ=255)
                                       gggTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGAT‑CATCTCTTTTTTACCAAC‑AATAATTCCTGCC‑TGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTC‑ACCAGACTGACGCCCGCCGC‑AATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCg                                                                                                                                                                                                >  1:45182/1‑251 (MQ=255)
                                              tAATTTTTCACTCAGAGCATCAGGGTGTAACTAAAGACGCAACGTGGATTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGAGGCTTTGCAGGCGGGCGGT‑CAACGCTTTTTTACCAAC‑AATAATTCCTGCC‑AGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACCCCACCACA‑TGACGCCCGCCCC‑AAACCACGCCCGCGGCATAGGATCTTTCGCCaaccccttctggc                                                                                                                                                                                         >  2:384240/1‑240 (MQ=255)
                                                aTTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGAT‑CATCTCTTTTTTACCAAC‑AATAATTCCTGCC‑TGCGGCCCGCCTAACAAATTGTCGCCGGAGAAAATC‑ACCAGACTGACGCCCGCCGC‑AATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCt                                                                                                                                                                                        >  1:108453/1‑250 (MQ=255)
                                                 tttttCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGAT‑CATCTCTTTTTTACCAAC‑AATAATTCCTGCC‑TGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTC‑ACCAGACTGACGCCCGCCGC‑AATCAACTCCTGCGGCATTGGCTCTTTTGGCAAACCGTACTGGCTa                                                                                                                                                                                       >  1:9771/1‑250 (MQ=255)
                                                    ttccctaagagCTTCAGGGTGAAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGGGGGTCATTTTACCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGAT‑CATCTCTTTTTTCCCAAC‑AATAATTCCTGCC‑TGCGGCCCGCCTAACACCTTGTCGCCGGAGAAACTC‑ACCAGACTGACGCCCGCCGC‑AATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGa                                                                                                                                                                                    <  1:270413/243‑1 (MQ=255)
                                                         ccAGAGCTTCAGGGTGAAAATAAAGACGCAACGTGGCTCCCAGCGCCGCGAGGGTCTTTTTATCCGCGCGTAATGCACGCTTCAGCGGGGGGCTTTGCGGGCGGGCGAT‑CATCTCTTTTTTACCAAC‑AATAATCCCTGCC‑TGCGGCCCGCCTAACACCTGGTCCCCGGAGAAACTC‑CCCAGACTGACCCCCGCCGC‑AATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGAcc                                                                                                                                                                              <  1:434493/250‑1 (MQ=255)
                                                                                     cAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGAT‑CATCTCTTTTTTACCAAC‑AATAATTCCTGCC‑TGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTC‑CCCAGACTGACGCCCCCCGC‑AATCAACTCCTGCGGCATTGGCTCTTTCTGCCAACCGTACTGGCCAAGATCGACCCGCGAGCCACTGCCCAAATCAGTCAcc                                                                                                                                                   >  2:127829/1‑249 (MQ=255)
                                                                                         gTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGAT‑CCTCTCCTTTTTACCAAC‑AATAATTCCTGCC‑TGGGGCCCGCCTAAAAACGTGTCGCCGGGGCAAATC‑ACCAGACTTACGCCCGCCGCCAAT‑AACTCCTGCCGCATTGGGCCTTTCGGCCAACCCTCCTGGCCACGATAGACCCGCGGGCCAcctcctcaatccaaccacacggg                                                                                                                                              >  2:234518/1‑229 (MQ=255)
                                                                                           ggCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGTT‑CATCTCTTTTTTACCAAC‑AATAATTCCTTGC‑TGCGGCCCGCCTAACAAATTTTCGCCGGAGAAAATC‑ACCAGACTGACTCCCGCCCC‑AATCCACTCCTTCGGAATTTGCCCTTTCGGCAAACCGTAACGGCTCAGAACGACCACCGACCCACCGCCtcacaccaccactcatgga                                                                                                                                             >  2:405753/1‑232 (MQ=255)
                                                                                                  aGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGAT‑CATCTCTTTTTTAACAAC‑AATAATTCCTGCC‑TGCGGCCCGCCTAACAAATTGTCGCCGGAGAAACTC‑ACCAGACTGACGCCCGCCGC‑AATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATcc                                                                                                                                       >  1:167592/1‑249 (MQ=255)
                                                                                                              gtctttttACCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGAT‑CATCTCTTTTTTACCAAC‑AATAATTCCTGCC‑TGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTC‑CCCAGACTGACGCCCGCCGC‑AATCAACTCCTGCgg                                                                                                                                                                                                                      <  1:357700/154‑1 (MQ=255)
                                                                                                                 aTTTAATCCGCGTTTAATGCACGCTTCATGGGGTTGCTTTGCGGGCGGG‑GATTCTTCTTTTTTTCACCAAC‑AATAATTCCTGCT‑TGGGGCCCGCCTAACATCTTTCCGCCGAAGAACCTC‑ACCAGACTGCCGCCCGCCGC‑AATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCTCTACGGGAACATCCAGCTCTTTGCCGAgcg                                                                                                                       <  1:162147/250‑1 (MQ=255)
                                                                                                                      tcccGCGCGTAGTGCACCTTTCGGGGGGTGTCTTTGCGGGCGGTCGAT‑CTTCTCTTTTTCACCAAC‑AATATTCCCTCCC‑TGCGCCCCCCCTAACACCTTGCCGCCGAAGAACCCC‑CCCAGACTGCCCCCCGCCGC‑AATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGccacc                                                                                                                  <  1:85100/248‑1 (MQ=255)
                                                                                                                              gaaaggcacgtctCAGGGGGTGGCTTTGCGGGCGGGCGAT‑CTTCT‑TTTTTTACCAACCAATAATTCCTCCC‑TGCGGCCCCCCTACCAACTTGTCGCCGGAGAACCCC‑CCCAGACTGACGCCCGCCGC‑AATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCt                                                                                                         <  1:66173/239‑1 (MQ=255)
                                                                                                                              gTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGAT‑CAACTCTTTTTTACCAAC‑AATAATTCCTGCC‑TGCGGCCCGCCTAACAACTTGTCGCCGGAGAAAATC‑ACCAGACTGACGCCCGCCGC‑AATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATACAGCTCTTTGCCGAGCGCCACCAGTTCCGCt                                                                                                         >  1:520240/1‑251 (MQ=255)
                                                                                                                                 aTGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGAT‑CATCTCTTTTTTACCAAA‑AATAATTCCTGCC‑TGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTC‑ACCAGACTGACGCCCGCCGC‑AATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGCGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCCGTTCCGCTTcc                                                                                                      >  1:313545/1‑250 (MQ=255)
                                                                                                                                     aCGCTTCAGCGGGTGGCTTTGCAGGCGGGCGGT‑CATCTCTTTTTTTCCAAC‑AATAATTCCTGCC‑TGCGGCCCGCCTAAAAAATTGTCGCCGGAGAAACTC‑ACCCGACTGGCGCCCGCCCC‑AATCAACTCCTGCGGGATTGGCTCTTTCGGCAAACCGTACTGGCGAAGATCGCCCAACGAACCACTGCCTAAATCAGTCACTACGGGGACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCt                                                                                                   >  1:264590/1‑250 (MQ=255)
                                                                                                                                            aGCGGGTGGCTTTGCAGGCGGGCGAT‑CATCTCTTTTTTACCAAC‑AATAATTCCTGCC‑TGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTC‑ACCCGACTGTCGCCCGCCGC‑AATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGGACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGGACATCCAGCTCTTTTCCGAGCGCCCCCAGTTCCGCCTCATCTATCGCttt                                                                                           >  2:115658/1‑251 (MQ=255)
                                                                                                                                                          cAGGCGGGCGAT‑CATCTCTTTTTTACCAAC‑AATAATTCCTGCC‑TGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTC‑ACCAGACTGACGCCCGCCGC‑AATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGa                                                                              >  1:200354/1‑250 (MQ=255)
                                                                                                                                                            ggctggcgAT‑CATCTATTTTTTACCAAG‑AATAATTCCTGCC‑TGCGGCCCGCCTAACAACTTGTCGCCGGAGACACTC‑ACCAGACTGACGCCCGCCGC‑AATCAACTCCTGCGGCATTGCCTCTTTCGGCACACCGTACTGGATAAGATCGCCCAACGAGCCACTGCCTAAATCAGGCACTACGGGACCAGCCAGCTCTTTGCCGGGCGCACCCCGTTCCGCTTCATCTATCGCATTGGTGAACCCCTGactt                                                                           >  2:434446/5‑248 (MQ=255)
                                                                                                                                                               gggcgCT‑CATCTTTTTTTTACCAAC‑AATAATCCCTCCC‑TGCGGCCCCCCTAACAATTTGTCGCCGGAGAAACCC‑CCCAGCCTGACGCCCGCCGC‑AATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTg                                                                        <  1:438601/251‑1 (MQ=255)
                                                                                                                                                                   gAT‑CATCTCTTTTTTACCAAC‑AATAATTCCTGCC‑TGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTC‑ACCAGACTGACGCCCGCCGC‑AATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGt                                                                    >  1:523510/1‑251 (MQ=255)
                                                                                                                                                                        atcTCTTTTTTACCAAC‑AATAATTCCTGCC‑TGCGGCCCGCCTAACAACTTTTCGCCGGAGAAACTC‑ACCAGACTGACGCCCGCCGC‑AATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACAACCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACt                                                                >  2:108047/1‑251 (MQ=255)
                                                                                                                                                                              tttttccaaac‑aaTAATCCCTGCC‑TGCGGCCCGCCTAACAACTTGTCGCCGGAGAACCCC‑CCCAGCCGGACGCCCGCCGC‑AATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTAt                                                           <  1:315912/242‑1 (MQ=255)
                                                                                                                                                                                  tACCAAC‑AATAATTCCTGCC‑TGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTC‑ACCAGACTGACGCCCGCCGC‑AATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGATTGCTGTAGTTACTGGTATGTACt                                                      >  1:7584/1‑251 (MQ=255)
                                                                                                                                                                                    ccaaaaaataatCCCTGCC‑TGCGCCCCGCCTAACACCTTGCCCCCGGAAAAACCC‑ACCAGACTGCCGCCCGCCGC‑AATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACtt                                                     <  1:354476/245‑1 (MQ=255)
                                                                                                                                                                                           ataatTCCTGCC‑TGCGGCCCCCCTAACACCTGGCCGCCGGAGAACCTC‑CCCAGACTGACGCCCGCCGC‑AATCAACTCCTGCGGCATCGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCaa                                              <  1:333342/251‑1 (MQ=255)
                                                                                                                                                                                             aatTCCTGCC‑TGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTC‑ACCAGACTGACGCCCGCCGC‑AATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAAc                                             >  2:133014/1‑250 (MQ=255)
                                                                                                                                                                                                  ctgcc‑tgcGGCCCGCCTAACAACTTGTCGCCGGAGAAACTC‑ACCAGACTGACGCCCGCCGC‑AATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCg                                       >  1:47963/1‑251 (MQ=255)
                                                                                                                                                                                                     cc‑tgcGGCCCGCCTAACAACTTGTCGCCGGAGAAACTC‑ACCAGACTGACGCCCGCCGC‑AATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGt                                     >  1:396665/1‑250 (MQ=255)
                                                                                                                                                                                                         gcGGCCCGCCTAACACCTTGCCGCCGGAAAACCTC‑CCCAGACTGACGCCCGCCGC‑AATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTAt                                   <  1:78687/249‑1 (MQ=255)
                                                                                                                                                                                                                               gCCGGAGAAACTC‑CCCAGACTGACGCCCGCCGC‑AATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGCCTGACGATAAt           <  1:231809/251‑1 (MQ=255)
                                                                                                                                                                                                                                  ggAGAACCTC‑ACCAGACTGACGCCCGCCGC‑AATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGCCTGACGATAATCAt        <  1:127829/251‑1 (MQ=255)
                                                                                                                                                                                                                                   gagaAACCC‑CCCAGACTGACGCCCGCCGC‑AATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGCCTGACGATAATCAtt       <  1:133014/251‑1 (MQ=255)
                                                                                                                                                                                                                                       aaCTC‑ACCAGACTGACGCCCGCCGC‑AATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGCCTGACGATAATCATTcgcg   >  2:97197/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aCTC‑ACCAGACTGACGCCCGCCGC‑AATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTAACTACGGGAACATCCAGCTATTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAAAAGTGCGGTATTTTCATTCACCGCCTGACGATAATCATTAgcgt  >  1:297259/1‑251 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                       
GGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGAT‑CATCTCTTTTTTACCAAC‑AATAATTCCTGCC‑TGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTC‑ACCAGACTGACGCCCGCCGC‑AATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGCCTGACGATAATCATTCGCGT  >  NC_000913/3760151‑3760630

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