Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	17,219	G→T	intergenic (+209/‑270)	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,219	0	G	T	92.1%	106.2 / ‑0.1	39	intergenic (+209/‑270)	sokC/nhaA	small regulatory RNA antitoxin SokC/Na(+):H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base T (22/13);  total (23/16)
Fisher's exact test for biased strand distribution p-value = 5.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

GCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGGCCTAACAAGCATCGCCGACTGACAACAAATTAATTATTACTTTTCCTAATTAATCCCTCAGGAATCCTCACCTTAAGCTATG‑ATTATCTAGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑A‑C‑GCATCTCACCGCTGATGGCGCAAATTCTTCAATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCT  >  NC_000913/16999‑17464                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                           gCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAaat                                                                                                                                                                                                                              >  1:151532/1‑251 (MQ=255)  cTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAaatt                                                                                                                                                                                                                             >  1:401456/1‑251 (MQ=255)                   ttCAGCTAATGCCTGAGACAGATAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCt                                                                                                                                                                                                             >  2:48024/1‑250 (MQ=255)                                            ccccaaccccccGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCt                                                                                                                                                                                   <  1:100885/242‑1 (MQ=255)                                             cccAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCtt                                                                                                                                                                                  <  1:50099/249‑1 (MQ=255)                                                aaGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTACTTTTTTAATGATTTCATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAGCTGGCTTATTTTTGAACTATTGGCCTAACAAGCATCGCCGCCTGACAACAAATTAATTATTACTTTTCCTAATTAATCCCTCAGGAATCCtaccctta                                                                                                                                                                                 >  2:301432/1‑242 (MQ=255)                                                 aGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGc                                                                                                                                                                              >  1:427373/1‑251 (MQ=255)                                                          cgcTATTATATCGCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTAt                                                                                                                                                                    <  1:203224/251‑1 (MQ=255)                                                                      gCTCTCTTTAACCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAggg                                                                                                                                                        >  1:60793/1‑251 (MQ=255)                                                                             ttAACCCATTCTGTTTTATCGATTCTAATCCTTAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATTATTTAATTATTTAACTTGAATTTATCTCATTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGGCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑AATATCTGGGCTTAGGGTCAc                                                                                                                                                    <  2:396938/248‑1 (MQ=255)                                                                               aaCCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGacaaga                                                                                                                                                                                                                                 <  2:389162/170‑3 (MQ=255)                                                                               aaCCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGacaaga                                                                                                                                                                                                                                 >  1:389162/1‑168 (MQ=255)                                                                               aaCCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTg                                                                                                                                               >  1:417908/1‑251 (MQ=255)                                                                                aCCCATTCTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGa                                                                                                                                              <  1:79782/251‑1 (MQ=255)                                                                                    attctGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAgcgc                                                                                                                                          >  1:37814/5‑251 (MQ=255)                                                                                           tttATCGATTCTACTCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGc                                                                                                                                   >  2:330748/1‑251 (MQ=255)                                                                                           tttATCGATTCTAATCCTGACGCCGCCTCGCATTTTTGTGGCGTAATTTTTTACTGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGATGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTATTTAATCCCTCAGAAATCCTAACCTTAAGCTATG‑ATTATCTTGCCTTAGGGTCACTCGTGAGCGCTTAccgc                                                                                                                                   >  2:91218/1‑248 (MQ=255)                                                                                             tATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTGTCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCg                                                                                                                                 >  2:406917/1‑251 (MQ=255)                                                                                                              aaGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCaa‑aa‑aactcatctc                                                                                                               >  1:505923/1‑241 (MQ=255)                                                                                                                   gCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAAGCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑AAC‑TCATCTCACCGCt                                                                                                        <  2:417908/251‑1 (MQ=255)                                                                                                                                          tttttAATTATTTAATTATTTAACTTTAATTTAACTCTTTAACGTAATTATTGACGACAATCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGATTGACAATAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCTCCTTAATCTATG‑ATTATCTGGGCTTAGGTTCACTCGTGAGCGTTTACAGCCGTCAA‑AA‑AAC‑TCATCCCACCGTTGATGTCGCAAATTCTTCAAtgg                                                                                  <  2:494893/250‑3 (MQ=255)                                                                                                                                                aTGATTTAATTATTTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTATAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑AAC‑TCATCTCACCGATGATGGCGCACATTATTCAATGGATCGTa                                                                            >  2:397388/1‑250 (MQ=255)                                                                                                                                                           attTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTTGATTATTTTTGAAATATTGTCCAAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCGTAAGCTATG‑ATTATCGGGGCTTAGGGTCACTCGTGAGCGCTGTCAGCCGTCAA‑AA‑AAC‑ACATCTCTCCGCGGATGGCGGAAATTCTTCAATGGATCGTAAAAAACGagtt                                                                 <  2:308009/250‑4 (MQ=255)                                                                                                                                                           attTAACTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAACCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGagtt                                                                 <  2:151532/250‑4 (MQ=255)                                                                                                                                                                aCTTTAATTTATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTACTTAATCCCTCAGAAATCCTCACCTTAAGCTATA‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AACA‑C‑TAATATCACCGCTGATGGCACAAATTCTTCAATGGATCGTAAAAAACgtgtaattcct                                                           >  2:278499/1‑241 (MQ=255)                                                                                                                                                                        ttATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGTAATCTTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑AAC‑TCATCTCACCGCTGATGACGCAAATTCTTCAATGGCTCGTAAAAAACGAGTTATTCATACACtat                                                    >  2:426814/1‑250 (MQ=255)                                                                                                                                                                            ctctTTATCGCACATATTGACGACAAGCTGGATTATTTTTGAAATATTGTCATAACAAATAACGCCGCCTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCCCCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTAGTGAGCGCTTACAGCCGTCAA‑AA‑A‑CCTCATATCACCGCTGATGGCGCAAATTCTTCAATGGTTCGAAAACAAAGAGTTATTCCTACACTATCATCt                                               >  2:378540/1‑251 (MQ=255)                                                                                                                                                                                  agcGCAAATATTGACGACAAGCTGGATTATTTGTGTAAGATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGTGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGAtgt                                          <  2:204937/248‑3 (MQ=255)                                                                                                                                                                                  aTCGCAATTATTGACGACAAGATGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGAtttt                                         >  2:323167/1‑251 (MQ=255)                                                                                                                                                                                         ttattGACGACAAGCTGGATTATTGTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGACATCCTCACATTAAGCTATG‑ATTATCTGGACTTAGGGTCACTCGTGAGCGCTTACATCCGTCAA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACTAGTTATTCCTACACTATAATCTGATTTTAACgatg                                  >  2:378945/1‑251 (MQ=255)                                                                                                                                                                                                   cAAGCTGGATTATTTTTGAAATATTGCCCTACCAATTACCGCCGCCTGCCAAGAAATTAAT‑ATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCggg                       <  1:223063/251‑1 (MQ=255)                                                                                                                                                                                                   aaaGCTGGATTATTTTGTAAAAATTTTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTGCCTAATTAATCCCTCAGAAATCCTCACCTTAAGC‑ATGGATTATCTGGGCTTAGGGGCACTCGTGAGCGCTTACAGTCGTCAA‑AA‑AAC‑GCATCTCCCCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCgg                        <  2:303611/250‑1 (MQ=255)                                                                                                                                                                                                        tGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGag                                                                   >  1:56401/1‑202 (MQ=255)                                                                                                                                                                                                        tGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAAGCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGag                                                                   <  2:56401/203‑2 (MQ=255)                                                                                                                                                                                                           attattTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAAt                <  1:341175/251‑1 (MQ=255)                                                                                                                                                                                                            ttattTTTGAAATCTTGTCCTAACAAATAACGCCGAATGAAAAGAACTTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATA‑ATTATCTGGGATTAGGGTCACTCGTGAGCGCTTACAGCCGTCAACAA‑A‑C‑TCATCTCACCGCTGATGGCGCAACTTCTTCAATGGATCGTAAAAAACGTGTTATTCCTACACTATAACATGATTTTAACgatgctttattcggggaacaat                >  1:365850/1‑232 (MQ=255)                                                                                                                                                                                                                 tttGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGTTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGAAAAATAGTaa           >  1:293822/1‑250 (MQ=255)                                                                                                                                                                                                                   tGAAATATTGGCCTAACAAGCATCGCCGACTGACAACAAATTAATTATTACTTTTCCTAATTAATCCCTCAGGAATCCTCACCTTAAGCTATG‑ATTATCTAGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑A‑C‑GCATCTCACCGCTGATGGCGCAAATTCTTCAATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTaaaaa        <  1:190485/250‑1 (MQ=255)                                                                                                                                                                                                                         atTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGACACTCGTGAGCGCTTACAGCCGTCAA‑AATA‑C‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGATATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAGACTAGTAAACACGATCt  >  1:244306/1‑251 (MQ=255)                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                           GCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGATTTAATTATTTAACTTTAATTTATCTCTTCATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGGCCTAACAAGCATCGCCGACTGACAACAAATTAATTATTACTTTTCCTAATTAATCCCTCAGGAATCCTCACCTTAAGCTATG‑ATTATCTAGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAA‑AA‑A‑C‑GCATCTCACCGCTGATGGCGCAAATTCTTCAATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCT  >  NC_000913/16999‑17464

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