Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 17,436 C→T intergenic (+426/‑53) sokC → / → nhaA small regulatory RNA antitoxin SokC/Na(+):H(+) antiporter NhaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091317,4360CT100.0% 120.1 / NA 39intergenic (+426/‑53)sokC/nhaAsmall regulatory RNA antitoxin SokC/Na(+):H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/17);  total (22/17)

TTTTTGAAATATTGGCCTAACAAGCATCGCCGACTGACAACAAATTAATTATTACTTTTCCTAATTAATCCCTCAGGAATCCTCACCTTAAGCTATGATTATCTAGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CAAAAA‑CGCATCTCACCGCTGATGGCGCAAATTCTTCAATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATG  >  NC_000913/17205‑17681
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     
tttttGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGt                                                                                                                                                                                                                                      >  2:420499/1‑250 (MQ=255)
     gAAATATTGTCCTAACAAATAACGCCGAATGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTaaaat                                                                                                                                                                                                                                 <  2:205211/250‑2 (MQ=255)
        atatTGTCCTAACAAATAACGCCGACTGACAAGCAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGAt                                                                                                                                                                                                                             >  1:497876/1‑251 (MQ=255)
                       ataacGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGaa                                                                                                                                                                                                                >  1:39926/5‑249 (MQ=255)
                       ataacGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCACTCAGAAATCCTCACCTTAAGCTCTGTTTATCTCGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGCTTTGCGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGaa                                                                                                                                                                                                                >  1:414704/5‑249 (MQ=255)
                               gACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATaa                                                                                                                                                                                                      <  1:380938/251‑1 (MQ=255)
                                      aagaaaTTAATTATTTCTTTTCCTAATTAATCCCTCAGAAAGCCGCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGA‑CAAAAAACTCATCTCTCCGCTGATTGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACATTATAAGCTGATTTTAACGATGATTTGGGCGGTGTAAAATAGTAAAAACGATCTAGTCGCCTGAAAGAGAAAGAAAAAGTGa                                                                                                                                                                                               <  2:564695/248‑1 (MQ=255)
                                        gaaaTTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGaa                                                                                                                                                                                              >  1:408485/2‑250 (MQ=255)
                                              aattaTTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCt                                                                                                                                                                                        <  1:624850/250‑1 (MQ=255)
                                                    taCTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCg                                                                                                                                                                                  >  1:386136/1‑250 (MQ=255)
                                                           ccTAATGAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCATAGTGTCAATCGTGAGCGCTTACAGCCGT‑CAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGGGCGGGGTAAAATAGTAAAAACGATCGATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGTTTCTTTa                                                                                                                                                                          <  2:51532/251‑1 (MQ=255)
                                                            cTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCttt                                                                                                                                                                           >  2:160993/1‑249 (MQ=255)
                                                                 taatCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTg                                                                                                                                                                    >  2:6187/1‑251 (MQ=255)
                                                                        tCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTc                                                                                                                                                             >  2:374203/1‑251 (MQ=255)
                                                                         cagaaaTCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CAAAAAACTCATCTCACCGCTGATGTCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACTCTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATATTAAAAACGATCTATTCCCCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTc                                                                                                                                                             <  2:554746/246‑1 (MQ=255)
                                                                         cagaaaTCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCg                                                                                                                                                            >  1:141432/5‑251 (MQ=255)
                                                                                        gaaGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCTGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCtt                                                                                                                                             <  1:595755/250‑1 (MQ=255)
                                                                                                          gCTTAGGGTCACTCGTGAGCGCTTACAGCCGTTAAAAAA‑CTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTg                                                                                                                           <  1:41290/251‑1 (MQ=255)
                                                                                                                    aCTCGTGAGCGCTTACAGCCGT‑CAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATACCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGAtt                                                                                                                  >  1:443913/1‑250 (MQ=255)
                                                                                                                       cGTGAGCGCTTACAGCCGT‑CAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTAt                                                                                                                <  1:428952/249‑1 (MQ=255)
                                                                                                                                  aCAGCCGT‑CAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATGTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGTGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAgcggc                                                                                                   <  2:291636/251‑1 (MQ=255)
                                                                                                                                             aaaa‑CTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGCGTTATTCCTACACTATAATCTGATTTTAACGATGCTTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGACATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTgg                                                                                        >  2:612077/1‑251 (MQ=255)
                                                                                                                                                          aCCGCTGATGGCGCAAATTCTTCAATGGCTCGTCAAAAACGAGTTGTTCGTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGGCAAATAGTAACAACGATCTATTCACCTGAAAGAGAAATACAAAATGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCCACCAGTGGATGGTATCACgc                                                                            >  2:34682/1‑250 (MQ=255)
                                                                                                                                                                        aaaTTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACATGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGGGACGc                                                              >  1:510549/1‑251 (MQ=255)
                                                                                                                                                                                    ggatCGTAAAAAACGAGTTAATCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTAGTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGAGTCTTTAGCAGTGATGCCTCGGGAGGCAGTATTCTTATCATTGCCGCTATCCTGGCGATGATGATGTCCAATAGCGGCGCAACCAGTGTATGGTATTACGTCTTTCTGGATACGcgcgtgcagctc                                                   <  2:95923/247‑12 (MQ=255)
                                                                                                                                                                                      atCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCgg                                                <  1:613947/250‑1 (MQ=255)
                                                                                                                                                                                        cGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCggg                                               <  2:386136/250‑1 (MQ=255)
                                                                                                                                                                                                       attCCTACACTATAATCTGAGTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACTCCGGTTCAGCTCCGGGTTGGTTCACTCGaaa                                <  1:476368/250‑1 (MQ=255)
                                                                                                                                                                                                             acacTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGCTTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATcaaca                          >  2:497519/1‑250 (MQ=255)
                                                                                                                                                                                                             acacTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATcaacaa                         <  1:124606/251‑1 (MQ=255)
                                                                                                                                                                                                             acacTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATcaaca                          >  1:158800/1‑250 (MQ=255)
                                                                                                                                                                                                              cacTATCATATGATTTTCACGATGAATTGTGAGCGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGCGAAATCAAAAGTGCAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGCTGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGCCGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATcaacaa                         >  2:566583/1‑250 (MQ=255)
                                                                                                                                                                                                              cacTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGACCCATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCATTACCGGCGCACCCAGTGGCTGTTATCACGCCTTTCTGAACCCGCCGGTTCAGCTCCGGGTTGTTTCACTCgcattcaccaa                         >  1:79715/1‑240 (MQ=255)
                                                                                                                                                                                                                  ataATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAATAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACaaaaa                      <  2:491399/249‑1 (MQ=255)
                                                                                                                                                                                                                                 cGATGATTTGTGCGGGGTACAATAGTAAAAACGATCTATTCACCTGAAAGTGAAATAAAAAGTGAAACATCTGCATAGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAAATTGCTGTTATGGATa      >  1:304943/1‑250 (MQ=255)
                                                                                                                                                                                                                                   atgatTTGTGCGGGGTAAAATAGTAAGAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAAt   >  2:116922/1‑251 (MQ=255)
                                                                                                                                                                                                                                   atgatTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAAt   <  1:347926/251‑1 (MQ=255)
                                                                                                                                                                                                                                     gatTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGACATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGCAATCAACAAAAACATGCTGTTATGGATAAATg  >  2:486971/1‑250 (MQ=255)
                                                                                                                                                                                                                                     gatTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATg  >  1:219454/1‑250 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     
TTTTTGAAATATTGGCCTAACAAGCATCGCCGACTGACAACAAATTAATTATTACTTTTCCTAATTAATCCCTCAGGAATCCTCACCTTAAGCTATGATTATCTAGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CAAAAA‑CGCATCTCACCGCTGATGGCGCAAATTCTTCAATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAAATG  >  NC_000913/17205‑17681

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