Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	17,348	(A)5→6	intergenic (+338/‑141)	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_000913	17,348	1	.	A	94.6%	111.8 / 0.6	37	intergenic (+338/‑141)	sokC/nhaA	small regulatory RNA antitoxin SokC/Na(+):H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base A (19/16);  total (20/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

AAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGGCCTAACAAGCATCGCCGACTGACAACAAATTAATTATTACTTTTCCTAATTAATCCCTCAGGAATCCTCACCTTAAGCTATGATTATCTAGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CA‑AAAA‑CGCATCTCACCGCTGATGGCGCAAATTCTTCAATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGG  >  NC_000913/17108‑17595                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                        aaGACGCCTCGCATTTTTGTGGAGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑Ca‑aaaaactcatctc                                                                                                                                                                                                                                                 >  2:428952/1‑241 (MQ=255)       ccTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTGTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCGCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGGCACTCGTGAGCGCTTACAGCCGT‑CA‑AAAAACTCATCTCACCGCTg                                                                                                                                                                                                                                          <  1:269619/251‑1 (MQ=255)                  tgtgGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTCTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CA‑AAAAACTCATCTCACCGCTGATGGCGCAAAt                                                                                                                                                                                                                               >  2:613947/1‑251 (MQ=255)                   gtgGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAATTATTGTCCTAACAAATAACGCCGCCTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTAACCTTAAGCTATGATTTTCTGTGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CA‑AAAAACTCATCTCACCACTGATGGCGCAAAtt                                                                                                                                                                                                                              >  2:380938/1‑251 (MQ=255)                           aTTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAAtg                                                                                                                                                                                                                      >  1:119527/1‑250 (MQ=255)                           aTTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAAtg                                                                                                                                                                                                                      <  2:209526/251‑2 (MQ=255)                                 tAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGAAATTATTGAAGACAATCTGGTTTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAATAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAATCTATGATTATCTGTGATTAGGGTCACTCGTGCACGCTTACAGCTGT‑CA‑AAAAACTCTTCTCACCGCTGATGTCGCCATTTCTTCAAtggatcgt                                                                                                                                                                                                                >  2:356044/1‑244 (MQ=255)                                        ttAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTTAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCATAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CA‑AAAAACTAATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTaaaaaa                                                                                                                                                                                                          <  2:388717/250‑1 (MQ=255)                                                      ttAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CA‑AAAAACTCATCTCACCGCt                                                                                                                                                                                                                                           >  1:219522/1‑203 (MQ=255)                                                      ttAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CA‑AAAAACTCATCTCACCGCt                                                                                                                                                                                                                                           <  2:219522/203‑1 (MQ=255)                                                      ttAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTac                                                                                                                                                                                            <  2:254263/250‑1 (MQ=255)                                                                    taTCGCAATTATTGACGACAAGATGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGAttt                                                                                                                                                                             >  1:491399/2‑251 (MQ=255)                                                                     aTCGCAATTAATGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGAttt                                                                                                                                                                             <  2:216223/250‑1 (MQ=255)                                                                          aaTTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACgt                                                                                                                                                                       >  2:347926/1‑250 (MQ=255)                                                                             tattGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAAAGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACgatg                                                                                                                                                                     >  1:78231/1‑250 (MQ=255)                                                                              attGACGACAAGCTGGATTATTTTTGAAATATTGTCCTATCAAATAACGCCGACTGACAAGAAATTTATTATTACTTTTCCTAATGAATCACTCAGAAATCCTCACCTTAAGCTATGATTATCTGGACTTAGGTTCACTCGTGAGCGCTTACAGCCGT‑CA‑AAAAACTCATCTCACCGCTGATGACTCAAATTCTTCAATTTATCGTAAAAAACGAGTTATTCCTCCACTATAATCTGATTTTAACgatgat                                                                                                                                                                   <  2:342178/251‑1 (MQ=255)                                                                                       acgCTGGCTTATTTTTGAAATATAGACCTACCACTTAAAGCCGCCTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCATCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGACACTCGTGAGCGCTTACAGCCGT‑CACAAAC‑CTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAACAAACGCGTTATTCCTACACTATAATCTGATTTTAACgatgctttgtgaggg                                                                                                                                                          >  2:570014/3‑240 (MQ=255)                                                                                                  tttttGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGt                                                                                                                                                >  2:420499/1‑250 (MQ=255)                                                                                                       gAAATATTGTCCTAACAAATAACGCCGAATGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTaaaat                                                                                                                                           <  2:205211/250‑2 (MQ=255)                                                                                                          atatTGTCCTAACAAATAACGCCGACTGACAAGCAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGAt                                                                                                                                       >  1:497876/1‑251 (MQ=255)                                                                                                                         ataacGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGaa                                                                                                                          >  1:39926/5‑249 (MQ=255)                                                                                                                         ataacGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCACTCAGAAATCCTCACCTTAAGCTCTGTTTATCTCGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGCTTTGCGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGaa                                                                                                                          >  1:414704/5‑249 (MQ=255)                                                                                                                                 gACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATaa                                                                                                                <  1:380938/251‑1 (MQ=255)                                                                                                                                        aagaaaTTAATTATTTCTTTTCCTAATTAATCCCTCAGAAAGCCGCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGA‑CA‑AAAAACTCATCTCTCCGCTGATTGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACATTATAAGCTGATTTTAACGATGATTTGGGCGGTGTAAAATAGTAAAAACGATCTAGTCGCCTGAAAGAGAAAGAAAAAGTGa                                                                                                         <  2:564695/248‑1 (MQ=255)                                                                                                                                          gaaaTTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGaa                                                                                                        >  1:408485/2‑250 (MQ=255)                                                                                                                                                aattaTTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCt                                                                                                  <  1:624850/250‑1 (MQ=255)                                                                                                                                                      taCTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCg                                                                                            >  1:386136/1‑250 (MQ=255)                                                                                                                                                             ccTAATGAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCATAGTGTCAATCGTGAGCGCTTACAGCCGT‑CA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGGGCGGGGTAAAATAGTAAAAACGATCGATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGTTTCTTTa                                                                                    <  2:51532/251‑1 (MQ=255)                                                                                                                                                              cTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCttt                                                                                     >  2:160993/1‑249 (MQ=255)                                                                                                                                                                   taatCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTg                                                                              >  2:6187/1‑251 (MQ=255)                                                                                                                                                                          tCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTc                                                                       >  2:374203/1‑251 (MQ=255)                                                                                                                                                                           cagaaaTCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CA‑AAAAACTCATCTCACCGCTGATGTCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACTCTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATATTAAAAACGATCTATTCCCCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTc                                                                       <  2:554746/246‑1 (MQ=255)                                                                                                                                                                           cagaaaTCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCg                                                                      >  1:141432/5‑251 (MQ=255)                                                                                                                                                                                          gaaGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCTGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCtt                                                       <  1:595755/250‑1 (MQ=255)                                                                                                                                                                                                            gCTTAGGGTCACTCGTGAGCGCTTACAGCCGTTAA‑AAAA‑CTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTg                                     <  1:41290/251‑1 (MQ=255)                                                                                                                                                                                                                      aCTCGTGAGCGCTTACAGCCGT‑CA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATACCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGAtt                            >  1:443913/1‑250 (MQ=255)                                                                                                                                                                                                                         cGTGAGCGCTTACAGCCGT‑CA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTAt                          <  1:428952/249‑1 (MQ=255)                                                                                                                                                                                                                                    aCAGCCGT‑CA‑AAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATGTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGTGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAgcggc             <  2:291636/251‑1 (MQ=255)                                                                                                                                                                                                                                                aaaa‑CTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGCGTTATTCCTACACTATAATCTGATTTTAACGATGCTTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGACATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTgg  >  2:612077/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                        AAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGGCCTAACAAGCATCGCCGACTGACAACAAATTAATTATTACTTTTCCTAATTAATCCCTCAGGAATCCTCACCTTAAGCTATGATTATCTAGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGT‑CA‑AAAA‑CGCATCTCACCGCTGATGGCGCAAATTCTTCAATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGG  >  NC_000913/17108‑17595

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