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%	147.4 / NA	45	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 (22/23);  total (22/23)

TTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTG‑CCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGG‑CAGCGCCGCACCGCCATGATGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAA‑ATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCT  >  NC_000913/3560211‑3560679                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                  ttttACAATCAACGAGTTAGCAACGTTGAGATTGTTGGCGACAATGCGCAAATTGCTGTGATTGAGCAGTGAGTGCGCTACCGCTTCCGGCGAGGTG‑CCGATATCGAAAAACAACGTCGAGCCATTGGGGATGGGCTCGG‑CCACGTTGCGGACGATGCGCTCTTTTTATTCGGCATGTGTGGCCTTGCGATCGTACCAAGGCGTGTTAACCGAACTGG‑GAGG‑CAGCGCCGCACCCACATGATGGcgcaggg                                                                                                                                                                                                                              >  2:258242/1‑245 (MQ=255)              gtgtTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTG‑CCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCG‑CCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGG‑CAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTc                                                                                                                                                                                                                 >  1:438611/1‑251 (MQ=255)                gtTAGCAACTTTGATATTGTTGGTGACAATGCGCAAATTGCTGTGATTGGGCAGTGCGTGCGCTACCGCTTCGGGCG‑GGTGGCCGATATCGATAAACGGCTTCGAGCAATTGGTGTTTTGCTCCG‑CCACTTTGCGGACGATGCGCGCTTTTTCTTCGTTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCTAACTGG‑AAGG‑CAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCg                                                                                                                                                                                                               <  1:199071/251‑1 (MQ=255)                      aaCGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTG‑CCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCG‑CCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGG‑CAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCt                                                                                                                                                                                                         >  1:34261/1‑251 (MQ=255)                        cGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTG‑CCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCG‑CCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGG‑CAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTc                                                                                                                                                                                                        >  1:403270/1‑250 (MQ=255)                             agaTTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTG‑CCGCTATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCG‑CCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGCTCGTGCCACGGCGTGTTAACCGAACTGG‑CAGG‑CAGCGCCGCACCGCCATGATGGCGCCGGATCAGGTTTTTCTCCGCCCGCTCATTg                                                                                                                                                                                                    >  2:85535/1‑249 (MQ=255)                                    tgGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTG‑CCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCG‑CCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGG‑CAGCGCCGCACCGCCCTGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcgg                                                                                                                                                                                           >  2:445859/1‑251 (MQ=255)                                       tGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTG‑CCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCG‑CCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGG‑CAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggcg                                                                                                                                                                                         >  1:212766/1‑250 (MQ=255)                                          cAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTG‑CCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCG‑CCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGG‑CAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAAt                                                                                                                                                                                      >  2:85024/1‑250 (MQ=255)                                            aTGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTG‑CCGATCTCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCG‑CCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGG‑CAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGt                                                                                                                                                                                   >  1:205599/1‑251 (MQ=255)                                                            gATTGAGCAGTGCGTGAGCTACCGCTTCCGGCGTGGTG‑CCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCG‑CCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGG‑CAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGGGGTAGCGGCGACTAGTCTGCGGGCTGACgg                                                                                                                                                                     >  2:102939/1‑249 (MQ=255)                                                                       gcgtgcgCTCCCGCTCCCGGCTGGGTC‑CCGATATCAATAACCAGCGTCGGCCCATTGGGGGTTTGCTCCG‑CCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTTTTAACCGAACTGG‑AAGG‑CAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAA‑ATGCTCTAc                                                                                                                                                       <  1:409653/251‑1 (MQ=255)                                                                             gcTACCGCTTCCGGCGTGGTG‑CCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCG‑CCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑CAGG‑CAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCCGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAT‑ATGCTCTACAAGctc                                                                                                                                                 >  2:318568/1‑251 (MQ=255)                                                                                       ccgggcgggggacagATATGGAAAAACAGGGTCGA‑CCATTGGGGGTTTGGTCCGGCCACTTTGCGGACGATGCGCTCTTTTCCTTCGGTGTGGGTGGCCTTGCGATCGTGCCACGGCGTTTTACCCGAACGGG‑AAGG‑CGGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAA‑ATGCTCTACCAGCTCTTCGGTa                                                                                                                                          <  2:236679/235‑1 (MQ=255)                                                                                       ccGGCGGGGTG‑CATATATCAAAAAACAGCGTCGGGCCATTGGGGTTTTGCCCCG‑CCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGG‑CAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAA‑ATGCTCTACCAGCTCTTCGGTACt                                                                                                                                        <  1:479921/250‑1 (MQ=255)                                                                                       ccGGCGGGGGG‑CCAATATCGATAAACGGCGCCGGGCAATGGGGGTTTTGCTCCG‑CCACTTGGCGGACGATGCGCTCTTTTTCTTCGGTCGGGGTGGCCTGGCGATCGTGCCACGGCGTGTTACCCGAACTGG‑AAGG‑CAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAA‑ATGCTCTACCAGCTCTTCGGTACTg                                                                                                                                       <  1:245431/251‑1 (MQ=255)                                                                                                tg‑CCGATAACGCTAAACAGCGTCAAGCCATTGGGGATTTGCTCCG‑CCACTTTGCGGACGAAGCGCTCTTTTTCTTCGGTCTGGGTGGCCTAGCGATCGTGCCACGGCGCGTTAACCGAACTGG‑AAGG‑CATCGCGGCACCACCAAGATGGCGCAGGATCAGGATTTGCTCCGCCAGCTCATTGAGGGCACGACGAATCATCTGCGGGCTGACGGAGAA‑ATGCTCTAACAAGTCTTCGGTAcggccatcaccc                                                                                                                              >  2:436504/1‑240 (MQ=255)                                                                                                 ggctgattacgaaaagcggCGTCGAGCCATTGGGGTTTTCCCCCG‑CCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCGGGGGGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGG‑CAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAA‑ATGCTCTACCAGCTCTTCGGTACTGACATAAccc                                                                                                                              <  1:445859/233‑1 (MQ=255)                                                                                                   ccGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCG‑CCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAACGAACTGG‑AAGG‑CAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCCGCTCATTGAGGTCGCGGCGAATAGACTGCGGGCTGACGGAGAA‑ATGCTCTACCAGCTCTTCGGTAATGACATAACCctt                                                                                                                            >  1:285036/1‑250 (MQ=255)                                                                                                      atatCGAAAAACAGCTTCGAGCCATTGGGGATTTGCTCCG‑CCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGG‑CAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAA‑ATGCTCTACCAGCTCTTCGGTACTGACATAACCctgctg                                                                                                                         <  1:27026/251‑1 (MQ=255)                                                                                                           gataaccggctgcgggCCATGGGGGTTTTGCTCCG‑CCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGG‑CAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAA‑ATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTaa                                                                                                                    <  1:411853/236‑1 (MQ=255)                                                                                                           gATAAACCGAGTCGAGCCATTGGGGATTTGCTCCG‑CCACTTTGCGGGAGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACAGAACTGG‑AAGG‑CAGCGCCGCACCGCCATGATGGCGCCGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGGCTGCGGGCGGACGGGGAACATG‑TATACCAGCTCTTCGGTACTGAAATAACCCTGCTGTTTa                                                                                                                     >  2:56022/1‑250 (MQ=255)                                                                                                           gATAAACAGCGTCGAGCCATTGGGGATTTGCTCCG‑CCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGG‑CAGCGCCGCACCGCCATGATGGCGCCGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGCCGGAGAA‑ATGCTCTACCAGCTCTTCGGTACTGACATAACCCTTCTGTTTa                                                                                                                     >  2:122443/1‑250 (MQ=255)                                                                                                                       cGAGCCATTGGGGATTTGCTCCG‑CCACTTTGCGGACGCTGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGACCTGG‑AAGG‑CAGCGCCGCACCGCCATGGTGGCGAAGGACCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGTCGACTATTCTCCGGGCTGCCGTGGAA‑ATGCACTACCAACTCTTCGTTACTGACCTAACCCTGCTGTTTAACCAGTTCGata                                                                                                         >  2:247380/1‑250 (MQ=255)                                                                                                                                        gccccG‑CCACTTTGCGGACGAGGCTCTCTTTTTCTGCGGTCTGGGTGGCCTGGCGATCGTGCCACGGCGTGTTAACCGAACGGG‑AAGG‑CAGCGCCGAACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAA‑ATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtg                                                                                        <  1:122443/247‑1 (MQ=255)                                                                                                                                            ccgtcACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTGCGATGGTGCCACGGCGTGTTAACCGAACTGG‑AAGG‑CAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAA‑ATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtgtg                                                                                      <  1:385351/245‑1 (MQ=255)                                                                                                                                              gccctttgtcggtcgaggcgcTCTTTTTCTTCGTTCTGGGTGGCCTGGCGATCGTGCCCCGGGGTGTTAACCGAACTGG‑AAGG‑CAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAA‑ATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtttg                                                                                  <  2:109438/234‑1 (MQ=255)                                                                                                                                              ccctcTTTGCGGACGATGCGCTTTTTTTCGTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTACCCGAACTGG‑AAGG‑CAGCGCCGCCCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAA‑ATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTTACGTTGtgtttg                                                                                  <  1:425714/247‑1 (MQ=255)                                                                                                                                               ccccTTTGCGGCCGAGGCGCTCTTTCGCTTCGGTCTGGGTGTCCTGGCGATCGGGCCACGGCGTGTTAACCGAACTGG‑AAGG‑CAGGCCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAA‑ATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgttt                                                                                   <  1:203499/246‑1 (MQ=255)                                                                                                                                                 atttgtcggaggatgcctcctTTTTCTGCGGTCTGGGTGCCCTTGCGACCGTGCCCCTGCTTGTTAACCGAACTGG‑AAGG‑CAGCGCCGCCCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAA‑ATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATCCCGTTGTGACGTTGTgtttgttt                                                                               <  1:271872/232‑1 (MQ=255)                                                                                                                                                    tggCGGCCGATGCGCTCTTTTTCTTGGGGTGGGGTGGCCTGGCGAGCGTGCCCCGGGGTGTTACCCGAACTGG‑AAGG‑CTGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAA‑ATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTGCGATAATACCGTTGTGACGTTGTGTTTGTTTCa                                                                             <  2:205599/248‑1 (MQ=255)                                                                                                                                                       cggacGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTTTTAACCGAACTGG‑AAGG‑CAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAA‑ATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAcc                                                                                                    <  1:357555/220‑1 (MQ=255)                                                                                                                                                       cggacGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGG‑CAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAA‑ATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAcc                                                                                                    >  2:357555/5‑224 (MQ=255)                                                                                                                                                       cggacGATGAGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGG‑CAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGATGACGGAGAA‑ATGCTCTACCAGCTCTTCGGGACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCAttt                                                                          >  2:239280/5‑250 (MQ=255)                                                                                                                                                             aTGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGG‑CAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAA‑ATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATc                                                                   >  2:237301/1‑251 (MQ=255)                                                                                                                                                              ggcgctcgttcccttggggttGGGTGGCCGTGCGCCGGTGCCGCGGCGAGTTACCCGA‑CTGGTAAAG‑CACCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTGTCGCCTCCAGCTCAGTGAGGTCGCGGCGAATAGTCTGCGGGCTGTCGGAGAA‑ATGCGCTACCAGCTCTTCGGTATTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATc                                                                   <  2:289998/230‑1 (MQ=255)                                                                                                                                                               gcgcTCTTTTTCTTCGGTCTGGGTGGCCTGGCGATCGGGCCACGGCGTGTTAACCGAACTGG‑AAGG‑CAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAA‑ATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATccc                                                                 <  1:74422/251‑1 (MQ=255)                                                                                                                                                                 gcTCTTTTTCTTCGGTTGGGGTGCCGTTGCGATCGTGCCACGGCGTGTTTACCGA‑CCAG‑AAGGTGAGCGCCGCACCGCCATGATGGCGCAGGATCATGTTTTGCTCCGCCAGCTCATTGAGGTCGCGTCGAATAGTCTGCGGGCTGACGGAGAA‑ATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATTCCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTg                                                               <  1:318568/251‑1 (MQ=255)                                                                                                                                                                   tctTTTTCTTCGGTCGGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGG‑CAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAA‑ATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCCGTTCGCTAATACCGTTGTGCCGTTGTGTTTGTTTCATTTATAAATCCCTggg                                                             >  2:183715/1‑250 (MQ=255)                                                                                                                                                                         tcttcGGGCGGGGTGGCCTTGCGATCGGGCCACGGCGGGTTAACCGAACGGG‑AAGG‑CAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAA‑ATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAattat                                                        <  2:92261/250‑1 (MQ=255)                                                                                                                                                                                       ggCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGG‑CAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAA‑ATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCgcgc                                          >  2:246277/1‑250 (MQ=255)                                                                                                                                                                                                  cGTGACACGGCGTGTTAACCGACCTGT‑AAGG‑CAGCGCCTCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGTGAATAGTCTGCGGGCTGACGGAGAA‑ATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAAt                               <  1:371767/250‑1 (MQ=255)                                                                                                                                                                                                   gtgcccgggcttgttaccggAACTGG‑AAGG‑AAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAA‑ATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATca                             <  1:212086/232‑1 (MQ=255)                                                                                                                                                                                                        aCGGGGTGTTAACCGAATGGG‑AAGG‑CGGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAA‑ATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACaaa                        <  1:85024/251‑1 (MQ=255)                                                                                                                                                                                                                         cTGG‑AAGG‑CAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAA‑ATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATAAACAAAAGCCATCGCTAAAccc        >  1:240197/1‑250 (MQ=255)                                                                                                                                                                                                                                gg‑CAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAA‑ATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCt  >  1:355573/1‑250 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                  TTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTG‑CCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGG‑AAGG‑CAGCGCCGCACCGCCATGATGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAA‑ATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCT  >  NC_000913/3560211‑3560679

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