Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	11,861	C→T	intergenic (‑75/‑302)	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,861	0	C	T	95.1%	111.6 / NA	42	intergenic (‑75/‑302)	yaaI/dnaK	DUF2541 domain‑containing protein YaaI/chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (18/21);  minor base A (2/0);  total (21/21)
Fisher's exact test for biased strand distribution p-value = 2.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGGGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGA‑TTCAGGCACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATATGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGTTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACG  >  NC_000913/11614‑12101                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                 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)                                                                                                                                                                                                                                               aaggttCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGAATTGGCTGCGATTCATTCTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACg  >  2:577177/5‑251 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                 GAGAGTTGTATGCGTTTCACAATGCGACAAACAGGAAGTTTCAGCGCCAGATCGTTGGTTTCGTTACGCGGCATTGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGGGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGA‑TTCAGGCACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATTCTTTATATGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGTTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACG  >  NC_000913/11614‑12101

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