Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,638,567	A→G	E130E (GAA→GAG)	creD →	putative inner membrane protein CreD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,638,567	0	A	G	98.0%	156.5 / NA	51	E130E (GAA→GAG)	creD	putative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (31/19);  minor base T (1/0);  total (32/19)
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.99e-01

AGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTAC  >  NC_000913/4638325‑4638794                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                     aGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGCTGTTTCGCGTCTTAGCGAGCTCAAcgc                                                                                                                                                                                                                              >  2:331689/1‑251 (MQ=255)    gggCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGcc                                                                                                                                                                                                                           >  1:452112/1‑251 (MQ=255)        cGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAat                                                                                                                                                                                                                       >  1:323601/1‑251 (MQ=255)              aaCTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCatc                                                                                                                                                                                                                 >  2:332330/1‑249 (MQ=255)                   gTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTAAACGCGCCAAATATCATCTTAgg                                                                                                                                                                                                            >  2:41693/1‑251 (MQ=255)                            cttaTCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGccctt                                                                                                                                                                                                   >  2:270342/4‑248 (MQ=255)                                        ccTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGA‑GAACGCAAGATAGGGATTTATACCTGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGc                                                                                                                                                                                      <  1:171384/251‑1 (MQ=255)                                            tgACCGAGCTTTATACGGTGCAGGAAGAGGATACAACCGTGGAGCGGAAACGAAGTTTTCTCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAGTCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCAATTTCGAGGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGCCAAGCCCTTTATTGTGATTAGCGTc                                                                                                                                                                                   >  1:435018/1‑251 (MQ=255)                                                   gCTTTATACTGTGAAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCGTTATTGTGATTAGCGTAGGGGATg                                                                                                                                                                            <  2:31723/251‑1 (MQ=255)                                                            ggTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGCAGTTTTATCCATTTTTGGTTACCTGAGTCATTGCTGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGCTGCGCG‑GGGt                                                                                                                                                                   >  2:466985/1‑250 (MQ=255)                                                              ggCAGTAAGAGGATAAAGCGGTGGATCGGAAACGAAGTTGTATCCAGTTTTGGTTAGCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGAAAGATAGGGATTTATCCCGGTCAGGTCTGGCACAGGGATTTAACGCTAAAAGCTGATTTCGATGTTTCGCGTCGTAGCGAGCTCAACGCGGCAAATAGCATCTTAGGCAATCCCTTTATTGTGATTAGCGTCGGGGAGGCGCT‑TGGTagt                                                                                                                                                                <  2:412619/250‑3 (MQ=255)                                                                cAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATtgg                                                                                                                                                              >  1:122133/1‑251 (MQ=255)                                                                  ggAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGCGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATtggtg                                                                                                                                                            >  2:3361/1‑251 (MQ=255)                                                                       agGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATtggtgtggtg                                                                                                                                                       <  1:362077/251‑1 (MQ=255)                                                                             aaaCCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAgcg                                                                                                                                                 <  1:339553/251‑1 (MQ=255)                                                                              aaCCGTGGAGCGGTAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGATGAAGGCAATCAGAACGTAGCAGAACGCAAGATAGGGATTTATAACGGTCAGGTCTGGCAAAGTGATTTACCGCTAAAAGCCAATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATAATCTTAGTCAAGACCTTTTTTGTGATTAGCGTCGGCGATGCTCG‑TGATTTTGGTGAGGAGAATgcgc                                                                                                                                                >  2:266832/1‑251 (MQ=255)                                                                                         ggAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGGGTGGTGGAAGGGCCTggaggtaa                                                                                                                                      >  2:305282/1‑243 (MQ=255)                                                                                          gAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGAt                                                                                                                                                                             <  1:212420/211‑1 (MQ=255)                                                                                          gAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGAt                                                                                                                                                                             >  2:212420/1‑211 (MQ=255)                                                                                            aaggAAGTTTTATCCAATTTTGGTTACATGAGTCATTGATGGTTGATGGCAATCAGAACGTGGTAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTGACGCTAAAAGCCGATTTCGATGTTTCGCGTCTAAGCGAGCTCAACGCGCCAACTATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGAGTTGGGGTGGTGAAGGCGCCTGCAGTTaacgg                                                                                                                                   >  2:28315/4‑250 (MQ=255)                                                                                             acggaGTTTTATACATTTTTGGTTACCGGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGCACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCAAAGTAATTTAAAGATAAAGGCCTATTTCGATGTTTAGTGTCTAAGAGTGCTCAACGCGACAAATATCATCTTAGGAAAGACCTTTATTGTGCTTAGAGTCGGGGGTGCGCGGTGG‑ATTGGTGAGGTAAAAgagacgagatttagctga                                                                                                                                  >  2:261412/5‑233 (MQ=255)                                                                                               gAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTaacggaacg                                                                                                                               <  1:332330/251‑1 (MQ=255)                                                                                                     ttATCAATTTTTGGTTACATGAGTCATTGATGGTTGATGGCAATAAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGGAAGAGCCTGAAGTTACCGGAACGGCGCTg                                                                                                                         >  2:238391/1‑251 (MQ=255)                                                                                                      tATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTgg                                                                                                                        >  1:183988/1‑250 (MQ=255)                                                                                                              tttGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGAGTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGa                                                                                                                 <  2:452112/250‑1 (MQ=255)                                                                                                                   tcaCCTGAGACATTGATGGTTGATGGCAATCAGAACGTGGTATAACGAAAGATAGGGATTTATACCGGTCAGGTCTGTCACAGTGGTTTAAAGATTAAAGACGATTTCGATGTTTAGCGTCTTAGCGATCTCAACGCGCAAAATATCATCATAGGCAAGCCCTTTATTGAGATTAGCGTCGGGGGTGAGAG‑TGGTATTGGTGTGGTGAAAGCGACGGGAGTTaacggagcagcgcataccatttaaacggg                                                                                                           >  2:33505/3‑228 (MQ=255)                                                                                                                    tACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGGGGTGAAAGCGCCTGAAGTTAACGGAACGGCCCTGACCATTGGACCCGGc                                                                                                          >  2:91092/1‑251 (MQ=255)                                                                                                                      ccTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGTAGTTAACGGAACGGCGCTGACCATTGAACCCggaac                                                                                                        >  2:2537/1‑248 (MQ=255)                                                                                                                             cATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑GGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGtt                                                                                                  >  2:474511/1‑250 (MQ=255)                                                                                                                               ttgatggttgatggCAATCAGCACGTGGAAGCACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGAGATTTAACTCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGCTTAGCGTCGGGGAATCGCG‑TGGGATTGGTGTGGTGAAAGAGCCTGAAGTTAACGGGACGGCGCTGGCCCATGAGCACGGCACCGGGtgtgg                                                                                               >  2:302763/1‑247 (MQ=255)                                                                                                                                tgatggttgatggCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTag                                                                                                >  1:190724/1‑249 (MQ=255)                                                                                                                                 gatggttgatggCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAg                                                                                                                                                                                       >  1:127159/1‑162 (MQ=255)                                                                                                                                 gatggttgatggCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAg                                                                                                                                                                                       <  2:127159/162‑1 (MQ=255)                                                                                                                                          tggCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCgggg                                                                                    >  2:115586/1‑250 (MQ=255)                                                                                                                                           ggCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCa                                                                                   <  1:270342/251‑1 (MQ=255)                                                                                                                                                       gTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATGGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCTCCGGGTTAGAGCAAGGCGGGCAGGGCGTGCatat                                                                       <  2:435018/251‑1 (MQ=255)                                                                                                                                                           aagaaCGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTTGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTGGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGCATTGGTGTGGTGCAAGCGCCTGAAGTTAACGGACCGGCGCTGACCATTGACCCCGGCACCGGGTTAGAGCAAGGCGGGCCGGGCgggccaatccct                                                                   >  1:354183/1‑240 (MQ=255)                                                                                                                                                            agaaCGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCtt                                                                  <  2:236827/251‑1 (MQ=255)                                                                                                                                                                     gATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAg                                                                                                                                                   <  2:228257/161‑1 (MQ=255)                                                                                                                                                                     gATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAg                                                                                                                                                   >  1:228257/1‑161 (MQ=255)                                                                                                                                                                      aTAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAgg                                                        >  1:77590/1‑251 (MQ=255)                                                                                                                                                                      aTAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATACGCG‑TAGCCTTGGTGTGGTGACAGCCGCTGCAGTTAACGGCACGGCGCTGACCATTGAACCCGGCCACGGGGTACAGCAAGGTGGGGACGGCGTGGTTATCCCTTTAccccaa                                                          >  2:146644/1‑245 (MQ=255)                                                                                                                                                                                ataCCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCgg                                              <  2:277960/251‑1 (MQ=255)                                                                                                                                                                                ataCCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCATCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCgg                                              <  1:331689/251‑1 (MQ=255)                                                                                                                                                                                      gTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCAAATCCCTTTACCCGAAGGGGACTGGAGGAAGCa                                         >  1:360391/1‑250 (MQ=255)                                                                                                                                                                                        cAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCAGTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGaa                                      <  1:302763/251‑1 (MQ=255)                                                                                                                                                                                                gCACAGTGAGTTAACGCTAAAAGCCGAGTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCGTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGc                              <  2:485858/251‑1 (MQ=255)                                                                                                                                                                                                            aaCGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGccc                  >  1:374090/1‑251 (MQ=255)                                                                                                                                                                                                                       ccGAGTTCGATGTTGCGCGGCTTAGCGCGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGGTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGTTGAATATGGCCCTGAATTTacgc       <  2:259949/251‑4 (MQ=255)                                                                                                                                                                                                                          aTTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGt    <  1:186780/251‑1 (MQ=255)                                                                                                                                                                                                                            ttCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGGGCATATCCCTTTCCCCGAAGGGGACTGGCGGGAGCCGAACCTGAAGCTGGAAATGGCCCTGAGTTTAAACggaac  >  2:38607/1‑248 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                     AGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCG‑TGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTAC  >  NC_000913/4638325‑4638794

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