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%	257.8 / NA	70	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 (35/35);  total (35/35)

AACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGA‑TAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGT  >  NC_000913/3560221‑3560694                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                aaCGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGAGATTGAGCAGTGCGTGCGCTACCGCTTCCAGCGTGGTGCCGATATCGA‑TAAACAGCGTCGCGACATTGGGGATTTG‑CTCCGCCACATTGCGGGCGATGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTGCGCTCGTGCCACGGCGTGTTAACCGGACTGGACGGCAGCGCCGC‑GCCGCCATTATGGcacaagaccaagtttgg                                                                                                                                                                                                                                  >  2:233020/1‑235 (MQ=255)    gtgtTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGA‑TAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑AACGCCATGATGGCGCAGGATCAGGTTTTGCTc                                                                                                                                                                                                                               >  2:329191/1‑251 (MQ=255)                 tgggttgtttgGGGACAATGCGCAAATTGTTGTGATTGGGCAGGGGGGGCGCTCCCGCTTCCGGGGTGGTGCCGATATCGA‑TAAACAGCGCCGGGCCATGGGGGATTTG‑CCCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTGCGATCGTGCCACGGCGTGTTATCCGAACTGGAAGGCAGCGCGGC‑ACCGCGATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCAtt                                                                                                                                                                                                                   <  1:164458/242‑1 (MQ=255)                  gagaTTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGA‑TAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTg                                                                                                                                                                                                                  >  2:167904/1‑250 (MQ=255)                    ggttgttgGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGA‑TAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCCTTGAgg                                                                                                                                                                                                               >  2:182701/3‑251 (MQ=255)                                        aaaTTGATGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGA‑TAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCgg                                                                                                                                                                                           >  1:352399/1‑251 (MQ=255)                                                  gATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGA‑TAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGGTTTGCTCCGCCAGCTCATTGCGGTCGCGGCGAATAGTCTGCGGGCTGACgg                                                                                                                                                                                   >  2:130184/1‑249 (MQ=255)                                                   aTTGAGCAGTGCGTGCGCTCCCGATTCCGGCGTGGAGCCGATATCGA‑TAAACCGCGCCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGCTGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATCGTCTGCGGGCGGActgggg                                                                                                                                                                                >  2:50795/1‑246 (MQ=255)                                                         cAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGA‑TAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGct                                                                                                                                                                          >  2:228886/1‑251 (MQ=255)                                                          aGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGA‑TAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGctc                                                                                                                                                                         >  2:123883/1‑251 (MQ=255)                                                          aGTGCGGGCGCTTCCGCTCCCGGGGTGGGCCCGAT‑GCGAATAAACGGCGTCGACCCGTTGGGAATTTC‑CTCCGCCCCTTTGCGGGGGGTGCGCTCTTTTTCTTCGGGCTGGGGGGCCTGGCGATCGTGCCACGGCTTGTTAACCGAACGGGAAGGCAGCGCCGC‑GCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGctc                                                                                                                                                                         <  1:418051/251‑1 (MQ=255)                                                           gtgcgtgcGCTACCGCTTCCGGCGTGGTGCCGATATCGA‑TAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGTGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGACGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGCCGGAGAACTGctct                                                                                                                                                                        >  2:3024/1‑251 (MQ=255)                                                               gtgcgCTACCGCTTCCGGCGTGGTGCCGATATCGA‑TAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAccc                                                                                                                                                                    >  2:142344/1‑250 (MQ=255)                                                                                   tggtgCCGATATCGA‑TAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTCCCAGCTCTTCGGTACTGACATa                                                                                                                                                 >  2:395197/1‑250 (MQ=255)                                                                                   ggggcccgattccGA‑TAAACAGCGTGGATCCTTGGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTGCTCTCGGGCCACGGCGGGTTAACCGAACTGGAAGTCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATaa                                                                                                                                                <  1:68286/239‑1 (MQ=255)                                                                                    ggtgCCGATACCGA‑TAAACAGCGTCGGGCTATTGGGTTTTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTGCTTCGTTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGACCTGGAAGGCAGCGCCGC‑CCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATaa                                                                                                                                                <  1:402727/250‑1 (MQ=255)                                                                                        ccGATATCGA‑TAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctt                                                                                                                                           >  2:20499/1‑250 (MQ=255)                                                                                          gATATCGA‑TAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATCACCctgct                                                                                                                                         >  1:34916/1‑251 (MQ=255)                                                                                               tcaaaaaaCAGCGCCGAGC‑ATTGGGGTTTTGGCCCCGCCACTTGGGGGGGGTTGCTCTTTTTTTTTGCGGTCTGGGTGGCCTTGCGAGCGGGCCACGGCGTGTTAACGGAACGGGAGGGCAGGGCCGC‑CCCGCCATGATGGCGCAGGATCAGGTTTTGCGCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGtt                                                                                                                                      <  2:288012/245‑1 (MQ=255)                                                                                               cGA‑TAAACAGCGTCGGGCCATTGGGTATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAcc                                                                                                                                              <  1:448221/241‑1 (MQ=255)                                                                                               cGA‑TAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGACTAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAcc                                                                                                                                              >  2:448221/1‑241 (MQ=255)                                                                                                gA‑TAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCAGTACTGACATAACCCTGCTGTTTa                                                                                                                                    >  1:177103/1‑250 (MQ=255)                                                                                                   tAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAcc                                                                                                                                 >  2:248302/1‑251 (MQ=255)                                                                                                        ggcgtctaccctttGGGGTTTTG‑CTCCGCCCCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTTTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTc                                                                                                                            <  1:297176/239‑1 (MQ=255)                                                                                                        aGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTc                                                                                                                            >  1:68072/1‑251 (MQ=255)                                                                                                         gcggggacccATTTGGTATTTG‑CCCCGCCCCTTTGCGGGCGATGCGCTCTTTTTCTTCGTTCTGGTTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCg                                                                                                                           <  1:30026/243‑1 (MQ=255)                                                                                                             cGAGCCTTTGGGGATTTG‑CCCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTACGataa                                                                                                                       <  1:123883/251‑1 (MQ=255)                                                                                                             cGAGCCATTGGGGCTTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGata                                                                                                                        >  2:374396/1‑250 (MQ=255)                                                                                                                  ccatttgttgttctcTCCGCCACTTAGCGGTCTATTCCCGCTTTTTCTTCGGTCGGGGGGTTCTTGCGATCGGTCCACGGCGTGTTAACCGAACTGGAAGCCTGCCCCGC‑GCCGCCTTGATGGCGCAGGATCATGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTTCCAGCTCTGCTGTACTGTCATAACCCTGCTGTTTAACCAGTTCGATAATAcc                                                                                                                   <  1:231705/237‑1 (MQ=255)                                                                                                                      gggggttttcccccGCCACTTTGCGGGCGATGCGCTCTTTTTCTCCGGTCTGGGTGGCCTGGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTtg                                                                                                               <  1:351167/239‑1 (MQ=255)                                                                                                                            ttG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGAAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGCGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGtt                                                                                                        >  1:330433/1‑251 (MQ=255)                                                                                                                            ttG‑CCCCGCCTCTTTGCGGGCGATGGGCTCTTTTTCTTCGGTCGGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACGGGAAGGCAGCGCCGC‑TCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGt                                                                                                         <  2:68072/250‑1 (MQ=255)                                                                                                                                cTCCGCCCCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGCGTGGACTTGCGATCGTGCCACGGCGTGTTAACCGAAATGGAAGGCCGCGCCGC‑ACCGACATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtgtg                                                                                                     >  1:270335/1‑251 (MQ=255)                                                                                                                                cTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtgtg                                                                                                     >  1:88003/1‑251 (MQ=255)                                                                                                                                       aCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGACGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGCAATGCTCTACCAGCTCTTCGGTACTg                                                                                                                                                      >  2:437594/1‑195 (MQ=255)                                                                                                                                       aCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACGGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTg                                                                                                                                                      <  1:437594/195‑1 (MQ=255)                                                                                                                                          ttGCGGGCGATGGGCTCTTTTTCTGCGGTCTGGGTGGCCTTGCGTTCGTCCCGCGGCGTGTTAACCGAACTGGAAGGCGGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCa                                                                                            <  2:363670/250‑1 (MQ=255)                                                                                                                                              gggCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTACCCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTtat                                                                                       <  1:443176/251‑1 (MQ=255)                                                                                                                                                 cGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATaa                                                                                     >  1:419223/1‑250 (MQ=255)                                                                                                                                                   aTGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGACGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTAGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATc                                                                                  >  1:371533/1‑251 (MQ=255)                                                                                                                                                         tctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑TCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAcc                                                                                                                   <  2:255116/212‑1 (MQ=255)                                                                                                                                                         tctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAcc                                                                                                                   >  1:255116/1‑212 (MQ=255)                                                                                                                                                          ttttttCTTCGGTCGGGGTGGCCTTGCGATCGTGCCGGGGGGTGTTAACCGAGCTGGAGGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGa                                                                            <  2:448812/249‑1 (MQ=255)                                                                                                                                                              ttcttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgc                                                                                                                                          <  2:408903/184‑1 (MQ=255)                                                                                                                                                              ttcttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgc                                                                                                                                          >  1:408903/1‑184 (MQ=255)                                                                                                                                                              ttcttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAattat                                                                       >  1:420768/1‑251 (MQ=255)                                                                                                                                                               tctccggtgggggtggcttggggATCGCGCCACGGGGTGTTAACCGAACTGGAAGGCAGCTCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGGCGTTGTGTTTGTTTCATGTATAAATCCCTGGAattat                                                                       <  2:412769/228‑1 (MQ=255)                                                                                                                                                                   cGGTCTGGGTGGCCTTTCGATCGTGCCACGGCGTGTTACCCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCTGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcg                                                                  <  1:50795/251‑1 (MQ=255)                                                                                                                                                                     gTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgtt                                                                <  1:350609/251‑1 (MQ=255)                                                                                                                                                                             ggCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGACTTATTTTCGTTTTCgcgc                                                         >  1:64131/1‑250 (MQ=255)                                                                                                                                                                              gCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCAt                                                       <  1:248302/251‑1 (MQ=255)                                                                                                                                                                                      aTCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtgtg                                                                                                     >  1:87470/1‑197 (MQ=255)                                                                                                                                                                                        cgcgcctggtcGTGTAAACGGAACTGCAGGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTGTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTGTTCGCGCATTGAGCGAATc                                             <  2:448009/241‑1 (MQ=255)                                                                                                                                                                                         gTGCCACGGCGTCTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGACTAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCCTTTATAATTCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATc                                             >  2:253078/1‑250 (MQ=255)                                                                                                                                                                                            ccACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATcaaca                                         <  1:142344/251‑1 (MQ=255)                                                                                                                                                                                               gggCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACaaaa                                      <  1:142262/250‑1 (MQ=255)                                                                                                                                                                                               cggggggtaaaCCGAACGGGAAGCCAGCCCCGC‑ACCGCCATGATGGCGTAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAAGCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACaaaa                                      <  2:419223/242‑1 (MQ=255)                                                                                                                                                                                                  ggtgtTAACCAACCTGGAAGGCGGCGCCGC‑CCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGcc                                   <  2:385759/250‑1 (MQ=255)                                                                                                                                                                                                   gtgtTACCCGAACTGGAAGGCAGCGCCGA‑ACAGCCATGATGGCGCAGGATCAGGATTTGCTCCGCCAGCTCATTGAGGGCGCGGCGAATAGTCTGCGGGATGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATACTACCGTTGAGACGTTGTGTTTGTTTCCTTTATAAATCCATGGAATTATTTTCGTTTTCGCGCCTTGAGCGAATCAACAAAAGccc                                  >  1:403783/1‑250 (MQ=255)                                                                                                                                                                                                   gtgtTAACCGAACTGGAGGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCa                                  <  1:15598/251‑1 (MQ=255)                                                                                                                                                                                                     gtTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCAt                                 <  1:374396/250‑1 (MQ=255)                                                                                                                                                                                                        aacctaccgggatggCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGTCATCGCt                             <  2:160657/238‑1 (MQ=255)                                                                                                                                                                                                            gAACTGGAAGGCAGCCCCGC‑GCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAAc                         <  2:88003/251‑1 (MQ=255)                                                                                                                                                                                                            aaaCTGGAAGGCAGC‑CCGCGTCCGCCAGTATGGCGCAGGAGCGGGGTTGGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGGAATGCTCTTCCGGCTTTTCTGTACTGACATATCCCTGCTGTTTAACCAGTTCGATAATACCGTTGTTTCGTTTTTTATTTTTCATTTATAAAGCCCTGGAATTATTTTCGTTTTCGTGCATTGAGCGAATCAGCAGAAGCCATCGCTAAAc                         <  2:270335/250‑1 (MQ=255)                                                                                                                                                                                                             aaCTGGAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAAcc                        <  1:228886/251‑1 (MQ=255)                                                                                                                                                                                                                 ggAAGGCCGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAGACCCAc                     <  1:130184/250‑1 (MQ=255)                                                                                                                                                                                                                  gAAGGCAGCGCCGC‑ACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACgg                   >  1:436265/1‑251 (MQ=255)                                                                                                                                                                                                                   caGGCAGCGCCGC‑ACCGCCATGATGACGCAGGATCACGTTTTGCACCCCCAGCCCATTGAGGTAGAGCCGACTAGTCTGCGGGATGACGCAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGAAAATACCGTTGTGACGTTGTGTTTGCTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGCATCAACAAAAGCCATCGCTAAACCCACGGc                  >  2:445829/2‑251 (MQ=255)                                                                                                                                                                                                                   aagtcggcgcCGC‑ACCGCCATGATGGCGCAGGAGCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGTCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGc                  <  1:33259/245‑1 (MQ=255)                                                                                                                                                                                                                                    gCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCCTTGAGGTCGCGGCGAATAGTCTGCGTGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGAtgtgt  >  2:342994/1‑251 (MQ=255)                                                                                                                                                                                                                                    gCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGAtgtgt  >  2:202951/1‑251 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                AACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGA‑TAAACAGCGTCGAGCCATTGGGGATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGC‑ACCGCCATGATGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGT  >  NC_000913/3560221‑3560694

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