Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,894,839	C→T	intergenic (+68/+80)	ygcP → / ← ygcQ	putative anti‑terminator regulatory protein/putative flavoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,894,839	0	C	T	100.0%	142.2 / NA	44	intergenic (+68/+80)	ygcP/ygcQ	putative anti‑terminator regulatory protein/putative flavoprotein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/24);  total (20/24)

ATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTAACGCTGCTGC  >  NC_000913/2894592‑2895080                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  aTCCTGCAAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGAGCGTAATGCGATTAACGCGGGTGACGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTg                                                                                                                                                                                                                                                >  1:182215/1‑251 (MQ=255)                aTGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGAcg                                                                                                                                                                                                                                 <  1:396829/251‑1 (MQ=255)                                   gggTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACaa                                                                                                                                                                                                              <  1:28050/251‑1 (MQ=255)                                    ggTGACAGAGAAAATCCGCCCACCGCTGAGTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAg                                                                                                                                                                                                             <  2:212268/251‑1 (MQ=255)                                              aaaaTCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGATTAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTTGCGTTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTTATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGcc                                                                                                                                                                                                   <  1:72007/251‑1 (MQ=255)                                                     gCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCATACAAGTCCCCAGGCCAGAACCg                                                                                                                                                                                            >  2:243766/1‑251 (MQ=255)                                                         aCCGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAgg                                                                                                                                                                                        <  2:258069/251‑1 (MQ=255)                                                                              ggTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGt                                                                                                                                                                   >  1:241342/1‑251 (MQ=255)                                                                              ggTGTGCGAGAAAGAAGAGGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAAAACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGCCGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACg                                                                                                                                                                    <  2:45107/250‑1 (MQ=255)                                                                                           gaagaTGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCCCGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTg                                                                                                                                                      >  1:398524/1‑251 (MQ=255)                                                                                             agaTGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTCGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCg                                                                                                                                                    >  2:438337/1‑251 (MQ=255)                                                                                                     gTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGtt                                                                                                                                            <  1:99764/251‑1 (MQ=255)                                                                                                       aaTGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTg                                                                                                                                          <  2:161336/251‑1 (MQ=255)                                                                                                       aaTGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTGTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTg                                                                                                                                          <  2:336471/251‑1 (MQ=255)                                                                                                         ggCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCTTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTc                                                                                                                                        <  1:150041/250‑1 (MQ=255)                                                                                                                aaCGCGGGAGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTg                                                                                                                                 >  1:136896/1‑251 (MQ=255)                                                                                                                                  cTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCCCGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGAcgcg                                                                                                                                                                                                                               >  2:421184/1‑139 (MQ=255)                                                                                                                                  cTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGAcgcg                                                                                                                                                                                                                               <  1:421184/139‑1 (MQ=255)                                                                                                                                                    ggggTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCa                                                                                              >  1:34362/1‑250 (MQ=255)                                                                                                                                                            gACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGACCCGTAGGCAGGATCATCCGGCCACCCGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGcc                                                                                     >  1:168443/1‑251 (MQ=255)                                                                                                                                                                        aaaaTTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATAcc                                                                         <  2:352789/251‑1 (MQ=255)                                                                                                                                                                            ttACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCAGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGc                                                                       >  1:325186/1‑249 (MQ=255)                                                                                                                                                                                        gATAATCGTAACCAATTGCATTTGGTTTGCTTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGCTGCGACGCTTGACGCGTCTTCTCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCAt                                                         >  2:197332/1‑251 (MQ=255)                                                                                                                                                                                        gATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCAt                                                         <  2:74191/251‑1 (MQ=255)                                                                                                                                                                                        gATAACCGTACCCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCAt                                                         <  2:241342/251‑1 (MQ=255)                                                                                                                                                                                          tAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAacac                                                                                                                                                                     >  1:373623/1‑141 (MQ=255)                                                                                                                                                                                          tAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAacac                                                                                                                                                                     <  2:373623/141‑1 (MQ=255)                                                                                                                                                                                            aTCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTg                                                                                       >  2:7732/1‑217 (MQ=255)                                                                                                                                                                                            aTCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTg                                                                                       <  1:7732/217‑1 (MQ=255)                                                                                                                                                                                               gTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCGGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCg                                                                                                                                                                            >  2:325524/1‑129 (MQ=255)                                                                                                                                                                                               gTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCGGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCg                                                                                                                                                                            <  1:325524/129‑1 (MQ=255)                                                                                                                                                                                                   ccAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGcc                                              >  2:352771/1‑251 (MQ=255)                                                                                                                                                                                                   ccAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGcc                                              <  1:412639/251‑1 (MQ=255)                                                                                                                                                                                                       ttGAATTTGGTTGGCTGTTTAGGCCGCTCGCCACATCCGACATTCAGCTTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAAGCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCGGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGccatc                                           <  2:136896/250‑3 (MQ=255)                                                                                                                                                                                                         gAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCcaccac                                         >  1:304328/1‑250 (MQ=255)                                                                                                                                                                                                             ttGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAAc                                     <  1:285389/250‑1 (MQ=255)                                                                                                                                                                                                              gggtttgtttgttgggCCCCACGCCACATCCGACATGCAGCGTCTGATGCGACGCTTGAAGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTAGTGGCAATCTGCGTGAATGGTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCCGTGAAAATACCGCAGCGCTTGCGTCATGATTAATCGCCCCCACAAACt                                    <  2:262448/236‑1 (MQ=255)                                                                                                                                                                                                                       ttGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGc                           <  1:438337/250‑1 (MQ=255)                                                                                                                                                                                                                           aGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTAcg                      >  1:84352/1‑251 (MQ=255)                                                                                                                                                                                                                               cGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTAcgcccg                  >  1:185775/1‑251 (MQ=255)                                                                                                                                                                                                                                  aCGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCa               <  1:408396/251‑1 (MQ=255)                                                                                                                                                                                                                                               cATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACAAGTTATTGACAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTAACGCtgctgc  >  1:221996/1‑251 (MQ=255)                                                                                                                                                                                                                                                    aGCGTCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAGTCCCCAGGCCCGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGcc                                                                                               >  2:436105/1‑153 (MQ=255)                                                                                                                                                                                                                                                    aGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGcc                                                                                               <  1:436105/153‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  ATCCTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTAACGCTGCTGC  >  NC_000913/2894592‑2895080

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