Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,760,3840AC32.8% 110.8 / 10.1 65V289G (GTG→GGG) selAselenocysteine synthase
Reads supporting (aligned to +/- strand):  ref base A (40/3);  new base C (3/18);  total (44/21)
Fisher's exact test for biased strand distribution p-value = 4.04e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑‑AACAATAATTCCTGCCTGCGGCCCGCCT‑AACAACTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCG  >  NC_000913/3760147‑3760609
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                  
gCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑‑AACAATAATTCCTGCCTTCGGCCCGCCT‑AACAACTTGTCGCC‑GGAGAAACTCa                                                                                                                                                                                                                                    >  1:2144923/1‑238 (MQ=255)
       tttGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑‑AAAAATAATTCCTGCCTGCGGCCCGCCT‑AACAAATTTTCGCC‑GGGGAAACTCACCAg                                                                                                                                                                                                                                >  1:1581169/1‑235 (MQ=255)
                   ctgcgctgcgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGGTCATCTCTTTTTTACC‑‑AACAATAATTCCTGCCTGCGGTCCGCCT‑AACAACATTTCGCC‑CGGGTAAATCACCCGCATTAcgccccgccc                                                                                                                                                                                                                 >  2:923236/1‑233 (MQ=255)
                     gcgctgcgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGGTCATCTCTTTTTTACC‑‑AACAAAAATTCCTGCCCGCGGCCCGCCC‑CAAAACTTTTCGCC‑GGAGAAACACACCAGACCGAcgcccgcccccat                                                                                                                                                                                                              >  2:2239666/1‑234 (MQ=255)
                           cgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGGTCATCTCTTTTTTTCC‑‑AACAATAATTCCTGCCTGCGGCCCGCCT‑AACAACTTGTCGCC‑GGAGAAACTCACCCGACTCACGCCCGCCGCAATCCACt                                                                                                                                                                                                         >  2:97170/1‑238 (MQ=255)
                           cgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTATC‑‑CACAATAATTCCTGCCTGCGGCCCGCCT‑AACAACTTTTTGCC‑GGAGAAAATCACCAGACTGACGCCCGCCGCAATCAACt                                                                                                                                                                                                         >  1:2235754/1‑238 (MQ=255)
                           cgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑‑AACAATAATTCCTGCCTGCGGCCCGCCT‑AACAACTTTTCGCC‑GGAGAAAATCACCAGACTGACGCCCGCCGCAATCAACt                                                                                                                                                                                                         >  1:1656371/1‑238 (MQ=255)
                           cgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑‑AACAATAATTCCTGCCTGCGGCCCGCCT‑AACAAATTTTTCCC‑GGAGAAACTCACCAGACTGACCCCCGCCGCAATCAACt                                                                                                                                                                                                         >  1:93133/1‑238 (MQ=255)
                                        gcAGGGGCGGTAATTTTTCACTCAGAGCTTCAGGGAGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGTGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGAGGATTTGCAGGCGGGCGATCCTCCCTTTTTTTCC‑‑AAGCACAATTCCTCCCTGGGGCCCACCT‑AACAACTTTCCCCC‑GTCCGAACTCTACAGAACACCACCCGCCGCAAACAACTCCTGCGGGCTTgg                                                                                                                                                                                            >  2:227082/1‑238 (MQ=255)
                                                tgTATTTTTTCACTAAGAGCTCCAGGGTGAAAATAAAGCCGCAAGGTGGCTCCCAGCGCGGGGAGGTTCTTTTTACCCGCTCTTAAGCCACGCTTCGGGGGGTGGTTTTGCGGGCGGGCGACCATCTTTTTTTTACA‑‑AAAAATAATCCCTGCCTGGGCCCCGCCA‑AACACCTTGCCGCG‑GGAAAACCCCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCg                                                                                                                                                                                    <  2:614654/237‑1 (MQ=255)
                                                      ttttCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTTCC‑‑AACAATAATTCCTGCCTGCGGCCCGCCC‑AACAACCTTTTGCC‑GGAGCAACTCACCAGACTGACGCCCGCCCCCATCAACTCCTGCGGCATCGGCCCTTTCGGCcaaac                                                                                                                                                                             >  2:1858669/1‑237 (MQ=255)
                                                      ttttCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑‑AACAATAATTCCTTCCTGCGGCCCGCCT‑AACAACTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCCACTCCTGCGGCATTGGCTCTTTCGGCCAAcc                                                                                                                                                                             >  1:1831394/1‑239 (MQ=255)
                                                        ttCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑‑AACAATAATTCCTGCCTGCGGCCCGCCT‑AACAACTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCCCTTTCGGCAAAcccg                                                                                                                                                                           >  1:614654/1‑237 (MQ=255)
                                                         tcccTCAGAGCTTCAGGGTTAAAAAAAAGCCGCAAGGTGGCTCCCGGCGCCGGGGGGGTCTTTTTATCCGCGGTTAATGCCCGCTTCAGCGGGTGGTTTTGCGGGCGGGCGTTCATCTTTTTTTTCCC‑‑AAAAAAAATCCCTGCCGGCGCCCCCCCT‑AACAACTTGTCGCC‑GAAAAAACCCCCCAGCCTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑a                                                                                                                                                                         <  2:1041159/236‑1 (MQ=255)
                                                            cTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑‑AACAATAATTCCTGCCTGCGGCCCGCCT‑AACAACTTGTCGCC‑GGAGAAACTCACCAGACCGACGCCCGCCCCAATCAACTCCTCCGGCCTTGGCTCTTTCGGCCAAcccg                                                                                                                                                                           >  1:1986586/1‑233 (MQ=255)
                                                                           tgtAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGTTCACCTCTTTTTTACC‑‑AACAAAAATTTCTTCCTGCGGCCCGCCT‑AACAACTTTTCGCC‑GGAGAAACTCAACACACTTACGCCCGCCGCCATCAACTCCTGCCGGATTGGCTCTTTCTGCCAAACGT‑ACTCGCCAAGATCGACCa                                                                                                                                                        >  2:552925/1‑238 (MQ=255)
                                                                              aaataaaGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑‑AACAATAATTCCTGCCTGCGGCCCGCCT‑AACAACTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGATCGACCAGCg                                                                                                                                                     >  1:1428176/1‑238 (MQ=255)
                                                                                       cGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGGGGCTTTGCAGGCGGGCGCTCATCTCTTTTTTACC‑‑AACCAAAATTCCTGCCCGGGGGCCGCCG‑AACAACTTTGCCCC‑GGAGAAAATCACCACAACGACTCCCGCCGCAAACACATCCCCCGGCAATGGCTCTTTTGGCAAAACGTAACTGCCCAA‑AACCACCACCGAGCCAcaccc                                                                                                                                           >  2:1328576/1‑235 (MQ=255)
                                                                                          aaCGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCCTCTCTTTTTTACC‑‑AACAATAATTCCTGCCTGCGGCCCGCCT‑AACAACTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCCTTGGCCCTTTCGGCAAACCGT‑CCTGGCTAAGa                                                                                                                                                               >  1:1744624/1‑216 (MQ=255)
                                                                                          aaCGTGGATTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGGGATCATCTCCTTTTTGCC‑‑AAAAATAATTCCCGCCTGCCGCCCGCCT‑CAAAACATGTCGCC‑GGAGAAACTCCCCAGACTGACGCCCCCCCCAATCCACTCCTCCGGCCGTCGCCCTTTCGGCAAAACGA‑ACAGGCCAAAACAGACCCCCCAGCCAcctcctta                                                                                                                                        >  2:1569492/1‑232 (MQ=255)
                                                                                            gggcgtttCCGGCCCCGCGGGGTTTTTTTTACCCGGCAGTAAGGCCCGTTTCGGGGGGGGGTTTGGCGGGGGGCCAATCTTCT‑TTTTTTACC‑‑AAAAAATATTCCCGGCCGGGGCCCCCCTAAACAACTTGCCGCC‑GGAAAACCCCCCCAGCCTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGATCGACCAGCGAGCCa                                                                                                                                                <  2:1889496/223‑1 (MQ=255)
                                                                                                  ttCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑‑CACAATAATTCCTTCCTTCGGCCCGCCT‑AACAACTTGTCGCC‑GGGGAAACTCCCCAGACCGGCGCCCGCCGCAATCAACCCCTGCGGCATTGGCCCTTTTCGCTAACCGG‑ACTGGCTAAAAACGGCCAGCGACCCACCGCCTCaatcaatc                                                                                                                                 >  2:1220990/1‑238 (MQ=255)
                                                                                                        cgccgcGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑‑AACAATAATTCCTGCCTGCGGCCCGCCT‑AACAACTTTTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCCAACCGT‑ACTGGCTAAGATCGACCCGCGAGCCACTGCCTAAAt                                                                                                                                      >  1:1889496/1‑227 (MQ=255)
                                                                                                           cgcgAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCCGTCATCTCTTTTTTTCC‑‑AACAATAATTCCTGCCTGCGGCCCGCCC‑AACAACTTGTCGCC‑GGAGAAACTCACCCGACTGACCCCCGCCCCAATCCACCCCTGCCGCCTTTGCTCTTTCGGCCAACCGC‑ACTGGCTAAGAACGAACACCGAGCCACCGCCTCAATCAGTCACCACgggga                                                                                                                       >  2:2127908/1‑237 (MQ=255)
                                                                                                                    cATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGGTCATCTCTTTTTTTCC‑‑AACAATAATTTCCGCCCTGGGCCCGCCC‑CACAACTTGGTCCC‑GGAGAAAATCCACCGACCGACGCCCGCCGCAATCAAATCCTGCGGGGTTGGCTCCTTCGGCAAAACGT‑AATGGCTAAAATCGACCCACGAGCCCCAGCCCATAACAGTCACTACAGGGAAATCCCGc                                                                                                               >  2:205839/1‑238 (MQ=255)
                                                                                                                           tCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGGTCACCTCTTTTTTTAC‑‑CACAAAAATTTCTGCCTGCGGCCCGCCT‑AACAAACTGTCGCC‑TGAGAAAATCACCAGACTGACGCCCCCCGCAATCAAACCCTGCGGCATTGGCTCTTTTGGCAAACCGT‑ACTGGCTTACATAGACAAACGAGACCCCGCCGAAATAAATCACCAATGGGACATACATATCTTTTcc                                                                                                       >  2:1183347/1‑239 (MQ=255)
                                                                                                                                 cgTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑‑AACAATAATTCCTGCCTGCGGCCCGCCT‑AACAACTTTTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAgcg                                                                                                  >  1:1819086/1‑238 (MQ=255)
                                                                                                                                        cccGCTCCCGCGGGTGCTTTTGCGGGGGGCCATTCCTTTTTTTTTCACC‑‑AAAAAAATTCCCTGCCGCCGCCCCCCAT‑AACCACTTGCCGCC‑GAAAAACCCCCCCAAACGACCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAgc                                                                                                   <  2:1693542/228‑1 (MQ=255)
                                                                                                                                              tCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑‑AACAATAATTCCTGCCTGCGGCCCGCCT‑AACAACTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGATCGAc                                                                                                                                                          >  1:135892/1‑169 (MQ=255)
                                                                                                                                                aGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTTCC‑‑CACAATAATTCCTGCCTTCGGCCCGCCT‑AACAACTTGTCGCC‑GGAGAAACTCACCAGACCGACGCCCGCCGCAATCAACTCCTCCGGCATTGGCTCTTTCGGCCAAACCT‑ACTGGCTAAGATCGGCCAGCCAGCCCCCGCCCAAATCAGCTACCAAGGGAACCACCAGCCCTTTTCCCAGCGCCCCCACTTTCGcct                                                                                   >  2:1689931/1‑236 (MQ=255)
                                                                                                                                                  cGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTTCC‑‑AACAATAATTCCTGCCTGCGGCCCGCCT‑AACAACTTTTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCCGCAAACCGT‑ACTGGCTAAGAACGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACAACCAGACCTTTGTCCAGCGCCCCCAGATcc                                                                                       >  2:1802090/1‑232 (MQ=255)
                                                                                                                                                  cGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑‑AACAATAATTCCTTCCTGCGGCCCGCCT‑AACAACTTTTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCCGtt                                                                                         >  1:1693542/1‑230 (MQ=255)
                                                                                                                                                       ggCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑‑AACAATAATTCCTGCCTGCGGCCCGCCT‑AACAACTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATc                                                                           >  1:1816410/1‑239 (MQ=255)
                                                                                                                                                            tGCCGGGGGGCCATCATCTTTTTTTTACC‑‑AAAAAACATCCCCCCCCGCCGCCCCCCT‑AAACACTGTGCGCCGGGAAAACCCCACCAGCCGACCCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCtt                                                                       <  2:505964/239‑1 (MQ=255)
                                                                                                                                                                 gggggcgcttctcttTTTTTTCCA‑‑AAAAATAACCCCCCCCCGCGGCCCCCAA‑AACATTTTGCCGCG‑GGAAAACCCCCCCACACTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTg                                                                  <  2:2159718/227‑1 (MQ=255)
                                                                                                                                                                 gcgggcgACCTTCTCTTTTTACCC‑‑AAAAATAATTCCTCCCGGCGGCCCCCCT‑AACAACTTGCCGCG‑GGAAAAACTCCCCAGACTGACCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTg                                                                  <  1:2239666/238‑1 (MQ=255)
                                                                                                                                                                 gcgggcgAACTTTTTTTTTAAACA‑‑AACAATAATCCCGGCCTGCGCCCCGCAA‑ACCAACTTGCCGCG‑GAAAAACCCCCCCAGCCGGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTg                                                                  <  2:753329/238‑1 (MQ=255)
                                                                                                                                                                      cgATCATCTCTTTTTTACC‑‑AACAATAATTCCTGCCTGCGGCCCGCCT‑AACAACTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCCATCAACCCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGATCTACCAGCGAGCCACCGCCTAAATCAGTCACTACCGGAACATCCAGCCCTTTGCCGACCGCCCCCAGTTCCGCTTCATCTATCGCCTTGGTGAAcc                                                              >  2:1081619/1‑237 (MQ=255)
                                                                                                                                                                         caatcTTTTTTTCACC‑‑AAAAAAATTCCCTGCCGGCGGCCCCCCT‑AACAACTAGCCGCC‑GGAAAACCCCACCAGACGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGa                                                         <  1:2100868/237‑1 (MQ=255)
                                                                                                                                                                           tcctctTTTTTACCCAAACAATATTCCCGCCCGGCGCCCCGCCA‑AAAACCTGGTCCCG‑GGAAAAACTCCCCAGCCGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGa                                                         <  2:1986586/237‑1 (MQ=255)
                                                                                                                                                                              ttttttttCCA‑‑AAAAATAATCCCTGCCTGCGGCCCCCCT‑AACAACTTGGCGCG‑GGAGAAACTCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCt                                                    <  1:944979/237‑1 (MQ=255)
                                                                                                                                                                                ttttttACC‑‑AACAATAATTCCTGCCTGCGGCCCGCCT‑AACAACTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAAt                                                                                                                                      >  1:1547747/1‑155 (MQ=255)
                                                                                                                                                                                   tttACC‑‑AACAATAATTCCTGCCTGCGGCCCGCCT‑AACAACTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGt                                               >  1:1766583/1‑239 (MQ=255)
                                                                                                                                                                                       tccaaaaaattatCCCTCCCGGCGCCCCCCCT‑AACATCTTGTCGCG‑GGAGAACCCCCCCAACCTCCCCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTa                                             <  1:205839/232‑1 (MQ=255)
                                                                                                                                                                                            acaaTATTCCCTGCGTGCGGCCCGCCA‑AACATCTTGGCCCC‑GGAAAACCCCCCCAGACGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGt                                        <  1:97170/239‑1 (MQ=255)
                                                                                                                                                                                            acaaTAATTCCTGCCTGCGGCCCGCCT‑AACAACTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGt                                        >  1:1808670/1‑239 (MQ=255)
                                                                                                                                                                                            acaaTAATTCCTGCCTGCGGCCCGCCT‑AACAACTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGCTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGt                                        >  1:424369/1‑239 (MQ=255)
                                                                                                                                                                                            aaaataatCCCTCCCGGCGGCCCGCCT‑AACACCTTGTCGCG‑GGAGAAACCCCCCAGGCGGGCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGt                                        <  1:1569492/237‑1 (MQ=255)
                                                                                                                                                                                              aataatTCCTGCCTGCGGCCCGCCT‑AACAACTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTAt                                      >  1:1817395/1‑239 (MQ=255)
                                                                                                                                                                                               ataatCCCTGCCGGCCCCCCGCAT‑AACAATTGGGCGCC‑GGAAAACCGCACAATGCGGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATg                                     <  2:1674832/239‑1 (MQ=255)
                                                                                                                                                                                                 attcccggccggcGGCCCCCCA‑AACAACTTGCCGCG‑GGAAAACCCCCCAAGGCTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                   <  1:923236/232‑1 (MQ=255)
                                                                                                                                                                                                   tccctgcctgCGCCCCCCCA‑AACACCTTGCCCCC‑GGAGAACCCCCCCACCCGAACCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑a                                                                                                                                                                         <  2:787776/101‑1 (MQ=255)
                                                                                                                                                                                                    tCCTGCCTGCGGCCCGCCT‑AACAACTTGTCGCC‑GGAGAAACTCa                                                                                                                                                                                                                                    >  2:101784/1‑44 (MQ=255)
                                                                                                                                                                                                    tCCTGCCTGCGGCCCGCCT‑AACAACTTGTCGCC‑GGAGAAACTCa                                                                                                                                                                                                                                    <  1:101784/44‑1 (MQ=255)
                                                                                                                                                                                                    ccctgcctgCGCCCCCCCA‑AACAACTTGTCGCC‑GGAAAACCCCCCCAGCCGGACCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑a                                                                                                                                                                         <  2:1420018/101‑1 (MQ=255)
                                                                                                                                                                                                     cctgcctgCGGCCCGCCT‑AACAACTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGt                                    >  1:1240894/1‑234 (MQ=255)
                                                                                                                                                                                                        gcctgcGGCCCGCCT‑AACAACTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTtca                             >  2:1647666/1‑238 (MQ=255)
                                                                                                                                                                                                        gcctgcGGCCCGCCT‑AACAACTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTtca                             >  2:1495516/1‑238 (MQ=255)
                                                                                                                                                                                                             cGGCCCGCCT‑AACAACTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGa                                                                                                                                                               >  1:787776/1‑103 (MQ=255)
                                                                                                                                                                                                              ggCCCGCCT‑AACAACTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGa                                                                                                                                                               >  1:1420018/1‑102 (MQ=255)
                                                                                                                                                                                                                     cT‑AACAACTTGCCGCC‑GGAAAACCCCCCCAGACGGCCCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑a                                                                                                                                                                         <  2:627024/85‑1 (MQ=255)
                                                                                                                                                                                                                          ccactttgcGCC‑GGAAAAACTCCCCAGCCTGGCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTAtttt           <  1:1245359/231‑1 (MQ=255)
                                                                                                                                                                                                                           tcctggTCGCC‑GAAAAACCCCCCAAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                   <  2:1897480/209‑1 (MQ=255)
                                                                                                                                                                                                                           aacTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGAATTttc          >  1:1032628/1‑239 (MQ=255)
                                                                                                                                                                                                                                gTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGa                                                                                                                                                               >  1:627024/1‑85 (MQ=255)
                                                                                                                                                                                                                                 tCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGAGGTATTttcattca     >  1:596656/1‑238 (MQ=255)
                                                                                                                                                                                                                                     c‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCaa                         >  1:1897480/1‑214 (MQ=255)
                                                                                                                                                                                                                                     c‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCg  >  1:1276281/1‑237 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                  
GCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑‑AACAATAATTCCTGCCTGCGGCCCGCCT‑AACAACTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGT‑ACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCG  >  NC_000913/3760147‑3760609

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