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%	205.0 / NA	54	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 (24/30);  total (24/30)

TCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGT‑GATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAG‑CCATTGGGGATTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCG‑AACTGGAAGGCAGCGCCGCACCGCCATGATGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGA  >  NC_000913/3560209‑3560689                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                           tcttTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGT‑GATTGAGCAGGGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAACCAGCGTCGAG‑CCATTGGGGATTG‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCG‑AACTGGAAGGCAGCGCCGCACCGCCATGATGGcg                                                                                                                                                                                                                                            <  1:70313/248‑2 (MQ=255)            cAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGT‑GATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAG‑CCATTGGGGATTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGCGTTAACCG‑AACTGGAAGGCAGCGCCGCCCCGCCCTGATGGCGCAGGATCAGGtttt                                                                                                                                                                                                                              >  2:416652/1‑251 (MQ=255)             aaCGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGT‑GATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAG‑CCATTGGGGATTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCG‑AACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTg                                                                                                                                                                                                                             >  2:367349/1‑251 (MQ=255)              aCGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGT‑GATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAG‑CCATTGGGGATTT‑GCTCAGCCACTTTGCGGACGATGCCCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCG‑AACTGGAAGGAAGCGCCGCACCGCCATGATGTCGCAGGATCAGGTTTTGc                                                                                                                                                                                                                            >  2:229121/1‑251 (MQ=255)                     aGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGT‑GATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAG‑CCATTGGGGATTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCG‑AACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCgccgccc                                                                                                                                                                                                                     >  2:193936/1‑250 (MQ=255)                     aGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGT‑GATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAG‑CCATTGGGGATTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCG‑AACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGcc                                                                                                                                                                                                                      >  1:265029/1‑250 (MQ=255)                        accttagagatgtttgGTCACAATGCGCAAATTG‑TGTTGATGGACCAGTGCGTGCGCTACCGCTTCCGGCTTGGTCCCATTATC‑GATAAACAGCGTCGAC‑CCATGGGGAATTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCG‑AACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCt                                                                                                                                                                                                                  <  1:69394/238‑1 (MQ=255)                         ccGTTGAGATTGTGGTTGACAATGCGCAAATTGTTGT‑GTTTGGGCAGGGCGTCCGCTACCGCTCCGGGCGTGGTCCCAATATG‑GATAACCAGCGACGAG‑CCATGGGGGATTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCG‑AACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTc                                                                                                                                                                                                                 <  1:204430/250‑1 (MQ=255)                           gTTGAGATTGTTGGTGACAATGCGCAAATTGCTGT‑GATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAG‑CCATTGGGGATTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCG‑AACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCAt                                                                                                                                                                                                               >  1:458760/1‑251 (MQ=255)                                      tgGTGACAATGCGCAAATTGCTGT‑GATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAG‑CCATTGGGGATTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCG‑AACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcgg                                                                                                                                                                                                    >  1:297928/1‑251 (MQ=255)                                           aCAATGCGCAATTTGCTGT‑TATTGAGCAGTGCGTGCGCTCCCGCTTCCGGCGGGGGGCCAATATC‑GATAACCAGCGTCGAG‑CCTTTGGGGATTT‑CCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGGGCCACGGCGTGTTAACCG‑AACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAAt                                                                                                                                                                                               <  1:428302/251‑1 (MQ=255)                                           aCAATGCGCAAATTGCTGT‑GATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAG‑CCATTGGGGATTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCG‑AACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAAt                                                                                                                                                                                               >  1:23125/1‑251 (MQ=255)                                           aCAATGCGCAAATTGCTGT‑GATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAG‑CCATTGGGGATTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCG‑AACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAAt                                                                                                                                                                                               >  2:321189/1‑251 (MQ=255)                                             aaagccgcaattgttggttaTTGAGCGGTGCTGGCGCTCCCGCTTGCGGCGTGTTGCCGAAACC‑GA‑AAACCGGGGCGAGGCCATTGGGGTTTT‑GCCCCGCCACTTTGCGGGCGGTGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTTTGAACCG‑AACTGGAGGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATa                                                                                                                                                                                              <  2:246985/232‑1 (MQ=255)                                                 gccaaaattgctgttattgccgggtgcgtgcGCCCCCGCTTCCGGCGTGTCGCAGATATG‑TATAAACAGCGTGGCG‑CCGTTGGGGATTT‑CCTCCGCCGCTTTGCGGGCGATGCGCTCGTTTTCTTCGTTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCG‑AACTGGAAGGCAGCGCGGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTc                                                                                                                                                                                           <  1:171136/227‑1 (MQ=255)                                                        tggtgt‑gATTGAGCAGTGCGTGCGCTCCCGTTTCCGGCGGGGTGCCGATATC‑GATAAACAGCGTCGAG‑CCTTTGGGGTTTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCG‑AACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCt                                                                                                                                                                                   <  1:352269/247‑1 (MQ=255)                                                               gATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAG‑CCATTGGGGATTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCG‑AACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACgg                                                                                                                                                                              >  1:362667/1‑249 (MQ=255)                                                                            gtgcgCTCCCGCTCCCGGCGTGGT‑CCCATTTCCGATAAACAGCGTCGAG‑CCTTTGGGTTTTT‑GCTCCGCCCCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCG‑AACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCa                                                                                                                                                               <  1:127261/251‑1 (MQ=255)                                                                                   aCCGCTCCCGGCGTGGTGCTGATATC‑GTTAAACAGCTTCGAG‑CCAT‑GGGGATTTTGCTCCGCCACTTTGCGGGCGTTGCTCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTAAACCG‑AACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCtt                                                                                                                                                         <  1:352229/250‑1 (MQ=255)                                                                                              ggtggtcCCGATATC‑GATAAACGGCGTCGAG‑CCATTGGGGTTTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCCCGGCGTGTTAACCG‑AACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACAt                                                                                                                                             <  1:152699/250‑1 (MQ=255)                                                                                                   cccggATATT‑GATAAACAGCGTCGGG‑CTTTTGGGTATTT‑GCCCCTCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGGGATCGTGCCCCGGGGTGTTAACCG‑AACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAccc                                                                                                                                        <  1:369707/247‑1 (MQ=255)                                                                                                               gttaaaaggcggcgggctttgggggTTTT‑GTTCCGCCACTGTGCGGTCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCG‑AACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTaa                                                                                                                              <  1:137921/230‑1 (MQ=255)                                                                                                                        gTCGTG‑CCATGGGGGTTTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCG‑AACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGat                                                                                                                    <  1:348655/251‑1 (MQ=255)                                                                                                                           cctagccattgtgtttttccTCCGCCACTTTGCGGGCGATGCTCTCTTTTTCTTCGGTCGGGGTGGCCTTGCTATCGTGCCACGGCGTGTTACCCG‑ACCTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCGCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGat                                                                                                                    <  1:389822/231‑1 (MQ=255)                                                                                                                            aG‑CCATTGGGGATTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCG‑AACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataat                                                                                                                 >  1:368726/1‑250 (MQ=255)                                                                                                                            aG‑CCATTGGGGATTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTATGGGAGGCCTTGCGATCGTGCCACGGCGTGTTAACCG‑AACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAACTGCTCTACCAGCTCTTCGGTCCTGACATAACCCTGCTGTTTAACCCGTTCGataat                                                                                                                 >  2:388926/1‑250 (MQ=255)                                                                                                                                  tttggggtttggtcccgccctttggCGGGCGATGCGCTTTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGGGCCCGGGCGTGTTAACCG‑AACTGGAAGGAAGGGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCg                                                                                                             <  2:458760/226‑1 (MQ=255)                                                                                                                                  ttGGGTATTT‑CCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCTTGTTAACCG‑AACTGGAAGGCAGTGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCg                                                                                                             <  1:190711/249‑1 (MQ=255)                                                                                                                                                    cACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCG‑AACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtttgtt                                                                                          >  2:373698/1‑251 (MQ=255)                                                                                                                                                      cTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGGGCCACGGCGTGTTAACCG‑AACGGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtttgtt                                                                                          <  2:8670/249‑1 (MQ=255)                                                                                                                                                       tttGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCG‑AACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCATGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGACTAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCCGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTcc                                                                                       >  2:117825/1‑250 (MQ=255)                                                                                                                                                         ggCGGGGGATGCGCTCTTTTTATTCTGTCTGGGTGGACTTGCGATCGTGACACGGGGTGTTA‑CCGTACCTGGAAGGCAGCGCCGCCCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATCCCGTTGTGACGTTGTGTTTGTTTCAt                                                                                      <  2:210629/249‑1 (MQ=255)                                                                                                                                                          gCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCG‑AACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCAttt                                                                                    <  1:112042/251‑1 (MQ=255)                                                                                                                                                              gggATGCGTTCTCTTTCTTCGCTCGGGTTGGCCTGCCGATCGGGCCCCGGCGGTTTACCCG‑AACTGGAAGGCCGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATaa                                                                                <  2:124230/249‑1 (MQ=255)                                                                                                                                                                 aTGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCG‑AACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATc                                                                             <  1:416652/251‑1 (MQ=255)                                                                                                                                                                  tGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCG‑AACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATcc                                                                            >  1:226585/1‑251 (MQ=255)                                                                                                                                                                   gcgcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCG‑AACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTCACCAGTTCGATAATACCGTTGTGAAGTTGTGTTTGTTTCATTTATAAATccc                                                                           >  2:159444/1‑251 (MQ=255)                                                                                                                                                                           tttCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGCGTTAACCG‑ACCTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCCGGTTTTGCTCCGCAAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGactt                                                                    >  2:353268/1‑247 (MQ=255)                                                                                                                                                                                              cTTGCGATCGTGCCACGGCGTGTTAACCG‑AACTGGAAGGCAGCGCCGCACCGCCATGCTGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGGCCCGGCGAATAGTCTGCGGGCTGACGGAGAAATTCTCTACCAGCTCTTCGGTACTGACATAACACTGCTGTTTAACCAGTTCGCTAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCACTGGCATTATTTTCGTTTTcgcgatttg                                                >  2:429391/1‑246 (MQ=255)                                                                                                                                                                                                ggcgctCGTGCCGCGGGGTGTTCACCG‑AACGGGAAGTCAGCTCGGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAACTATTTTCGTGTTCGCGCATTGa                                               <  2:392483/245‑1 (MQ=255)                                                                                                                                                                                                         gCCACGGCGTGTTAACCG‑AACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATcaac                                     >  1:302667/1‑251 (MQ=255)                                                                                                                                                                                                          cctcggggtgtTAACCG‑AACTGGAAGGCAGCGCCGCCCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATcaaca                                    <  2:36811/244‑1 (MQ=255)                                                                                                                                                                                                            ccGGCGTTTTAACCG‑AACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACaaa                                  <  1:321189/250‑1 (MQ=255)                                                                                                                                                                                                             cGGCGTGTTAACCG‑AACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACaaa                                  >  2:15106/1‑250 (MQ=255)                                                                                                                                                                                                                   gtTAACGG‑GACTGGAAGGCCGCGCCGCCCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCAt                            <  2:451376/250‑1 (MQ=255)                                                                                                                                                                                                                     tAACCG‑AACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCg                          <  1:367349/250‑1 (MQ=255)                                                                                                                                                                                                                               tggcaGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGCCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCAc                >  1:357342/5‑251 (MQ=255)                                                                                                                                                                                                                                ggAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACg               <  1:120692/251‑1 (MQ=255)                                                                                                                                                                                                                                  gcgggcgcgccgcgccGCCATGATGGCGCAGGATGAGGTTTTGCTCAGCCAGCTCATGGAGGTCGCGGCGAATAGGCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACAGAACCCTGCTGTTTAACCAGTTCTATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATGATTTTCGTTGTCGCGCAGGGAGCGAATGAACAAAAGCCGTCGCTAAACCCACGGc             <  2:158675/245‑1 (MQ=255)                                                                                                                                                                                                                                      cAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCt                        <  1:269172/236‑1 (MQ=255)                                                                                                                                                                                                                                      cAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCt                        >  2:269172/1‑236 (MQ=255)                                                                                                                                                                                                                                         cgccgcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCg      <  2:302667/251‑1 (MQ=255)                                                                                                                                                                                                                                             gcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTCACCCGGCg   >  2:437344/1‑250 (MQ=255)                                                                                                                                                                                                                                             gcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCg   >  2:28789/1‑250 (MQ=255)                                                                                                                                                                                                                                             gcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGa  <  2:449796/251‑1 (MQ=255)                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                           TCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGT‑GATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAG‑CCATTGGGGATTT‑GCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCG‑AACTGGAAGGCAGCGCCGCACCGCCATGATGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGA  >  NC_000913/3560209‑3560689

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