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‑‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCG  >  NC_000913/3760147‑3760609
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                  
gCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑‑AACAATAATTCCTGCCTTCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCa                                                                                                                                                                                                                                    >  1:2144923/1‑238 (MQ=255)
       tttGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑‑AAAAATAATTCCTGCCTGCGGCCCGCCTAACAAATTTTCGCCGGGGAAACTCACCAg                                                                                                                                                                                                                                >  1:1581169/1‑235 (MQ=255)
                   ctgcgctgcgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGGTCATCTCTTTTTTACC‑‑AACAATAATTCCTGCCTGCGGTCCGCCTAACAACATTTCGCCCGGGTAAATCACCCGCATTAcgccccgccc                                                                                                                                                                                                                 >  2:923236/1‑233 (MQ=255)
                     gcgctgcgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGGTCATCTCTTTTTTACC‑‑AACAAAAATTCCTGCCCGCGGCCCGCCCCAAAACTTTTCGCCGGAGAAACACACCAGACCGAcgcccgcccccat                                                                                                                                                                                                              >  2:2239666/1‑234 (MQ=255)
                           cgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGGTCATCTCTTTTTTTCC‑‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCCGACTCACGCCCGCCGCAATCCACt                                                                                                                                                                                                         >  2:97170/1‑238 (MQ=255)
                           cgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTATC‑‑CACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTTTTGCCGGAGAAAATCACCAGACTGACGCCCGCCGCAATCAACt                                                                                                                                                                                                         >  1:2235754/1‑238 (MQ=255)
                           cgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTTTCGCCGGAGAAAATCACCAGACTGACGCCCGCCGCAATCAACt                                                                                                                                                                                                         >  1:1656371/1‑238 (MQ=255)
                           cgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAAATTTTTCCCGGAGAAACTCACCAGACTGACCCCCGCCGCAATCAACt                                                                                                                                                                                                         >  1:93133/1‑238 (MQ=255)
                                        gcAGGGGCGGTAATTTTTCACTCAGAGCTTCAGGGAGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGTGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGAGGATTTGCAGGCGGGCGATCCTCCCTTTTTTTCC‑‑AAGCACAATTCCTCCCTGGGGCCCACCTAACAACTTTCCCCCGTCCGAACTCTACAGAACACCACCCGCCGCAAACAACTCCTGCGGGCTTgg                                                                                                                                                                                            >  2:227082/1‑238 (MQ=255)
                                                tgTATTTTTTCACTAAGAGCTCCAGGGTGAAAATAAAGCCGCAAGGTGGCTCCCAGCGCGGGGAGGTTCTTTTTACCCGCTCTTAAGCCACGCTTCGGGGGGTGGTTTTGCGGGCGGGCGACCATCTTTTTTTTACA‑‑AAAAATAATCCCTGCCTGGGCCCCGCCAAACACCTTGCCGCGGGAAAACCCCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCg                                                                                                                                                                                    <  2:614654/237‑1 (MQ=255)
                                                      ttttCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTTCC‑‑AACAATAATTCCTGCCTGCGGCCCGCCCAACAACCTTTTGCCGGAGCAACTCACCAGACTGACGCCCGCCCCCATCAACTCCTGCGGCATCGGCCCTTTCGGCcaaac                                                                                                                                                                             >  2:1858669/1‑237 (MQ=255)
                                                      ttttCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑‑AACAATAATTCCTTCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCCACTCCTGCGGCATTGGCTCTTTCGGCCAAcc                                                                                                                                                                             >  1:1831394/1‑239 (MQ=255)
                                                        ttCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCCCTTTCGGCAAAcccg                                                                                                                                                                           >  1:614654/1‑237 (MQ=255)
                                                         tcccTCAGAGCTTCAGGGTTAAAAAAAAGCCGCAAGGTGGCTCCCGGCGCCGGGGGGGTCTTTTTATCCGCGGTTAATGCCCGCTTCAGCGGGTGGTTTTGCGGGCGGGCGTTCATCTTTTTTTTCCC‑‑AAAAAAAATCCCTGCCGGCGCCCCCCCTAACAACTTGTCGCCGAAAAAACCCCCCAGCCTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                         <  2:1041159/236‑1 (MQ=255)
                                                            cTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACCGACGCCCGCCCCAATCAACTCCTCCGGCCTTGGCTCTTTCGGCCAAcccg                                                                                                                                                                           >  1:1986586/1‑233 (MQ=255)
                                                                           tgtAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGTTCACCTCTTTTTTACC‑‑AACAAAAATTTCTTCCTGCGGCCCGCCTAACAACTTTTCGCCGGAGAAACTCAACACACTTACGCCCGCCGCCATCAACTCCTGCCGGATTGGCTCTTTCTGCCAAACGTACTCGCCAAGATCGACCa                                                                                                                                                        >  2:552925/1‑238 (MQ=255)
                                                                              aaataaaGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCg                                                                                                                                                     >  1:1428176/1‑238 (MQ=255)
                                                                                       cGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGGGGCTTTGCAGGCGGGCGCTCATCTCTTTTTTACC‑‑AACCAAAATTCCTGCCCGGGGGCCGCCGAACAACTTTGCCCCGGAGAAAATCACCACAACGACTCCCGCCGCAAACACATCCCCCGGCAATGGCTCTTTTGGCAAAACGTAACTGCCCAAAACCACCACCGAGCCAcaccc                                                                                                                                           >  2:1328576/1‑235 (MQ=255)
                                                                                          aaCGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCCTCTCTTTTTTACC‑‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCCTTGGCCCTTTCGGCAAACCGTCCTGGCTAAGa                                                                                                                                                               >  1:1744624/1‑216 (MQ=255)
                                                                                          aaCGTGGATTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGGGATCATCTCCTTTTTGCC‑‑AAAAATAATTCCCGCCTGCCGCCCGCCTCAAAACATGTCGCCGGAGAAACTCCCCAGACTGACGCCCCCCCCAATCCACTCCTCCGGCCGTCGCCCTTTCGGCAAAACGAACAGGCCAAAACAGACCCCCCAGCCAcctcctta                                                                                                                                        >  2:1569492/1‑232 (MQ=255)
                                                                                            gggcgtttCCGGCCCCGCGGGGTTTTTTTTACCCGGCAGTAAGGCCCGTTTCGGGGGGGGGTTTGGCGGGGGGCCAATCTTCTTTTTTTACC‑‑AAAAAATATTCCCGGCCGGGGCCCCCCTAAACAACTTGCCGCCGGAAAACCCCCCCAGCCTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCa                                                                                                                                                <  2:1889496/223‑1 (MQ=255)
                                                                                                  ttCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑‑CACAATAATTCCTTCCTTCGGCCCGCCTAACAACTTGTCGCCGGGGAAACTCCCCAGACCGGCGCCCGCCGCAATCAACCCCTGCGGCATTGGCCCTTTTCGCTAACCGGACTGGCTAAAAACGGCCAGCGACCCACCGCCTCaatcaatc                                                                                                                                 >  2:1220990/1‑238 (MQ=255)
                                                                                                        cgccgcGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTTTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCCAACCGTACTGGCTAAGATCGACCCGCGAGCCACTGCCTAAAt                                                                                                                                      >  1:1889496/1‑227 (MQ=255)
                                                                                                           cgcgAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCCGTCATCTCTTTTTTTCC‑‑AACAATAATTCCTGCCTGCGGCCCGCCCAACAACTTGTCGCCGGAGAAACTCACCCGACTGACCCCCGCCCCAATCCACCCCTGCCGCCTTTGCTCTTTCGGCCAACCGCACTGGCTAAGAACGAACACCGAGCCACCGCCTCAATCAGTCACCACgggga                                                                                                                       >  2:2127908/1‑237 (MQ=255)
                                                                                                                    cATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGGTCATCTCTTTTTTTCC‑‑AACAATAATTTCCGCCCTGGGCCCGCCCCACAACTTGGTCCCGGAGAAAATCCACCGACCGACGCCCGCCGCAATCAAATCCTGCGGGGTTGGCTCCTTCGGCAAAACGTAATGGCTAAAATCGACCCACGAGCCCCAGCCCATAACAGTCACTACAGGGAAATCCCGc                                                                                                               >  2:205839/1‑238 (MQ=255)
                                                                                                                           tCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGGTCACCTCTTTTTTTAC‑‑CACAAAAATTTCTGCCTGCGGCCCGCCTAACAAACTGTCGCCTGAGAAAATCACCAGACTGACGCCCCCCGCAATCAAACCCTGCGGCATTGGCTCTTTTGGCAAACCGTACTGGCTTACATAGACAAACGAGACCCCGCCGAAATAAATCACCAATGGGACATACATATCTTTTcc                                                                                                       >  2:1183347/1‑239 (MQ=255)
                                                                                                                                 cgTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTTTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAgcg                                                                                                  >  1:1819086/1‑238 (MQ=255)
                                                                                                                                        cccGCTCCCGCGGGTGCTTTTGCGGGGGGCCATTCCTTTTTTTTTCACC‑‑AAAAAAATTCCCTGCCGCCGCCCCCCATAACCACTTGCCGCCGAAAAACCCCCCCAAACGACCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAgc                                                                                                   <  2:1693542/228‑1 (MQ=255)
                                                                                                                                              tCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGAc                                                                                                                                                          >  1:135892/1‑169 (MQ=255)
                                                                                                                                                aGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTTCC‑‑CACAATAATTCCTGCCTTCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACCGACGCCCGCCGCAATCAACTCCTCCGGCATTGGCTCTTTCGGCCAAACCTACTGGCTAAGATCGGCCAGCCAGCCCCCGCCCAAATCAGCTACCAAGGGAACCACCAGCCCTTTTCCCAGCGCCCCCACTTTCGcct                                                                                   >  2:1689931/1‑236 (MQ=255)
                                                                                                                                                  cGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTTCC‑‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTTTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCCGCAAACCGTACTGGCTAAGAACGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACAACCAGACCTTTGTCCAGCGCCCCCAGATcc                                                                                       >  2:1802090/1‑232 (MQ=255)
                                                                                                                                                  cGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑‑AACAATAATTCCTTCCTGCGGCCCGCCTAACAACTTTTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCCGtt                                                                                         >  1:1693542/1‑230 (MQ=255)
                                                                                                                                                       ggCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATc                                                                           >  1:1816410/1‑239 (MQ=255)
                                                                                                                                                            tGCCGGGGGGCCATCATCTTTTTTTTACC‑‑AAAAAACATCCCCCCCCGCCGCCCCCCTAAACACTGTGCGCCGGGAAAACCCCACCAGCCGACCCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCtt                                                                       <  2:505964/239‑1 (MQ=255)
                                                                                                                                                                 gggggcgcttctcttTTTTTTCCA‑‑AAAAATAACCCCCCCCCGCGGCCCCCAAAACATTTTGCCGCGGGAAAACCCCCCCACACTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTg                                                                  <  2:2159718/227‑1 (MQ=255)
                                                                                                                                                                 gcgggcgACCTTCTCTTTTTACCC‑‑AAAAATAATTCCTCCCGGCGGCCCCCCTAACAACTTGCCGCGGGAAAAACTCCCCAGACTGACCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTg                                                                  <  1:2239666/238‑1 (MQ=255)
                                                                                                                                                                 gcgggcgAACTTTTTTTTTAAACA‑‑AACAATAATCCCGGCCTGCGCCCCGCAAACCAACTTGCCGCGGAAAAACCCCCCCAGCCGGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTg                                                                  <  2:753329/238‑1 (MQ=255)
                                                                                                                                                                      cgATCATCTCTTTTTTACC‑‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCCATCAACCCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCTACCAGCGAGCCACCGCCTAAATCAGTCACTACCGGAACATCCAGCCCTTTGCCGACCGCCCCCAGTTCCGCTTCATCTATCGCCTTGGTGAAcc                                                              >  2:1081619/1‑237 (MQ=255)
                                                                                                                                                                         caatcTTTTTTTCACC‑‑AAAAAAATTCCCTGCCGGCGGCCCCCCTAACAACTAGCCGCCGGAAAACCCCACCAGACGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGa                                                         <  1:2100868/237‑1 (MQ=255)
                                                                                                                                                                           tcctctTTTTTACCCAAACAATATTCCCGCCCGGCGCCCCGCCAAAAACCTGGTCCCGGGAAAAACTCCCCAGCCGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGa                                                         <  2:1986586/237‑1 (MQ=255)
                                                                                                                                                                              ttttttttCCA‑‑AAAAATAATCCCTGCCTGCGGCCCCCCTAACAACTTGGCGCGGGAGAAACTCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCt                                                    <  1:944979/237‑1 (MQ=255)
                                                                                                                                                                                ttttttACC‑‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAAt                                                                                                                                      >  1:1547747/1‑155 (MQ=255)
                                                                                                                                                                                   tttACC‑‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGt                                               >  1:1766583/1‑239 (MQ=255)
                                                                                                                                                                                       tccaaaaaattatCCCTCCCGGCGCCCCCCCTAACATCTTGTCGCGGGAGAACCCCCCCAACCTCCCCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTa                                             <  1:205839/232‑1 (MQ=255)
                                                                                                                                                                                            acaaTATTCCCTGCGTGCGGCCCGCCAAACATCTTGGCCCCGGAAAACCCCCCCAGACGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGt                                        <  1:97170/239‑1 (MQ=255)
                                                                                                                                                                                            acaaTAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGt                                        >  1:1808670/1‑239 (MQ=255)
                                                                                                                                                                                            acaaTAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGCTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGt                                        >  1:424369/1‑239 (MQ=255)
                                                                                                                                                                                            aaaataatCCCTCCCGGCGGCCCGCCTAACACCTTGTCGCGGGAGAAACCCCCCAGGCGGGCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGt                                        <  1:1569492/237‑1 (MQ=255)
                                                                                                                                                                                              aataatTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTAt                                      >  1:1817395/1‑239 (MQ=255)
                                                                                                                                                                                               ataatCCCTGCCGGCCCCCCGCATAACAATTGGGCGCCGGAAAACCGCACAATGCGGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATg                                     <  2:1674832/239‑1 (MQ=255)
                                                                                                                                                                                                 attcccggccggcGGCCCCCCAAACAACTTGCCGCGGGAAAACCCCCCAAGGCTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                   <  1:923236/232‑1 (MQ=255)
                                                                                                                                                                                                   tccctgcctgCGCCCCCCCAAACACCTTGCCCCCGGAGAACCCCCCCACCCGAACCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                         <  2:787776/101‑1 (MQ=255)
                                                                                                                                                                                                    tCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCa                                                                                                                                                                                                                                    >  2:101784/1‑44 (MQ=255)
                                                                                                                                                                                                    tCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCa                                                                                                                                                                                                                                    <  1:101784/44‑1 (MQ=255)
                                                                                                                                                                                                    ccctgcctgCGCCCCCCCAAACAACTTGTCGCCGGAAAACCCCCCCAGCCGGACCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                         <  2:1420018/101‑1 (MQ=255)
                                                                                                                                                                                                     cctgcctgCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGt                                    >  1:1240894/1‑234 (MQ=255)
                                                                                                                                                                                                        gcctgcGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTtca                             >  2:1647666/1‑238 (MQ=255)
                                                                                                                                                                                                        gcctgcGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTtca                             >  2:1495516/1‑238 (MQ=255)
                                                                                                                                                                                                             cGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGa                                                                                                                                                               >  1:787776/1‑103 (MQ=255)
                                                                                                                                                                                                              ggCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGa                                                                                                                                                               >  1:1420018/1‑102 (MQ=255)
                                                                                                                                                                                                                     cTAACAACTTGCCGCCGGAAAACCCCCCCAGACGGCCCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                         <  2:627024/85‑1 (MQ=255)
                                                                                                                                                                                                                          ccactttgcGCCGGAAAAACTCCCCAGCCTGGCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTAtttt           <  1:1245359/231‑1 (MQ=255)
                                                                                                                                                                                                                           tcctggTCGCCGAAAAACCCCCCAAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                   <  2:1897480/209‑1 (MQ=255)
                                                                                                                                                                                                                           aacTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGAATTttc          >  1:1032628/1‑239 (MQ=255)
                                                                                                                                                                                                                                gTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGa                                                                                                                                                               >  1:627024/1‑85 (MQ=255)
                                                                                                                                                                                                                                 tCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGAGGTATTttcattca     >  1:596656/1‑238 (MQ=255)
                                                                                                                                                                                                                                     cGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCaa                         >  1:1897480/1‑214 (MQ=255)
                                                                                                                                                                                                                                     cGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCg  >  1:1276281/1‑237 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                  
GCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACC‑‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCG  >  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: