Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,078,954	A→G	intergenic (+23/+78)	yoeF → / ← tmaR	putative uncharacterized protein YoeF/PTS system regulator TmaR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,078,954	0	A	G	100.0%	162.8 / NA	41	intergenic (+23/+78)	yoeF/tmaR	putative uncharacterized protein YoeF/PTS system regulator TmaR
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/19);  total (22/19)

GGTAAAAACGCCCTACCGGGCCATGAATCACCCGGTCATTGGAGTAGTAACCAAAGCAGATTTAGCCAGCATGGAACAAATATCTTTGGTGAAAAGTTGGTTACGGGAAGCTGGAGCGCACAACGTGTTAGTAACCAGTGCAGTTAATAACAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGTCACCTTTCATGCTGGCGATGATGCCCTGGATTGCTTCAACGCTC  >  NC_000913/2078727‑2079194                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 ccgacaaaatgcgttattgacagcggtttctgtggattgacacatttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGAt                                                                                                                                                                                                                            >  1:446219‑M2/204‑250 (MQ=255)                   gacagcggtttctgtggattgacacatttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATcaatgagcaatg                                                                                                                                                                                                          <  2:169221‑M2/65‑1 (MQ=255)                            ttctgtggattgacacatttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGcc                                                                                                                                                                                                 <  2:251910‑M2/74‑1 (MQ=255)                                  ggattgacacatttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCaa                                                                                                                                                                                           <  1:476995‑M2/80‑1 (MQ=255)                                           catttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGa                                                                                                                                                                                 <  1:51064‑M2/90‑1 (MQ=255)                                                  ttcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAAc                                                                                                                                                                          >  1:160445‑M2/155‑251 (MQ=255)                                                   tcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACg                                                                                                                                                                         >  2:16618‑M2/154‑251 (MQ=255)                                                                              gaacaaccagcaaacggtcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaataacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACCCAGATAAGAATCAATGAGCATTGCGTGTTGCCTGGCAACCCTCATTGAATAAAACTGGAATTACGTCGCTTCGCCGGTTTTCa                                                                                                                                              <  2:139141‑M2/125‑1 (MQ=255)                                                                                         aaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCAt                                                                                                                                    >  2:386005‑M2/116‑250 (MQ=255)                                                                                         aaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATa                                                                                                                                   <  1:331773‑M2/136‑1 (MQ=255)                                                                                                cattgtgttctgaaatcaaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTc                                                                                                                             >  1:128831‑M2/109‑250 (MQ=255)                                                                                                                  aggctattgagtcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAAGGAGCAATGAGGGTTGCCGGGCAATCCTCATTGAATAAAACGGGAATTTCTTCGCTTCGCCGGTTTTCATTTTGCCCATTGCTTTCAGCTTTGGGGAGAGATcgc                                                                                                          <  2:156734‑M2/161‑1 (MQ=255)                                                                                                                       attgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGt                                                                                                     >  1:30747‑M2/86‑251 (MQ=255)                                                                                                                            ttcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCttt                                                                                                >  1:110580‑M2/81‑251 (MQ=255)                                                                                                                               aacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTgg                                                                                              >  2:427373‑M2/78‑250 (MQ=255)                                                                                                                                   gcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATAACGCGGCGTTCTTTGGAGAGc                                                                                         <  2:271963‑M2/178‑1 (MQ=255)                                                                                                                                         acagagctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGGTGCCGGGCAACCCTCAGTGAATAAAACGGGAATTTCTTCGCTTCGCCGTTTTTCATTTCGCCTATAGTTTTCAGCTTTTGGGAGATATCGCTGCGGTCTTTGGAGAGCTCGGCa                                                                                   <  2:138627‑M2/184‑1 (MQ=255)                                                                                                                                            gatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCAttt                                                                                >  1:353034‑M2/65‑251 (MQ=255)                                                                                                                                                 gcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTtga                                                                            <  1:443693‑M2/191‑1 (MQ=255)                                                                                                                                                      gtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTGTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTGgatgat                                                                        <  2:454333‑M2/195‑1 (MQ=255)                                                                                                                                                                 tgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCaa                                                             <  2:491439‑M2/206‑1 (MQ=255)                                                                                                                                                                    aaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAAcgcg                                                         >  2:403916‑M2/41‑250 (MQ=255)                                                                                                                                                                         aacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTc                                                   >  2:42612‑M2/36‑251 (MQ=255)                                                                                                                                                                             taatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCAt                                                 >  1:182155‑M2/32‑249 (MQ=255)                                                                                                                                                                              aatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGt                                              >  1:505934‑M2/31‑251 (MQ=255)                                                                                                                                                                               atagaattatccggaatggtggcgactatGATTTCGATTGGTTATAGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGGTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTAGACGTTTTTCATTTCTCACATAGCTTTCAGCTTTTTGGATATATCGCGGCGTTCTTTGGAGGGCTCAGCATTTTTGATGATGTGATAGTCAACGCGATCTTCATAGTc                                             >  2:258457‑M2/30‑251 (MQ=255)                                                                                                                                                                                   aattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAATCGGGAATGACTTCGCGGCGCCGTTGTTCATTTCGCCCATAGCTTGCATCTTTTTGGAGATATCGCGGCGTTCTTGTGAGAGCTCGGCATTGTTGAGGAGGTAAGCGTCTACGCGATCTTCATAGTCaccg                                         <  2:403376‑M2/226‑2 (MQ=255)                                                                                                                                                                                     ttatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGTCACCttt                                       >  1:65339‑M2/24‑251 (MQ=255)                                                                                                                                                                                     ttatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGTCACCttt                                       >  1:482968‑M2/24‑251 (MQ=255)                                                                                                                                                                                     ttatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGTCACCttt                                       >  1:434904‑M2/24‑251 (MQ=255)                                                                                                                                                                                      tatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGTCACCttt                                       <  1:342976‑M2/228‑1 (MQ=255)                                                                                                                                                                                             aatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGTCACCTTTCATGCTgg                               >  2:420175‑M2/16‑251 (MQ=255)                                                                                                                                                                                             aatggaggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGTCACCTTTCATGCTgg                               >  1:512448‑M2/16‑251 (MQ=255)                                                                                                                                                                                                       actatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTTGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGTCACCTTTCATGCTGGCGATGATGcc                     <  2:471807‑M2/246‑1 (MQ=255)                                                                                                                                                                                                       actatAAGTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGTCACCTTTCATGCTGGCGATGATGcc                     <  2:255772‑M2/246‑1 (MQ=255)                                                                                                                                                                                                               ttCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGTCACCTTTCATGCTGGCGATGATGCCCTGGATTg             >  1:251271/1‑251 (MQ=255)                                                                                                                                                                                                                 cGATTTGTTCTCGGATGATGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACTGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGTCACCTTTCATGCTGGCGATGATGCCCTGGATTGCt           <  2:353034/251‑1 (MQ=255)                                                                                                                                                                                                                   aTTGGTTCGCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGTCACCTTTCATGCTGGCGATTATGCCCTGGATTGCTTc         <  1:420175/251‑1 (MQ=255)                                                                                                                                                                                                                      ggTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCa                                                                                                                                              >  2:436280/1‑115 (MQ=255)                                                                                                                                                                                                                      ggTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCa                                                                                                                                              <  1:436280/115‑1 (MQ=255)                                                                                                                                                                                                                           tcGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGTCACCTTTCATGCTGGCGATGATGCCCTGGATTGCTTCAACGCTc  >  2:23920/1‑250 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 GGTAAAAACGCCCTACCGGGCCATGAATCACCCGGTCATTGGAGTAGTAACCAAAGCAGATTTAGCCAGCATGGAACAAATATCTTTGGTGAAAAGTTGGTTACGGGAAGCTGGAGCGCACAACGTGTTAGTAACCAGTGCAGTTAATAACAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGTCACCTTTCATGCTGGCGATGATGCCCTGGATTGCTTCAACGCTC  >  NC_000913/2078727‑2079194

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