Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,560,443	C→G	intergenic (+577/+179)	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,443	0	C	G	100.0%	120.7 / NA	37	intergenic (+577/+179)	rtcR/glpG	DNA‑binding transcriptional activator RtcR/rhomboid protease GlpG
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (10/27);  total (10/27)

TGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTA  >  NC_000913/3560198‑3560680                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               tgCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGgcagcgccgcagcgcc                                                                                                                                                                                                                                            >  2:98301/1‑250 (MQ=255)                 tACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATc                                                                                                                                                                                                                            <  1:131321/250‑1 (MQ=255)                 tACCATCAACGTGTTAGCAACGTTGAGATGGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCATTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATc                                                                                                                                                                                                                            <  2:31546/250‑1 (MQ=255)                   ccATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAgg                                                                                                                                                                                                                         >  2:243482/1‑251 (MQ=255)                    cATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGt                                                                                                                                                                                                                        >  1:276407/1‑251 (MQ=255)                     aTCAACGTGTTAGCAACTTTGGGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGtt                                                                                                                                                                                                                       <  2:84537/251‑1 (MQ=255)                        aaCGTTTTAGCATCGTTGAGATTGTTGGTGACAATGCACAAATTGCTGTGATTGAGCGGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCTATAGCGATAAGCAGCGGCGAGCCTTTGGGGATTTTCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACGGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTg                                                                                                                                                                                                                    <  2:259796/251‑1 (MQ=255)                        aaCGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTg                                                                                                                                                                                                                    <  1:158888/251‑1 (MQ=255)                            ggtTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATGGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTcc                                                                                                                                                                                                                <  2:92204/250‑1 (MQ=255)                                aGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCa                                                                                                                                                                                                            <  1:98301/251‑1 (MQ=255)                                   aaCGTTGAGATTGTTGGTGCCAAGGCGCAATTTGCTGTGTGGGAGGATGGCGTGCGCTACCGCTTCCGGCGTGGTGCCGTTATCGATAAACAGCGTCGGGCCATTGGGGATTTTCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCt                                                                                                                                                                                                         <  2:230567/251‑1 (MQ=255)                                         gagaTTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTg                                                                                                                                                                                                    >  1:201321/1‑250 (MQ=255)                                               gtgggTGACAATGCGCAAATTGGTGTGATTGAGCAGTGCGTGCTCTCCCGCTTCCGGCGGGTTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcg                                                                                                                                                                                              <  2:64463/247‑1 (MQ=255)                                                 tggggaaaaTGCGCAAATTGCTGTGATTGAGCAGGGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcgcg                                                                                                                                                                                            <  2:201321/243‑1 (MQ=255)                                                                       gggagtgaggggtgcgtgcGCTACCGATTCAGGAGTGGGGCTGATTTCGATAAACAGCGTCGAGCCAGTGGGGATGTGGTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCAGTGAGGTCGCGGCGAATAGTCTGCGGGCTGACgg                                                                                                                                                                     <  2:201839/240‑1 (MQ=255)                                                                        gggtggAGCAGGGCGTGCGCTACCGCTGCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTGGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTATTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGga                                                                                                                                                                    <  2:225185/246‑1 (MQ=255)                                                                                                          gggtgacgatttcgcaaaaCAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTGGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATaa                                                                                                                                  <  2:7672/235‑1 (MQ=255)                                                                                                             ggCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAccc                                                                                                                               <  1:103798/250‑1 (MQ=255)                                                                                                                    atCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGtt                                                                                                                        <  1:184516/251‑1 (MQ=255)                                                                                                                              ggCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTc                                                                                                              <  2:41264/250‑1 (MQ=255)                                                                                                                                 gTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAAt                                                                                                                                                                                      >  1:262317/1‑176 (MQ=255)                                                                                                                                 gTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAAt                                                                                                                                                                                      <  2:262317/176‑1 (MQ=255)                                                                                                                                       cctttGGGGATTTGCGCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAcc                                                                                                     <  2:272362/248‑1 (MQ=255)                                                                                                                                         aTTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGt                                                                                                   >  2:111634/1‑251 (MQ=255)                                                                                                                                               gATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGAc                                                                                             >  1:132691/1‑251 (MQ=255)                                                                                                                                                        cGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgttt                                                                                    >  1:255326/1‑251 (MQ=255)                                                                                                                                                        cGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgttt                                                                                    <  2:276407/251‑1 (MQ=255)                                                                                                                                                             cTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTc                                                                               >  1:40833/1‑251 (MQ=255)                                                                                                                                                                  cGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGAGGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTta                                                                          <  1:232989/251‑1 (MQ=255)                                                                                                                                                                   gggCGATGCTCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGACATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCAACCAGCTCTGCGGTACTGACATAACCCTGCTGTTTAACCAGTACGATAATACCGTTGTGACGTTGTGTTTGTTTCATTtat                                                                         <  1:186586/251‑1 (MQ=255)                                                                                                                                                                                              gggTGGTCTTGCGATC‑TGTCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGGCTCGGCGAATAGTCGGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATATCCCTGCTGTTTAACCAGCTCGATAATACCGTTGTGACGGTGTGTTTGTTTCATTTATAACTCCCTGGAATTAttttccttttc                                               <  2:232117/250‑1 (MQ=255)                                                                                                                                                                                              gggTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgttttcg                                              <  1:228260/251‑1 (MQ=255)                                                                                                                                                                                                gTGGCCTTGCGATCGTG‑CCACGGCGTTTTAACCGAACTGGAAGGCAGCGCCGCCGCGCCATGATGGCGCAGGATCAGGTTTTGCTCAGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTa                                                                                                                      >  2:129632/1‑177 (MQ=255)                                                                                                                                                                                                gTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTGTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTa                                                                                                                      <  1:129632/177‑1 (MQ=255)                                                                                                                                                                                                       tgagagcGTG‑ACACGGCGTGTTAACCGAACTGTAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTGGCTCTGCCAGCTCATGGAGTTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCATCTCGTCGGTACTGACATAACTCTGCTGTTTAACCAGTTCGATAATACCGTTCTGACGTTGTGTTTGTTTCACTTATAGTTCCCTGGAATGATGTTCGTGTGCGCGCATTGa                                      <  2:133664/244‑1 (MQ=255)                                                                                                                                                                                                              gTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATc                               >  1:240846/1‑250 (MQ=255)                                                                                                                                                                                                                     cGGCGTGGTAACCGAACTGGAAGGCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACaaaa                        <  1:152993/251‑1 (MQ=255)                                                                                                                                                                                                                                           ggCAGCGCCGCAGCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTa  <  1:243482/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               TGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTG‑CCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTA  >  NC_000913/3560198‑3560680

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