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	97.7%	146.6 / NA	44	A156V (GCC→GTC)	ygfS	putative electron transport protein YgfS
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (24/19);  minor base C (1/0);  total (25/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.91e-01

ATGGCTGGCCTGGTGTGTGCCGTGTTCATCGTTCTGACTTGTAATATCGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGGCAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAATACAACGGGCAGAATTGGTCTGCACCACCTGCT  >  NC_000913/3028303‑3028786                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 aTGGCTGGCCTGGTGTGTGCCGTGTTCATCGTTCTGACTTGTAATATCGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATCTTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGCCAACAAt                                                                                                                                                                                                                                           >  2:328340/1‑251 (MQ=255)     ctggcCTGGTGTGTGCCGTGTTCATCGTTCTGACTTGTAATATCGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACgc                                                                                                                                                                                                                                       >  1:315138/1‑251 (MQ=255)       ggcCTGGTGTGTGCCGTGTTCATCGTTCTGACTTGTAATATCGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACgctg                                                                                                                                                                                                                                     >  2:401833/1‑251 (MQ=255)                                 tCTGACTTCTATTATAGTTACCCGTATTATGCTGAGATTTGCGACACTGGTCGTAGGCCGCGCCTTTGCACGCGAACGGCCAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTCGTCGCATTTTACGCGGGGGGTTGGGCAAACTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTATTTCAGAGAGGATTCTCACCGCTGACACCAATTCGCTGCTGGCGGACTCGCCTTAGTTCTCttt                                                                                                                                                                                                           >  2:362052/1‑251 (MQ=255)                                   tGACTTGTAATATCGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTTAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTACTGCCGGACTCGCCTGCGTTCTCTTTcc                                                                                                                                                                                                         >  1:399141/1‑250 (MQ=255)                                         gTAATATCGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAg                                                                                                                                                                                                    <  2:75645/250‑1 (MQ=255)                                         gTAATATCGTTACCGGTATTATGCTGGGATTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGc                                                                                                                                                                                                   >  2:153183/1‑251 (MQ=255)                                            atatCGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAAc                                                                                                                                                                                                >  1:172981/1‑251 (MQ=255)                                                 gTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTGTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTGTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAg                                                                                                                                                                                           <  1:74513/251‑1 (MQ=255)                                                         tattatGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCtgtg                                                                                                                                                                                   >  1:148163/1‑251 (MQ=255)                                                             atGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTg                                                                                                                                                                                <  2:148163/250‑1 (MQ=255)                                                                cTGGGCTTTGTGACACTGGACGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCATTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGc                                                                                                                                                                             >  2:342353/1‑250 (MQ=255)                                                                      tttGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTc                                                                                                                                                                      >  2:382735/1‑251 (MQ=255)                                                                           gACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACGGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATGCTCACCGCTGACAACAATACGCTGCTGGCGGCCTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACg                                                                                                                                                                  <  1:68115/250‑1 (MQ=255)                                                                            acacTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAGTATTGAAACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCATTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGCATCGCCTGAGTTCGCTTTCAGTAAGCAACTGCAGAGCCTGTGTGGGGCAAGATTCAACGcc                                                                                                                                                                >  2:554970/1‑250 (MQ=255)                                                                             cacTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGc                                                                                                                                                                 >  1:228579/1‑249 (MQ=255)                                                                                             cgTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTccc                                                                                                                                               <  1:528743/251‑1 (MQ=255)                                                                                              gTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCg                                                                                                                                              <  2:42849/251‑1 (MQ=255)                                                                                                 tttGAACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACGGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCGCCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCt                                                                                                                                            <  2:245693/250‑1 (MQ=255)                                                                                                  ttGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTg                                                                                                                                           <  1:242170/250‑1 (MQ=255)                                                                                                           aaTGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACag                                                                                                                                    >  1:532555/1‑248 (MQ=255)                                                                                                                  aaaGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCAtt                                                                                                                          >  1:177569/1‑251 (MQ=255)                                                                                                                      cTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTcac                                                                                                                      >  1:63674/1‑251 (MQ=255)                                                                                                                       tGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTcac                                                                                                                      >  2:266523/1‑250 (MQ=255)                                                                                                                                                      aCTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACt                                                                                      <  1:382735/251‑1 (MQ=255)                                                                                                                                                         ggTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGa                                                                                    <  1:220132/250‑1 (MQ=255)                                                                                                                                                                ttttACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGAt                                                                            >  2:43719/1‑251 (MQ=255)                                                                                                                                                                     cgcgGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTc                                                                       <  2:460334/251‑1 (MQ=255)                                                                                                                                                                      gcgGGTGGTTGGGCAATATAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGAAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCACGGCCACGACTGAATGGTGAGCATTcc                                                                      >  2:533811/1‑251 (MQ=255)                                                                                                                                                                       cgGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCa                                                                     <  1:401833/251‑1 (MQ=255)                                                                                                                                                                                   cAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGGCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGgcacgca                                                         <  1:62609/251‑1 (MQ=255)                                                                                                                                                                                      tCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTg                                                      >  1:201448/1‑251 (MQ=255)                                                                                                                                                                                                tctcCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAActc                                             <  1:535248/250‑1 (MQ=255)                                                                                                                                                                                                   ccGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGAGCATTCCAAACGGGCACGCACTGACGCAACTCTGAc                                         <  2:64915/251‑1 (MQ=255)                                                                                                                                                                                                       gcCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACGGACGCAACTCTGACAGcc                                     >  2:155166/1‑251 (MQ=255)                                                                                                                                                                                                         cTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACACTACGCTGCTGGCGGAATCGCCTGAGTTCTCTTTCAGTAAACAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCGCGCTGTCAACAGAGATCGCATTTCACAATTTGTTGTCGGGTATACCAAGGCAACGCCTGAATGGTGATCATTCCAAACGCGCACGCACTGCCGCAACTATGACAGCCa                                    >  1:351612/1‑250 (MQ=255)                                                                                                                                                                                                             taGCGCGGAGATTTCTTTCAAAGAGGATTCTCCCCGCCGACAACAATTCGCTGCTGGGGGGCTCGCCTTAGTTCTCTTTCAATAAGCAACTGCAGCGCCTGTGTTGGGGAAGATTCAACCCAGGCTGGCCCCTCTTCCCCCTGTTCACAGAGAGCGCATTTCACAATTTGTTGGCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAATAc                                <  2:397604/249‑1 (MQ=255)                                                                                                                                                                                                                cgcgGAGATTTCTTTCAGAGAGGCTTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAAt                                                                                  >  1:147096/1‑197 (MQ=255)                                                                                                                                                                                                                cgcgGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCGGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTGCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAAt                                                                                  <  2:147096/197‑1 (MQ=255)                                                                                                                                                                                                                  cgGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAATACAACggg                          >  1:182575/1‑251 (MQ=255)                                                                                                                                                                                                                        tttctttcAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAATACAACGGGCAGAAt                    >  2:161982/1‑251 (MQ=255)                                                                                                                                                                                                                               cAGAGAGGATTCGCACCGCTGACACCAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTGCAGTAAGCAACTGCAGCGCCTGTTTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATGTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACCCTGACAGCCAATACAACTGGCAGAATTGGTCTg             <  2:114377/251‑1 (MQ=255)                                                                                                                                                                                                                                      gATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAATACAACGGGCAGAATTGGTCTGCaccacc      <  1:161982/251‑1 (MQ=255)                                                                                                                                                                                                                                          ctcACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAATACAACGGGCAGAATTGGTCTGCACCACCTGCt  >  1:501905/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 ATGGCTGGCCTGGTGTGTGCCGTGTTCATCGTTCTGACTTGTAATATCGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGGCAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAATACAACGGGCAGAATTGGTCTGCACCACCTGCT  >  NC_000913/3028303‑3028786

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