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

ATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGG‑CGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGC  >  NC_000913/3560205‑3560696                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                  atttttttttACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGAAAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTCCCGCTCCCGG‑CGTGGTGCCGATATCGATAAACAGCGTCGAGCCATGGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGc                                                                                                                                                                                                                                                    <  1:53098/245‑1 (MQ=255)     tcttTTACCACCAACGTGTTAGAACCGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGGGTGCGCTACCGCTTCCGG‑CGTGGTCCCGATATCGATAAACGGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGc                                                                                                                                                                                                                                                    <  1:103106/247‑1 (MQ=255)        tttACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGG‑CGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGAt                                                                                                                                                                                                                                             >  1:87654/1‑251 (MQ=255)           aCCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGAAGTGCGTGCGCTACCGCTTCCGG‑CGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATccg                                                                                                                                                                                                                                          >  1:457545/1‑249 (MQ=255)                   ggtgttagaaagttggagaTGGTGGGTGACAAGGCGCAAATTGTTGTGATGGGGCAGTGCGTGCGCTACGGCTTCGGG‑GGTGGGGCCGATATGGAGAAACGGCGTCGCGCCATTGGGGTTTTCCACCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCAGGGCGGGTTAACCGAACTGGAAGGCAGAGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCt                                                                                                                                                                                                                                  <  1:524468/236‑1 (MQ=255)                   cGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGG‑CGTGGTGCGGATAGCGATAAACGGCGTCGAGCCTTTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCt                                                                                                                                                                                                                                  <  1:143597/251‑1 (MQ=255)                            aatgtggaggttgttgGGAACAATGCGCAAGTGGCTGGGATGGGGCTGTGGGTGCGC‑ACCCGTTCCGGGCGGGGTGCGGATATCGATAAACAGCGGGGGGCCTTTGGGGATTGGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGGGCCCCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑TCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCt                                                                                                                                                                                                                         <  2:116278/241‑1 (MQ=255)                               gTTGAGATTGTTGGTGACAATGCCCACATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGG‑CGTGGTGCCGATATCGACAAACAGCGTCGAGCCATTGGGGATTTGATACGCCACTTTGCGGCCGCTGCGCTCTTTTTCTTCGGTCTGGGCGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCTCCGC‑CACGCCATGAGGGCGCAGGATACGGTTTTGCTCCGCCAGCTCAt                                                                                                                                                                                                                      >  2:525203/1‑251 (MQ=255)                                            gTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGG‑CGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggcg                                                                                                                                                                                                         >  2:494266/1‑251 (MQ=255)                                                    ccgaaatttttgttattggAGCGGTGCGTGCGCTCCCGTTTCCGG‑TGTGGTCCCGTTCTTCATAAACAGCGTCGAGCCATTGGGTATTTGCGCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCGGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGTTCGCGGCGAATAGTCt                                                                                                                                                                                                 <  1:439088/233‑1 (MQ=255)                                                     cgcAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGG‑CGTGGTGCCGATATCGATAAACAGCGTCGAGACATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGAGGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCGCCGCCAGTTCATTGAGGTCGCGGCGAATCGTCTg                                                                                                                                                                                                >  2:18921/1‑251 (MQ=255)                                                         aaTTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGG‑CGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCggg                                                                                                                                                                                            >  2:552214/1‑251 (MQ=255)                                                                  gATTGAGCAGTGCGTGCGCTACCGCTTCCGG‑CGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACgg                                                                                                                                                                                     >  1:148554/1‑249 (MQ=255)                                                                       aGCAGTGCGTGCGCTACCGCTTCCGG‑CGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATg                                                                                                                                                                              >  1:287962/1‑251 (MQ=255)                                                                               gtgcgCTACCGCTTCCGG‑CGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAcc                                                                                                                                                                       >  1:350645/1‑250 (MQ=255)                                                                                        cGCTTCCGG‑GGTGGTGCCGATATCGATAAACAGCGTCGAGCCTTTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCgg                                                                                                                                                             <  1:313497/251‑1 (MQ=255)                                                                                         gCTTCCGG‑CGTGGTGCCGATATCGATAAACAGCGTCGAGCCATGGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGt                                                                                                                                                            <  1:102940/251‑1 (MQ=255)                                                                                               gg‑CGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGAc                                                                                                                                                      >  1:520602/1‑251 (MQ=255)                                                                                                       ggCCGATATCGATAAACAGCGTCGAGCCTTTGGGGATTTGCCCCGCCACTTTGCGGGGGATGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTGCGATCGTGCCACGGCTTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAccc                                                                                                                                               <  1:232042/250‑1 (MQ=255)                                                                                                         ccGATATCGTTAAACAGCGCCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctg                                                                                                                                             <  1:525585/251‑1 (MQ=255)                                                                                                                ggaaaaaCAGCGTCGGGCCATGGGGGATTTGCCCCGCCCCTTGGCGGGAGATGCGCTCTTTTTCTTCGGTCGGGGGGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGttt                                                                                                                                       <  2:385579/246‑1 (MQ=255)                                                                                                                 gATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTa                                                                                                                                      >  1:73489/1‑250 (MQ=255)                                                                                                                      ccAGCGTCGAGCCTTTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGt                                                                                                                                <  1:294074/250‑1 (MQ=255)                                                                                                                      aCAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGt                                                                                                                                >  1:380781/1‑251 (MQ=255)                                                                                                                        ggCGTCGGGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTc                                                                                                                              <  1:373753/250‑1 (MQ=255)                                                                                                                              tctcccttGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTTTGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGACTTCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataat                                                                                                                        >  1:429360/7‑251 (MQ=255)                                                                                                                                   aTTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGCTAATACCGt                                                                                                                   >  2:305556/1‑251 (MQ=255)                                                                                                                                             tGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGtt                                                                                                          <  1:155708/250‑1 (MQ=255)                                                                                                                                             tGCTCCGCCACTTGGCGTGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCCCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGtt                                                                                                          <  1:18921/250‑1 (MQ=255)                                                                                                                                               cTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGt                                                                                                                                         >  1:168385/1‑217 (MQ=255)                                                                                                                                                cccgacgcttggCGGGCGAGGCGCTCTGTTTTTGCGGTGTGGGTGGCATTGCGAGGGTGCCACGGCGTGTTAACCGACCTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATATCCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTgtgt                                                                                                      <  2:429151/240‑1 (MQ=255)                                                                                                                                                        tttGCGGGCGTTGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTc                                                                                               >  1:297270/1‑250 (MQ=255)                                                                                                                                                               gCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTtata                                                                                        >  1:223872/1‑250 (MQ=255)                                                                                                                                                                                 tcgggttGGGTGGCCTGGCGATAGTGCCCCGGCGTGTAAACAGGACTGGGAGGCGGCGC‑GCGACGGCCAGGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAtttt                                                                      <  2:224645/244‑1 (MQ=255)                                                                                                                                                                                 tcGGTCTGGGTGTCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAtttt                                                                      <  1:233438/250‑1 (MQ=255)                                                                                                                                                                                  gggTCTGGGGGGCCTGGCGACCGTGCCACGGCGTGTTAACCGAACGGGAGGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcg                                                                    <  2:520602/250‑1 (MQ=255)                                                                                                                                                                                  cGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGt                                                                                                                                                                                                   <  1:189097/124‑1 (MQ=255)                                                                                                                                                                                  cGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGt                                                                                                                                                                                                   >  2:189097/1‑124 (MQ=255)                                                                                                                                                                                    gtcggggTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCCCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgtt                                                                  <  2:350645/247‑1 (MQ=255)                                                                                                                                                                                        ggggggCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgtttt                                                                <  1:547431/245‑1 (MQ=255)                                                                                                                                                                                        gggTGGCTTTGCGATCGTGCCACGGCGTGTTAACCGAACGGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgttttcg                                                              <  1:37560/251‑1 (MQ=255)                                                                                                                                                                                              cgtggCGATCGTGCCACGGCGTGTTAACAGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCAtt                                                        <  2:80155/247‑1 (MQ=255)                                                                                                                                                                                                               ggcttgtTACCCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAg                                       <  1:552214/247‑1 (MQ=255)                                                                                                                                                                                                                               gggAAGGCGGCCCCGC‑CCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAACAGTTCGATAATGCCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATGGAGAGAATCAACAAAAGCCATCGCTaaaaaaac                       <  2:429360/250‑6 (MQ=255)                                                                                                                                                                                                                                       aGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAAcc               <  1:223986/251‑1 (MQ=255)                                                                                                                                                                                                                                        gcgcCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAAccc              <  1:302281/251‑1 (MQ=255)                                                                                                                                                                                                                                         cgccgc‑CCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTGCCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCg             <  2:389775/251‑1 (MQ=255)                                                                                                                                                                                                                                           ccgc‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGc           >  1:560597/1‑251 (MQ=255)                                                                                                                                                                                                                                                 ccGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGAtgt      >  2:440187/1‑251 (MQ=255)                                                                                                                                                                                                                                                   gcgatgatgGCGCAGGATCAGGGTTTGCCCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCGCTGTAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCCTCGCTAAACCCACGGCTAAACCGGCGAtgtgt    <  2:457545/248‑1 (MQ=255)                                                                                                                                                                                                                                                    ccATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGAtgtgtg   >  1:197614/1‑251 (MQ=255)                                                                                                                                                                                                                                                     cATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGc  <  2:223872/251‑1 (MQ=255)                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                  ATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGG‑CGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGC  >  NC_000913/3560205‑3560696

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