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%	232.5 / NA	64	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 (33/31);  total (33/31)

TTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATT‑GAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGA‑TATCG‑ATAAACAGCGTCGAGCCATTGGGGATTTGCT‑CCGCCACTTTGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGT  >  NC_000913/3560208‑3560694                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                ttcttTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATT‑GAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGA‑TATCG‑ATAAACAGCGTCGAGCCATTGGGGATTTGCT‑CCGCCACTTTGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGcgc                                                                                                                                                                                                                                                >  1:255063/1‑248 (MQ=255)        aCCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATT‑GAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGA‑TATCG‑ATAAACAGCGTCGAGCCATTGGGGATTTGCT‑CCGCCACTTTGCGGGCGATGCGC‑TCTTTTTCTTCGGCCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATc                                                                                                                                                                                                                                          >  2:360934/1‑249 (MQ=255)            tCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATT‑GAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGA‑TATCG‑ATAAACAGCGTCGAGCCATTGGGGATTTGCT‑CCGCCACTTTGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGttt                                                                                                                                                                                                                                    >  2:501710/1‑251 (MQ=255)                   gtTAGCAACGTTGAGATGGTTGGTGACAATGCGCAATTTGCTGTGATT‑GAGCAGTGCGTGCGCTACCGTTTCCGGCGTGGTGCCGA‑TAGCG‑ATAAACAGCGTCGAGCCATTGGGGATTTGCT‑CCGCCACTTTGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTcc                                                                                                                                                                                                                              <  1:111096/250‑1 (MQ=255)                      agaaaCGTTGAGATTGTTGGTGACAATGCGAAATTTGCTGTGATT‑GAGCAGGGCGTGCGCTCCCTCTTCCGGGGTGGTGCCGA‑TATCG‑ATAAACAGCGTCGAGCCATGGGGGATTTGCT‑CCGCCACTTTGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGcc                                                                                                                                                                                                                           <  1:381069/247‑1 (MQ=255)                             ttGAGATTTTTGGGGACAATGCGAAAATTGCTGTGATT‑GAGCAGTGCGTGCGCTCCCGCTTCCGGCGTGGTGCCGA‑TATCG‑ATAACCAGCGTCGAGCCATTGGGGATTTGCT‑CCGCCACTTTGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCa                                                                                                                                                                                                                     <  1:68967/249‑1 (MQ=255)                                 gaTTGTTGGTGACAATGCGCAAATTGCTGTGATT‑GCGCAGTGCGTGCGATACCGCTTCCGGAGTGGTGCCGC‑TATCG‑ATAAACAGCGTCGCGCCATTGGCGATTTGCT‑CCGCCACTTTGCGGGCGCTGCGC‑TCTTTTTCGTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGACGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCGCCGCCAGCTCATTGAgg                                                                                                                                                                                                               >  2:375265/1‑251 (MQ=255)                                   ttgttgGTAACACTGCGCAAATTGCTGTGATT‑GAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGA‑TATCG‑ATAAACAGCGTCGAGCCATTGGGGATTTGCT‑CCGCCACTTTGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCACCTCATTGAGGTc                                                                                                                                                                                                             >  2:292958/1‑251 (MQ=255)                                    tggtgGTGACAATGCGCAAATGGCTG‑GATTTGAGAAGTGCGTGCGCTCCCGCTTCCGGC‑TGGTGGCGAATATGG‑ATAAACAGCGTCTAGCCATTGGGGATTTGCT‑CCGCCCCTTTGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGGGGC‑CTTGCGATCGTGCCACGGCGTTTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTc                                                                                                                                                                                                             <  1:13819/247‑1 (MQ=255)                                           gACAATGCGCAAATTGCTGTGATT‑GAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGA‑TATCG‑ATAAACAGCGTCGAGCCATTGGGGATTTGCT‑CCGCCACTTTGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGa                                                                                                                                                                                                      >  1:359881/1‑250 (MQ=255)                                                    cAAATTGCTGTGATT‑GAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGA‑TATCG‑ATAAACAGCGTCGAGCCATTGGGGATTTGCT‑CCGCCACTTTGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGCACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGCGGTCGCGGCGAATAgcatgcg                                                                                                                                                                                            >  2:324481/1‑245 (MQ=255)                                                        ttGCTGTGATT‑GAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGA‑TATCG‑ATAAACAGCGTCGAGCCATTGGGGATTTGCT‑CCGCCACTTTGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCt                                                                                                                                                                                        >  1:249555/1‑251 (MQ=255)                                                         ggCTGTGTTT‑GAGCGGTGGGTGCGCTACCGCTTCCGGGGTGGGGCCGA‑AATCG‑ATAAACAGCGCCGAGCCATTGGGGATTTGCT‑CCGCCACTTTGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCt                                                                                                                                                                                        <  1:273541/249‑1 (MQ=255)                                                               gATT‑GAGCAGTGCGGGCGCTACCGCTTCCGGCGTGGTGACGA‑TATCG‑ATAAACAGCGTCGAGCCCTTGGGGATTTGCT‑CCGCAACTTTGCGGGCGATGAGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGAACAGGTTTTGCTCCGCCAGCTAATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACgg                                                                                                                                                                                   >  2:442111/1‑249 (MQ=255)                                                                     ggcggtgcgtgcGCTACCGCTCCGGGGGTGGTGCCGA‑TATCG‑ATAAACAGCGTCGAGCCTTTGGGGATTTGCC‑CCGCCACTTTGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATg                                                                                                                                                                            <  1:344470/247‑1 (MQ=255)                                                                     aGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGA‑TATCG‑ATAAACAGCGTCGAGCCATTGGGGATTTGCT‑CCGCCACTTTGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATg                                                                                                                                                                            >  1:331021/1‑251 (MQ=255)                                                                      gCAGTGCGTGCGCGCCCGCTCCCGGCGTGGGGCCGA‑T‑TCGGATAAACAGCTTGGAGCCATTGGGGATTGGCT‑CCGCCACTTTGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGc                                                                                                                                                                           <  1:136938/251‑1 (MQ=255)                                                                       cAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCTA‑TATCG‑ATAAACAGCGTCGAGCCATTGGGGATTTGCT‑CCGCCACTTTGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGct                                                                                                                                                                          <  1:392201/251‑1 (MQ=255)                                                                             gtgcgCTACCGCTTCCGGCGTGGTGCCGA‑TATCG‑ATAAACAGCGTCGAGCCGTTGGGGATTTGCT‑CCGCCACTTTGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAcc                                                                                                                                                                     >  1:314096/1‑250 (MQ=255)                                                                             gtgcgCTACCGCTTCCGGCGTGGTGCCGA‑TATCG‑ATAAACAGCGTCGAGCCATTGGGGATTTGCT‑CCGCCACTTTGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAcc                                                                                                                                                                     >  1:444483/1‑250 (MQ=255)                                                                              tgcgCTACCGCTTCCGGCGTGGTGCCGA‑TATCG‑ATAAACAGCGTCGAGCCATTGGGGATTTGCT‑CCGCCACTTTGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAcc                                                                                                                                                                     >  1:161858/1‑249 (MQ=255)                                                                               gcgcTCCCGCTTCCGGCTGGGTGCCGA‑TATCA‑ATAAACAGCGTCGAGCCATTGGTGATTTGCT‑CCGCCACTTTGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGACCTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGc                                                                                                                                                                  <  1:233048/251‑1 (MQ=255)                                                                               gcgcTACTCCTTCCGGCGAGGTGCCGA‑TATCG‑ATAAACAGAGTCGAGCAATGCGGGCTTTGCC‑CTGACCCTTCGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGTACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGCCGAACCGTCTGCGGGATGACggcggactcatctacccac                                                                                                                                                                  >  2:360328/1‑234 (MQ=255)                                                                                 gcgcccGCTCCCGGCTGGGTCCCGA‑TATCG‑ATAAACAGGGTCGCGCCATTGGGGATTTGCT‑CCGCCACTTTGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGctc                                                                                                                                                                <  1:218794/247‑1 (MQ=255)                                                                                             ggCGTGGTGCGGA‑TATCG‑AAAAACAGCGTCGAGCCATTGGGGATTTGCT‑CCGCCACTTTGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGAc                                                                                                                                                    <  1:360934/251‑1 (MQ=255)                                                                                               ggtggtgCGGA‑TATCG‑ATAACCAGCGTCGAGCCATTGGGGTTTTGCT‑CCGCCACTTTGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACAt                                                                                                                                                  <  1:443710/250‑1 (MQ=255)                                                                                                  ggtgCCGA‑TATCG‑ATAAACAGCGTCGAGCCATTGGGGATTTGCT‑CCGCCACTTTGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAc                                                                                                                                               >  1:362549/1‑251 (MQ=255)                                                                                                    ggtgccaaattcaaaaaaCAGCTTCGTGC‑ATTGGGGGTTTTCTGCCGCCTCTTGGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAc                                                                                                                                               <  1:262493/236‑1 (MQ=255)                                                                                                     gCCGA‑TATCG‑ATAAACAGCGTCGAGCCATTGGGGATTTGCT‑CCGCCACTTTGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAccc                                                                                                                                             <  1:503925/250‑1 (MQ=255)                                                                                                        gA‑TATCG‑ATAACCGGCGTGGAGCCATTGGGGTTTTGCT‑CCGCCACTTTGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgct                                                                                                                                         <  1:11582/251‑1 (MQ=255)                                                                                                        gA‑TATCG‑ATAAACAGCGTCGAGCCATTGGGTATTTGCT‑CCGCCACTTTGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGCGAAATGCTCTCCCAGCTCTTCGGTACTGACATAACCctgct                                                                                                                                         >  2:172890/1‑251 (MQ=255)                                                                                                          gatatagataacgcggttcgccccATTGGGGTTTTGCC‑CCGCCCCTTTGCGGGCGATGCGC‑TCTTTTTCTTGGGTCTGGGGGGC‑CTGGCGATCGTGCCACGGCGGGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgct                                                                                                                                         <  2:331021/229‑1 (MQ=255)                                                                                                             tCG‑ATAAACAGCGTCGAGCCATTGGGGATTTGCT‑CCGCCACTTTGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGttt                                                                                                                                     >  2:343490/1‑251 (MQ=255)                                                                                                                   aaaCAGCGTCGAGCCATTGGGGATTTGCT‑CCGCCACTTTGCGAGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCCGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTATGTAATGACATAAACCTGCGGTTTAAccc                                                                                                                                >  2:21973/1‑250 (MQ=255)                                                                                                                       aGCGTCGAGCCATTGGGGATTTGCT‑CCGCCACTTTGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTc                                                                                                                            >  1:321481/1‑251 (MQ=255)                                                                                                                            cGAGCCATTGGGGATTTGCT‑CCGCCACTTTGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGTCGAATAGTCTGCTGTCTGACGGAGAAATGCTCTACCCGCTCTTCTGTACTGCCATAACCCTGCTGTTTAACCAGTTCGata                                                                                                                        >  2:243179/1‑250 (MQ=255)                                                                                                                             gAGCCATGGGGGTTTTGCT‑CCGCCACTTTGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataat                                                                                                                      <  1:501710/251‑1 (MQ=255)                                                                                                                                 cATTGGGGATTTGCT‑CCGCCACTTTGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCg                                                                                                                  <  1:340086/251‑1 (MQ=255)                                                                                                                                           ttGCT‑CCGCCACTTTGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGt                                                                                                         <  1:375265/250‑1 (MQ=255)                                                                                                                                             gCT‑CCGCCACTTGGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTACCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtgt                                                                                                      <  1:284011/251‑1 (MQ=255)                                                                                                                                                cccGCCTCTTGGCGGCCGATGCGC‑TCTTTTTCTTCGGTCGGGGTGGC‑CTTGCGATCGTGCCCCGGTGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTTCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtgt                                                                                                      <  2:365870/248‑1 (MQ=255)                                                                                                                                                   ccctcTTTGCGGGCG‑TGGGCCTCTTTTTCTTCGGTGTGGGTGGCTCTT‑CGACCGTGCCCGGGCGTTTTACCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtttg                                                                                                 <  2:321481/247‑1 (MQ=255)                                                                                                                                                    ccACTTTGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGGGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtttgt                                                                                                <  1:384776/251‑1 (MQ=255)                                                                                                                                                      aCTTTGCGCGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCCTGGTGGCGCAGGATCAGGATTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGGAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTgtttgttt                                                                                              >  1:342065/1‑251 (MQ=255)                                                                                                                                                              ggCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTtat                                                                                       >  1:33766/1‑250 (MQ=255)                                                                                                                                                              ggCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGGCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTtat                                                                                       >  2:96420/1‑250 (MQ=255)                                                                                                                                                                      tc‑tctTTTTCTTCGGTCTGGGTGGC‑CTTGCGATAGTGCAACGGCGTGTTAACCGCACTGCAAGGCAGAGCCGCACCGCCATGAAGGCGCAGGATCAGGTTTTGCTCCGCCCGCTCATTGAGATCGCGGCGAATAGCTTGCGGGCCGGAGGCGAAATGCTCTACCAGCTCTTCGGTACTGACATAAACCTGCTGCTTAAAACGTTCGATAATACCGTTGAGCCGTTGTGTTTGTTTCCTTTATAACTCCCTg                                                                              >  2:170344/2‑251 (MQ=255)                                                                                                                                                                      gc‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTg                                                                              >  1:97032/1‑251 (MQ=255)                                                                                                                                                                       gtttttttcctgcGGTCTGGGTGGC‑GTTGCGATGGGGCCCCGGCGTGTTAACCGACCGGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTg                                                                              <  2:128675/239‑1 (MQ=255)                                                                                                                                                                                cttaggTCTGGGCGGC‑AATGCGATACTGCCACGCAGTGTTAAACGAACTGAAAGGCAGCGCAGCACAGCCATGATGGCGCAGGATCAGGTTGTACTACGCAAGCCCATTGAGGTCGCGGCAAAAAGTCAGCAGGATGACGGAGAAATGATCTACACGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCCGATCAATCATACCGTTGTGACGTACAGTGTGTTTAATTTATAACTCCATGGAATTAttt                                                                     >  2:85894/5‑251 (MQ=255)                                                                                                                                                                                cggcGGTCGGGGGGGC‑CTTGCGCGCGTGCCACGGGGTGTTAACCGCACTGGGAGGCGGCGCCGCTCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttt                                                                     <  2:88999/248‑1 (MQ=255)                                                                                                                                                                                      tCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgttt                                                               >  1:237782/1‑251 (MQ=255)                                                                                                                                                                                         gggTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgttttcg                                                            <  1:428622/251‑1 (MQ=255)                                                                                                                                                                                         gggTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgtttt                                                              <  1:96420/249‑1 (MQ=255)                                                                                                                                                                                           gTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTcgcg                                                          <  1:442111/251‑1 (MQ=255)                                                                                                                                                                                               c‑cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCAtt                                                      <  1:403439/251‑1 (MQ=255)                                                                                                                                                                                                               aCGGCGTGTTCACCGAACTGGAAGGCAGCGCCGCACCGCCATGAAGGCGCAGGACCAGGCTTTGCTCCGCCAGCTCCTTGTGGCCGCGGCGAATAGTCTGCGGGCTGACGGAGCAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAAAACCGTTGTGACGTTGTGTTTGTTTTATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATtaacaa                                        >  2:474525/1‑250 (MQ=255)                                                                                                                                                                                                                              aaCTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAAcc                        >  1:365197/1‑251 (MQ=255)                                                                                                                                                                                                                                 tGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCAc                     >  1:314179/1‑251 (MQ=255)                                                                                                                                                                                                                                  ggAAGGCAGCGCGGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACg                    <  1:315526/251‑1 (MQ=255)                                                                                                                                                                                                                                     aGGCAGAGCCGCACCGCCATCATAGCGCAGGAACAGGTTTTGCTCCGCCCGCTCATTGAGGTCGCGGCGAATAGACTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGAGACGTTGTGTTTGTTTCATTTATAAATCTCTGGAATTATTTTCGTTTTAGCGAATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGcc                 >  2:330790/1‑250 (MQ=255)                                                                                                                                                                                                                                            gccgcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAAAAATGCTCTACCAGCTCTTCGGGACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGGTTCATTTATAAATCACTGGAAATATATTCGATTTAGCGCAGTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCgg          >  2:53864/1‑251 (MQ=255)                                                                                                                                                                                                                                                cACCTCCATGATGGCGCAGGAACAGGTTTTGCTCCGCCAGCTCATTGCGGTTGCGGCGACTAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTGCGGGACTGACAGAACCCTGCTGCTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATACATCCGTGGAATTATTTTCGTCTTAGCGCCTTGAGCGAATCAACAAAAGCCATCGCTAAGCCCACGGCTAACCAGGCg        >  2:304871/1‑249 (MQ=255)                                                                                                                                                                                                                                                    gCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGAtgtgt  <  1:324481/251‑1 (MQ=255)                                                                                                                                                                                                                                                    gCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGAtgtgt  >  2:90209/1‑251 (MQ=255)                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                TTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATT‑GAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGA‑TATCG‑ATAAACAGCGTCGAGCCATTGGGGATTTGCT‑CCGCCACTTTGCGGGCGATGCGC‑TCTTTTTCTTCGGTCTGGGTGGC‑CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGT  >  NC_000913/3560208‑3560694

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