Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	11,854	A→G	intergenic (‑68/‑309)	yaaI ← / → dnaK	DUF2541 domain‑containing protein YaaI/chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	11,854	0	A	G	97.6%	126.1 / ‑6.4	42	intergenic (‑68/‑309)	yaaI/dnaK	DUF2541 domain‑containing protein YaaI/chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (19/22);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGGGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGA‑TTCAGGCACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATATGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGTTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACC  >  NC_000913/11610‑12096                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    ggCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGG‑TTCAgg                                                                                                                                                                                                                                              <  1:430755/251‑1 (MQ=255)     gagagTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGG‑TTCAGGTACg                                                                                                                                                                                                                                          >  1:590683/1‑251 (MQ=255)        agTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGG‑TTCAGGTACGTaa                                                                                                                                                                                                                                       <  1:582369/251‑1 (MQ=255)              aTGCGTTTCTCAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTGCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATTCCCGCTAGTTTACTGCTGATAAAGAGAAGG‑TTCAGGTACGTAAtctttt                                                                                                                                                                                                                                 <  1:601987/251‑1 (MQ=255)               tGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGCGTTTTTGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAACCGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTAATGCTGATAAAGAGA‑GAGTTCAGGTACGTAAtcttttt                                                                                                                                                                                                                                >  1:35444/1‑250 (MQ=255)                   tttCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGG‑TTCAGGTACGTAATcttttctttt                                                                                                                                                                                                                            >  2:137067/1‑251 (MQ=255)                    ttCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGG‑TTCAGGTACGTAATCTTTTCttttt                                                                                                                                                                                                                           >  2:298707/1‑251 (MQ=255)                     tCACAATGCGCCAAACAGGAAGTTTGTGCACCAAATCGATGGCTTCGTTACGCAGCCATGCAATAATGACGAGGAGTTTATGGTCGGTTGCCTGCGCCGTGCAGAACAGCATCAGGCTAATCGCACGGCTGGCGGAACTCGTAACAACGGATTTACTAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAACGAGAAGG‑TTCAGGATCGTACTCTTTTCttttt                                                                                                                                                                                                                           >  2:587811/1‑250 (MQ=255)                        caATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGATAAGG‑TTCAGGTACGTAATCTTTTCTTTttatta                                                                                                                                                                                                                       <  2:110387/251‑1 (MQ=255)                            gCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGG‑TTCAGGTACGTAATCTTTTCTTTTTATTACAAt                                                                                                                                                                                                                   <  2:19042/251‑1 (MQ=255)                                  aaCAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGG‑TTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTtg                                                                                                                                                                                                             <  2:431747/251‑1 (MQ=255)                                    cagtaaGTTTCAGCGCTATATCGTGTGTTACGTTACGCGGCATTGCAATGACACCTATGAGTTTATGGGCGTTTGCCTTCGCCGTGCAGCGCAGCATCAGGCTAATCGCCAGGCTGGCGGAAATAGTAAAAACGGATTTCATAAGGATTCTCTAAGTGGGAAGAGGTAGAGGGATGATTTCCCACTAGTTTACTGCTGATAAATAGAAGG‑TGCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTtgat                                                                                                                                                                                                           <  2:454584/247‑1 (MQ=255)                                                  gcCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGG‑TTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTg                                                                                                                                                                                             >  1:591364/1‑251 (MQ=255)                                                      gATCGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGG‑TTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGat                                                                                                                                                                                         <  2:177710/251‑1 (MQ=255)                                                         cGTTGGTTTCGTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGG‑TTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGattca                                                                                                                                                                                      <  2:591364/251‑1 (MQ=255)                                                                 tCGTTGCGCGGCAGTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCATAGCATCAGGCTAATCGCCAGTCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGG‑TTCAGGTACGTAATGTTTTCTTTTTATTACGATTGTTTGATGAATGCCTTGGCTGCGAGTCATTCTTta                                                                                                                                                                               <  2:390326/250‑1 (MQ=255)                                                                   gTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGG‑TTCAGGTACGTCATCTTTTCTTTTTCTTACAATTTTTTGATGAATGCCTTGGCTGCGattcatatatcgga                                                                                                                                                                             >  2:78159/1‑242 (MQ=255)                                                                   gTTACGCGGCATTGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGG‑TTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGattcat                                                                                                                                                                                     <  1:78159/242‑1 (MQ=255)                                                                                 cAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCATTATTTTACTGCTGATAAAGAGAAGG‑TTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTTATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTg                                                                                                                                                              <  1:326101/251‑1 (MQ=255)                                                                                     gacgcCGAGGAGTTTATGGACGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGCGAAGG‑TTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGTCa                                                                                                                                                           >  2:628077/3‑250 (MQ=255)                                                                                           gaggagTTTATGGACGCTTGCCTCCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGCAATCGTAAAAACGGCTTTCATAAGGATTCTCTTAGTGGGAAGAGGAAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGG‑TTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGTCAATTTTAc                                                                                                                                                    >  2:41792/1‑251 (MQ=255)                                                                                             ggagTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGG‑TTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGTCAATTTTACGt                                                                                                                                                  <  2:159421/251‑1 (MQ=255)                                                                                               agTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTAATGCTGATAAAGAGAAGG‑TTCAGGTACGTAATCTTTTATTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGTCAATTTTACGTCt                                                                                                                                                >  1:343239/1‑251 (MQ=255)                                                                                                      ggTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGG‑TTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCt                                                                                                                                          <  1:44105/250‑1 (MQ=255)                                                                                                      ggTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAAGCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGG‑TTCAGTTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCt                                                                                                                                          <  1:305290/250‑1 (MQ=255)                                                                                                       gTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAACCGGATTTAATAAGGATTATCTTAGTGGGAAGAGGGAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGG‑ATCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGTTGAATGCCTTGGCTGAGATTCATTCTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGc                                                                                                                                        >  2:581148/1‑251 (MQ=255)                                                                                                              gCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGG‑TTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATcg                                                                                                                                 >  2:298081/1‑251 (MQ=255)                                                                                                                 tGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGATGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTAATGCTGATAAAGAGAAGG‑TTCAGGGACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATcgcg                                                                                                                               >  1:17739/1‑250 (MQ=255)                                                                                                                  gcgcCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGG‑TTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGTCAATTTTACGGCTTGTCCTGCCATATCGCGaa                                                                                                                             <  2:213372/251‑1 (MQ=255)                                                                                                                      cGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTACAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGG‑TTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGTCAATTTTCCGTCTTGTCCTGCCATATCGCGAAAttt                                                                                                                         >  2:433122/1‑251 (MQ=255)                                                                                                                          aagcacagcaGAAGGCTAAGCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTTCTGCTGATAAAGAGAAGG‑TTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGAGGAATGCCTTGGCTGCGATTCATTCTTTATCTGAAGAAAATTGCTGTCAATTTTACGTCTTGGCCTGCCATATCGCGAAATTTCTgc                                                                                                                     <  2:132169/250‑1 (MQ=255)                                                                                                                                 gcaTCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGG‑TTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAg                                                                                                              <  1:137067/251‑1 (MQ=255)                                                                                                                                               gCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGG‑TTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAAttttt                                                                                                >  1:414064/1‑251 (MQ=255)                                                                                                                                                         cGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGG‑TTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTcccc                                                                                      <  1:160355/251‑1 (MQ=255)                                                                                                                                                                cGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGG‑TTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTtgatg                                                                               >  2:478993/1‑251 (MQ=255)                                                                                                                                                                   aaaaaCGGATTTCATAAGGATTCTCTTAGTGGAAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGCTAAAGAGAAGG‑TTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACCAAAAATTTTTGCATCTCCACCTTGATGACg                                                                            >  1:534804/1‑251 (MQ=255)                                                                                                                                                                       aCGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGG‑TTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTggg                                                                        <  1:298707/251‑2 (MQ=255)                                                                                                                                                                              tCATACGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGG‑TTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGAc                                                                 >  2:412237/1‑251 (MQ=255)                                                                                                                                                                                              aGTGGGAAGAGGTAGAGGGATGAATACCAACTAGTTTACTGCTGATAAAGAGAAGG‑TTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGCTGACGGGGGTTACGACCCCATTTAGTAGTCa                                                  >  2:372422/1‑250 (MQ=255)                                                                                                                                                                                                 gggAAGAGGTAGAGGGTTGAATACCCTCTAGTTTACTTCTGATAAAGAGAAGG‑TTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTGTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAAcc                                               <  2:332354/250‑1 (MQ=255)                                                                                                                                                                                                       agGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGG‑TTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTGCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCagtga                                        <  1:269426/251‑1 (MQ=255)                                                                                                                                                                                                                                             aaGAGAAGG‑TTCAGGAACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGACTTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAATTTGGGCAGTTGAAAcc  >  2:73563/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    GGCAGAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGGGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGA‑TTCAGGCACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATATGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGTTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACC  >  NC_000913/11610‑12096

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