Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,566,642	G→A	S202L (TCG→TTG)	yjiN ←	DUF445 domain‑containing protein YjiN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,566,642	0	G	A	97.8%	154.1 / ‑3.9	46	S202L (TCG→TTG)	yjiN	DUF445 domain‑containing protein YjiN
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (28/17);  total (28/18)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.13e-01

CGGCTCGCGCTGCCATTTCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGC  >  NC_000913/4566400‑4566884                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   cGGCTCGCGCTGCCATTTCCGGATCGATTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCAAGATCGCGACTAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGAGGGGAGGATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGTCGATAAACTTGCGCAATTTATcg                                                                                                                                                                                                                                            >  1:381461/1‑251 (MQ=255)        agcTGCCATTTCCGGATCGTTTTTCAGGTTGTCGATCAGGGCGAAGGTCGCGCGCTCAAACGCATGGCGGATCTGATGCGCAAGCTCGCGACTAATATCATCAAGCAAAGAGTTCACCGCGTAAGAAACCAACTCCGCGCTATGTTCGCACAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGCTGCTCGCTCTCCAGCCAGCGAAACATTTGCTGGGCGATAAACTTGCGCAATTTATCGCGCTgg                                                                                                                                                                                                                                      >  2:319139/2‑250 (MQ=255)                    cGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCAATTTATCGCGCTGGAGAAGGGCGATc                                                                                                                                                                                                                          >  1:304591/1‑250 (MQ=255)                     ggATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCa                                                                                                                                                                                                                         >  1:43699/1‑250 (MQ=255)                     ggATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCa                                                                                                                                                                                                                         >  1:333937/1‑250 (MQ=255)                            ttttCAGCTTGTCGATCAGGGCGAAGGTGGCGCAATCACACGCATGGCGGATCTGCTGCGCACGATCGCGACTAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCAATTTATCGCGCTGGATAAGGGCGCTCAACTGTGc                                                                                                                                                                                                                  <  2:84287/250‑1 (MQ=255)                                cAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATc                                                                                                                                                                                                              <  1:37116/250‑1 (MQ=255)                                     tGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGt                                                                                                                                                                                                         <  2:253703/250‑1 (MQ=255)                                            cAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCcagc                                                                                                                                                                                                  >  1:403174/1‑250 (MQ=255)                                              gggCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCagcagc                                                                                                                                                                                               >  1:438238/1‑251 (MQ=255)                                                  gAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCt                                                                                                                                                                                           >  1:77525/1‑251 (MQ=255)                                                        ggCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAATCTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCCCCTGATGAc                                                                                                                                                                                     <  2:158387/251‑1 (MQ=255)                                                            cgATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATATCATCAAGCAAAGAATTCACCGAGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATc                                                                                                                                                                                 >  1:114406/1‑251 (MQ=255)                                                                cAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGttt                                                                                                                                                                             >  1:386775/1‑251 (MQ=255)                                                                       aTGGCGGATCTGATGCGCACGATCGCGACTAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAAATCCGCGCTATGTTCGCCCAACCATTCGGGGGGGAGCATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGGCAGCGAACAATTTGCTGGGCGATAAACTTGCGCAATTTATCGCGCTGGAGAAGAGCGATCAACTGTGCGCTCAGCGTATCCCGCAGCACCTGATGACGATCGTTTTTGGTc                                                                                                                                                                       >  2:394189/1‑250 (MQ=255)                                                                                      cgcACGATAGCGCCTAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCAAGCATc                                                                                                                                                        >  2:45854/1‑250 (MQ=255)                                                                                            aTCGCGACTAATATCATCAAGCAAAGAATGCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCAc                                                                                                                                                 <  2:22359/251‑1 (MQ=255)                                                                                                   cTAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGcc                                                                                                                                          >  1:514144/1‑251 (MQ=255)                                                                                                    tAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCa                                                                                                                                                                                                             <  2:242864/183‑1 (MQ=255)                                                                                                    tAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCa                                                                                                                                                                                                             >  1:242864/1‑183 (MQ=255)                                                                                                      atatCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGCTGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCgg                                                                                                                                        >  1:317605/1‑250 (MQ=255)                                                                                                            acaAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATc                                                                                                                                 <  2:53798/250‑1 (MQ=255)                                                                                                              aaGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGa                                                                                                                               <  2:317605/251‑1 (MQ=255)                                                                                                                 cAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCt                                                                                                                            >  1:320975/1‑251 (MQ=255)                                                                                                                      aaTTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATc                                                                                                                        >  1:264679/1‑250 (MQ=255)                                                                                                                            ccGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGcc                                                                                                                 <  2:438238/251‑1 (MQ=255)                                                                                                                                gTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGAt                                                                                                             >  1:217650/1‑251 (MQ=255)                                                                                                                                gTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATCTGCTGGGCGATAAACTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGAt                                                                                                             <  2:90455/251‑1 (MQ=255)                                                                                                                                     aaaCCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGAcc                                                                                                        >  1:185106/1‑251 (MQ=255)                                                                                                                                                 cgcTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAGc                                                                                            >  1:270289/1‑251 (MQ=255)                                                                                                                                                         tCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACGCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAGCAGGCGCt                                                                                     <  2:114406/250‑1 (MQ=255)                                                                                                                                                                  cATTCGGTGGGGAGACTTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGCGCGATAAACTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCCACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCACCGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAGCAGGCGCTGAATACgcgc                                                                           >  2:250890/1‑251 (MQ=255)                                                                                                                                                                                      ggCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGaa                                                       <  2:270289/251‑1 (MQ=255)                                                                                                                                                                                               tGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTc                                               >  2:384839/1‑250 (MQ=255)                                                                                                                                                                                                  aTGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTGTGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAGCCGCTCATGa                                           <  2:191265/251‑1 (MQ=255)                                                                                                                                                                                                       cgctcTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGc                                      <  1:394189/251‑1 (MQ=255)                                                                                                                                                                                                           ctcCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAAc                                   >  1:186390/1‑250 (MQ=255)                                                                                                                                                                                                             ccagccagCGAACAATTTGCTGGGCGATAAACGTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAg                                 <  2:430355/250‑1 (MQ=255)                                                                                                                                                                                                                      gAACAATTTGCTGGGCGATAAACTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGAcc                        >  1:381367/1‑250 (MQ=255)                                                                                                                                                                                                                          aaTTTGCTGGGCGATAAACTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAAcgc                   <  2:403174/251‑1 (MQ=255)                                                                                                                                                                                                                             ttGCTGGGCGATAAACTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACgcggc                >  2:252932/1‑251 (MQ=255)                                                                                                                                                                                                                                      gATAAACTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGtttt       >  1:164607/1‑251 (MQ=255)                                                                                                                                                                                                                                      gATAAACTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGtttt       >  1:444956/1‑251 (MQ=255)                                                                                                                                                                                                                                          aaCTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTa                                                                                                                          >  2:9156/1‑132 (MQ=255)                                                                                                                                                                                                                                          aaCTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTa                                                                                                                          <  1:9156/132‑1 (MQ=255)                                                                                                                                                                                                                                           aCTTGCGCAATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTctggc  <  2:380854/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   CGGCTCGCGCTGCCATTTCCGGATCGTTTTTCAGCTTGTCGATCAGGGCGAAGGTGGCGCGATCAAACGCATGGCGGATCTGATGCGCACGATCGCGACTAATATCATCAAGCAAAGAATTCACCGCGTCAGAAACCAACTCCGCGCTATGTTCGCCCAACCATTCGGTGGGGAGAATTTTGGCTTTCAGTGGATGCTCGCTCTCCAGCCAGCGAACAATTTGCTGGGCGATAAACTTGCGCGATTTATCGCGCTGGAGAAGGGCGATCAACTGTGCGATCAGCGTATCCAGCAGCACCTGATGACGATCGTTTTTGGTCATACTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGC  >  NC_000913/4566400‑4566884

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