Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,588,910	C→T	G57D (GGC→GAC)	yjiX ←	PF04328 family protein YjiX

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,588,910	0	C	T	97.8%	141.5 / NA	45	G57D (GGC→GAC)	yjiX	PF04328 family protein YjiX
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (15/29);  minor base G (0/1);  total (15/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GAACAACAGATGCAGCCGTTGGTCAGCGTTTTGATCTGCGTGGCGCGGTCGCCAATCAGTTGATCATCAACAGAGACTTCGCCGAATTCGTTTTCAATCACGGCAATCTTGTAGCCGTGCTGTTCGTTAAGAATATGGCGCAGCAGGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAATCGGGTTCATGCTCTCTCCATTAACAACAGCGCATACCGCCTTTACCATCGCCGCCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATAGCTCATGTACGGCTTGTCGGGATGGTTGGTCTTCATATGCTCGACGTAGTTGTCGTAGTCTGGAATGCCAATCAACATCTTCGCCGCCTGGCCGAGATATTTTTTTGCCTGTCCTAAGTTACCAAACATAGTTCACTCTGATAAGAACAAAGCCCCGCCGAAGCGGGGCTAAACACGGTTAGTGGTGCGAAGAGATCTTCACGCCG  >  NC_000913/4588663‑4589156                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       gAACAACAGATGCAGCCGTTGGTCAGCGTTTTGATCTGCGTGGCGCGGTCGCCAATCAGTTGATCATCAACAGAGACTTCGCCGAATTCGTTTTCCATCACGGCAATCTTGTAGCCGAGCTGTTCGTTAAGAATATGGCGCAGCAGGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGCAACTGCAATCGGGTTCATGCTCTCTCCATTAACAACAGCGCATACCGCCTTTACCATCGCcgtcgt                                                                                                                                                                                                                                                     >  1:255778/1‑251 (MQ=255)   acaacaGATGCAGCCTGTGGTCAGCTTTTTGATCTGCGTGGCGCGGTCGCCAATCAGGTGATCATCAACAGAGACTTCGCCGAATTCGTTTTCAATCACGGCAATCTTTAAGCCGTGCTGTTCGTTAAGAATATGGCGCAGCAGGTTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAAGCGGGTTCATGCTCGCTCCATTAACAACAGCGCATACCGCCTTTCCCATCTCCGTCGTAg                                                                                                                                                                                                                                                   <  2:283590/251‑1 (MQ=255)                      gTCAGCGTTTTGATCTGCGTGGCGCGGTCGCCAATCAGTTGATCATCAACAGAGACTTCGCCGAATTCGTTTTCAATCACGGCAATCTTGTAGCCGTGCTGTTCGTTAAGAATATGGCGCAGCAGGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAATCGGGTTCATGCTCTCTCCATTAACAACAGCGCATACCGCCTTTACCATCGCCGTCGTAGCGCGCGTTCTGGCGTTcgc                                                                                                                                                                                                                                <  2:306418/251‑1 (MQ=255)                          gCGTTTTGATCTGCGTGGCGCGGTCGCCAATCAGTTGATCATCAACAGAGACTTCGCCGAATTCGTTTTCAATCACGGCAATCTTGTAGCCGTGCTGTTCGTTAAGAATATGGCGCAGCAGGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGTGTAACTGCAATCGGGTTCATGCTCTCTCCATTAACAACAGCGCATACCGCCTTTACCATCGCCGTCGTAGCGCGCGTTCTGGCGTTCGCGgaa                                                                                                                                                                                                                            <  1:191788/251‑1 (MQ=255)                             ttttGCTCTGCGTGGCGCGGTCGCCAATCAGTTGATCATCAACAGAGACTTCGCCGAATTCGTTTTCAATCACGGCAATCTTGTAGCCGTGCTGTTCGTTAAGAATATGGCGCAGCAGGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAATCGGGTTCATGCTCTCTCCATTAACAACAGCGCATACCGCCTTTACCATCGCCGTCGTAGCGCGCGTTCTGGCGTTCGCGgaaga                                                                                                                                                                                                                          >  2:250762/1‑250 (MQ=255)                                 gATCTGCGTGGCGCGGTCGCCAATCAGTTGATCATCAACAGAGACTTCGCCGAATTCGTTTTCAATCACGGCAATCTTGTAGCCGTGCTGTTCGTTAAGAATATGGCGCAGCAGGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAATCGGGTTCATGCTCTCTCCATTAACAACAGCGCATACCGCCTTTACCATCGCCGTCGTAGCGCGCGTTCTGGCGTGCGCGTAAGAAttc                                                                                                                                                                                                                      <  1:357931/250‑1 (MQ=255)                                  aTCTGCGTGGCGCGGTCGCCAATCAGTTGATCATCAACAGAGACTTCGCCGAATTCGTTTTCAATCACGGCAATCTTGTAGCCGTGCTGTTCGTTAAGAATATGGCGCAGCAGGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAATCGGGTTCATGCTCTCTCCATTAACAACAGCGCATACCGCCTTTACCATCGCCGTCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAAttctt                                                                                                                                                                                                                    >  1:38235/1‑251 (MQ=255)                                     tGCGTGGCGCGGTCGCCAATCAGTTGATCATCAACAGAGACTTCGCCGAATTCGTTTTCAATCACGGCAATCTTGTAGCCGTGCTGTTCGTTAAGAATATGGCGCAGCAGGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAATCGGGTTCATGCTCTCTCCATTAACAACAGCGCATACCGCCTTTACCATCGCCGTCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCAt                                                                                                                                                                                                                 >  2:171687/1‑251 (MQ=255)                                       cGTGGCGCGGTCGCCAATCAGTTGATCATCAACAGAGACTTCGCCGAATTCGTTTTCAATCACGGCAATCTTGTAGCCGTGCTGTTCGTTAAGAATATGGCGCAGCAGGGTGGTTTTTCCGGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAACCGGGTTCATGCTCGCTCCATGAACAACAGCGCATACCGCCTTTACCATCGCCGGCGTAGCGCGCGTTCTGGTGTTCGCGGAAGAATTCGTCATAg                                                                                                                                                                                                               <  2:184495/251‑1 (MQ=255)                                        gTGGCGCGGTCGCCAATCAGTTGATCATCAACAGAGACTTCGCCGAATTCGTTTTCAATCACGGCAATCTTGTAGCCGTGCTGTTCGTTAAGAATATGGCGCAGCAGGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAATCGGGTTCATGCTCTCTCCATTAACAACAGCGCATACCGCCTTTACCATCGCCGTCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATAGc                                                                                                                                                                                                              >  2:9847/1‑251 (MQ=255)                                            cgcgGTCGCCAATCAGTTGATCATCAACAGAGACTGCGCCGAATTCGTTTGCAATCACGGCAATCTTGTAGCCGTGCTGTTCGTTAAGAATATGGCGCAGCAGGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAATCGGGTTCATGCTCTCGCCATGAACGACAGCGCATTCCGCCTTTACCATCGCCGTCGTCGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATAGCTCAt                                                                                                                                                                                                          <  1:32559/251‑1 (MQ=255)                                             gcgGTCGCCAATCAGTTGATCATCAACAGAGACTTCGCCGAATTCGTTTTCAATCACGGCAATCTTGTAGCCGTGCTGTTCGTTAAGAATATGGCGCAGCAGGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAATCGGGTTCATGCTCTCTCCATTAACAACAGCGCATACCGCCTTTACCATCGCCGTCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATAGCTCATg                                                                                                                                                                                                         <  1:228720/251‑1 (MQ=255)                                               ggTCGCCAATCAGTTGATCATCAACAGAGACTTCGCCGAATTCGTTTTCAATCACGGCAATCTTGTAGCCGTGCTGTTCGTTAAGAATATGGCGCAGCAGGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAATCGGGTTCATGCTCTCTCCATTAACAACAGCGCATACCGCCTTTACCATCGCCGTCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATAGCTCATGTa                                                                                                                                                                                                       <  1:336958/251‑1 (MQ=255)                                                gTCGCCAATCAGTTGATCATCAACAGAGACTTCGCCGAATTCGTTTTCAATCACGGCAATCTTGTAGCCGTGCTGTTCGTTAAGAATATGGCGCAGCAGGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAATCGGGTTCATGCTCTCTCCATTAACAACAGCGCATACCGCCTTTACCATCGCCGTCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATAGCTCATGTAc                                                                                                                                                                                                      >  1:278798/1‑251 (MQ=255)                                                      aaTCAGTTTATCATCAACAGAGACTTCGCCGAATTCGTTTTCAATCACGGCAATCTTGTAGCCGTGCTGTTCGTTAAGAATATGGCGCAGCAGGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAATCGGGTTCATGCTCTCTCCATTAACAACAGCGCATACCGCCTTTACCATCGCCGTCGTAGCGCGCTTTCTGGCGTTCGCGGAAGAATTCTTCATAGCTCATGTACGGCTTg                                                                                                                                                                                                <  1:64497/251‑1 (MQ=255)                                                             tGATCATCAACAGAGACTTCGCCGAATTCGTTTTCAATCACGGCAATCTTGTAGCCGTGCTGTTCGTTAAGAATATGGCGCAGCAGGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAATCGGGTTCATGCTCTCTCCATTAACAACAGCGCATACCGCCTTTACCATCGCCGTCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATAGCTCATGTACGGCTTGTCGGGAt                                                                                                                                                                                         >  1:427720/1‑251 (MQ=255)                                                              gATCATCAACAGAGACTTCGCCGAATTCGTTTTCAATCACGGCAATCTTGTAGCCGTGCTGTTCGTTAAGAATATGGCGCAGCAGGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAATCGGGTTCATGCTCTCTCCATTAACAACAGCGCATACCGCCTTTACCATCGCCGTCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATAGCTCATGTACGGCTTGTCGGGAtg                                                                                                                                                                                        >  2:271483/1‑251 (MQ=255)                                                                    caACAGAGACTTCGCCGAATTCGTTTTCAATCACGGCAATCTTGTAGCCGTGCTGTTCGTTAAGAATATGGCGCAGCAGGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAATCGGGTTCATGCTCTCTCCATTAACAACAGCGCATACCGCCTTTACCATCGCCGTCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATAGCTCATGTACGGCTTGTCGGGAtggttgg                                                                                                                                                                                   <  1:196620/250‑1 (MQ=255)                                                                     aaCAGAGACTTCGCCGAATTCGTTTTCAATCACGGCAATCTTGTAGCCGTGCTGTTCGTTAAGAATATGGCGCAGCAGGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAATCGGGTTCATGCTCTCTCCATTAACAACAGCGCATACCGCCTTTACCATCGCCGTCGTAGCGAGCGTTCTGGCGTTCGCGGAAGCATTCTTCATAGCTc                                                                                                                                                                                                            >  2:70087/1‑224 (MQ=255)                                                                     aaCAGAGACTGCGCCGAATTCGTTTTCAATCACGGCAATCTTGTAGCCGTGCTGTTCGTTAAGAATATGGCGCAGCAGGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAATCGGGTTCATGCTCTCTCCATTAACAACAGCGCATACCGCCTTTACCATCGCCGTCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATAGCTc                                                                                                                                                                                                            <  1:70087/224‑1 (MQ=255)                                                                        agagACTTCGCCGAATTCGTTTTCAACCCCGGCAATCTTGTAGCCGTGCTGTTCGTGAAGAATATGGCGCAGCAGGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAATCGGGTTCATGCTCTCTCCATTAACAACAGCGTATACCGCCTTTCCCATCGCCGTCGTAGCGCGCTTTCTGGCGTTCGCGGAAGAATTCTTCATAGCTCCTGTACGGCTTGTCGGGATGGTTGGTCTTc                                                                                                                                                                              <  2:96257/251‑1 (MQ=255)                                                                            aCTTCGCCGAATTCGTTTTCAATCACGGCAATCTTGTAGCCGTGCTGTTCGTTAAGAATATGGCGCAGCAGGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAATCGGGTTCATGCTCTCTCCATTAACAACAGCGCATACCGCCTTTACCATCGCCGTCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATAGCTCATGTACGGCTTGTCGGGATGGTTGGTCTTCatat                                                                                                                                                                          >  1:117463/1‑251 (MQ=255)                                                                             cTTCGCCGAATTCGTTTTCAATCACGGCAATCTTGTAGCCGTGCTGTTCGTTAAGAATATGGCGCAGCAGGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAATCGGGTTCATGCTCTCTCCATTAACAACAGCGCATACCGCCTTTACCATCGCCGTCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATAGCTCATGTACGGCTTGTCGGGATGGTTGGTCTTCATATg                                                                                                                                                                         >  1:312414/1‑251 (MQ=255)                                                                                cGCCGAATTCTTTTTCAATCACGGCAATCTTGTAGCCGTGCTGTTCGTTAAGAATATGGCGCAGCAGGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAATCGGGTTCATGCTCTCTCCATTAACAACAGCGCATACCGCCTTTACCATCGCCGTCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATAGCTCATGTACGGCTTGTCGGGATGTTTGGTCTTCATATGCTc                                                                                                                                                                      <  1:339173/251‑1 (MQ=255)                                                                                      aTTCGTTTTCAATCACGGCAATCTTGTAGCCGTGCTGTTCGTTAAGAATATGGCGCAGCAGGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAATCGGGTTCATGCTCTCTCCATTAACAACAGCGCATACCGCCTTTACCATCGCCGTCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATAGCTCATGTACGGCTTGTCGGGATGGTTGGTCTTCATATGCTCGACGTa                                                                                                                                                                <  2:117463/251‑1 (MQ=255)                                                                                                     cGGCAATCTTGTAGCCGTGCTGTGCGTTAAGAATATGGCGCAGCAGGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAAGCGGGTTCATGCTCTCTCCATTAACAACAGCGCATACCGCCTTTACCATCGCCGTCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATAGCTCATGTACGGCTTGTCGGGATGGTTGGTCTTCATATGCTCGACGTAGTTGTCGTAGTCTgg                                                                                                                                                 <  1:397009/251‑1 (MQ=255)                                                                                                        cAATCTTGTAGCCGTGCTGTTCGTTAAGAATATGGCGCAGCAGGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAATCGGGTTCATGCTCTCTCCATTAACAACAGCGCATACCGCCTTTACCATCGCCGTCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATAGCTCATGTACGGCTTGTCGGGATGGTTGGTCTTCATATGCTCGACGTAGTTGTCGTAGTCTGGAAt                                                                                                                                              <  1:35512/251‑1 (MQ=255)                                                                                                              tGTAGCCGTGCTGTTCGTTAAGAATATGGCGTAGCAGGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAACCGGTTTCATGCTCTCTCCATTAACACCAGCGCATCCCGCCTTTACCATCGCCGTCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATAGTTCATGTACGGCTTGTCGGGATGGTTGGTCTTCATATGCTCGACGTAGTTGTCGTAGTCTGGAATGCCa                                                                                                                                          <  1:73052/249‑1 (MQ=255)                                                                                                                       gCTGTTCGTTAAGAATATGGCGCAGCAGGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAATCGGGTTCATGCTCTCTCCATTAACAACAGCGCATACCGCCTTTACCATCGCCGTCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATAGCTCATGTACGGCTTGTCGGGATGGTTGGTCTTCATATGCTCGACGTAGTTGTCGTAGTCTGGAATGCCAATCAACATCt                                                                                                                                <  1:9847/250‑1 (MQ=255)                                                                                                                          gTTCGTTAAGAATATGGCGCAGCAGGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAATCGGGTTCATGCTCTCTCCATTAACAACAGCGCATACCGCCTTTACCATCGCCGTCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATAGCTCATGTACGGCTTGTCGGGATGGTTGGTCTTCATATGCTCGACGTAGTTGTCGTAGTCTGGAATGCCAATCAACATCTTc                                                                                                                              <  2:369674/249‑1 (MQ=255)                                                                                                                              gTTAAGAATATGGCGCAGCAGGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAATCGGGTTCATGCTCTCTCCATTAACAACAGCGCATACCGCCTTTACCATCGCCGTCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATa                                                                                                                                                                                                                <  1:7297/163‑1 (MQ=255)                                                                                                                              gTTAAGAATATGGCGCAGCAGGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAATCGGGTTCATGCTCTCTCCATTAACAACAGCGCATACCGCCTTTACCATCGCCGTCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATa                                                                                                                                                                                                                >  2:7297/1‑163 (MQ=255)                                                                                                                                tAAGAATATGGCGCAGCAGGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGTTAACTGCACTCGGGTTCATGCTCTCTCCATTACCAACAGCGCATACCGCCTTTACCATCGCCGTCGTCGCGCGCGGTCTGGCGTTCGCGGCAGAATTCTTCATAGCTCATGTACgg                                                                                                                                                                                                    >  2:101818/1‑173 (MQ=255)                                                                                                                                tAAGAATATGGCGCAGCAGGGTGGGTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAATCGGGTTCATGCTCTCTCCATTAACAACAGCGTATTCCGCCTTTACCATCGCCGTCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATAGCTCATGTACgg                                                                                                                                                                                                    <  1:101818/173‑1 (MQ=255)                                                                                                                                    aaTATGGCGCAGCAGGGGGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAATCGGGTTCATGCTCTCTCCATTAACAACAGCGCATACCGCCTTTACCATCGCCGTCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATAGCTCATGTACGGCTTGTCGGGATGGTTGGTCTTCATATGCTCGACGTAGTTGTCGTAGTCTGGAATGCCAATCAACATCTTCGCCGCCTGGCCg                                                                                                                  <  2:38235/251‑1 (MQ=255)                                                                                                                                     atatGGCGCAGCAGGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAATCGGGTTCATGCTCTCTCCATTAACAACAGCGCATACCGCCTTTACCATCGCCGTCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATAGCTCATGTACGGCTTGTCGGGATGGTTGGTCTTCATATGCTCGACGTAGTTGTCGTAGTCTGGAATGCCAATCAACATCTTCGCCGCCTGGCCga                                                                                                                 <  1:266589/251‑1 (MQ=255)                                                                                                                                            gcagcaGGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAATCGGGTTCATGCTCTCTCCATTAACAACAGCGCATACCGCCTTTACCATCGCCGTCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATAGCGCATGTACGGCTTGTCGGGATGGTTGGTCTTCATATGCTCGACGTAGTTGTCGTATTCTGGAATGCCAATCAACATCTTCGCCTCCTGGCCGAGATAttt                                                                                                          <  2:427720/251‑1 (MQ=255)                                                                                                                                                cagGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAATCGGGTTCATGCTCTCTCCATTAACAACAGCGCATACCGCCTTTACCATCGCCGTCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATAGCTCATGTACGGCTTGTCGGGATGGTTGGTCTTCATATGCTCGACGTAGTTGTCGTAGTCTGGAATGCCAATCAACATCTTCGCCGCATGGCCGAGATAttttttt                                                                                                      <  1:130568/251‑1 (MQ=255)                                                                                                                                                 agGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAATCGGGTTCATGCTCTCTCCATTAACAACAGCGCATACCGCCTTTACCATCGCCGTCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATAGCTCATGTACGGCTTGTCGGGATGGTTGGTCTTCAGATGCTCGACGTAGTTGTCGTAGTCTGGAATGCCAATCAACATCTTCGCCGCCTGGCCGAGATATTTTTTTg                                                                                                     <  1:178834/251‑1 (MQ=255)                                                                                                                                                                                           aCTGCAATCGGGTTCATGCTCTCTCCATTAACAACAGCGCATACCGCCTTTACCATCGCCGTCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATAGCTCATGTACGGCTTGTCGGGATGGTTGGTCTTCATATGCTCGACGTAGTTGTCGTAGTCTGGAATGCCAATCAACATCTTCGCCGCCTGGCCGAGATATTTTTTTGCCTGTCCTAAGTTACCAAACATAGTTCACTCTGATAAGAACa                                                           <  1:126525/251‑1 (MQ=255)                                                                                                                                                                                                       ttCATGCTCTCTCCATTAACAACAGAGCATACCGCATTTACAATCGCACTCGTAGCGCGCGTTCTTGCGTTCGCGGCAGAATTCTTCATAGCTCATGTACGGCTTGTAGGGATGGTTGGTCTTCATATGCTCGTCGTAGTTGTCGTAGTCTGGAATGCCAATCAACATCTTTGCCGCCTGGCCGAGATTTTTTTTTGCCTGTCCTAAGTTACCAACCATCACTCACTCTGATCAGAACAAAGCCctgccg                                                >  2:22543/1‑250 (MQ=255)                                                                                                                                                                                                                 ctcCATTAACAACAGCGCATACCGCCTTTACCATCGCCGTCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATAGCTCATGTACGGCTTGTCGGGATGGTTGGTCTTCATATGCTCGACGTAGTTGTCGTAGTCTGGAATGCCAATCAACATCTTCGCCGCCTGGCCGAGATATTTTTTTGCCTGTCCTAAGTTACCAAACATAGTTCACTCTGATAAGAACAAAGCCCCGCCGAAGCGGGGCt                                      >  1:164114/1‑250 (MQ=255)                                                                                                                                                                                                                     aTTAACAACAGCGCATACCGCCTTTACCATCGCCGTCGGAGCGCGCGTTCTGGCGTTCACGGAAGAATTCTTCATAGCTCATGTACGACTTGTCGGGAGGGTTGGTCTTCATATGCTCGACGTAGTTGTCGTAGTCTGGAATGCCAATCAACATCTTCGCCGCCTGGCCGAGATATTTTTTTGACTGTCCTAAATTACCAAAAATAGATCACTCTGATAAGAACAAAGCCCCGCCGAATCGGGGCTAAac                                  >  1:255670/1‑250 (MQ=255)                                                                                                                                                                                                                             caGCGCATACCGCCTTTACCATCGCCGTCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATAGCTCATGTACGGCTTGTCGGGATGGTTGGTCTTCATATGCTCGACGTAGTTGTCGTAGTCTGGAATGCCAATCAACATCTTCGCCGCCTGGCCGAGATATTTTTTTGCCTGTCCTAAGTTACCAAACATAGTTCACTCTGATAAGAACAAAGCCCCGCCGAAGCGGGGCTAAACACGGTTAgt                         <  1:271483/251‑1 (MQ=255)                                                                                                                                                                                                                                         ccTTTACCATCGCCGTCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATAGCTCATGTACGGCTTGTCGGGATGGTTGGTCTTCATATGCGCGCCGTAGTTGTCGTAGTCTGGAATGCCAATCAACATCTTCGCCGCCTGGCCGAGATATTTTTTTGCCTGTCCTAAGTTACCAAACATAGTTCACTCTGATAAGAAAAAAGCCCCGCCGAAGCGGGGCTAAACACGGTTAGTGGAGCGAAgaga             <  2:134887/251‑1 (MQ=255)                                                                                                                                                                                                                                               ccATCGCCGTCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATAGCTCATGTACGGCTTGTCGGGAGGGATGGTCTTCATATGCTCGACGTAGTTGTCGTAGTCTGGAATGCCAAGCAACATCTTCGCCGCCTGGCCGAGATATTTTTTTGCCTGTCCTAAGTTACCAAACATAGTTCACTCTGATAAGAACAAAGCCCCGCCGAAGCGGGGCTAAACACGGTTAGTGGTGCGAAGAGATCTTCa       <  1:67857/251‑1 (MQ=255)                                                                                                                                                                                                                                                    gccgtcGTAGCGCGCGGTCTGGCGTTCGCGGAAGAATTCTTCATAGCTCATGTACGGCTTGTCGGGATGGTTGGTCTTCATATGCTCGACGTAGTGGTCGGAGTCTGGAATGCCAATCAACATCTTCGCCGCCTGGCCGAGATATTTTTTTGCCGGTCCTAAGTTACCAAACATAGTTCACTCTGATAAGAACAAAGCCCCGCCGAAGCGGGGCTAAACACGGTTAGTGGTGCGAAGAGATCTTCAcgccg  <  2:334062/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       GAACAACAGATGCAGCCGTTGGTCAGCGTTTTGATCTGCGTGGCGCGGTCGCCAATCAGTTGATCATCAACAGAGACTTCGCCGAATTCGTTTTCAATCACGGCAATCTTGTAGCCGTGCTGTTCGTTAAGAATATGGCGCAGCAGGGTGGTTTTTCCTGCGCCTAAAAAACCGGTGAGTAGGGTAACTGCAATCGGGTTCATGCTCTCTCCATTAACAACAGCGCATACCGCCTTTACCATCGCCGCCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATAGCTCATGTACGGCTTGTCGGGATGGTTGGTCTTCATATGCTCGACGTAGTTGTCGTAGTCTGGAATGCCAATCAACATCTTCGCCGCCTGGCCGAGATATTTTTTTGCCTGTCCTAAGTTACCAAACATAGTTCACTCTGATAAGAACAAAGCCCCGCCGAAGCGGGGCTAAACACGGTTAGTGGTGCGAAGAGATCTTCACGCCG  >  NC_000913/4588663‑4589156

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