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	100.0%	154.9 / NA	42	intergenic (+589/+167)	rtcR/glpG	DNA‑binding transcriptional activator RtcR/rhomboid protease GlpG
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (21/21);  total (21/21)

ATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCG‑GTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGGCAGCGCCGCACCGCCATGATGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTC‑GGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTG  >  NC_000913/3560205‑3560695                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                  aTCTTCTTTTACCATCAACGCGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCG‑GTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGGCAGCGCCGCACCACCATGATGGc                                                                                                                                                                                                                                                    >  2:75132/1‑251 (MQ=255)     tcttTTCCCAGCAGCGTGTGACCAATGTGGAGATTGTGGGTGAAAGTGCG‑AAATTGGTTGTGTTGGAGCGGGGCGTGCTCCACCGCTCCCGGCGGGGGGCCGATTGAGATAAACGGGTTCGAGCCTTGGGGGTTTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCG‑GTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGGCAGCGCCGCACCGCCATGATGGcg                                                                                                                                                                                                                                                   <  2:332649/248‑2 (MQ=255)        tttACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGGCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTGTTCG‑GTCTGGGGGGCCTTGCGATCGTGCCCCGGCGTGTTAACCGAACTGG‑AAGGCAGCGCCGCACCGCCATGATGGCGCAGGa                                                                                                                                                                                                                                              >  1:43776/1‑250 (MQ=255)                    gtgtTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATAG‑CTGTGACTGAGCAGAGCGTGCGCTACCGCTTACGGAGTGGTGACGATATCGATAAACAGTGTCGAGACATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCG‑GTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTc                                                                                                                                                                                                                                 >  1:337216/1‑251 (MQ=255)                            aaCGTTGAGATTGTTGGTGACAAGGCGTAATTGG‑CTGTGATGGAGCAGTGCGTGCGCTCCCGCGTCCGGGGGGGGGCTGATTTCGATAAACAGGGTGGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCG‑GTCTGGGGGGCCGTGCGATCGTGCCACGGCGTGTTAACCTAACTGG‑AATGCAGCGCGGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCt                                                                                                                                                                                                                         <  1:457412/251‑1 (MQ=255)                               gTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATCTCGATAAACAGCGTCGAGCCCTTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCG‑GTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCCGCTCAt                                                                                                                                                                                                                      >  2:187272/1‑251 (MQ=255)                                          tgGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCG‑GTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcgg                                                                                                                                                                                                           >  1:3665/1‑251 (MQ=255)                                            ggtgaaaaggcgaaaattcttgggATTGAGCAGGGCGGGCGCTCCCGCTCCCGGCGTGGTGCCGATACCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCG‑GTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggc                                                                                                                                                                                                          <  1:158305/228‑1 (MQ=255)                                                    gcgaaaaTTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGGGTGGGGCCGATATCGATAAGCAGCGTCGAGCAATGGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCG‑GTCTGGGTGGCCTTGCGATCGGGCCACGGCGTGTTAACCGAACTGG‑AAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCt                                                                                                                                                                                                 <  1:348465/247‑1 (MQ=255)                                                     cgcAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCG‑GTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTg                                                                                                                                                                                                >  1:73928/1‑251 (MQ=255)                                                                   gATTGAGCGGGGCGTGCGCTGCCGCTTCGGGGGGGGGGCCGATTGCGAAAAACAGCGGCGGGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGAGGCGCTCTTTTTCTTCG‑GTCTGGGGGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACGGG‑AAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGgag                                                                                                                                                                                   <  1:399956/251‑1 (MQ=255)                                                                    attgagcagtgagacagcTACCGCTTCCGGTGTGGTGACGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCG‑GTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGGCAGCGCCGCACCGCCATGCTGGCGCAGGATCCGAGTTTGCTCCGCCCGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACgg                                                                                                                                                                                     >  2:69161/17‑248 (MQ=255)                                                                        aGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCG‑GTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGtttt                                                                                                                                                                                                                                     >  1:277343/1‑196 (MQ=255)                                                                               cgtgcgCTACCGCTTCCGGCGTGGTGCCGATCTCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCG‑GTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGATCCGCCAGCTCATTGCGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAcc                                                                                                                                                                       >  2:80961/1‑251 (MQ=255)                                                                                              ccGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCG‑GTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGATCTTC‑GGTACTg                                                                                                                                                       >  1:253393/1‑251 (MQ=255)                                                                                                                    aaaCAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCG‑GTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTC‑GGTACTGACATAACCCTGCTGTTTAAcc                                                                                                                                  >  1:196179/1‑250 (MQ=255)                                                                                                                      actgctcggcgcatttGGGGTTTTGCCCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCG‑GTCTGGGTGGCCTTGCGATCGTCCCACGGCGTGTTAACCGAACTGA‑AAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTC‑GGTACTGACATAACCCTGCTGTTTAACCAGt                                                                                                                               <  1:75132/237‑1 (MQ=255)                                                                                                                      aCAGCGTAGAGCCATTGGGGATTTGCGTCGCCACTTTGCGGGCGTTGCGCTCTTTTTCTTCG‑GTCTGGGTGGCCTTGCGATCGTGCCGCGGCGTGTTAACTGAACTGG‑CAGGCGGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGCATAGTCTGAAGGGTGACGGAGAATTGCTATACCAGCTCGTC‑GGTTCTGACATaaacttgatggttgcacagt                                                                                                                               >  2:284142/1‑233 (MQ=255)                                                                                                                               gtgtgtgggggttcgcccccgcccttttGCGGTGGATGCCCTCTTTTTCTTCGTGT‑TGGGTGGCCTTGCATTCGCGCCACGGCGTGTTAACCGAACGGG‑AAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTTCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACTGAGAAATGCTCTACCAGCTCTTC‑GGTACTGACATAACCCTGCTGTGTAACCAGGTCGataata                                                                                                                      <  2:277622/226‑1 (MQ=255)                                                                                                                                   gtggggatgtggcaccGCCACTTTGCGGGGGTTGCGCTCTTTTTCTGCG‑TGCTGGGGGGCCTTGCGCTCGTGCCACGGCGTGTTAACCGAACTGG‑AAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTC‑GGTACTGACATAACCCTTCTGTTTAACCAGTTCGATAATAcc                                                                                                                    <  2:173085/235‑1 (MQ=255)                                                                                                                                    ttGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCG‑GTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGCACTGG‑AAGGCAGCGCCGCACCGCCATGATGGCGCCGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTC‑GGTACTGACATAACCCTGCTGTTTAACCCGCTCGATAATACCGtt                                                                                                                 >  2:310866/1‑251 (MQ=255)                                                                                                                                     ggggattttcccccGCCACTTGGCGGGCGTGGCGCTCTTTTTCTTCG‑GTCTGGGTGGCCTTGCGATCGGGCCCGGGCGTGTTACCCGAACTGG‑AAGGCGGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTC‑GGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTtg                                                                                                                <  2:120510/239‑1 (MQ=255)                                                                                                                                         gATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCG‑GTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTC‑GGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGAc                                                                                                            <  1:187272/251‑1 (MQ=255)                                                                                                                                          tttcgtaccGCCACTTGGCGGGCGGTGCGACCTTGTTCTTCG‑GTCTGGTTGGCCTTGATGTCGTGCCCTGGGTTCGTAACAGAA‑TGGTAAGGCAGGGCGGCGCCGCCATGATGGCGCAGGATCAGGTTTTGCGCCGCCTGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCGGACGGAGAAATTCTCTACCAGCTCTTCTGGT‑TTGATAAAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGAc                                                                                                            <  1:419395/243‑1 (MQ=255)                                                                                                                                            ttGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCG‑GTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGAGGGCTGACGGAGAAATGCTCTACCAGCTCTTC‑GGTACTGACATAACCCTGCTGTTTAACCCGTTCGCTAATACCGTTGTGACGtt                                                                                                         >  2:96279/1‑251 (MQ=255)                                                                                                                                               cgccgccgtgtggggggCGATGCGCTCTTTTTCTTCG‑GTCTGGGTGGCCTGGCGCTCGTGCCACGGCGTGTTAACCGAACTGG‑AAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTC‑GGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtgt                                                                                                       <  2:75869/236‑1 (MQ=255)                                                                                                                                                               ggcgtggtctctTTTTCTTCG‑GTCTGGGTGGCCTTGCGACCGGCCCACGGCGTGTTAACCAACCTGG‑AAGGCAGCGCCGCCCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTC‑GGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATaa                                                                                      <  2:93166/243‑1 (MQ=255)                                                                                                                                                               gCGATGCGCTCTTTTTCTTCG‑GTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTC‑GGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTtata                                                                                       >  1:194621/1‑250 (MQ=255)                                                                                                                                                                   ggcgcTCTTTTTCTCCG‑GTTGGGGGGGCCTGGCGCTCGTCCCACGGCGTGTTAACCGAACGGG‑AAGGCACCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTC‑GGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATcc                                                                                  <  2:73928/250‑1 (MQ=255)                                                                                                                                                                     cgcTCTTTTTCTTCG‑GGCTGGGTGGCCTTGCGATCGTGCCACGGCGGTTTACCCGACCGGG‑AGGGCGGCGCCGCTCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCTGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTC‑GGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCt                                                                                <  1:45997/251‑1 (MQ=255)                                                                                                                                                                     cgcTCTTTTTCTGCG‑GTCGGGGTCGCCTTGCGCTGGTGCCCCGGCGTGTTAACCGAACTGG‑AAGGCAGCGCCCCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCTGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTC‑GGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCt                                                                                <  1:284142/251‑1 (MQ=255)                                                                                                                                                                                  cG‑GTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTC‑GGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcg                                                                   >  1:44317/1‑251 (MQ=255)                                                                                                                                                                                          ggTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGAGGGCTGACGGAGAAATGCTCTACCAGCTCTTC‑GGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTcgc                                                            <  1:96279/251‑1 (MQ=255)                                                                                                                                                                                              gCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTC‑GGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCAt                                                        <  1:105101/251‑1 (MQ=255)                                                                                                                                                                                               ccTTGCGATCGTGCCACGGCGTGTTAACCGAACGGG‑AAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTC‑GGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCAtt                                                       <  1:45156/251‑1 (MQ=255)                                                                                                                                                                                                           gCCACGGCGTGTTAACCGAACTGG‑AAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTC‑GGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCCTTGAGCGAATcaac                                           >  2:9257/1‑251 (MQ=255)                                                                                                                                                                                                                    tgtTAACCGAACTGG‑AAGGCAGCGCCGCACCGCCATGCTGGCGCAGGATCAGGTTTTGCTCCGCCAGCCCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTC‑GGTACTGACATAACCCTGCTGTTTCACCCGTTCGATAATACCGTTGTGACGTTGTGTTTGTCTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGACTCCACAAAAGCCAt                                  >  2:460497/1‑251 (MQ=255)                                                                                                                                                                                                                       gaaCCGAGCGGG‑AAGGCCGCCCCGCCCAGCCAGGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATGGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTC‑GGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCg                                <  2:253393/249‑1 (MQ=255)                                                                                                                                                                                                                                    aagtcgggcccgcACCGCCATGAGGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTGC‑GGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGGTTCATTTATAAATCCCTGGAAGGATTTTCGTTTTCGCGCATGGAGCGAAGCAACAAAAGCCATCGCGAATCCCTCGGc                   <  2:53656/242‑1 (MQ=255)                                                                                                                                                                                                                                      ggCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCCGCTAATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTC‑GGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCt                  >  1:295897/1‑250 (MQ=255)                                                                                                                                                                                                                                        cAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGAAAGCTCATTGAGGACGCGGCGAGTAGTCTGCGGGATGACGGAGAACTGCACTACCAGCACTTC‑GGTACTGACGTAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGCCTTATTTTCGTTTTCGCGCATTGAGCGAATCACCAAAAGCGATCGCGAAACACACGGCTaa                >  1:409440/1‑250 (MQ=255)                                                                                                                                                                                                                                           cgccgcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTC‑GGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCg            <  1:18708/251‑1 (MQ=255)                                                                                                                                                                                                                                                  ccGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTC‑GGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGAtgt     >  1:216068/1‑251 (MQ=255)                                                                                                                                                                                                                                                     ccATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTC‑GGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGAtgtgtg  >  1:465173/1‑251 (MQ=255)                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                  ATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCG‑GTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGGCAGCGCCGCACCGCCATGATGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTC‑GGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTG  >  NC_000913/3560205‑3560695

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