Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,560,455	+G	intergenic (+589/+167)	rtcR → / ← glpG	DNA‑binding transcriptional activator RtcR/rhomboid protease GlpG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,560,455	1	.	G	97.2%	127.7 / ‑4.3	36	intergenic (+589/+167)	rtcR/glpG	DNA‑binding transcriptional activator RtcR/rhomboid protease GlpG
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (20/15);  total (20/16)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑C‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGG  >  NC_000913/3560206‑3560686                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                        tcttctTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑cg                                                                                                                                                                                                                                         >  1:189749/1‑250 (MQ=255) ccttctttacccagaaaCGTGTTAGCATCGTTGGGATTGTTGGTGTCAATGCGCAAAGTGCTTGGATTGAGAAGTGCTTGCCCTCCCGCTTCCGGCGTGGGGCCGATATCGATAACCAGCTCCGGGCCATTGGGAATTTGCTCCGCCCCTTTGCGGGCGATGCGCTCTTTTTCTCCGTTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑c                                                                                                                                                                                                                                          <  1:3981/235‑1 (MQ=255)       tttACCATAAACGTTTTAGCAACGTTGAGATTGTTGGTGACAAGGCGCAAATGGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGGGGGGTGCCGATATCTATAAACAGCGTGGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGGGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGa                                                                                                                                                                                                                                    <  1:177538/250‑1 (MQ=255)                  cGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGCGCCAGTGGGGATTTGCTCAGAGACTTTGCGGGCGATCCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTACCACGGCGTGTTAACCACACTGGACTGCAGCGCCGCACCGCCATGATGG‑CGCAGGCTCAGAtttttct                                                                                                                                                                                                                        >  2:168182/1‑248 (MQ=255)                  cGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTAACGCTTCCGGCGTGGTGCCGCTATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGCGTTCACCGAACTGGAAGGCAGCGCCGCACAGCCATAATGG‑CGCAGGATCAGGTTTTGCt                                                                                                                                                                                                                        >  2:264075/1‑251 (MQ=255)                           aaCGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCt                                                                                                                                                                                                               >  1:32571/1‑251 (MQ=255)                            atttggagtttgtggGTGACATGGCCCAAATTGTTGTGGTGGAGCGGTGCGGGCGCTCCCCCTTCCGGCGGGGTGCCGATACCGAAAAACGGCTTCGAGCCTTTGGGGATTTGCTCCGCCGCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTc                                                                                                                                                                                                              <  1:64682/242‑1 (MQ=255)                                   gaTTGTTGTGGACAATGCGAAAATTGCGGAGATTGAGCAGTGCGGGCGCTCCCGCTCCCGGCGTGGGACCTATATCGATAAACAGGGGCGGCCCATGGGGGTTTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTGGCGATCGTGCCACGGCGTGTTAACCGAACGGGAAGGCAGCGCCGCGCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAgg                                                                                                                                                                                                       <  1:411186/251‑1 (MQ=255)                                    aTTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTCCCGCTTCGGGCGGGGTGCCGATATCGATAAACAGCGTCGAGCCTTTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGt                                                                                                                                                                                                      <  1:265436/251‑1 (MQ=255)                                    aTTGTTGGTGAAAATGCGCAAATGGGTGTGATTGAGCAGGGCGTGCGCTACCGCTCCCGGCGGGGTGCCGATACCGAAAAACAGCGTCGGGCCATGGGGTTTTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTGGCGAGCGGGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGt                                                                                                                                                                                                      <  1:181229/251‑1 (MQ=255)                                      tgttgGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcg                                                                                                                                                                                                    >  1:79334/1‑251 (MQ=255)                                           gTGACAATGCGCAACTTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGGGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGCTGG‑CGCAGGATCAGGTTTTGCTCCCCCAGCTCATTGAGGTCtcggcg                                                                                                                                                                                               >  2:100143/1‑251 (MQ=255)                                             gACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCCTTGAGGTCGCGGCGa                                                                                                                                                                                              >  1:360788/1‑250 (MQ=255)                                              tcAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGACCTGGAAGGCAGCGCAGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTAATTGAGGTCGCGGCGAAt                                                                                                                                                                                            >  2:108293/2‑251 (MQ=255)                                                              tgtgATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGCTTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTc                                                                                                                                                                                                                       >  2:403466/1‑208 (MQ=255)                                                              tgtgATTGAGCAGGGCGGGCGCTCCCGCTTCCGGCGTGGTGCCGATACCGATAAACAGCGTCGAGCCATTGGGGATTGGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGTCCTTGCGATCGGGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTc                                                                                                                                                                                                                       <  1:403466/208‑1 (MQ=255)                                                                     gAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGt                                                                                                                                                                                         >  1:352827/1‑231 (MQ=255)                                                                     gAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAAt                                                                                                                                                                     >  1:50554/1‑251 (MQ=255)                                                                                     cccGCTTCGGGCGTGGCGCGGTTCTCGATAACCGGCGTCTGGCCATTGGGTTTTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTc                                                                                                                                                     <  1:34254/250‑1 (MQ=255)                                                                                        gCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCATTGCGATCGTGCCACGGCGTGTTAACCGTACTGGAAGGCAGCGCCGCACAGCCATGATGA‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATATCCTGCGGGCTGACGGAGAAATGCTCTACCAGGTCTTCGGt                                                                                                                                                  >  2:9731/1‑251 (MQ=255)                                                                                                      gcctatatCGATAAACAGCGTCGAGCCATGGGGTATTTGCCCCGCCCCTTTGGGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTGGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCct                                                                                                                                    <  1:217050/247‑1 (MQ=255)                                                                                                              cGATAAACGGGGTCGGGCTATTGGGGTTTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTGGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGttt                                                                                                                             <  1:339798/250‑1 (MQ=255)                                                                                                                             aGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataat                                                                                                              >  1:77184/1‑250 (MQ=255)                                                                                                                               cctttgtggttgtgccccGCCACTTTGCGGGGGATGCGCTTTTTTTCTTCGGTCTGGGTGGCTTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataatgcc                                                                                                           <  2:79334/235‑4 (MQ=255)                                                                                                                                    ggggATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTtgt                                                                                                      >  1:346646/1‑251 (MQ=255)                                                                                                                                           tGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGAGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGCCGGAGACATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGttt                                                                                               >  2:170139/1‑250 (MQ=255)                                                                                                                                                   cgttttggggcggATGCTCTCTTTTTCTGCGGGCGGGGTGTCCTTGCGATCGGGCCACGGCGTGTTAACCGAACGGGAAGGCGGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtttgtt                                                                                       <  2:165625/239‑1 (MQ=255)                                                                                                                                                   cattgggCGGGCGATGCGCTTTTTTCCTTCGGTCGGGGTGGCCGTGCGATGGGGCCACGGCGGGTGAACCGAACTGGAGGGCAGAGCGGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtttgtt                                                                                       <  2:301747/245‑1 (MQ=255)                                                                                                                                                         gCGGGCGAGGCGCTTTTTTTTTTCGGTCTGGGTGGCCTTGCGATCGTGCCCCGGGGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCAtt                                                                                  <  1:24760/250‑1 (MQ=255)                                                                                                                                                                   cgcTCTTTTTCTTCGGTCTGGGTGGCCTTGCAATAGTGACACGGCGTGTTAACCGAACTGGAAGGCCGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCACATTGAGGCCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTACTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCt                                                                       >  1:276687/1‑251 (MQ=255)                                                                                                                                                                         ttttCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAatt                                                                 >  1:255982/1‑251 (MQ=255)                                                                                                                                                                               tcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAtttt                                                            <  1:102124/250‑1 (MQ=255)                                                                                                                                                                                   tCTGGGTGGCCTTGCGGTCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgttt                                                       >  1:162820/1‑251 (MQ=255)                                                                                                                                                                                                     cGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGCTGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATc                                     >  2:219306/1‑251 (MQ=255)                                                                                                                                                                                                            aggCGTGTTACCCGAACTGGAATGCAGCGCCGCACCGCCATGATGGAG‑CAGGATCAGTTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCTGTACTGACATAACCCTGCTGTTTAACCAGTTCTATTATCCCGTTGTGACGTTGTGTTTGTTTCATTTATAAATTCCTTGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACaaaa                              <  1:287669/250‑1 (MQ=255)                                                                                                                                                                                                                               gAAGGCAGCGCAGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAATCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTAATGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCCACAAAAGCCATCGCTAAACCCACgg           >  2:220290/1‑251 (MQ=255)                                                                                                                                                                                                                                       cgccgcACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCg   <  1:416726/251‑1 (MQ=255)                                                                                                                                                                                                                                        gccgcACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCgg  >  1:304686/1‑251 (MQ=255)                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                        TCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑C‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGG  >  NC_000913/3560206‑3560686

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