Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,658,1630CT29.2% 68.8 / 18.0 48S32F (TCT→TTT) ynfEputative selenate reductase YnfE
Reads supporting (aligned to +/- strand):  ref base C (10/24);  new base T (14/0);  total (24/24)
Fisher's exact test for biased strand distribution p-value = 8.13e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGG‑TTTTTCTTTGCCGTTTA‑CCCTGCGC‑AATGCAGCAGCAGCGG‑TACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCA‑CTTCGTCTACA‑TGTTAAAGATAATGAAGTG‑ACCTGGGTGGAAACTG‑ACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCTTGAATTACC  >  NC_000913/1657927‑1658411
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                               
aaTCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAAGCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGG‑TTTTTCTTTGCCGTTTA‑cc                                                                                                                                                                                                                                                   <  1:46611/250‑1 (MQ=255)
           attaCTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGG‑TTTTTCTTTGCCGTTTT‑CCCTGCGC‑AATgcc                                                                                                                                                                                                                                      >  1:192186/1‑250 (MQ=255)
              cctcctACCCTGCTATTGATATTGATTAATAATTCTTAACCCAAAATGGTTAGACTCCCTCTAGGGTGAGCGCGCTAAATATTGAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTTAGTTATGTCCAAAAATGAACGAATGTTGGTCATCAGCCGTAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCGGGGTGCAGTCGC‑TTGTTCTTTGCCTTGTA‑CCCTGTGC‑AATgcagca                                                                                                                                                                                                                                   <  2:653506/246‑1 (MQ=255)
                  aTACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGG‑TTTTTTTTTTCCGTTTA‑CCCTGCGC‑AAAgcagaagagc                                                                                                                                                                                                                               >  1:551720/1‑245 (MQ=255)
                    aCCCTGCTATGGATATGGATTAATAAGTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATAGATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCGCTGGCGCTGGCTGCAGGCGG‑TTTTTCTTTGCCGTTTA‑CCCTGCGC‑AATGCAGCAGCAGCg                                                                                                                                                                                                                             <  2:329060/251‑1 (MQ=255)
                      ccTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCTATAGGTTCTCTGGCGCTGGCTGCAGGCGG‑TTTTTCTTTGCCGTTTA‑CCCTGCGC‑AATGCAGCAGCAGCGG‑t                                                                                                                                                                                                                          <  1:455007/251‑1 (MQ=255)
                                               tCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGG‑TTTTTCTTTTCCGGTTA‑CCCTGCGC‑AATGCAGCAGCAGCGG‑GACAAAAGGGCCGCGCAAAAGTgggc                                                                                                                                                                                                 >  1:154412/1‑249 (MQ=255)
                                                      ccAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGG‑TTTTTCTTTGCCGTTTA‑CCCTGCGC‑AATGCAGCAGCAGCGG‑Tacaaca                                                                                                                                                                                                                    <  2:187536/225‑1 (MQ=255)
                                                          aaTGGGTAGACTCCCTCGATTGTTAACGCGCAAAATATTAAATACATAAACTGTTAAATAACGATAAAGAACACGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCGGCAGTCGG‑TTTTTCGTTGCCGTGCA‑CCCTGCGC‑ATGGCAGCAGCAGCGG‑TACAACAGGCCCGCGAAAAATTGGTCTGGGTTGCCTg                                                                                                                                                                                      <  2:305734/251‑1 (MQ=255)
                                                            tGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGG‑TTTTTCTTTTCCGGTTA‑CCCTGCGC‑CATGCAGCAGCCGCGG‑TACAACAGGCCCGCGCAAAAATGGTCTGGGGTGCCCGtt                                                                                                                                                                                    >  1:568273/1‑251 (MQ=255)
                                                             gggTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTTAGTTATGTCCAAAAATGAACGAATGGGGGGCATCAGCCGCAGAACACTCTTTAAATCTACCGCGATAGGTTCTCTGGCTCTGGCTGCAGGCGG‑TTTTTCTTTGCCGTTTA‑CCCTGCGC‑AATGCAGCAGCAGCGG‑TACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGtt                                                                                                                                                                                    <  2:623226/250‑1 (MQ=255)
                                                                 aGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGG‑TTTTTTTTTTCCCTTTT‑CCCTGCGC‑CATTCAGCCGCCACGG‑TACAACAAGCCCGCGACAACATGGTCTGGGGTGCCTTTTTCGTc                                                                                                                                                                               >  1:425173/1‑251 (MQ=255)
                                                                 aGACTCCCTCTATTGGTAGCGCGCTAAATAATCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGG‑TTTTTCGTTGCCGTTTA‑CCCTGCGC‑AATGCAGCAGCAGCGG‑TACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGt                                                                                                                                                                                <  2:430289/250‑1 (MQ=255)
                                                                  gACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGG‑TTTTTCTTTTCCGTTTAACCCT‑CGC‑AATGCAGCAGCAGCGG‑GACAAAAGGGCCCCGGAAAAATGGGCTGGGGTGGCTGTTTCGTcc                                                                                                                                                                              >  1:707736/1‑250 (MQ=255)
                                                                                  aGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGG‑TTTTTCTTTGCCGTTTA‑CCCTGCGC‑AATGCAGCAGCAGCGG‑TACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCt                                                                                                                                                               <  1:272162/250‑1 (MQ=255)
                                                                                  aGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGG‑TTTTTCTTTGCCGTTTA‑CCCTGCGC‑AATGCAGCAGCAGCGG‑TACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCt                                                                                                                                                               <  1:532405/250‑1 (MQ=255)
                                                                                     gcgcTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGGATTTTTTTTTTGCTTTTA‑CC‑TGCGC‑AATGCAGAAGCAGCGG‑GACAACAAGGCCGAGAAAAAATGGGTTGGGGGGCCTGTTcctccaacatcggtagccgccgtgc                                                                                                                                                           >  2:417371/1‑228 (MQ=255)
                                                                                        cTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGG‑TTTTTTTTTGCCGTTTT‑CCCCGCGC‑AATGCAGGAACAGCGG‑GACAAACGGCCCGCGGAAAAGGGGTCTGGGGTGCCTGTTTCGGCCACTGTGGTAGCCGCTGTGCA‑cc                                                                                                                                                       >  1:86888/1‑250 (MQ=255)
                                                                                         tAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGG‑TTTTTCTTTTCCGTTTT‑CCCCTCGC‑CATGCAGCAGCCGCGG‑GAGCACAAGCCCGCGACAAAAAGGTCTTGGGGGCCCGTTCCGTCaaacttgggaggcgccgtggactt                                                                                                                                                       >  2:268466/1‑229 (MQ=255)
                                                                                              atTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGG‑TTTTTCTTTGCCGTTTA‑CCCTGCGC‑AATGCAGCAGCAGCGG‑TACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCA‑CTTCGTc                                                                                                                                                  <  1:234221/250‑1 (MQ=255)
                                                                                                       ataAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGG‑ATTTTTTTTTCCGTTTT‑CCCTGCGC‑CATGCAGCAGCCGCGG‑TACAACAGGGCCCAGAAAAAGAGGTCTGGGGGGGCCGTTCCGGCAAATGGGGGAGCCGCTGTGGG‑CTTCGtgtcaatgttt                                                                                                                                         >  1:105703/1‑240 (MQ=255)
                                                                                                       agaaaCTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATTGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGG‑TTTTTCTTTGCCGTTTA‑CCCTGCGC‑AATGCAGCAGCAGCGG‑TACAACAGGCCCGCGAAAAAGTGGTCTTGGGTTCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCA‑CTTCGTCTACA‑TGTTaa                                                                                                                                       <  2:78093/249‑1 (MQ=255)
                                                                                                                tatgtaacgtgaaaaaacgggtatggagtTATGTCCAAAAATGCACGAAGGGTGGGCATCATCCGCAGAACACTCGTTAAATCTACCGCGCGAGGTTTTCTGGCGCTAGCGGCAGGCGG‑GTTATCTTTGCCGTTTA‑CCCTGCGC‑AATGCAGCAGCAGTGG‑TACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTTTTCCGACAACTGTGGTAGCCGCTGTGCA‑CTTCGTCTCCA‑TGTTAAAAATAATg                                                                                                                               <  2:195079/225‑1 (MQ=255)
                                                                                                                       cGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGG‑TTTTTCTTTGCCGTTTA‑CCCTGCGC‑AATGCAGCAGCAGCGG‑TACAACAGTCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAg                                                                                                                                                                    <  2:681363/208‑1 (MQ=255)
                                                                                                                         aTAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGG‑TTTTTTTTTTCCGTTTA‑ACCCGCGC‑AATGCAGCAGCAGCGG‑TACAACAGGCCCCCGAAAAAGTGGTCTTGGGTGGCTGTTCCTGCCACTGTGGTAGACGCCTTGCA‑CTTCGGCT‑CACTGTTTAAAATAATGAAGAG‑Accgg                                                                                                                    >  1:177406/1‑249 (MQ=255)
                                                                                                                          tAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGG‑TTTTTCTTTGACGTTTA‑ACCTCCGC‑CATTCACCAGCAGCGGATA‑AAAAGGCCCGCCAAAAAGGGGTCTGGGGTGCCTGTTCAGTCAACTTTGGAAGCCGCTGTGCAACTT‑GTCTACA‑TGTTAAAGATAaagaagggtcatgg                                                                                                                    >  1:567884/1‑243 (MQ=255)
                                                                                                                             aGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGG‑TTTTTTTTTTCCGGTTA‑CCCTGCGC‑CATGCAGCAGCAGCGG‑TACAACAGGCCCGCGGAAAAGGGGGCTGGGGTGCCTGTTTCGTCAACTGTGGTAGCCGCTGTGGA‑CTTCGTCTACA‑TGTTTAAGATAATGAAGTG‑ACCTggggg                                                                                                                >  1:209991/1‑249 (MQ=255)
                                                                                                                             aGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGG‑TTTTTCTTTGCCGTTTA‑CCCTGCGC‑AATGCAGCAGCAGCGG‑TACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCA‑CTTCGTCTACA‑TGTTAAAGATAATGAAGTG‑ACCTGGGTg                                                                                                                <  1:457767/251‑1 (MQ=255)
                                                                                                                                                 tCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGG‑TTTTTCTTTGCCGTTTA‑CCCTGCGC‑AATGCAGCAGCAGCGG‑TACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGtt                                                                                                                                                                                    <  1:73918/166‑1 (MQ=255)
                                                                                                                                                   cAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGG‑TTTTTCTTTGCCGTTTA‑ACCTGCGC‑AATGCAGCAGCAGCGG‑GAAAAAAGGCCCGCGGAAAAATGGGCTGGGGGGCCTGTTTCGTCCACTGTGGGAGCCGCTGTGCA‑CTTCGTCTACA‑TGTTAAAAATAAATAAGTG‑GCCTGGGGGGGAAATG‑ACAAAAACGGCaacg                                                                                         >  1:477202/1‑248 (MQ=255)
                                                                                                                                                           aaCGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGG‑TTTTTCTTTGCCGTTTA‑CCCTGCGC‑AATGCAGCAGCAGCGG‑TACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCA‑CTTCGTCTACA‑TGTTAAAGATAATGAAGTG‑ACCTGGGTGGAAACTG‑ACAATACCGGCAGCGATGAGTAc                                                                                 <  1:179048/251‑1 (MQ=255)
                                                                                                                                                                          tCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGG‑TTTTTTTTTTCCGTTTA‑CCCTTCGC‑AATGCAGCAGCAGCGG‑GACAACAAGCCCGCGGAAAAGGGGGCTGGGGGGCCTGTTTCGTCAACTGTGGTAGCCGCTGTGCA‑CTTCGTCTACA‑TGTTAAAGATAATGAAGTG‑GCCTGGGTGGAAAATG‑AAAATACCGGCAGCGGTGAGTACGGcgaacctcaggg                                                                   >  1:97609/1‑239 (MQ=255)
                                                                                                                                                                                gCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGG‑TTTTTTTTTTCCGTTTA‑CCCTGCGC‑AATGCCGCAGCAGCGG‑GACAACAGGCCCCCGAAAAAGGGGGCTGGGGGGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCA‑CTTCGTCTACA‑TGTTAAAGTTAATGAAGTG‑ACCTGGGGGGGAACTG‑GCAAAACCGGCAGCGTTGAGGACGGCCACCATCAGGGACGCGcc                                                            >  1:403926/1‑251 (MQ=255)
                                                                                                                                                                                    aaCACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGG‑TTTTTCTTTGCCGTTTA‑CCCTGCGC‑AATGCAGCAGCAGCGG‑TACAACAGGCCCGCGAAAAAGTGGTCTGGGTTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCA‑CTTCGTCTACA‑TGGTAAAGATAATGAAGTG‑ACCTGGGTGGAAACTG‑ACAATACCGGCAGCTATGAGTACGGCAACCATCAGGTACGCGCCTGtt                                                        <  2:707736/251‑1 (MQ=255)
                                                                                                                                                                                          cGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGG‑TTTTTTTTTTCCGTTTA‑CCCCGCGC‑AATGGAGCAACAGCGG‑TACAACAGGCCCGCGAAAAAAGGGTCTGGGGGGCCCGTTTCGGCAACTGTGGTAGCCGCTGGGCC‑CTTCGTCTACA‑TGTTTAAGATTATGAAGGG‑AGCTGGGGGGAAAATG‑ACAAAAACGGCAGCGGTGGGTACGGGAACCATCAGGTACGCGCCAGTTTGtgcg                                                  >  1:615615/1‑251 (MQ=255)
                                                                                                                                                                                            ttAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGG‑CTTTTTTTTGTCGTTTT‑CCCTGCGCCAGTG‑AGCAGAAGCGG‑GACAACCGGCCCGCGAAAAAGTGGTCTGGGGGGCCGGTTCCGTCAACTGTGGGAGGCGCTGTGCG‑GTTCGTCTACA‑AATTTAAGAAAATGGAGGG‑AGATGGGGGGAAACTG‑CCAATACCGGCAGCGATGAAGAAGGCGCCCATCAGGGAAGAGCCTTTTTGCGCGGt                                                >  1:260394/1‑251 (MQ=255)
                                                                                                                                                                                             tAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGG‑TTTTTTTTTGCCGGTTA‑CCCTGCGC‑AATGCAGCAGCAGCGG‑TACAACAGGCCCGCGGAAAAGGGGGCTGGGGGGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCA‑CTTCGTCTACA‑TGTTAAAGATAATGAAGTG‑GACTGGGGGGAAAATG‑ACAATACCGGCAGCGATGAGTACGGGAACCATCAGGTACGCGCCTGTTTGCGCGGTc                                               >  1:525198/1‑251 (MQ=255)
                                                                                                                                                                                              aaaTCTACCTCGATAGGTTCTCTGGCGCTGGCTGCAGGCGG‑TTTTTCTTTGCCGTTTA‑CCCTGCGC‑AATGCAGCAGCAGCGG‑TACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCA‑CTTCGTCTACA‑TGTTAAAGATAATGAAGTG‑ACCTGGGTGGAAACTG‑ACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTc                                               <  2:403926/250‑1 (MQ=255)
                                                                                                                                                                                                    aCCGCGATAGGTTCTCTGGCGCTGGCGGCAGGCGG‑TTTTTCTTTGCCGTTTT‑CCCTGCGC‑AATGTAGCAGCAGCGG‑TACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCA‑CTTCGTCTTCA‑TGTTAAAGATAATGAAGTG‑ACCTGGGTGGAAACTG‑ACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTcc                                          <  1:370838/249‑1 (MQ=255)
                                                                                                                                                                                                      cgcgATAGGTTCTCTGGCGCTGGCTGCAGGCGG‑TTTTTTTTTGCCGGTTA‑ACCTGCGC‑AATTCAGCAGCAGCGG‑GACAACAGGCCCGCGAAAAAGGGGTCTGGGGGGGCTGTTTCGTCCACTGTGGGAGCCCCTGTGGA‑CTTCGTCTACA‑TGTTAAAGATAATGAAGTGGCCCTGGG‑GGGAACTGGACAA‑ACCGGCAGCGAGGGGGACGGCAACCCTCAGGGACGCGCCTGTTTTCGCGGgcgccccctcc                                      >  1:39051/1‑244 (MQ=255)
                                                                                                                                                                                                         gATAGGTTCTCTGGCGCTGGCTGCAGGCGG‑TTTTTCTTTGCCGTTTA‑CCCTGCGC‑AATGCAGCAGCAGCGG‑TACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCA‑CTTCGTCTACA‑TGTTAAAGATAATGAAGTG‑ACCTGGGTGGAAACTG‑ACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATccg                                     <  1:344016/249‑1 (MQ=255)
                                                                                                                                                                                                              gTTCTCTGGCGCTGGCTGCAGGCGG‑TTTTTCTTTGCCGTTTA‑CCCTGCGC‑AATGCAGCAGCAGCGG‑TACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCA‑CTTCGTCTACA‑TGTTAAAGATAATGAAGTG‑ACCTGGGTGGAAACTG‑ACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGt                              <  1:641505/251‑1 (MQ=255)
                                                                                                                                                                                                               ttCTCTGGCGCTGGCTGCAGGCGG‑TTTTTCTTTGCCGTTTA‑CCCTGCGC‑AATGCAGCAGCAGCGG‑TACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCA‑CTTCGTCTACA‑TGTTAAAGATAATGAAGTG‑ACCTGGGTGGAAACTG‑ACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCt                                                           <  1:254719/221‑1 (MQ=255)
                                                                                                                                                                                                                  tctGGCGCTGGCTGCAGGCGG‑TTTTTTTTTTACGTTTAACCCT‑CGC‑AATGCAGGAGCCACCG‑TACAACAAGCCCGCGGAAAAGGGGGCTGGGGTGCCTGTTCCGTCAACTTGGGTGGCCGCCGTGCA‑CTTCCTCTTCC‑TTTTAAAGATAATGAAGTG‑ACCTGGGTGGAAACTT‑ACAATAACGGCCACGCTGAGCACGGCCACACTCAGGTAAGCGCCTGGTTTCGCGGGCGCTCCCTCCCCCGGCGtcttt                          >  1:441030/1‑247 (MQ=255)
                                                                                                                                                                                                                   ctGGCGCTGGCTGCAGGCGG‑TTTTTCTTTGCCGTTTA‑CCCTGCGC‑AATGCAGCAGCAGCGG‑TACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCA‑CTTCGTCTACA‑TGTTAAAGATAATGAAGTG‑ACCTGGGTGGAAACTG‑ACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTAtt                           <  1:320322/249‑1 (MQ=255)
                                                                                                                                                                                                                                      g‑TTTTTCTTTGCCGTTTA‑CCCTGCGC‑AATGCAGCAGCAGCGG‑TACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCA‑CTTCGTCTACA‑TGTTAAAGATAATGAAGTG‑ACCTGGGTGGAAACTG‑ACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCTTGAAt      >  1:693797/1‑251 (MQ=255)
                                                                                                                                                                                                                                          tttCTTTGCCGTTTA‑CCCTGCGC‑AATGCAGCAGCAGCGG‑TACAACTGGCCCGCGAAAACGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCA‑CTTCGTCTACA‑TGTTAAAGATAATGAAGTG‑ACCTGGGTGGAACCTG‑ACAATACCGGCAGCGATGAGTACGGCAACCTGCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTACTCATCCCGATCGCTTGAGTTAc   >  2:265305/1‑251 (MQ=255)
                                                                                                                                                                                                                                           ttCTTTGCCGTTTA‑CCCTGCGC‑AATGCAGCAGCAGCGG‑TACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCA‑CTTCGTCTACA‑TGTTAAAGATAATGAAGTG‑ACCTGGGTGGAAACTG‑ACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCTTGAATTAcc  >  1:243099/1‑251 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                               
AATCCCACCTTATTACTCATACCCTTCTATTGATATGGATTAATAATTCTTAACCCAAAATGGGTAGACTCCCTCTATTGTTAGCGCGCTAAATATTCAATATATAAACTTTTATATAACGATAAAGAACAGGGAGTGAGTTATGTCCAAAAATGAACGAATGGTGGGCATCAGCCGCAGAACACTCGTTAAATCTACCGCGATAGGTTCTCTGGCGCTGGCTGCAGGCGG‑TTTTTCTTTGCCGTTTA‑CCCTGCGC‑AATGCAGCAGCAGCGG‑TACAACAGGCCCGCGAAAAAGTGGTCTGGGGTGCCTGTTCCGTCAACTGTGGTAGCCGCTGTGCA‑CTTCGTCTACA‑TGTTAAAGATAATGAAGTG‑ACCTGGGTGGAAACTG‑ACAATACCGGCAGCGATGAGTACGGCAACCATCAGGTACGCGCCTGTTTGCGCGGTCGCTCCATCCGCCGGCGTATTAATCATCCCGATCGCTTGAATTACC  >  NC_000913/1657927‑1658411

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