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	81.1%	91.0 / 14.0	37	intergenic (+338/‑141)	sokC/nhaA	small regulatory RNA antitoxin SokC/Na(+):H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base . (6/1);  new base A (13/17);  total (19/18)
Fisher's exact test for biased strand distribution p-value = 8.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGGCCTAACAAGCATCGCCGACTGACAACAAATTAATTATTACTTTTCCTAATTAATCCCTCAGGAATCCTCACCTTAAGCTATG‑ATTATCTAGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑A‑C‑GCATCTCACCGCTGATGGCGCAAATTCTTCAATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAAC  >  NC_000913/17108‑17589                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   aaGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑aa‑aa‑aactcatctc                                                                                                                                                                                                                                            >  1:505923/1‑241 (MQ=255)      gCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAAGCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCt                                                                                                                                                                                                                                     <  2:417908/251‑1 (MQ=255)                             tttttAATTATTTAATTATTTAACTTTAATTTAACTCTTTAACGTAATTATTGACGACAATCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGATTGACAATAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCTCCTTAATCTATG‑ATTATCTGGGCTTAGGTTCACTCGTGAGCGTTTACAGCCGTC‑AA‑AA‑AAC‑TCATCCCACCGTTGATGTCGCAAATTCTTCAAtgg                                                                                                                                                                                                               <  2:494893/250‑3 (MQ=255)                                   aTGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTATAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGATGATGGCGCACATTATTCAATGGATCGTa                                                                                                                                                                                                         >  2:397388/1‑250 (MQ=255)                                              attTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTTGATTATTTTTGAAATATTGTCCAAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCGTAAGCTATG‑ATTATCGGGGCTTAGGGTCACTCGTGAGCGCTGTCAGCCGTC‑AA‑AA‑AAC‑ACATCTCTCCGCGGATGGCGGAAATTCTTCAATGGATCGTAAAAAACGagtt                                                                                                                                                                                              <  2:308009/250‑4 (MQ=255)                                              attTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAACCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGagtt                                                                                                                                                                                              <  2:151532/250‑4 (MQ=255)                                                   aCTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTACTTAATCCCTCAGAAATCCTCACCTTAAGCTATA‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AACA‑C‑TAATATCACCGCTGATGGCACAAATTCTTCAATGGATCGTAAAAAACgtgtaattcct                                                                                                                                                                                        >  2:278499/1‑241 (MQ=255)                                                           ttATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGTAATCTTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGACGCAAATTCTTCAATGGCTCGTAAAAAACGAGTTATTCATACACtat                                                                                                                                                                                 >  2:426814/1‑250 (MQ=255)                                                               ctctTTATCGCACATATTGACGACAAGCTGGATTATTTTTGAAATATTGTCATAACAAATAACGCCGCCTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCCCCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTAGTGAGCGCTTACAGCCGTC‑AA‑AA‑A‑CCTCATATCACCGCTGATGGCGCAAATTCTTCAATGGTTCGAAAACAAAGAGTTATTCCTACACTATCATCt                                                                                                                                                                            >  2:378540/1‑251 (MQ=255)                                                                     agcGCAAATATTGACGACAAGCTGGATTATTTGTGTAAGATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGTGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGAtgt                                                                                                                                                                       <  2:204937/248‑3 (MQ=255)                                                                     aTCGCAATTATTGACGACAAGATGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGAtttt                                                                                                                                                                      >  2:323167/1‑251 (MQ=255)                                                                            ttattGACGACAAGCTGGATTATTGTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGACATCCTCACATTAAGCTATG‑ATTATCTGGACTTAGGGTCACTCGTGAGCGCTTACATCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACTAGTTATTCCTACACTATAATCTGATTTTAACgatg                                                                                                                                                               >  2:378945/1‑251 (MQ=255)                                                                                      cAAGCTGGATTATTTTTGAAATATTGCCCTACCAATTACCGCCGCCTGCCAAGAAATTAAT‑ATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCggg                                                                                                                                                    <  1:223063/251‑1 (MQ=255)                                                                                      aaaGCTGGATTATTTTGTAAAAATTTTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTGCCTAATTAATCCCTCAGAAATCCTCACCTTAAGC‑ATGGATTATCTGGGCTTAGGGGCACTCGTGAGCGCTTACAGTCGTC‑AA‑AA‑AAC‑GCATCTCCCCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCgg                                                                                                                                                     <  2:303611/250‑1 (MQ=255)                                                                                           tGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGag                                                                                                                                                                                                >  1:56401/1‑202 (MQ=255)                                                                                           tGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAAGCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGag                                                                                                                                                                                                <  2:56401/203‑2 (MQ=255)                                                                                              attattTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAAt                                                                                                                                             <  1:341175/251‑1 (MQ=255)                                                                                               ttattTTTGAAATCTTGTCCTAACAAATAACGCCGAATGAAAAGAACTTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATA‑ATTATCTGGGATTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AACAA‑A‑C‑TCATCTCACCGCTGATGGCGCAACTTCTTCAATGGATCGTAAAAAACGTGTTATTCCTACACTATAACATGATTTTAACgatgctttattcggggaacaat                                                                                                                                             >  1:365850/1‑232 (MQ=255)                                                                                                    tttGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGTTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGAAAAATAGTaa                                                                                                                                        >  1:293822/1‑250 (MQ=255)                                                                                                      tGAAATATTGGCCTAACAAGCATCGCCGACTGACAACAAATTAATTATTACTTTTCCTAATTAATCCCTCAGGAATCCTCACCTTAAGCTATG‑ATTATCTAGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑A‑C‑GCATCTCACCGCTGATGGCGCAAATTCTTCAATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTaaaaa                                                                                                                                     <  1:190485/250‑1 (MQ=255)                                                                                                            atTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGACACTCGTGAGCGCTTACAGCCGTC‑AA‑AATA‑C‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGATATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAGACTAGTAAACACGATCt                                                                                                                               >  1:244306/1‑251 (MQ=255)                                                                                                                     acaaatcatgCCGACTGACAAGAAATTACTTATTACCCTTCCTAATTAACCCCTCACAAATACTCACCTTAAGCTATG‑ATTTTCTGGGCTTAGGCTCACTAGTGACCGCTTACAGCCGTC‑AA‑AA‑A‑CATCATCTCACCGCTGACGGCGCAAATTCTTCAATGGATCGTAAAAAAAGAATTATTCCTACACTATAATCTGATTTTAACAATGATTTATGCGGGGTAAACTAGTAAAAAAGCTCTATCCACCTg                                                                                                                      >  2:236474/10‑251 (MQ=255)                                                                                                                     aaaaataacGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCAACTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTg                                                                                                                      <  2:427193/243‑1 (MQ=255)                                                                                                                            acGCCGACTGACAAGAAAGTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAgaga                                                                                                               <  1:483827/250‑1 (MQ=255)                                                                                                                                        aagaaaTTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGa                                                                                                   <  2:37814/248‑1 (MQ=255)                                                                                                                                          aaaaTTAATTATTACTTTACCTAATTAATCCCTCAGAAAACCTCACCTTAAGCTATT‑ATTATCTGGGCTTAGGGTCACTCGTGATCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCATATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCTGGGTAAAATAATAAAAACGATCTATTCCCCTGAAAGAGAAATAAAAAGTGaaa                                                                                                 <  2:62093/250‑1 (MQ=255)                                                                                                                                                 attattACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGATATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAACCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTCACGATGATTTGTGCGGGGTAAATTAGTAAAAACGATCTTTTCACCTGAAAGAGAA‑Taaaaacgacaaacatatgc                                                                                         >  2:444916/1‑237 (MQ=255)                                                                                                                                                       aCTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATACAAAGTGAAACATCTGCATCGAt                                                                                    >  1:270606/1‑251 (MQ=255)                                                                                                                                                                aattaatCCCTCAAAAATCCTCACCTTAAGCTATG‑ATTAACTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCCCCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATGTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCa                                                                           <  2:99928/251‑1 (MQ=255)                                                                                                                                                                  ttaatCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGt                                                                         >  2:262957/1‑251 (MQ=255)                                                                                                                                                                    aatCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTATTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAACTAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACCTCTGCATCGATTCTTTAGCAGTgg                                                                       >  2:123799/1‑250 (MQ=255)                                                                                                                                                                           cagaaaTCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTc                                                                 <  1:406917/246‑1 (MQ=255)                                                                                                                                                                               aaTCCTTACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCTTGAGCGCTTACAGCCGTA‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCTCCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAg                                                            <  2:401456/251‑1 (MQ=255)                                                                                                                                                                               aaTCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGCGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCggg                                                              >  1:217932/1‑249 (MQ=255)                                                                                                                                                                                        cctAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCGATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCGCGGGAGGCATTATTc                                                   <  2:302391/249‑1 (MQ=255)                                                                                                                                                                                              cTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATc                                              >  2:367861/1‑250 (MQ=255)                                                                                                                                                                                                                 aGGGTCACTCGTGAGCGCTTACAGCCGTCTAA‑AA‑A‑C‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCg                            >  1:93143/1‑250 (MQ=255)                                                                                                                                                                                                                 aGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCgg                           >  1:103665/1‑250 (MQ=255)                                                                                                                                                                                                                                          cggcaaa‑aa‑a‑C‑TCATCGCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAATGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAAc  <  1:209338/246‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   AAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGGCCTAACAAGCATCGCCGACTGACAACAAATTAATTATTACTTTTCCTAATTAATCCCTCAGGAATCCTCACCTTAAGCTATG‑ATTATCTAGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑A‑C‑GCATCTCACCGCTGATGGCGCAAATTCTTCAATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAAC  >  NC_000913/17108‑17589

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