Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,340,953	C→T	G226R (GGG→AGG)	melR ←	DNA‑binding transcriptional dual regulator MelR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,340,953	0	C	T	97.8%	130.4 / ‑4.4	46	G226R (GGG→AGG)	melR	DNA‑binding transcriptional dual regulator MelR
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (22/23);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.00e-01

GATATCGTCTTTAGCCGGGAAACGTCTGGCGGCGCTGTTGGCTAAGTTTGCGGTATTGTTGCGGCGACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATG‑ACCCGCTGGAATATCCCCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGT‑CAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCCCAGCCAGAGAGGCTAAATCGCTTGAGCATCAGGCCAATTTCATCAATGGCGAGCTGGCGAATTTGCTCGTTCGGACT  >  NC_000913/4340710‑4341199                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        gataTCGTCTTTAGCCGGGAAACGTCTGGCGGCGCTGTTGGCTAAGTTTGCGGTATTGTTGCGGCGACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATG‑ACCCGCTGGAATATCCTCATTg                                                                                                                                                                                                                                                    <  2:324010/249‑1 (MQ=255) gataTCGTCTTTAGCCGGGAAACGTCTGGCGGCGCTGTTGGCTAAGTTTGCGGTATTGTTGCGGCGACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATG‑ACCCGCTGGAATATCCTCATTGCa                                                                                                                                                                                                                                                  <  1:244566/251‑1 (MQ=255)  atatCGTCTTTAGCCGGGAAACGTCTGGCGGCGCTGTTGGCTAAGTTTGCGGTATTGTTGCGGCGACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATG‑ACCCGCTGGAATATCCTCATTGCAt                                                                                                                                                                                                                                                 >  1:67200/1‑251 (MQ=255)   tatCGTCTTTAGCCGGGAAACGTCTGGCGGCGCTGTTGGCTAAGTTTGCGGTATTGTTGCGGCGACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATAATTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTTGTTCATCGTCAATTGCATG‑ACCCGCTGGAATATCCTCATTGCAt                                                                                                                                                                                                                                                 >  2:256409/1‑250 (MQ=255)           ttAGCCGGGAAACGTCTGGCGGCGCTGTTGGCTAAGTGTGCGGTATTGTTGCGGCGACATGCCGACATATTTGCCGAACGGGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATG‑ACCCGCTGGAATATCCTCATTGCATAGTTgg                                                                                                                                                                                                                                           <  2:173721/248‑1 (MQ=255)                                       tGGCTAAGTTTGCGGTATTGTTGCGGCGACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATG‑ACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGt                                                                                                                                                                                                             <  1:59616/250‑1 (MQ=255)                                                  gCGGTAGTGTTGCGGCGACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATG‑ACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGT‑CAgcg                                                                                                                                                                                                <  1:146915/251‑1 (MQ=255)                                                    ggTATTGTTGCGGCGACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATG‑ACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGT‑CAGCGcc                                                                                                                                                                                              <  1:307481/251‑1 (MQ=255)                                                     ggaTTGTTGCGGCGACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATG‑ACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTTATGGG‑CAGCGCCt                                                                                                                                                                                             <  2:253978/249‑1 (MQ=255)                                                       aTTGTTGCGGCGACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGTGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATG‑ACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGT‑CAGCGCCTGa                                                                                                                                                                                           <  2:391802/251‑1 (MQ=255)                                                                 cgACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATG‑ACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCCCATCGTTGCTGGT‑CAGCGCCTGATCATAGtttt                                                                                                                                                                                 >  2:421349/1‑251 (MQ=255)                                                                  gACATGCCGACTTATTTGCCGAAAGTGATGGAAAAACGACTCCTTGAACGCACGACTACCGTCAAGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGCACGAGGTTGCTGAGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCACG‑ACCCGCTGGAACATCCTCGTTGCATAGATGGCTTTAAGCTTGACTTGCATACCAACCTCGTTGATGGTACA‑CGCCTGTTCATAGTATTc                                                                                                                                                                                >  2:448312/1‑251 (MQ=255)                                                                           aCATATTTGCCGAACTGGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATG‑ACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGT‑CAGCGCCTGATCATAGTTTTCGGCAATaaa                                                                                                                                                                       <  1:431094/251‑1 (MQ=255)                                                                                 ttGCCGATCGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATG‑ACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGT‑CAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCa                                                                                                                                                                  <  2:121253/250‑1 (MQ=255)                                                                                   gCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCTGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCTCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATG‑ACCCGCTGGAATATCCTCATTTCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCTCATCGTTGATGGT‑CAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCa                                                                                                                                                               <  1:277789/251‑1 (MQ=255)                                                                                       aaCGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATG‑ACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGT‑CAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTg                                                                                                                                                           >  2:396341/1‑251 (MQ=255)                                                                                           tGCTGTAAAAACGACTACTGTAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATG‑ACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGT‑CAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCt                                                                                                                                                        <  1:336170/250‑1 (MQ=255)                                                                                                   aaaCGACTACTTGAAAGAAAGCCTGCCGGCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGAATCGCGGTAATGTACTGTTTCATCGTCAATTGCATG‑ACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGT‑CAGCGCCTGATCATAGTTGTCGTCAATAAAGCCCAGCATCTGGCTCACATaaa                                                                                                                                                >  2:382464/1‑250 (MQ=255)                                                                                                     aCGACTACTGGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATTCTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATG‑ACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGT‑CAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAAtt                                                                                                                                             <  2:326827/251‑1 (MQ=255)                                                                                                           cctTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTAGCGCTCAGTAACGCGCGAACGTGGTTGATGCGCAGCGCGGTAATGTACTGTTTCATCGTCAATTGCATG‑ACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGT‑CAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCAt                                                                                                                                       <  2:131009/250‑1 (MQ=255)                                                                                                            ctTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATG‑ACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGT‑CAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATg                                                                                                                                      >  1:127269/1‑251 (MQ=255)                                                                                                              tGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATG‑ACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGT‑CAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATg                                                                                                                                      <  1:410501/249‑1 (MQ=255)                                                                                                               gAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATG‑ACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGATGATGGT‑CAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGgcg                                                                                                                                   >  1:349789/1‑251 (MQ=255)                                                                                                               gAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACGTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGATGATGCGCATCGCGGTAATGTACTGTTTAATCGTCAATTGAATG‑ACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGGTTGAAGTGCTCAGTCACATCGTTGATGGT‑CAGCGCCTGATCATAGTTTTCGGGAATAAAGCCCAGCATCTGGTTAACATAAAATTGCGCATGgcg                                                                                                                                   >  2:16605/1‑251 (MQ=255)                                                                                                                      aGCCTGCCGTCAGAACAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATG‑ACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGT‑CAGCGCCTGCTCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCCTGGCGCGAGACg                                                                                                                            >  2:60540/1‑251 (MQ=255)                                                                                                                             cGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATG‑ACCCGCTGGAAAATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGT‑CAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGtttt                                                                                                                     >  1:487605/1‑251 (MQ=255)                                                                                                                                           cGAGAATACGTTTATCGGTATCGCTCAGTAACGCGCGATCGTGGTTGATGCTCATCGAGGTAATGTTCTGTTTCATCGTCAATTGCAGG‑ACCCGCGGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGGTGATGGT‑CAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGa                                                                                                        <  2:56539/250‑1 (MQ=255)                                                                                                                                               aaTACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATG‑ACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGT‑CAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGttt                                                                                                   >  1:66158/1‑251 (MQ=255)                                                                                                                                                aTACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATG‑ACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGT‑CAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATACAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGtttt                                                                                                  >  2:76255/1‑251 (MQ=255)                                                                                                                                                  ttttttATCGGTATCGCTCAGTTACGCGCGAACGTGGTTGTTGCGCATCGCGGTAATTTACTGTTTCATCTTCAATTGCATG‑ACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAATTTTGCCGTGTTCAGCCACATCGTTGATGGT‑CAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCTAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTa                                                                                                 <  2:490104/248‑1 (MQ=255)                                                                                                                                                  aCTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATG‑ACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGT‑CAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTa                                                                                                 <  1:339761/250‑1 (MQ=255)                                                                                                                                                              aTCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATG‑ACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGT‑CAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCg                                                                                    >  1:420936/1‑251 (MQ=255)                                                                                                                                                              aTCGCTAAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATG‑ACCCGCTGGAATATCCCCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGT‑CAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCg                                                                                    >  1:92425/1‑251 (MQ=255)                                                                                                                                                              aTCGCACAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATG‑ACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGT‑CACCGCCAGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAAATGCGCATGGCGCGAGAAGCTGGTTTTGTGTGTGCGCGAGGTTTTATTgccccgattcg                                                                                    >  2:147067/1‑241 (MQ=255)                                                                                                                                                                cGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATG‑ACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGT‑CAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGt                                                                                  >  2:27910/1‑251 (MQ=255)                                                                                                                                                                    cAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATG‑ACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGT‑CAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTc                                                                                <  1:375602/249‑1 (MQ=255)                                                                                                                                                                     aGTAACGCGCGAACTTGGTTGATTCGCATCGCGGTAATTTACTGTTTCAGCGGCAAGTG‑ATGGACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGAGGGT‑CAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAAGTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGATGTTTTATTGACCAGAATCGGTTCcca                                                                             <  2:326955/251‑1 (MQ=255)                                                                                                                                                                       tAACGCGCGAACGTGGTTGATGCGCATCGCGGTGATGTACTGTTTCATCGTCAATTGCATG‑ACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGT‑CAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGATTTATTGACCAGAATCGGTTCccagc                                                                           >  1:493001/1‑251 (MQ=255)                                                                                                                                                                         aCGCGCGAACGTGGTTGATGCGCAGCGCGGTAATGTACTGTTTCATCGTCAATTTCATG‑ACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGT‑CAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCccagcca                                                                         <  2:24083/251‑1 (MQ=255)                                                                                                                                                                                  cGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATG‑ACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGT‑CAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCCCAGCCGGAGAGGCt                                                                 >  2:88363/1‑250 (MQ=255)                                                                                                                                                                                                         attttattgttCATCGCCAAGTGTATG‑ACCCGCTGGCATATCCTCATTGTATATTTGGCTTTAAGTTTTCCGTTCTCAGCCCTTTCGTTGATGTT‑CAGCGACTGATCATAGTTTTCGGCAATAAAGCCCAGCATGTGGCTAATATAAAATTTCGCATGGCGCGAGACGCTGTGTTTGTGTGTGCTCTAGGTTTTATTGACAAGAATCGGTTCCCATCCAGAGAGGCTAAGTCGCGTGCGCATAAGGCCAAt                                         <  2:365995/242‑1 (MQ=255)                                                                                                                                                                                                         aggTACTGGTTCATCGTCAATTGCATG‑ACCCGCTGGCATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACCTCGTTGATGGT‑CAGCGCCTGATCATAGTTTTCGTCACTAAAGCCCAGCATCTGGCTAACATAAAATTGCGAATGGCGCGAGACGCTGTTTTGGTGTGTGAGCGAGGTTTTATTGCCAAGAATAGTGTCACAGCCAGAGAGGCTAAATCGCTTGAGCATCAGGCCact                                         >  2:382890/3‑249 (MQ=255)                                                                                                                                                                                                                     aTCGGCAATTGCAGG‑CCCCGCTGGAATATCCTCATTGCATAGGTGGCGTTAAGTTTGACGGGCTCAGCCACATCGTTGATGGT‑CAGCGCCTGATCATAGTTGTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCTCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCCCAGCCAGAGAGGCTAAATCTCTTGAGCATCAGGCCAATTTCATCAATGGc                             <  2:66158/251‑1 (MQ=255)                                                                                                                                                                                                                                cATG‑ACCCGCTGGAATATCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGT‑CAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCCCAGCCAGAGAGGCTAAATCGCTTGAGCATCAGGCCAATTTCATCAATGGCGAGCTGGCGa                   >  2:288512/1‑250 (MQ=255)                                                                                                                                                                                                                                         gCTGAAATATCCTCGTTGCATAGTTGGAGTTAATTTTGACATGCTCAGCAACATCGTTGATGGT‑CAGCGACTGATCATAGTTTTCGTCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGTGGTTTTCTTGACCAGACTCGGGTCCCAGCCAGAGAGGCTACATAGCTTGAGCATCAGGCCAATTTCATCAATGTCAAGCTGGCGAATTTGATCg          >  2:168391/1‑251 (MQ=255)                                                                                                                                                                                                                                                 atCCTCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGT‑CAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCCCAGCCAGAGAGGCTAAATCGCTTGAGCATCAGGCCAATTTCATCAATGGCGAGCTGGCGAATTTGCTCGTTCGGACt  >  1:59741/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        GATATCGTCTTTAGCCGGGAAACGTCTGGCGGCGCTGTTGGCTAAGTTTGCGGTATTGTTGCGGCGACATGCCGACATATTTGCCGAACGTGCTGTAAAAACGACTACTTGAACGAAAGCCTGCCGTCAGGGCAATATCGAGAATACTTTTATCGGTATCGCTCAGTAACGCGCGAACGTGGTTGATGCGCATCGCGGTAATGTACTGTTTCATCGTCAATTGCATG‑ACCCGCTGGAATATCCCCATTGCATAGTTGGCGTTAAGTTTGACGTGCTCAGCCACATCGTTGATGGT‑CAGCGCCTGATCATAGTTTTCGGCAATAAAGCCCAGCATCTGGCTAACATAAAATTGCGCATGGCGCGAGACGCTGTTTTTGTGTGTGCGCGAGGTTTTATTGACCAGAATCGGTTCCCAGCCAGAGAGGCTAAATCGCTTGAGCATCAGGCCAATTTCATCAATGGCGAGCTGGCGAATTTGCTCGTTCGGACT  >  NC_000913/4340710‑4341199

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