Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,571,280	C→T	E99K (GAG→AAG)	yjiR ←	fused putative DNA‑binding transcriptional regulator/putative aminotransferase YjiR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,571,280	0	C	T	100.0%	133.6 / NA	45	E99K (GAG→AAG)	yjiR	fused putative DNA‑binding transcriptional regulator/putative aminotransferase YjiR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/33);  total (12/33)

TGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGC‑ACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGTGACGATACAGCCCT  >  NC_000913/4571076‑4571524                                                                                                                                                                                                              |                                                                                                                                                                                                                                                     tGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGC‑ACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGgcgc                                                                                                                                                                                                        >  1:268221/1‑251 (MQ=255)  gACCGCCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCGCGCAATACTCGCTGACCGGCTAACAAGTCAAAACCGAGAACGGTTTGCAGATTATGCTGC‑ACCCCCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCt                                                                                                                                                                                                       <  1:301620/251‑1 (MQ=255)   accaccGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGC‑ACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTg                                                                                                                                                                                                      >  1:590550/1‑251 (MQ=255)    ccaccGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTCACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGC‑ACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCTGGCGCTGc                                                                                                                                                                                                     >  2:202088/1‑251 (MQ=255)                 gTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGCCCGGCTAACAAGTCATAACCGAAAACGGTTTGCAGATTATGCTGC‑ACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGCGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGCcgcctc                                                                                                                                                                                         >  2:454673/1‑247 (MQ=255)                      gCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGACATAACCGAGAACGGTTTGCAGATTATGCTGC‑ACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGCACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATcggc                                                                                                                                                                                   >  1:279011/1‑251 (MQ=255)                             gCGGGAAATCTGCTCCCGCATTAATCGCTGACCGGCGAACAGGCCATATCCGAGAACGTTTTGCAGATTATGCTTC‑ACCCCCCGTCTTAGCCCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCTCCTGATCCCACTGGGGAATGTCCACCGGGCGCTTCACCGGACGCGTCATcttctgtactg                                                                                                                                                                            <  2:564903/251‑11 (MQ=255)                              cGGGCAATCTGCTCTCGCAATACTCGCGGACCGGCTAAAAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGC‑ACCCCCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTgg                                                                                                                                                                           <  1:321383/251‑1 (MQ=255)                                   aaTCTGCCCGCGCAATACTCGCTGACCGGCTAACATGTCACAACCGGGAACGGCTTGTAGATTATGTTGC‑ACCACCCGCCTTAGGTCCCGCCTCAGCGGTTTCAGGCTGGGCGCTTCGACCTCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATCTCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTTCCTCGGACGCCTCATCGGCGGTACTGGCGGCt                                                                                                                                                                      <  1:425910/251‑1 (MQ=255)                                            ccGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGC‑CCCCCCCGGCTTAGCCCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTAc                                                                                                                                                             <  1:450021/250‑1 (MQ=255)                                             cgcacatcgtgctgaccgctcaaCAAGTCAGAACCGAGATGGGTTTGAAGATTATGCTGC‑ACCCCCGGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGTTTGGGCTTTAc                                                                                                                                                             <  2:516170/229‑1 (MQ=255)                                               caatccccGCTGACCGGCTAACAAGACAGAACCGAGAACGGTTTGCAGATTAGCCTGC‑ATCCCCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTTCTTTTGCTTAACGGAACAATGTAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCGCCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGtt                                                                                                                                                          <  1:213506/244‑1 (MQ=255)                                                aatccTCGCTGATCGCCTAACAAGCCATACCCGAGACTGGTTTGCAGATTAGGCTGC‑GCCCCCCGCCTTGGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTTGGCTGTACGTtg                                                                                                                                                         <  2:383052/247‑1 (MQ=255)                                                   actagCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCACATTATGCTGC‑ACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCAGGTGTGCTGTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACAGGACGCGTCATCGGCGGTTCTcgcggacgggggtatcgttcttc                                                                                                                                                      >  2:454819/5‑233 (MQ=255)                                                       ccggACCGGCTCACAGGGCATAGCAGAGAGCGTTTTGCAGATGAGGCTGCCACC‑CCCGGCTTACCCCCCGCCACAGCGGTGTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTGTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCCCCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACa                                                                                                                                                  <  2:484682/248‑1 (MQ=255)                                                        cTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGC‑ACCCCCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGGACTGGCGGCTGGGCttt                                                                                                                                                               <  1:481759/237‑1 (MQ=255)                                                            ccGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGC‑ACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTTCTTGCGGCTgggtcttatcttgtgtgacacaat                                                                                                                                              >  2:428484/1‑229 (MQ=255)                                                                cTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGC‑ACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTAATGGCGGCGGGGGttttcagtgttggggcaaaaaaacc                                                                                                                                         >  1:440183/1‑229 (MQ=255)                                                                  aaCAAGTCATAACCGAGACCGGTTTTCAGATTATGCTGC‑CCCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAAccc                                                                                                                                         <  1:399531/249‑1 (MQ=255)                                                                  aaCAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGC‑ACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAAccc                                                                                                                                         <  1:427704/249‑1 (MQ=255)                                                                    cAAGTCATACCCGAGAACGGTTTGCAAATTATGCTGC‑ACCACCCGGCTTATCCCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAAc                                                                                                                                     <  1:435478/251‑1 (MQ=255)                                                                          acaacagagagcggtttgagtaTTAGGCTGC‑CCCCCCTGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCTCTATGCGCCTTCAGCATATCCAGCCCCTGATCCCCCTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCt                                                                                                                               <  2:458356/230‑1 (MQ=255)                                                                            aaCCGAGAACGGTTTGCAGATTATGCTGC‑ACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGGTTTTCGCTGTGCGACaaaaaaaaccccacgccgcgcg                                                                                                                             >  1:306198/1‑233 (MQ=255)                                                                                               aTTATGCTGC‑CCCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCg                                                                                                          <  1:402612/251‑1 (MQ=255)                                                                                                ttATGCTGC‑ACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGGCTTACGTTGTGCGGCAAAATAACCCGAACGCGGCGGCGGGGtggttgggttccttgt                                                                                                         >  1:132556/1‑237 (MQ=255)                                                                                                      ccgccccacccgtacgaGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCCTCGACATCCTGTGTGCGTTTGCTTAACGGAACAATGGAACGGTCGCTATGCTCCTTCAGCATATCCAGCACCTGACCCCACTGGGTAATTTCCACCGGGCTCTGCACCTGACGCGGCATCGGCGTTACTTGCGGCTGGGCGTTACGTTTTGCGACAAAATAACCCGAACGCGGCTGGTGAGTGATGAGCATCAtcttctaca                                                                                                    <  2:551898/235‑4 (MQ=255)                                                                                                             ccccGGCTTGGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCGCCTTATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGTGTGATGAGCTTCATCGTCTCCAGCGTCTg                                                                                             <  2:434065/250‑1 (MQ=255)                                                                                                                      aGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATcgg                                                                                                                                                                                    >  1:101497/1‑155 (MQ=255)                                                                                                                      aGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATcgg                                                                                                                                                                                    <  2:101497/155‑1 (MQ=255)                                                                                                                        cccccGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGGCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGTCGCTGCACCGGACGCGTCAGCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCctgct                                                                                  <  1:393666/249‑1 (MQ=255)                                                                                                                                acaGCGGTTTCAGGCTGGTCGCTTCGACATCCGGTGTTCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCTCTGGGTAATTTCCACCGGGCGCTGCACCTGACGCGTCAGCTGCGGTACTGTCGGCTGGGCTTTACGGTGTGCGACAAAATaa                                                                                                                                            <  2:395484/185‑1 (MQ=255)                                                                                                                                      gTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATa                                                                    <  1:351704/251‑1 (MQ=255)                                                                                                                                             gCTGGGCGCTTCGACAGCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATAGCCAGCACCTGTTCCCTCTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCAACGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGc                                                             <  2:590550/251‑1 (MQ=255)                                                                                                                                                gggCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCg                                                           <  1:202088/250‑1 (MQ=255)                                                                                                                                                gggCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGt                                                          <  1:114147/251‑1 (MQ=255)                                                                                                                                                 ggCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTg                                                         <  1:546106/251‑1 (MQ=255)                                                                                                                                                    gcTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCGGATACTGACGCCGTGCTc                                                      <  1:175387/251‑1 (MQ=255)                                                                                                                                                                     gtttttGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCTCCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACTGTGCTGATACTGACGCCGTGCTACTGACTTAAGCTgcgc                                      <  2:393763/248‑1 (MQ=255)                                                                                                                                                                           gCTTAACGGAACAATGTAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGAc                                <  2:469205/250‑1 (MQ=255)                                                                                                                                                                            cTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGTGCGTTACGTTGTGCGACAAAAt                                                                                                                                              <  2:62013/139‑1 (MQ=255)                                                                                                                                                                                  ggAACAATGGAACTGTCGCTATGCGCCTTCAGCAGATCCAGCTCCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCTCCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAAttt                        <  2:173200/251‑1 (MQ=255)                                                                                                                                                                                                gTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGTGACg            <  1:136373/249‑1 (MQ=255)                                                                                                                                                                                                   gCTATGCGCCGTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCAGCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGATCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACt                                                                  <  2:418425/192‑1 (MQ=255)                                                                                                                                                                                                        gcgcCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGGCGGGCTTTTCGTTGTGCGGCAAAATAACCCCAACCCGGGGGCGGGGGGATGAGCGTCATCGTCTCCCACGGCTGATACGCCCGCCGCACGGGGCTGGTACTGCAGCCGGGCTCCTGGGTTAAGCGGCGCACCCACGGCACTTTTTccccctggccgtacacgccc  >  1:141734/1‑235 (MQ=255)                                                                                                                                                                                                        gcgcCTTCAGCATATCCAGCACCGGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATACCCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCGGCGCACCGACGTCAACTTCTCCCCGTGACGATACAGccc   <  2:497577/250‑1 (MQ=255)                                                                                                                                                                                                              |                                                                                                                                                                                                                                                     TGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGC‑ACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGTGACGATACAGCCCT  >  NC_000913/4571076‑4571524

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