Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091317,0000CT33.3% 48.7 / 17.5 42noncoding (49/59 nt)sokCsmall regulatory RNA antitoxin SokC
Reads supporting (aligned to +/- strand):  ref base C (20/8);  new base T (0/14);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 6.78e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.83e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGGATTCGTAAGCCGTGAAAACAGCAACCT‑CCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTC‑TTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGT‑GTTCAGCATATAGGAGGCCTCGGGTTG‑ATGGTAAAATATCACTCGGGGCTTTTC‑TCTATCTGCCGTTCAGCTAAT‑GCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGGCCTAACAAG  >  NC_000913/16756‑17228
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                      
cgaatgcttAAGCCGTAAAAACAGCAACCC‑CCTTTTGGCCAGTTCGGAGGTGAACCTCACAGGGGTC‑TTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGT‑GTTCAGCATATAGGAGGCCTCGGGTTG‑ATGGTAAAATATCACTCGGGGCtttt                                                                                                                                                                                                                                    <  2:60997/241‑1 (MQ=255)
cGGATTCGTAAGCCGTGAAAACAGCAACCT‑CCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTC‑TTTTCTCGTTCCCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGT‑GTTCAGCATATAGGAGGCCTCGGGTTG‑ATGGTAAAATATCACTCGGGGCTTTtc‑t                                                                                                                                                                                                                                 <  1:429761/251‑1 (MQ=255)
 ggttccGTAAGCCGAGAAAACAGCA‑CCTGCCGTCTGGCCATTTGGGATGTGAACCTCACAGAGGTC‑TTTTCTCTTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGT‑GTTCAGCATATAGGAGGCCTCGGGTTG‑ATGGTAAAATATCACTCGGGGCTTTtc‑tc                                                                                                                                                                                                                                <  2:545638/246‑1 (MQ=255)
 ggATTCGTAAGCCGTGAAAACAGCAACCT‑CCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTC‑TTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTTCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGT‑GTTCAGCATATAGGAGGCCTCGGGGTG‑ATGGTAAAATATCACTCGGGGCTTTtc‑tt                                                                                                                                                                                                                                >  2:544382/1‑250 (MQ=255)
         agaccGTGAAAACAGCAACCT‑CGTTCTGCCCGGTTGGGATTTGAACCCCTCAG‑TGTCATTTTTTCGTTCCCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGATGCACTCTACATGT‑GTGCAGCATATAGTAGGCCTCGGGTTG‑ATGGTAAAATATCACTCGGGGCTTGTC‑TCTATCTGcc                                                                                                                                                                                                                        <  2:501320/248‑1 (MQ=255)
         aaGCCGTGAAAACAGCAACCT‑CCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTC‑TTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATTACGCTCAAGGCGAAAATCCACGCCTTATGCTGCTTCAATGCCCTCTTCTCCTTGACCTTACCGTCAGTCAGAAGCACTCTACCTGT‑GTTCAGCATATAGGAGGCCTCCGGTTG‑ATTGTAAAATATCACTCTGGACTTTTc‑tctctcttcc                                                                                                                                                                                                                        >  2:130827/1‑248 (MQ=255)
                    aCAGCAACCT‑CCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTC‑TTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGT‑GTTCAGCATATAGGAGGCCTCGGGGTG‑ATGGTAAAATATCACTCGGGGCTTTTC‑TCTATCTGCCGTTCAGCTAAt                                                                                                                                                                                                             >  1:470617/1‑251 (MQ=255)
                     ccgCAACCC‑CCGTCGGGCCAGTTGGGGTGTGACCCTAACAGAGGTC‑TTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGT‑GTTCAGCATATAGGAGGCCTCGGGTTG‑ATGGTAAAATATCACTCGGGGCTTTTC‑TCTATCTGCCGTTCAGCTAAT‑g                                                                                                                                                                                                           <  2:325166/249‑1 (MQ=255)
                             ccccGTCGGGCCAGTCCGGATGTGACCCTCACAGAGGTC‑TTGTCTCGTTCCCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGT‑GTTCAGCATATAGGAGGCCTCGGGTTG‑ATGGTAAAATATCACTCGGGGCTTTTC‑TCTATCTGCCGTTCAGCTAAT‑GCCTgaga                                                                                                                                                                                                    <  2:505551/249‑1 (MQ=255)
                                    tGGCCAGTTCGGATGTGAACCTCACAGAGGTC‑TTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGTCGATCAGGGCGACAATCATCGCCTTATACTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGT‑GTTCAGCATATAGGAGGCCCCGGGTTG‑ATGGTAAAATATCACTCGGGGCTTTTC‑TCTATCTTCCGTTCCGCTCAT‑GCCTGagacagac                                                                                                                                                                                               >  2:466064/1‑249 (MQ=255)
                                          gTTCGGATGTGAACCTCACAGAGGTC‑TTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGCTAAAGATACCGAACAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCCTTCGGTCAGTAAAACGCACTCTACCTGT‑GTTCCGACTATAGGGGGCCTCGGGTTG‑ATGGGAAAATATCACTCGGGGCTTTTC‑TCTTTCTGCCGTTTAGCTAAT‑TCCTGAGACAGACAGCCTcc                                                                                                                                                                                        >  2:562752/1‑249 (MQ=255)
                                                aTGTGAACCTCACAGAGGTC‑TTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGCCAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGT‑GTTCAGCATATAGGAGGCCTCGGGTTG‑ATGGGAAAATATCACTCGGGGCTTTTC‑TCTATCTGCCGTTCAGCTAAT‑GCCCGGGGCAGACAGCCTACAGCAccc                                                                                                                                                                                 >  2:551213/1‑251 (MQ=255)
                                                 tgtgAACCTCACAGAGGTC‑TTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGGTACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACCCTACATGT‑GTTCAGCATATAGGAGGCCCCGGGTTG‑ATGGGAAAATATCACTCGGGGCTTTTC‑TCTATCTGCCGTTCAGCTAAT‑GCCTGAGACAGACAGCCTCCAGCccccc                                                                                                                                                                                >  2:573213/1‑250 (MQ=255)
                                                  gtgAACCTCACAGGGGTC‑TTTTCTCTTTCCCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGT‑GTTCAGCATATAGGAGGCCTCGGGATG‑ATGGTAAAATATCACTCGGGGCTTTTC‑TCTATCTGCCGTTCAGCTAAT‑GCCTGAGACAGACAGCCTCAAGCACccg                                                                                                                                                                                <  2:158514/250‑1 (MQ=255)
                                                     aaCCTCACCGAGGAC‑TTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGAGCCAGACGACGATCCGGGCCACAATCCTCGACCTATTCTTCTTCCTTGCTCCCTTCTCCTTGTCCTTTCAGTTAGTAAGAGGCCCCCTACATGC‑GTTCAGCCTATAGGAGGCCTCGGGGTG‑ATGGTAAAAAATAACTCTGGGCTTTTC‑TCTATCTGCCGTTCAGCTAAT‑TCCTGAGACAGACAGCCTCAAGCACCcgccgc                                                                                                                                                                            >  2:497089/1‑251 (MQ=255)
                                                               agGTC‑TTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCCTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGT‑GTTCAGCATATAGGAGGCCTCGGGTTG‑ATGGTAAAATATCACTCGGGGCTTTTC‑TCTATCTGCCGTTCAGCTAAT‑GCCTGAGACAGACAGCCTCAAGAACCCGCCGCTATTATATCg                                                                                                                                                                  >  2:124333/1‑251 (MQ=255)
                                                                                   cccccccccacgcggggttGATACAGATGACGATAGGGGCGAAAATCCTCGCCTTGTGCTTCTTCTTTCCTCTTTTCCCCTTGCCCTTAGGGGCAGAAAGAGCCACTCTAC‑TGTTGTCCAGCATAGAGGAGGCCTGGGTTTG‑ATGATAAAATACCCCTGGGGTCTTTTC‑TCCACCTCCCGTGCA‑CTAATTCCCTGAGACAGCCACCCTAAACCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTg                                                                                                                                               <  2:470617/233‑1 (MQ=255)
                                                                                    cccccccccTACGGCGGTAAAACAGATGACGATAAGGGCGACAACCACCGCCTTTTGCTGCTTCATTCCTCTCTTCCCCTTGCCCTTACGGTAAGTAAGGGCCCCTCAAAATTT‑GTTCAGCATATAGGGGCCCGCGGTTTG‑ATGGAAAAATATCACTCGGGTTTTTTC‑TCTATCTCCCTTTCAGCAAAG‑CCCTGAGACAGACCCCCTCAACCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGt                                                                                                                                              <  1:464235/243‑1 (MQ=255)
                                                                                     gccgccACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGT‑GTTCAGCATATAGGAGGCCTCGGGTTG‑ATGGTAAAATATCACTCGGGGCTTTTC‑TCTATCTGCCGTTCAGCTAAT‑GCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCCTTTTGtt                                                                                                                                             >  1:524475/1‑251 (MQ=255)
                                                                                                ggcggTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGT‑GTTCAGCATATAGGAGGCCTCGGGTTG‑ATGGTAAAATATCACTCGGGGCTTTTC‑TCTATCTGCCGTTCAGCTAAT‑GCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCt                                                                                                                                  >  1:135096/1‑251 (MQ=255)
                                                                                                       aTACAGATGACGACCGGGGCGACACTCACCTCTTTATGCTGCTTCTTTCCTCTCTTCCCTTTCTCCTTGGGGTAAGAAAGGGGCCTTCTAATTGT‑TTTCAGAATATAGGAGGCTTGGGGTTG‑AGGGTAAAATATCCCTGGGGGTTTTTT‑TCTTTCTGCCTTTAAGCAAAT‑GCTTGAGACAGACAGCCCCAAGCCCCCGCCGCTATTATATCGCTCTCTTGAACCCATTTTGTTTTATCGATTCTAATCCTg                                                                                                                           <  1:462434/251‑1 (MQ=255)
                                                                                                           aGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGT‑GTTCAGCATATAGGAGGCCTCGGGTTG‑ATGGTAAAATATCACTCGGGGCTTTTC‑TCTATCTGCCGTTCAGCTAAT‑GCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGa                                                                                                                          >  2:120995/1‑248 (MQ=255)
                                                                                                                        ggCGAAAATCATCGCTTTATGCTGCTTCTGTGCTCTCTCCCCCTTGCCCTTCCGGTCAGTAAGAGGCACTCT‑CATGTTGTGAAGAAAAGGGGAGGCCTGGGGTGG‑ATGTAAAAATATCACCGGGGGTTTTTC‑TCTATCTGCCGTCAAGTTATG‑GCCGGAGACAGCCAGCCCCAACCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCAtttt                                                                                                          <  1:201458/251‑1 (MQ=255)
                                                                                                                           gACAATCATCGCCTTAGGCTCCTTCATTGCTCTCTTCTCCTTGCCCTTACGGTCAGAAAGAGGCCCTCTACATGT‑GTTCAGCATATAGGAGGCCTGGGTTTG‑ATGGAAAAATACCCCTCGGGGTTTTTC‑TCTATCTGCCGTTCAGCTAAT‑GCCTGAGACAGACACCCTCAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTtgt                                                                                                       <  1:30799/251‑1 (MQ=255)
                                                                                                                               atcatcGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGT‑GTTCAGCATATAGGAGGCCTCGGGTTG‑ATGGTAAAATATCACTCGGGGCTTTTC‑TCTATCTGCCGTTCAGCTAAT‑GCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCg                                                                                                   >  2:583951/1‑251 (MQ=255)
                                                                                                                                   tctccTTTTGCTGCTTCATTCCTCTCTTCTCCTTGCCTTTACTGTAAGAAAGAGGCCCTCTATATGT‑GTTCATCAAAGAGGAGGCCT‑GGGTTGGATGATAAAATATCCCGCGGGTTTTTTC‑TCTTTCTGCCGTTCACATAAT‑GCCAGAAAAAGACAGCCTCAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTa                                                                                                 <  1:124333/246‑1 (MQ=255)
                                                                                                                                        ttATGCTGCTTCATTGCTCTCTTCTCCTTGCCCTTACGGGCAGAAAGGGCCACTCTACATGT‑GTACAGCATATAGGAGGCCGGGGGGTG‑ATGGAAAAATATCTCGGGGGTTTTTTC‑TCTATCTGCCGTGCAGTTATT‑GCGTGAAACAAACAGCCCAAACCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAAttttt                                                                                           <  1:514984/250‑1 (MQ=255)
                                                                                                                                               cctcctttcctctttcctcCTTGCCCTTAGGGACAGAAAGAGCCCCTCTACATGT‑GTTCAAAAAATGGGAGGCCTGGGGTGG‑ATGAAAAAATATCCCGGGGGTTTTTTC‑TCTTCCTCCCGTTCAGCAAAT‑GCCTGAGACAGACACCCTCAACCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGAtt                                                                                   <  1:299075/235‑1 (MQ=255)
                                                                                                                                                  tCATTGCTCTCTTCTCCTTGCTCTTACGGTCAGAAAGAGGCCCTCAACTTGT‑GTTCAGAATATAGGGGGCCTGGGTTTG‑ATGTAAAAAAATCCCGCGGGCTTTTTC‑TTTATCTGCCTTTCAGCTAAT‑GCCTAAGACAGACACCCTCAACCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGAttt                                                                                  <  1:554140/249‑1 (MQ=255)
                                                                                                                                                           tcttctCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGT‑GTTCAGCATATAGGAGGCCTCGGGTTG‑ATGGTAAAATATCACTCGGGGCTTTTC‑TCTATCTGCCGTTCAGCTAAT‑GCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTACTTATTTacc                                                                       >  2:406737/1‑249 (MQ=255)
                                                                                                                                                                cccTTGCCCTTACGGTCGATAAGAGCCCCTCTACTTGT‑GTTAAGCATAAAGGAGGCCTCGGGTTG‑ATGAAAAAATATCAC‑GGGGGTTTTTCCTCTTTCTGCCGTTCAGCTAAG‑TCCGGAGACAGACACCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTaa                                                                  <  1:369468/250‑1 (MQ=255)
                                                                                                                                                                          tACGGTAAGTAAGAGGCACTCTACTTGT‑TTTCAGCATATAGGGGGCCTCGGGTGG‑ATGTAAAAATATCCCGCGGGTTTTTTC‑TCTATCTGCCGTTCAGCAAAT‑CCCTGAGACAGACACCCTCAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTAtctc                                                          <  1:585598/249‑1 (MQ=255)
                                                                                                                                                                             ggTCAGTAAGAGGCACTCTACATGT‑GTTCAGCATATAGGAGGCCTCGGGTTG‑ATGGTAAAATATCACTCGGGGCTTTTC‑TCTATCTGCCGTTCAGCTAAT‑GCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCAt                                                     >  2:528943/1‑251 (MQ=255)
                                                                                                                                                                                                cattgt‑gTTGAGCATATGGGGGGCCTCGGGTGG‑ATGTAAAAATACCACGGGGGGTTTTTC‑TCTTTCTGCCGTTGACCTAAT‑CCCTGAGACAGACACCCCCAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACa                                    <  1:499867/246‑1 (MQ=255)
                                                                                                                                                                                                    gt‑gtTCAGCATATAGGAGGCCTCGGGTTG‑ATGGTAAAATATCACTCGGGGCTTTTC‑TCTATCTGCCGTTCAGCTAAT‑GCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAGCTgg                              >  1:588984/1‑251 (MQ=255)
                                                                                                                                                                                                       gtTCAGCATATAGGAGGCCTCGGGTTG‑ATGGTAAAATATCACTCGGGGCTTTTC‑TCTATCTGCCGTTCAGCTAAT‑GCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAGCTGGat                            >  1:124063/1‑251 (MQ=255)
                                                                                                                                                                                                          cagaatatgggggtcttggggTTG‑ATGTAAAAATATCCCTGGGGTTTTTTC‑TCTCTCTCCCGTCCAGCTAAT‑GCCTGAGACAGACACCCCCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAGCTGGattat                         <  1:296865/233‑1 (MQ=255)
                                                                                                                                                                                                          aagattaggggcggCCTGGGGTGG‑AGGGTAAAATTTCACTGGGGGTTTTTC‑TCTACCGGCTGTTCACCTAAG‑CCCGGAGACAGACAGCTTAAACCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAGCTGGattat                         <  1:271714/239‑1 (MQ=255)
                                                                                                                                                                                                                ataGGAGGCCTCGGGTTG‑ATGGTAAAATATCACTCGGGGCTTTTC‑TCTATCTGCCGTTCAGCTAAT‑GCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGCTTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAGCTGGATTATTTTTg                    >  2:438773/1‑250 (MQ=255)
                                                                                                                                                                                                                     aggCCTCGGGTTG‑ATGGTAAAATATCACTCGGGGCTTTTC‑TCTATCTGCCGTTCAGCTAAT‑GCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAatat              >  2:382882/1‑251 (MQ=255)
                                                                                                                                                                                                                           cGGGTTG‑ATGGTAAAATATCACTCGGGGCTTTTC‑TCTATCTGCCGTTCAGCTAAT‑GCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGGCCt        >  1:492867/1‑251 (MQ=255)
                                                                                                                                                                                                                                  gatgtaaaaaTATCACTCGGGGTTTTTC‑TCTATCGCCCGTTCAGCTAAT‑GCCTGAGACAGCCACCCCCAACCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGGCCTAACAAg  <  1:583951/245‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                      
CGGATTCGTAAGCCGTGAAAACAGCAACCT‑CCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTC‑TTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGT‑GTTCAGCATATAGGAGGCCTCGGGTTG‑ATGGTAAAATATCACTCGGGGCTTTTC‑TCTATCTGCCGTTCAGCTAAT‑GCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGGCCTAACAAG  >  NC_000913/16756‑17228

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