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%	145.2 / NA	40	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 (19/21);  total (19/21)

TTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGC‑GCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGG‑ATTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTG‑CGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCG‑GTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGT  >  NC_000913/3560208‑3560694                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                 ttcttTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGC‑GCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGG‑ATTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTG‑CGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGcgc                                                                                                                                                                                                                                                 >  1:82639/1‑248 (MQ=255)  tcttTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGC‑GCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGATCCATTGGGG‑ATTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTG‑CGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGcgc                                                                                                                                                                                                                                                 >  1:529135/1‑247 (MQ=255)    atttACCATCAACGTGTTCGCAACGGTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGC‑GCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGCCGAGCCATTGGGG‑ATTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTG‑CGATCGTGCCACGGCGTGTTACCCGAACTGGAAGGCAGCGCCCAACCGCACTGATGGcgca                                                                                                                                                                                                                                                >  2:85500/2‑245 (MQ=255)       tACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGC‑GCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGG‑ATTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTG‑CGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATc                                                                                                                                                                                                                                           >  1:209729/1‑250 (MQ=255)           aTCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGC‑GCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGG‑ATTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTG‑CGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCACGATGGCGCAGGATCAGGtt                                                                                                                                                                                                                                      >  1:475613/1‑251 (MQ=255)             cAACGTTTTAGCACCGTGGGGTTTGTTGGTGACATGGCGAAATTTG‑TGGTGATGGAGCAGTGCGTGC‑CCTCCCGTTCCCGGGGGGGGCCCGATATCGAAAAACAGGGCCGAGCCATGGGGG‑TTTT‑GCTCCGCCACTTGGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTG‑CGATCGTGCCAGGGCGTGTTAACCTAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGtttt                                                                                                                                                                                                                                    <  1:195725/251‑1 (MQ=255)                cGTGTTAGAAATGTGGAGATTTTGGGTGACAAT‑CACAAATTGTCTGTGATGGAGCAGTG‑GTGCCGCTCCCCTTTCCGGTGTGGTGCGGATATCGAAAAACGGCGGCGCGCCAT‑GGGG‑ATTTTGCTCCGCCACTTGGCGGGCGATGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTG‑CGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCt                                                                                                                                                                                                                                 <  2:278135/251‑1 (MQ=255)                             ttGAGAGTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGC‑GCAACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGG‑CTTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTG‑TGATCGTGCCACGGCGTGTTATCCGAACTGGAAGGCAGCGCCGCCCAGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCAtt                                                                                                                                                                                                                    >  2:529115/1‑251 (MQ=255)                               gagaTTGTTGGTGACAATGCGCAAATTG‑CTGTGCTTGAGCAGTGCGTGC‑GCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGG‑ATTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTG‑CGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGAAGCGCCGCACCGCCATGAGGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCCTTgc                                                                                                                                                                                                                  >  2:58939/1‑250 (MQ=255)                               gagaTTGTGGGTGACAAGGCGAAAATTG‑CTGTGATTGAGCAGTGCGGGC‑GCTACCGCTTGCGGCGTGGGGCCGATACCGAAAAACGGCGCCGGGCCATTGGGG‑TGTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGTTCTGGGTGGCCTTG‑CGATGGTGCCACGGCGTGTTACCCGAATTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGa                                                                                                                                                                                                                  <  1:198948/251‑1 (MQ=255)                                agaTTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGC‑GCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGG‑ATTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTG‑CGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTg                                                                                                                                                                                                                   >  2:537449/1‑249 (MQ=255)                                   ttgttgGTGACAATGCGCAAATTG‑CTGTGATTGAGCATTGCGTGC‑GCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGG‑CTTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTG‑CGATCGTGCCACGGCGTGTTAACCGCACTGGAAGGCAGCGCCGCACCGCCCTGATGGctcaggc                                                                                                                                                                                                                                             >  2:118363/1‑214 (MQ=255)                                        ggTGACAACGCGCAAATTT‑CTGTGTTTGAGCAGTGCGTGC‑CCTCCCGCTTCCGGCGTGGTGCCATTATCGATAAACAGCGTCGGGCCTTGGGGG‑TTTT‑GCTCCGCCACTTTGCGGGCGATGCGCGTTTTTTCGTCGGTCTGGGTGGCCTGG‑CGATCGTGCCACGGCGTGTTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggc                                                                                                                                                                                                         <  1:498632/251‑1 (MQ=255)                                         gTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGC‑GCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGG‑ATTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTG‑CGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggcg                                                                                                                                                                                                        >  1:350394/1‑251 (MQ=255)                                          ggACAATGCACAAATTG‑CTGTATTTGAGCAGTGCGGGC‑GCTACGGCTGCGGGGGTGTGGCCGATATCGAAAAACAGCGTCGGGCC‑TTGGGGTATTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTG‑CGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGa                                                                                                                                                                                                       <  1:130463/250‑1 (MQ=255)                                                       atttctttgtatggAGCAGTGCGGGC‑GCTCCCGTTTCCGGGGGGTTGCCGATACCGATAAACGGGTTCGAGCCTTGGGGG‑TTTT‑GTCCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTG‑CGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCggg                                                                                                                                                                                           <  1:531715/238‑1 (MQ=255)                                                                                  gcTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGG‑CTTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTG‑CGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAACGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCCTTGAGGTCGCGGCGAATACACTGCGGGCTGCCGGAGAACTGCTCTACCAGctc                                                                                                                                                                 >  2:61230/1‑251 (MQ=255)                                                                                  gcTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGG‑ATTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTG‑CGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGctc                                                                                                                                                                 >  2:398364/1‑251 (MQ=255)                                                                                     aCCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGG‑ATTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTG‑CGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTc                                                                                                                                                              >  1:44515/1‑251 (MQ=255)                                                                                      cgctttcgggCGTGGTGCCAATATCGAAAAACGGCGTGGGGCCTTGGGGG‑TTTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTG‑CGATCGTGCCGCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCg                                                                                                                                                             <  1:508847/243‑1 (MQ=255)                                                                                                              cGATAAACAGCGTCGAGCCATTGGGG‑ATTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTG‑CGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCG‑GTACTGACATAACCCTGCTGttt                                                                                                                                     >  1:473353/1‑250 (MQ=255)                                                                                                                              gctatggggg‑TTTT‑CCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTG‑CGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCG‑GTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAc                                                                                                                    <  1:162359/245‑1 (MQ=255)                                                                                                                              gCCATGGGGG‑GTTT‑GCTCCGCCTCTTGGCGGGAGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTGG‑CGATTGTGCCCCGGGGGGTGAACCGAACTGGAAGGCAGCGCCGAACCGCCATGATGGCGTAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGTTGACGGAGAAATGCTCTACCAGCTCTTCG‑GTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAc                                                                                                                    <  1:336885/251‑1 (MQ=255)                                                                                                                               ggcatgtggggttggtccccGCCCCTTTGGGGGGGAGGCGCTCTTTTTCTGCGGTGTGGGTGCCCTGG‑CGATCGGGCCACGGCGTGTTAACCGAGCTGGAGGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCG‑GTACTGACATAACCCTGCTGTTTAACCAGTTCGataata                                                                                                                     <  2:44515/233‑1 (MQ=255)                                                                                                                                   tGGGG‑TTTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTG‑CGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCG‑GTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGtt                                                                                                                <  1:13572/250‑1 (MQ=255)                                                                                                                                        gtttt‑GCCCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTG‑CGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCG‑GTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGAc                                                                                                           <  1:208310/249‑1 (MQ=255)                                                                                                                                               cTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGCGGCCTTG‑CGAACGTGCCACGGCGTGTTAACCGAACTGGACGGCAGCGCCGCACCGCCACGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGCGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAACTGCTCTACCAGCTCTTCG‑GTACTGACATAACCCTGCTGTTTNACCAGCTCGCTAATCCCGTTGTGACGTtgtg                                                                                                     >  2:271508/1‑251 (MQ=255)                                                                                                                                                  cGCCACTTGGGGGGCGTTGCGCTCTTTTTCTTCGGTGTGGGTGGCCTTG‑CGTTCGTGCCACGGCGGGTTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCG‑GTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgttt                                                                                                  <  2:126703/251‑1 (MQ=255)                                                                                                                                                           gCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTG‑CGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGCTCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCG‑GTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCCttt                                                                                         >  1:204094/1‑251 (MQ=255)                                                                                                                                                            cGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTG‑CGATCGTGCCACGGCGTGTTAACCGAACTGGACGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGACATGCTCTACCAGCTCTTCG‑GTACTGACATAACCCTGCTGTTTAACCCGTTCGATAATACCGATGTGACGTTGGGTTTGTTTCAtttt                                                                                        >  2:147130/1‑250 (MQ=255)                                                                                                                                                            cGGGCGATGCGCTCTTTTCCTTGGGTCTGGGTGGC‑TTGTAGAGCGGGCCACGGGGTTTTAACCAATCTGAAAGGGAGGGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCGCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCG‑GTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTta                                                                                        <  1:520368/251‑1 (MQ=255)                                                                                                                                                               gCGATGCGCTTTTTTTCTGCGGTCGGGGGGGCCGTG‑GGATCGTGCCACGGTGTGTTAACCGAACTGGAGGGCGGCGCCGCACCGCCATGATGGCGCAGGAGCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCG‑GTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATaa                                                                                     <  1:224496/251‑1 (MQ=255)                                                                                                                                                                 gATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTG‑CGACCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCGCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCG‑GTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAAt                                                                                   <  1:343031/251‑1 (MQ=255)                                                                                                                                                                   tGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTG‑CGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCG‑GTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGTCGTTGTGTTTGTTTCATTTATAAATcc                                                                                 >  1:337332/1‑251 (MQ=255)                                                                                                                                                                              tatacgggcTGGGTAGCCTTG‑CGATCGTGCCACGGCGTGTTAACCGGCCTTGACTGCCGAGCCGAACCGCTATGTTTGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATGGAGGACCCGGCGAATACTCCGACGGCTGAAAGAGAAAGTCTCAACCAGCTCTTCGCTTAC‑GCCGTCACCATGCTGTTTTACCCGTGCGCTAAGACCACTGTGACGATGAGTTTGTTTTATCTAACAATCCCTGGAattatt                                                                      >  2:334467/9‑251 (MQ=255)                                                                                                                                                                                   ggTCTGGGTGGCCTTG‑CGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCG‑GTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgt                                                                 >  1:256344/1‑251 (MQ=255)                                                                                                                                                                                        gggTGTCCTTG‑AGATCGTGCCCCGGGGTGTTAACCGACCTGGAGGGCAGCGCGGCACCGCCATTATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCG‑GTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTGCATTTATAAATCCCTGGAATTAttttcgttt                                                               <  2:54501/248‑1 (MQ=255)                                                                                                                                                                                             gCCTGG‑CGATCGTGCCACGGGGTGTTAACCGAACGGGAATGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCG‑GTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCAt                                                       <  1:81191/251‑1 (MQ=255)                                                                                                                                                                                               cTTG‑CGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAAAGCTCTACCAGCTCTTCG‑GTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGGGTTTGTTTCATTTCTAAATCCCTGGAAATAATTTCGTTTTCGCGCaatg                                                     >  2:498471/1‑248 (MQ=255)                                                                                                                                                                                                     gATCGTGCCACGGCGTGTTAACCGAACTGAAAGGCAGCGCCGCACCGCCATGATGGCGCCGGATCAGGCTTTGCTCAGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGCCGGAGAAATGCTCTACCAGCTCTTCG‑GTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATACATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCg                                                 >  2:541396/1‑250 (MQ=255)                                                                                                                                                                                                                           cGAACTGGAAGGCGGCGCCGCACCGCCACGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCG‑GTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTCTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTaaa                          <  2:475613/251‑1 (MQ=255)                                                                                                                                                                                                                                            ccgcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCG‑GTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGc         <  1:398364/251‑1 (MQ=255)                                                                                                                                                                                                                                                   gCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCG‑GTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGAtgtgt  >  1:114355/1‑251 (MQ=255)                                                                                                                                                                                                                                                   gCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGGCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCG‑GTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGAtgtgt  <  2:469478/251‑1 (MQ=255)                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                 TTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGC‑GCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGG‑ATTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTG‑CGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCG‑GTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGT  >  NC_000913/3560208‑3560694

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: ‑