Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,028,546	G→A	A156V (GCC→GTC)	ygfS ←	putative electron transport protein YgfS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,028,546	0	G	A	92.1%	117.1 / 2.0	39	A156V (GCC→GTC)	ygfS	putative electron transport protein YgfS
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (16/19);  total (18/21)
Fisher's exact test for biased strand distribution p-value = 5.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

CCTGGTGTGTGCCGTGTTCATCGTTCTGACTTGTAATATCGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGGCAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAATACAACGGGCAGAATTGGTCTGCACCAC  >  NC_000913/3028311‑3028781                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                            ccTGGTGTGTGCCGTGTTCATCGTTCTGACTTGTAATATCGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACgctgct                                                                                                                                                                                                                              >  1:102627/1‑251 (MQ=255) ccTGGTGTGTGCCGTGTTCATCGTTCTGACTTGTAATATCGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACgctgct                                                                                                                                                                                                                              <  1:271507/251‑1 (MQ=255)              gtgtTCATCGTTCTGACTTGTAATATCGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATTGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCt                                                                                                                                                                                                                 <  2:249279/251‑1 (MQ=255)                                 gTAATATCGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTCCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAg                                                                                                                                                                                               <  2:291199/250‑1 (MQ=255)                                 gTAATATCGTTACCGGTATTATGCTGAGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCTCGCCTTTCAGCGCGGAGATTTCTGTCAGAGAGGATTCTCACCGCTGACAACAATACGCGGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAg                                                                                                                                                                                               <  1:491556/250‑1 (MQ=255)                                        cGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGc                                                                                                                                                                                        >  2:91235/1‑250 (MQ=255)                                              cGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGccc                                                                                                                                                                                 >  2:317226/1‑250 (MQ=255)                                                  attatGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTgtgt                                                                                                                                                                             <  2:102627/251‑1 (MQ=255)                                                       gCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTggg                                                                                                                                                                         <  1:210648/250‑1 (MQ=255)                                                           ggCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGGCAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCa                                                                                                                                                                       >  1:66123/1‑248 (MQ=255)                                                                   gACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACg                                                                                                                                                             <  1:252805/250‑1 (MQ=255)                                                                    acacTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCAGTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCa                                                                                                                                                           <  1:374751/251‑1 (MQ=255)                                                                       cTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGc                                                                                                                                                        >  1:199391/1‑251 (MQ=255)                                                                               aGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATAGTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGcccc                                                                                                                                                 <  1:228553/250‑1 (MQ=255)                                                                                 gCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCctct                                                                                                                                              <  1:302053/251‑1 (MQ=255)                                                                                     ggTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTccc                                                                                                                                          <  2:192642/250‑1 (MQ=255)                                                                                            gCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGtt                                                                                                                                    <  2:262847/250‑1 (MQ=255)                                                                                                      ggCAAAAGCTGAATAGTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATc                                                                                                                          <  1:440086/250‑1 (MQ=255)                                                                                                                  atatTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAAtt                                                                                                             >  2:496050/1‑251 (MQ=255)                                                                                                                            ggTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCACCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCggg                                                                                                   >  1:541424/1‑251 (MQ=255)                                                                                                                                  cATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATccc                                                                                             >  1:133246/1‑251 (MQ=255)                                                                                                                                       cTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGa                                                                                                                                                                                                                        <  1:14135/123‑1 (MQ=255)                                                                                                                                       cTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGa                                                                                                                                                                                                                        >  2:14135/1‑123 (MQ=255)                                                                                                                                         gCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGc                                                                                        <  1:64265/249‑1 (MQ=255)                                                                                                                                           cGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGCGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACg                                                                                    >  2:369023/1‑251 (MQ=255)                                                                                                                                             cACTGGTCGCATTTTACGCGGGTTGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTGTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGAc                                                                                  <  1:317226/251‑1 (MQ=255)                                                                                                                                                             cgcgGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTc                                                                  >  2:500262/1‑251 (MQ=255)                                                                                                                                                             cgcgGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTc                                                                  <  1:297504/251‑1 (MQ=255)                                                                                                                                                                   tggtTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCACTGGCAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCACTGTTCACAGAGATAGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGGGATCATTCCAAACg                                                            >  2:265712/1‑251 (MQ=255)                                                                                                                                                                        gggCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGgcac                                                       >  2:190218/1‑251 (MQ=255)                                                                                                                                                                                 aaTAGATTCTCCGCGCCTTTCAGCGCGGTGATTTCTTTCAGAGAGGATTCTCACCGCTGCCAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCCCAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACt                                                  <  1:554017/247‑1 (MQ=255)                                                                                                                                                                                    aGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCa                                            >  1:391463/1‑250 (MQ=255)                                                                                                                                                                                     gATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAAc                                          <  1:201462/251‑1 (MQ=255)                                                                                                                                                                                                     cAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCACTGGCAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCACTGTTCACAGAGATAGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGatgctgtccagtcgctatacctt                          <  1:265712/251‑23 (MQ=255)                                                                                                                                                                                                           ggAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAATACAACGGGc                    >  1:518541/1‑251 (MQ=255)                                                                                                                                                                                                             agaTTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAATACAACGGGCAg                  <  1:91235/251‑1 (MQ=255)                                                                                                                                                                                                                       cAGAGAGGATTCTCACAGCTGACAACAATACGCTGCTGGCGGAGTCGCCTGAGTTATCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAATACAACGGGCAGAATTGGTCTg        >  1:510195/1‑251 (MQ=255)                                                                                                                                                                                                                           gagGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAATACAACGGGCAGAATTGGTCTGcacc    >  1:376697/1‑251 (MQ=255)                                                                                                                                                                                                                             ggATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAATACAACGGGCAGAATTGGTCTGcaccac  >  1:464124/1‑251 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                            CCTGGTGTGTGCCGTGTTCATCGTTCTGACTTGTAATATCGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGGCAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAATACAACGGGCAGAATTGGTCTGCACCAC  >  NC_000913/3028311‑3028781

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