Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,518,681	A→C	K52T (AAG→ACG)	insA7 →	IS1 family repressor protein InsA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,518,681	0	A	C	97.4%	91.9 / ‑3.7	38	K52T (AAG→ACG)	insA7	IS1 family repressor protein InsA
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (25/12);  total (26/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.55e-01

GACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCAC‑TGCCGGAC‑ATCAGCGCTATCTCTGCTCTCCTTGCCGTAAAACATGGCAACTA‑CAGT‑TC‑ACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGG  >  NC_000913/4518439‑4518915                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                           gACATCCTCACGCCCCAAAAGGATGGGGATGTCGGTAATGCTGCCAACTGCCTGATTTATTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATAAGATTTCCTTCCTTCTCCGCTACTTAAGCCGTGGTGCGTAACGGAAAAAGCAT‑TGCCGGAC‑ATCAGCGTTATCCCTCCTCTCATTGCCTTAAAACATGGCAACTA‑CGGT‑TC‑ACTTACACCGCCTCTCAGCCCGGTACGCACCAg                                                                                                                                                                                                                                     <  1:164657/250‑1 (MQ=255)                   aaGGATGGGGATTTCGGTAATGCTGACAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTACGCTACTGAAGGCGTGGTGCGTAACGGCAAACGCAC‑TGCCGGAC‑ATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTA‑CAGT‑TC‑AATTACACCGCCTCTCAGCCCGCGACGCACCAGAAAATCATTGATATGGcc                                                                                                                                                                                                                   >  1:305558/1‑250 (MQ=21)                    aGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGGGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGAAGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCAC‑TGCCGGAC‑ATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTA‑CAGT‑TC‑ACTTACACAGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCAt                                                                                                                                                                                                                 >  1:252988/1‑251 (MQ=255)                                 cGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCAC‑TGCCGGAC‑ATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTA‑CAGT‑TC‑ACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCgg                                                                                                                                                                                                     >  1:300932/1‑250 (MQ=255)                                                    cTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCAC‑TGCCGGAC‑ATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTA‑CAGT‑TC‑ACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTcgcg                                                                                                                                                                                             <  1:110552/239‑1 (MQ=255)                                                    cTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCAC‑TGCCGGAC‑ATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTA‑CAGT‑TC‑ACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTcgcg                                                                                                                                                                                             >  2:110552/1‑239 (MQ=255)                                                                 tgatgGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCAC‑TGCCGGAC‑ATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTA‑CAGT‑TC‑ACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTgg                                                                                                                                                                    >  1:87125/1‑251 (MQ=255)                                                                        gATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGGCCTTCCTGCTACGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCAC‑TGCCGGAC‑ATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGAAACTA‑CAGT‑TC‑ACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCaa                                                                                                                                                              >  1:33325/1‑250 (MQ=255)                                                                          tattAAAGTGCTTACGTGGCTTCCATTTCCCTCAGTTGTCCTTGCTGATCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCAC‑TTCAGTAC‑ATAATCCCTATCTCTGCTCTGATTGCAGTACACCATGGCACCTA‑CAGT‑TC‑ACTTACACCGCCTCTCAGCCCGGTACGCACCAGCAAATCATTGATATGGACATGAATGTCGTCGGATTACGCGCCAGTGCACGCATTATGGGcataggcatctaccc                                                                                                                                                           >  2:106845/3‑237 (MQ=255)                                                                            ttgtggggtttgtgTGGCTTCCATTTCCATCAGAGTTCCTTCCTGCTCCGCTACGGAAGGCGGGGGGCGTAACGCAAAAAGCAC‑TGCCGGAC‑ACCAGCGCTATCGCTGCTCTCATTCCCGTAAAATATGGCAACTA‑CAGT‑TC‑ATTTACGCCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACg                                                                                                                                                          <  1:345123/237‑1 (MQ=18)                                                                                    cTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCAC‑TGCCGGAC‑ATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTA‑CAGT‑TC‑ACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGt                                                                                                                                                 >  1:33105/1‑251 (MQ=255)                                                                                    cTTGCGTGGCTTCCATTTCCATAAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCAC‑TGCCGGAC‑TTAAGCGCTATCTCTGCTCTCAATGCCGTAAAACATGGCACCTA‑CAGT‑TC‑ACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATAATTGATATGGACATGCAGGGCCACTGATGTCGCGCCAATGCACGCATTATGGTCGTTGTCCACATCACGGCTTTActt                                                                                                                                                 >  2:119749/1‑249 (MQ=255)                                                                                                ccatttccatCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCAC‑TGCCGGAC‑ATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTA‑CAGT‑TC‑ACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAAc                                                                                                                                     >  1:292399/1‑251 (MQ=255)                                                                                                ccatttccatCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCAC‑TGCCGGAC‑ATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACAA‑CAGT‑TC‑ACTTACACCGCCTCTCAACCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCACGCCAGAGACCGAATTATGGGCGTTTGCCTCAACAAGGTGTTACGTCAATTaaaaaa                                                                                                                                      >  2:351027/1‑250 (MQ=17)                                                                                                  atttccatCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCAC‑TGCCGGAC‑ATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTA‑CAGT‑TC‑ACTTACACCGCCTCTCAGCCCGGTACGCACCCGAAAATCCTTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTc                                                                                                                                   >  1:147262/1‑251 (MQ=255)                                                                                                            gATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCAC‑TGCCGGAC‑ATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTA‑CAGT‑TC‑ACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCTCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGt                                                                                                                         >  1:89836/1‑251 (MQ=255)                                                                                                             aTGTCCTTCCTGCCCCGCTACTGAAGGGGTGGTGCGTAACGGCAAAAGCCC‑TGCCGGCC‑ATCAGCGCTATCTCTGCCCTCATTGCCGTAAAACAGGGCAACTA‑CAGT‑TC‑CCTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTc                                                                                                                        <  1:196850/251‑1 (MQ=17)                                                                                                              gttccttcctGCCCCGCTCCTGAAGGGGTGGTGGTTAACGGCAAAACCCC‑TGC‑GGACCACCAGCGCCTCCCCTCCTCTATTGGCCTTAAAACATGGCAACTA‑CCGT‑TC‑ACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCg                                                                                                                       <  2:87401/249‑1 (MQ=17)                                                                                                                                ggcgAAGGCGTGGTGCGTAGCGTCAAAATCTC‑GGCCTAAC‑ATCAGCGCTACCTCTGTTCCTTTTGCCGTAAAACATGGCAA‑TATCATT‑TC‑ACTTACACCGCCTCTCAGCCCGGTACGCGCCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACTCGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACaa                                                                                                     <  2:219962/248‑1 (MQ=17)                                                                                                                                aCTGAAGGCGTGGTGCGTAACGGCAAAAGCAC‑TGCCGGAC‑ATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTA‑CAGT‑TC‑ACTTCCACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTGGGCCTCAACACGGTTTTACGTCACTTCAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACa                                                                                                      >  1:175323/1‑250 (MQ=255)                                                                                                                                 cTGAAGGCGTGGTGCGTAACGGCAAAAGCAC‑TGCCGGAC‑ATCAGCGCTATCTCTGCTCTAATTGCCGTAAAACATGGCAACTA‑CAGT‑TC‑ACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACAAGGTTTTCCGTCACTTAAAAAACTCAGGCCGGAGCCGGTAACCTCGCGCATACAAc                                                                                                    >  2:328472/1‑251 (MQ=17)                                                                                                                                   acaGGCGTGGGGCGTAACGTCAAAA‑CTCGTGCCGGAC‑ATCGGCGCTATCCCTGCCCTCATTGCCTAAAAACTTGG‑AACCT‑CCGTGTC‑CCTGCCACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCg                                                                                                  <  1:105729/249‑1 (MQ=17)                                                                                                                                     aGGCGTGGTACGTAACGGCAAAAGCAC‑TGCCGGAC‑ATCAGCGCTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTA‑CAGT‑TC‑ACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGAACGCATTATGGGCGTTGGCCTCAACACGGATTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCggg                                                                                                >  2:103567/1‑251 (MQ=255)                                                                                                                                        cGTGGCGCGTAACGGCAAAAGCCA‑GGCCGTAC‑ATCAGCGCTATCTCTGCGCCCATGGCCGAAAAACAGGGCAACAA‑CAGT‑TC‑ATTGACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTGAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAg                                                                                             <  2:185920/251‑1 (MQ=17)                                                                                                                                                       cAAAAGCAC‑TGCCGGAC‑ATCAGCGCTATCTCTGCTCTCCTTGCCGTAAAACATGGCAACTA‑CAGT‑TC‑ACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTg                                                                              >  1:162458/1‑251 (MQ=255)                                                                                                                                                        aaaaGCAC‑TGCCGGAC‑ACCAGCGCTATCTCGGCTCTCATTGCCTTAAAACATGGCAACTA‑CAGT‑TC‑ACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGTCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACTGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTgc                                                                             <  1:16594/251‑1 (MQ=18)                                                                                                                                                                  gCCGGAC‑ATCAGCGCTACCTCTGCTCCCATTGCCGTAAAACAGGGCAACTA‑CAGT‑CC‑ACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGaaa                                                                      <  1:73632/249‑1 (MQ=18)                                                                                                                                                                                gcTATCTCTGCTCTCATTGCCGTAAAACATGGCAACTA‑CAGT‑TC‑ACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTggg                                                       >  1:281601/1‑251 (MQ=255)                                                                                                                                                                                            ctcattGCCGTAAAACATGGCAACTA‑CAGT‑TC‑ACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACCGTGGGGCTACGTCGGTg                                           >  1:319473/5‑251 (MQ=255)                                                                                                                                                                                                     atAAAACATGGCAACTC‑CAGC‑TCTATTTA‑GCCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGTCATGAATGGCGTCGGATTGCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCAc                                  <  2:119323/250‑1 (MQ=18)                                                                                                                                                                                                         aaCATGGCAACTA‑CAGT‑TC‑ACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATTGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGCACTGGACGCACAGTGGGGCTACGTCGGTGCTAAATCACGTc                               >  1:255578/1‑250 (MQ=255)                                                                                                                                                                                                              tgCAACTA‑CATT‑TC‑ACTTACACCGCCTCTCAGCCCGGTACGCACAAGAAAAACATTGATATGGCCCTGAATAGAGTCGGCTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTAAACCCGGTTTTACGTCACGTAAAAAACTCAGGCCGTAGTCGGTACCCTCGCGCATACAACCAGACAGTGATGTTATTGTCTGCGCTGAAATGGGCGAACCGTGGGGCTACGTCGGTGCTAAAGCACGTAAGCgct                         >  2:377061/2‑251 (MQ=255)                                                                                                                                                                                                                     a‑CAGT‑TC‑ACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAACGCTGGCTGtt                  >  1:215413/1‑251 (MQ=255)                                                                                                                                                                                                                       cAGT‑TC‑ATTTACCCCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCTCGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGttt                 <  1:280505/251‑1 (MQ=255)                                                                                                                                                                                                                            tC‑ACTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGCCCCCCCAGTGCCCGCCTTCTTGGCGTTGGCCTCAACACGGTTTTACGTCACTTACAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGCGCTAAATCACGTCAGCGCTGGCTGTTTTAcg             >  1:53999/1‑251 (MQ=17)                                                                                                                                                                                                                               aCTTACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCCTTGATATGGCCATGACTGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCATCAACACGGTTTTACGTCCCTTAAAAAACTCAGGCCGTAGACGGTAACCTCGCGCATACAACCGGGAAGTGATGTGCTTGTCTGCGCTGAAATGGACGAACAGTGGGGCGTCGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTAcggg           >  2:143875/1‑249 (MQ=255)                                                                                                                                                                                                                                 ttACACCGCCTCTCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGt          >  1:150524/1‑250 (MQ=255)                                                                                                                                                                                                                                          ctctCAGCCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGTAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGACATGGACGAACCGTAGGGCTACGTCGGGGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAgg  >  2:74220/1‑249 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                           GACATCCTCACTGCCCTAAAGGATGGGGATTTCGGTAATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTAAGGTGCTTGCGTGGCTTCCATTTCCATCAGATGTCCTTCCTGCTCCGCTACTGAAGGCGTGGTGCGTAACGGCAAAAGCAC‑TGCCGGAC‑ATCAGCGCTATCTCTGCTCTCCTTGCCGTAAAACATGGCAACTA‑CAGT‑TC‑ACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGG  >  NC_000913/4518439‑4518915

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