Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,759,2270AT20.9% 97.6 / 10.8 43V212E (GTA→GAA) selBselenocysteinyl‑tRNA‑specific translation factor
Reads supporting (aligned to +/- strand):  ref base A (14/20);  new base T (9/0);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 1.75e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGGCGTGGTGAATATGCAGCGGCTGCCACTGGGTCAGCGGTGTATGGGTTTGAAGCTCGACAATCACCCGTGTGAACGGCTCTGGCGGCACATCGGCAAGCAGCCAGTCGCCACGGTTAATCTGCTCTTTTTCCGCATCACCCGCGATGTTAAGCGCGATACGCTGCCCGGCATTGGCGGTTTCTGTTGGCTGGTTTTGCGCATGCAGCGCACGTACACGCATCGGTTTATTTACACCAGTCAGCCAGAGTGAATCGCCTACCTTCAC‑TTCCCCGCTTAACGCCGTACCGGTGACGACC‑A‑GCCCGGCACCTTTTACGGTAAATGCGCGGTCAATCGCGAGGCGGAAACTATGTTGGCTGGCGTGCTCGCGTTCCGGCAACTGAAGCAGATGCTCGCGCAGGGCATCCATTCCCCGACCTTCGGTTGCTGCGGTGATAAACAGTTTTGCCTCAGCAAAACCGTAT  >  NC_000913/3758994‑3759457
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                         
cggcgTGGTGAATATGCAGCGGCTGCCACTGGGTCAGCGGTGTATGGGTTTGAAGCTCGACAATCACCCGTGTGAACGGCTCTGGCGGCACATCGGCAAGCAGCCAGTCGCCACGGTTAATCTGCTCTTTTTCCGCATCACCCGCGATGTTAAGCGCGATACGCTGCCCGGCATTGGCGGTTTCTGTTGGCTGGTTTTGCGCATGCAGCGCACGTACACGCATCGGTTTATTTac                                                                                                                                                                                                                                          <  1:49649/235‑1 (MQ=255)
 ggcgTGGTGAATATGCAGCGGCTGCCACTGGGTCAGCGGTGTATGGGTTTGAAGCTCGACAATCACCCGTGTGAACGGCTCTGGCGGCACATCGGCAAGCAGCCAGTCGCCACGGTTAATCTGCTCTTTTTCCGCATCACCCGCGATGTTAAGCGCGATACGCTGCCCGGCATTGGCGGTTTCTGTTGGCTGGTTTTGCGCATGCAGCGCACGTACACGCATCGGTTTATTTacac                                                                                                                                                                                                                                        <  1:18375/236‑1 (MQ=255)
          aaTATGCAGCGGCTGCCACTGGGTCAGCGGTGTATGGGTTTGAAGCTCGACAATCCCCCGTGTGAACGGCTCTGGCGGCACATCGGCAAGCAGCCAGTCGCCACGGTTAATCTGCTCTTTTTCCGCATCACCCGCGATGTTAAGCGCGATACGCTGCCCGGCATTGGCGGTTTCTGTTGGCTGGTTTTGCGCATGCAGCGCACGTACACGCATCGGTTTATTTACACCa                                                                                                                                                                                                                                      <  1:371449/229‑1 (MQ=255)
                 aGCGGCTGCCACTGGGTCAGCGGTGTATGGGTTTGAAGCTCGACAATCACCCGTGTGAACGGCTCTGGCGGCACATCGGCAAGCAGCCAGTCGCCACGGTTAATCTGCTCTTTTTCCGCATCACCCGCGATGTTAAGCGCGATACGCTGCCCGGCATTGGCGGTTTTTGTTGGTTGGGTTTTGGCATGCAGCGCACGTACACGCATCGGTTTATTTTCACCAGTCAGCCACAGTGa                                                                                                                                                                                                                        >  1:505251/1‑236 (MQ=255)
                 aGCGGCTGCCACTGGGTCAGCGGTGTATGGGTTTGAAGCTCGACAATCACCCGTGTGAACGGCTCTGGCGGCACATCGGCAAGCAGCCAGTCGCCACGGTTAATCTGCTCTTTTTCCGCATCACCCGCGATGTTAAGCGCGATACGCTGCCCGGCATTGGCGGTTTCTGCTGGCTGGTTTTGTGCAGTCCGCGCCCGTGCACGCGTCGGGTTATTTTTACCAGTCAGccccagac                                                                                                                                                                                                                         >  2:371449/1‑229 (MQ=255)
                  gCGGCTGCCACTGGGTCAGCGGTTTATGGGTTTGAAGCTCGACAATCACCCGTGTGAACGGCTCTGGCGGCACATCGGCAAGCAGCCAGTCGCCACGGTTAATCTGCTCTTTTTCCGCATCACCCGCGATGTTAAGCGCGATACGCTGCCCGGCATTGGCGGTTTCTGTTGGCTGGTTTTGCGCATGCAGCGCACGTACACGCATCGGTTTATTTACACCAGTCAGCCAGAGTGaa                                                                                                                                                                                                                       <  1:91246/236‑1 (MQ=255)
                  gCGGCTGCCACTGGGTCAGCGGTGTATGGGTTTGAAGCTCGACAATCACCCGTGTGAACGGCTCTGGCGGCACATCGGCAAGCAGCCAGTCGCCACGGTTAATCTGCTCTTTTTCCGCATCACCCGCGATGTTAAGCGCGATACGCTGCCCGGCATTGGCGGTTTCTTTTTGCTGGTTTTGCGCCTGCACGCCACGTACACGCATCCGGTTATTTTCACCAGTCAGCCAGAGTGa                                                                                                                                                                                                                        >  1:12441/1‑235 (MQ=255)
                  gCGGCTGCCACTGGGTCAGCGGTGTATGGGTTTGAAGCTCGACAATCACCCGTGTGAACGGCTCTGGCGGCACATCGGCAAGCAGCCAGTCGCCACGGTTAATCTGCTCTTTTTCCGCATCACCCGCGATGTTAAGCGCGATACGCTGCCCGGCATTGGCGGTTTCTGTTGGCTGGTTTTTCGCATGCAGCGCACGTACACGCATCGGGTTATTTTCACCAGTCAGTCAGAGTGa                                                                                                                                                                                                                        >  1:62366/1‑235 (MQ=255)
                       tGCCACTGGGTCCGGGGTGTATGGGTTTGAAGCTCGACAATCCCCCGTGTGAACGGCTCTGGCGGCACATCGGCAAGCAGCCAGTCGCCACGGTTAATCTGCTCTTTTTCCGCATCACCCGCGATGTTAAGCGCGATACGCTGCCCGGCATTGGCGGTTTCTGTTGGCTGGTTTTGCGCATGCAGCGCACGTACACGCATCGGTTTATTTACACCAGTCAGCCAGAGTGAATCGcc                                                                                                                                                                                                                  <  1:231144/236‑1 (MQ=255)
                       tGCCACTGGGTCAGCGGTGTATGGGTTTGAAGCTCGACAATCACCCGTGTGAACGGCTCTGGCGGCACATCGGCAAGCAGCCAGTCGCCACGGTTAATCTGCTCTTTTTCCGCATCACCCGCGATGTTAAGCGCGATACGCTGCCCGGCATTGGCGGTTTCTGTTGGCTGGTTTTGCGCATGCAGCGCACGTACACGCATCGGTTTATTTACACCAGTCAGCCAGAGTGAATCGcc                                                                                                                                                                                                                  <  1:458428/236‑1 (MQ=255)
                         ccACTGGGTCAGCGGTTGTGGGGTTTGAAGCCCGACAATCCCCCGTGTGAACGGCTCTGGCGGCCAATCGGCAAGCAGCCAGTCGCCACGGTTAATCTGCTCTTTTTCCGCATCACCCGCGATGTTAAGCGCGATACGCTGCCCGGCATTGGCGGTTTCTGTTGGCTGGTTTTGCGCATGCAGCGCACGTACACGCATCGGTTTATTTACACCAGTCAGCCATAGTGAATCGCCt                                                                                                                                                                                                                 <  2:464950/235‑1 (MQ=255)
                                                  tGAAGCTCGACAATCACCCGTGTGAACGGCTCTGGCGGCACATCGGCAAGCAGCCAGTCGCCACGGTTAATCTGCTCTTTTTCCGCATCACCCGCGATGTTAAGCGCGATACGCTGCCCGGCATTGGCGGTTTCTGTTGGCTGGGTTTGCGCATGCAACGCACGTACACGCATCGGGTTATTTACACCAGTCAGCCAGAGTGAAACGCCTACCTTCAC‑TTCCCCCCTTAAccacgt                                                                                                                                                                                      >  1:487795/1‑231 (MQ=255)
                                                    aaGCTCGACAATCACCCGTGTGAACGGCTCTGGCGGCACATCGGCAAGCAGCCAGTCGCCACGGTTAATCTGCTCTTTTTCCGCATCACCCGCGATGTTAAGCGCGATACGCTGCCCGGCATTGGCGGTTTTTTTTGGGTGGGTTTTTGCACGCAGGCCACCTACCCCCACTGGGTTTTTTTTACCCGTCAGCCCCAGGGAAACGCCCACCTCTAC‑TTTCCCGCTTGacgacggaa                                                                                                                                                                                    >  2:311265/1‑233 (MQ=255)
                                                     aGCTCGACAATCACCCGTGTGAACGGCTCTGGCGGCACATCGGCAAGCAGCCAGTAGCCACGGTTAATCTGCTCTTTTTCCGCATCACCCGCGATGTTAAGCGCGATACGCTGCCCGGCATTGGCGGTTTCTGTTGGCTGGTTTTGCGCATGCAGCGCACGTACACGCATCGGTTTATTTACACCAGTCAGCCAGAGTGAATCGCCTACCTTCAC‑TTCCCCGCTTAACGCCGTAc                                                                                                                                                                                    <  1:362723/235‑1 (MQ=255)
                                                          gACAATCCCCCGTGTGAACGGCTCTGGCGGCACATCGGCAAGCAGCCAGTCGCCACGGTTAATCTGCTCTTTTTCCGCATCACCCGCGATGTTAAGCGCGATACGCTGCCCGGCATTGGCGGTTTCTGTTGGCTGGTTTTGCGCATGCAGCGCACGTACACGCATCGGTTTATTTACACCAGTCAGCCAGAGTGAATCGCCTACCTTCAC‑TTCCCCGCTTaa                                                                                                                                                                                            <  1:762592/222‑1 (MQ=255)
                                                                 aCCCGTGTGAACGGCTCTGGCGGCACATCGGCAAGCAGCCAGTCGCCACGGTTAATCTGCTCTTTTTCCGCATCACCCGCGATGTTAAGCGCGATACGCTGCCCGGCATTGGCGGTTTCTGTTTGCTGGGTTTTCGCATGCAGCGCACGTACACCCATCGGGTTATTTTACCCAATCAGCCCGAGGGAAAGTCCTCCCTTCCA‑TTTCCCGCCTAACAcccgac                                                                                                                                                                                    >  2:762592/1‑219 (MQ=255)
                                                                 aCCCGTGTGAACGGCTCTGGCGGCACATCGGCAAGCAGCCAGTCGCCACGGTTAATCTGCTCTTTTTCCGCATCACCCGCGATGTTAAGCGCGATACGCTGCCCGGCATTGGCGGTTTCTGTTGGCTGGTTTTGCGCATGCAGCGCACGTACACGCATCGGTTTATTTACACCAGTCAGCCAGAGTGAATCGCCTACCTTCAC‑TTcc                                                                                                                                                                                                    <  1:365081/207‑1 (MQ=255)
                                                                      tgtgAACGGCTCTGGCGGCACATCGGCAAGCAGCCAGTCGCCACGGTTAATCTGCTCTTTTTCCGCATCACCCGCGATGTTAAGCGCGATACGCTGCCCGGCATTGGCGGTTTTTGTTGGCGGGTTTTTCGCATGCCGCGCACGGACACGCCTCGGGTTTTTTTTACCAGGCAGCCCCAGTGAATCGCCTACATTCAC‑TTCCCCGCTTCATGCCCTACCAGTGCCGTAC‑A‑GCCCg                                                                                                                                                                 >  2:757115/1‑235 (MQ=255)
                                                                      tgtgAACGGCTCTGGCGGCACATCGGCAAGCAGCCAGTCGCCACGGTTAATCTGCTCTTTTTCCGCATCACCCGCGATGTTAAGCGCGATACGCTGCCCGGCATTGGCGGTTTCTGTTGGCTGGGTTTTCGCATGCAGCGCACGTACACGCCTCGGGTTTTTTTCAcc                                                                                                                                                                                                                                       >  2:278710/1‑168 (MQ=255)
                                                                             ggCTCTGGCGGCACATCGGCAAGCAGCCAGTCGCCACGGTTAATCTGCTCTTTTTCCGCATCACCCGCGATGTTAAGCGCGATACGCTGCCCGGCATTGGCGGTTTCTGTTGGCTGGGTTTGCGCATGGAGCGCCCCGACACCCCTCGGGTTATTTTAACCAGTCAACCCGAGTGAAACGCCCCACCTCCC‑TTTCCCGCCTCAAGCCCGAACGGCGGAGACC‑A‑GCCCCGcggccgt                                                                                                                                                         >  2:356325/1‑230 (MQ=255)
                                                                               ctctGGCGGCACATCGGCAAGCAGCCAGTCGCCACGGTTAATCTGCTCTTTTTCCGCATCACCCGCGATGTTAAGCGCGATACGCTGCCCGGCATTGGCGGTTTTTTTTGGCTGGTTTTGCGGCTGCAGCGCACCGAACCCCATCGGTTTTTTTAAACCAATCAGCCCGAGTGAAACGCCCACCTTCACTTTCCCCCCTTAA‑GCCGTAACAGTGACGACC‑ACGCCCgcgccc                                                                                                                                                            >  2:629379/1‑227 (MQ=255)
                                                                                         acaTCGGCAAGCAGCCAGTCGCCACGGTTAATCTGCTCTTTTTCCGCATCACCCGCGATGTTAAGCGCGATACGCTGCCCGGCATTGGCGGTTTCTGTTTGGTGGGTTTTCGCATGCAGCCCCCGTACACGCCTCGGGTTTTTTATACCAGTCCACCAGAGGGAATCGCCCACCCTCAC‑TTTCCCGCCTTACACCCTTCCGGTGACGACCAA‑CCCCGCCACTTTTTACgtgaaaac                                                                                                                                              >  2:585224/1‑229 (MQ=255)
                                                                                                         aggcGCCACGGTTAATCTGCTCTTTTTCCGCATCACCCGCGATGTTAAGCGCGATACGCTGCCCGGCATTGGCGGTTTCTGTTGGCTGGTTTTGCGCATGCAGCGCACGTACACGCATCGGTTTATTTACACCAGTCAGCCAGAGTGAATCGCCTACCTTCAC‑TTCCCCGCTTAACGCCGTACCGGTGACGACC‑A‑GCCCGGCACCTTTTACGGTAAATGCGCGGTCAATCGCGAgg                                                                                                                             <  1:585224/233‑1 (MQ=255)
                                                                                                               cACGGTTAATCTGCTCTTTTTCCGCATCACCCGCGATGTTAAGCGCGATACGCTGCCCGGCATTGGCGGTTTCTGTTGGCTGGTTTTGCGCATGCAGCGCACGTACACGCATCGGTTTATTTACACCAGTCAGCCAGAGTGAATCGCCTACCTTCAC‑TTCCCCGCTTAACGCCGTACCGGTGACGACC‑A‑GCCCGGCACCTTTTACGGTAAATGCGCGGTCAATCGCGAGGc                                                                                                                            <  1:484379/231‑1 (MQ=255)
                                                                                                               cACGGTTAATCTGCTCTTTTGCCGCATCACCCGCGATGTTAAGCGCGATACGCTGCCCGGCATTGGCGGTTTCTGTTGGCTGGTTTTGCGCATGCAGCGCACGTACACGCATCGGTTTATTTACACCAGTCAGCCAGAGTGAATCGCCTACCTTCAC‑TTCCCCGCTTAACGCCGTACCGGTGACGACC‑A‑GCCCGGCACCTTTTACGGTAAATGCGCGGTCAATCGCGAGGCGGaaa                                                                                                                       <  1:178390/236‑1 (MQ=255)
                                                                                                                        tCTGCTCTTTTTCCGCATCACCCGCGATGTTAAGCGCGATACGCTGCCCGGCATTGGCGGTTTCTGTTGGCTGGTTTTTCGCCTGCAGCGCCCGGACACCCCTCGGGTTTTTTACCCCAGTCCGCCAGAGTGGAACGCCCTCCTTCAC‑CTTCCCGCCTAACGCCGTACCCGTTGCGACC‑A‑GCCCGGCACCACTTTCGGTAAATTCCCGGGCAATCGCGCAGCGGCAACTAtagttg                                                                                                              >  2:156989/1‑231 (MQ=255)
                                                                                                                                            tctgagattatAAGCGCGATGCGCTGCCCGGCATTGGCGGTTTCTGTTGGCTGGTTTTGCGCATGCAGCGCACGTACACGCATCGGTTTATTTACACCAGTCAGCCAGAGTGAATCGCCTACCTTCAC‑TTCCCCGCTTAACGCCGTACCGGTGACGACC‑A‑GCCCGGCACCTTTTACGGTAAATGCGCGGTCAATCGCGAGGCGGAAACTATGTTGGCTGGCGTGCTCGCGTTCCgg                                                                                          <  2:669574/226‑1 (MQ=255)
                                                                                                                                                      tAAGCGCGATACGCTGCCCGGCATTGGCGGTTTCTGTTGGCTGGTTTTGCGCATGCAGCGCACGTACACGCATCGGTTTATTTACACCAGTCAGCCAGAGTGAATCGCCTACCTTCAC‑TTCCCCGCTTAACGCCGTACCGGTGACGACC‑A‑GCCCGGCACCTTTTACGGTAAATGCGCGGTCAATCGCGAGGCGGAAACTATGTTGGCTGGCGTGCTCGCGTTCCGGCAACTGAAg                                                                                 <  1:14305/235‑1 (MQ=255)
                                                                                                                                                               tACGCTGCCCGGCATTGGCGGTTTCTGTTGGCTGGTTTTTCGCATGCAGCGCACGGACAAGCCTCGGTTTTTTTACACCAGTCAGCCAGAGTGAATCGCCTACCTTCAC‑TTCCCCGCTTAACGCCGTACCGGGGACGACC‑A‑ACCCCGCACCCTTTTCGGTCAATGCGCGGGCAATCGCGAGGCGGGAACTATGTTTGGTTGGGGGCTTGCGGTCCGGGGACCGcagcaggagcgtg                                                                       >  2:533728/1‑229 (MQ=255)
                                                                                                                                                                      cccGGCATTGGCGGTTTCTGTTGGCTGGTTTTGCGCATGCAGCGCACGTACACGCATCGGTTTATTTACACCAGTCAGCCAGAGTGAATCGCCTACCTTCAC‑TTCCCCGCTTAACGCCGTACCGGTGACGACC‑A‑GCCCGGCACCTTTTACGGTAAATGCGCGGTCAATCGCGAGGCGGAAACTATGTTGGCTGGCGTGCTCGCGTTCCGGCAACTGAAGCAGATGCTCGCGCAggg                                                                <  1:113964/236‑1 (MQ=255)
                                                                                                                                                                       ccGGCATTGGCGGTTTCGGTTGGCTGGTTTTGCGCATGCAGCGCACGTACACGCATCGGTTTATTTACACCAGTCAGCCAGAGTGAATCGCCTACCTTCAC‑TTCCCCGCTTAACGCCGTACCGGTGACGACC‑A‑GCCCGGCACCTTTTACGGTAAATGCGCGGTCAATCGCGAGGCGGAAACTATGTTGGCTGGCGTGCTCGCGTTCCGGCAACTGAAGCAGATGCTCGCGCAggg                                                                <  1:629379/235‑1 (MQ=255)
                                                                                                                                                                          gCATTGGCGGTTTCTGTTGGCTGGTTTTGCGCATGCAGCGCACGTACACGCATCGGTTTATTTACACCAGTCAGCCAGAGTGAATCGCCTACCTTCAC‑TTCCCCGCTTAACGCCGTACCGGTGACGACC‑A‑GCCCGGCACCTTTTACGGTAAATGCGCGGTCAATCGCGAGGCGGAAACTATGTTGGCTGGCGTGCTCGCGTTCCGGCAACTGAAGCAGATGCTCGCGCAGGGCATc                                                            <  1:533728/236‑1 (MQ=255)
                                                                                                                                                                          gCATTGGCGGTTTCTGGTGGCTGGTTTTGCGCATGCAGCGCACGTACACGCATCGGTTTATTTACACCAGTCAGCCAGAGTGAATCGCCTACCTTCAC‑TTCCCCGCTTAACGCCGTACCGGTGACGACC‑A‑GCCCGGCACCTTTTACGGTAAATGCGCGGTCAATCGCGAGGCGGAAACTATGTTGGCTGGCGTGCTCGCGTTCCGGCAACTGAAGCAGATGCTCGCGCAGGGCATc                                                            <  1:504912/236‑1 (MQ=255)
                                                                                                                                                                                              cTGGTTTTGCGCATGCAGCGCACGTACACGCATCGGTTTATTTACACCAGTCAGCCAGAGTGAATCGCCTACCTTCAC‑TTCCCCGCTTAACGCCGTACCGGTGACGACC‑A‑GCCCGGCACCTTTTACGGTAAATGCGCGGTCAATCGCGAGGCGGAAACTATGTTGGCTGGCGTGCTCGCGTTCCGGCAACTGAAGCAGATGCTCGCGCAGGGCATCCA‑TCCCCGACCTTCGGTtgc                                       >  2:494855/1‑236 (MQ=255)
                                                                                                                                                                                              cTGGTTTTGCGCATGCAGCGCACGTACACGCATCGGTTTATTTACACCAGTCAGCCAGAGTGAATCGCCTACCTTCAC‑TTCCCCGCTTAACGCCGTACCGGTGACGACC‑A‑GCCCGGCACCTTTTACGGTAAATGCGCGGTCAATCGCGAGGCGGAAACTATGTTGGCTGGCGTGCTCGCGTTCCGGCAACTGAAGCAGATGCTCGCGCACGGCATCCATTCCCCGACCTTCGGCtg                                        >  1:562105/1‑236 (MQ=255)
                                                                                                                                                                                                         cATGCAGCGCACGTACACGCATCGGTTTATTTACACCAGTCAGCCAGAGTGAATCGCCTACCTTCAC‑TTCCCCGCTTAACGCCGTACCGGTGACGACC‑A‑GCCCGGCACCTTTTACGGTAAATGCGCGGTCAATCGCGAGGCGGAAACTATGTTGGCTGGCGTGCTCGCGTTCCGGCAACTGAAGCAGATGCTCGCGCAGGGCATCCA‑TCCCCGACCTTCGGTTGCTGCGGTGATaa                            <  1:494855/236‑1 (MQ=255)
                                                                                                                                                                                                                cgcACGTACACGCATCGGTTTATTTACACCAGTCAGCCAGAGTGAATCGCCTACCTTCAC‑TTCCCCGCTTAACGCCGTACCGGTGACGACC‑A‑GCCCGGCACCTTTTACGGTAAATGCGCgg                                                                                                                                         <  2:169980/121‑1 (MQ=255)
                                                                                                                                                                                                                cgcACGTACACGCATCGGTTTATTTACACCAGTCAGCCAGAGTGAATCGCCTACCTTCAC‑TTCCCCGCTTAACGCCGTACCGGTGACGACC‑A‑GCCCGGCACCTTTTACGGTAAATGCGCgg                                                                                                                                         >  1:169980/1‑121 (MQ=255)
                                                                                                                                                                                                                    cGTACACGCATCGGTTTATTTACACCAGTCAGCCAGAGTGAATCGCCTACCTTCAC‑TTCCCCGCTTAACGCCGTACCGGTGACGACC‑A‑GCCCGGCACCTTTTACGGTAAATGCGCGGTCAATCGCGAGGCGGAAACTATGTTGGCTGGCGTGCTCGCGTTCCGGCAACTGAAGCAGATGCTCGCGCAGGGCATCCATTCCCCGACCTTCGGTTGCTGCGGTGATAAACAGtt                      <  1:156989/232‑1 (MQ=255)
                                                                                                                                                                                                                            cATCGGTTTATTTACACCAGTCAGCCAGAGTGAATCGCCTACCTTCAC‑TTCCCCGCTTAACGCCGTACCGGTGACGACC‑A‑GCCCGGCACCTTTTACGGTAAATGCGCGGTCAATCGCGAGGCGGAAACTATGTTGGCTGGCGTGCTCGCGTTCCGGCAACTGAAGCAGATGCTCGCGCAGGGCATCCATTCCCCGACCTTCGGTTGCTGCGGTGATAAACAGTTTTGCCTCAGCa           >  2:488842/1‑235 (MQ=255)
                                                                                                                                                                                                                             aTCGGTTTATTTACACCAGTCAGCCAGAGTGAATCGCCTACCTTCAC‑TTCCCCGCTTAACGCCGTACCGGTGACGACC‑A‑GCCCGGCACCTTTTACGGTAAATGCGCGGTCAATCGCGAGGCGGAAACTATGTTGGCTGGCGTGCTCGCGTTCCGGCAACTGAAGCAGATGCTCGCGCAGGGCATCCATTCCCCGACCTTCGGTTGCTGCGGTGATAAACAGTTTTGCCTCAGCcaa         >  2:729769/1‑236 (MQ=255)
                                                                                                                                                                                                                             aTCGGTTTATTTACACCAGTCAGCCAGAGTGAATCGCCTACCTTCAC‑TTCCCCGCTTAACGCCGTACCGGTGACGACC‑A‑GCCCGGCACCTTTTACGGTAAATGCGCGGTCAATCGCGAGGCGGAAACTATGTTGGCTGGCGTGCTCGCGTTCCGGCAACTGAAGCAGATGCTCGCGCAGGGCATCCATTCCCCGACCTTCGGTTGCTGCGGTGATAAACAGTTTTGCCTCAGCaca         >  1:634890/1‑234 (MQ=255)
                                                                                                                                                                                                                               cGGTTTATTTACACCAGTCAGCCAGAGTGAATCGCCTACCTTCAC‑TTCCCCGCTTAACGCCGTACCGGTGACGACC‑A‑GCCCGGCACCTTTTACGGTAAATGCGCGGTCAATCGCGAGGCGGAAACTATGTTGGCTGGCGTGCTCGCGTTCCGGCAACTGAAGCAGATGCTCGCGCAGGGCATCCATTCCCCGACCTTCGGTTGCTGCGGTGATAAACAGTTTTGCCTCAGCTAAAc       >  1:146293/1‑236 (MQ=255)
                                                                                                                                                                                                                               cGGTTTATTTACACCAGTCAGCCAGAGTGAATCGCCTACCTTCAC‑TTCCCCGCTTAACGCCGTACCGGTGACGACC‑A‑GCCCGGCACCTTTTACGGTAAATGCGCGGTCAATCGCGAGGCGGAAACTATGTTGGCTGGCGTGCTCGCGTTCCGGCAACTGAAGCAGATGCTCGCGCAGGGCATCCATTCCCCGACCTTCGGTTGCTGCGGTGATAAACAGTTTTGCCTCAGCAAAAc       >  2:26928/1‑236 (MQ=255)
                                                                                                                                                                                                                               cGGTTTATTTACACCAGTCAGCCAGAGTGAATCGCCTACCTTCAC‑TTCCCCGCTTAACGCCGTACCGGTGACGACC‑A‑GCCCGGCACCTTTTACGGTAAATGCGCGGTCAATCGCGAGGCGGAAACTATGTTGGCTGGCGTGCTCGCGTTCCGGCAACTGAAGCAGATGCTCGCGCAGGGCATCCATTCCCCGACCTTCGGTTGCTGCGGTGATAAACAGTTTTGCCTCAGCAAAAc       >  1:330204/1‑236 (MQ=255)
                                                                                                                                                                                                                                     atttaCACCAGTCAGCCAGAGTGAATCGCCTACCTTCAC‑TTCCCCGCTTAACGCCGTACCGGTGACGACC‑A‑GCCCGGCACCTTTTACGGTAAATGCGCGGTCAATCGCGAGGCGGAAACTATGTTGGCTGGCGTGCTCGCGTTCCGGCAACTGAAGCAGATGCTCGCGCAGGGCATCCATTCCCCGACCTTCGGTTGCTGCGGTGATAAACAGTTTTGCCTCAGCAAAACCGAAt  >  1:411958/1‑235 (MQ=255)
                                                                                                                                                                                                                                       ttaCACCAGTCAGCCAGAGTGAATCGCCTACCTTCAC‑TTCCCCGCTTAACGCCGTACCGGTGACGACC‑A‑GCCCGGCACCTTTTACGGTAAATGCGCGGTCAATCGCGAGGCGGAAACTATGTTGGCTGGCGTGCTCGCGTTCCGGCAACTGAAGCAGATGCTCGCGCAGGGCATCCATTCCCCGACCTTCGGTtgctg                                     <  2:550071/198‑1 (MQ=255)
                                                                                                                                                                                                                                       ttaCACCAGTCAGCCAGAGTGAATCGCCTACCTTCAC‑TTCCCCGCTTAACGCCGTACCGGTGACGACC‑A‑GCCCGGCACCTTTTACGGTAAATGCGCGGTCAATCGCGAGGCGGAAACTATGTTGGCTGGCGTGCTCGCGTTCCGGCAACTGAAGCAGATGCTCGCGCAGGGCATCCATTCCCCGACCTTCGGTtgctg                                     >  1:550071/1‑198 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                         
CGGCGTGGTGAATATGCAGCGGCTGCCACTGGGTCAGCGGTGTATGGGTTTGAAGCTCGACAATCACCCGTGTGAACGGCTCTGGCGGCACATCGGCAAGCAGCCAGTCGCCACGGTTAATCTGCTCTTTTTCCGCATCACCCGCGATGTTAAGCGCGATACGCTGCCCGGCATTGGCGGTTTCTGTTGGCTGGTTTTGCGCATGCAGCGCACGTACACGCATCGGTTTATTTACACCAGTCAGCCAGAGTGAATCGCCTACCTTCAC‑TTCCCCGCTTAACGCCGTACCGGTGACGACC‑A‑GCCCGGCACCTTTTACGGTAAATGCGCGGTCAATCGCGAGGCGGAAACTATGTTGGCTGGCGTGCTCGCGTTCCGGCAACTGAAGCAGATGCTCGCGCAGGGCATCCATTCCCCGACCTTCGGTTGCTGCGGTGATAAACAGTTTTGCCTCAGCAAAACCGTAT  >  NC_000913/3758994‑3759457

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