Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	661,069	G→A	L104F (CTT→TTT)	ybeF ←	putative DNA‑binding transcriptional regulator YbeF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	661,069	0	G	A	97.4%	130.7 / NA	39	L104F (CTT→TTT)	ybeF	putative DNA‑binding transcriptional regulator YbeF
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (11/27);  minor base C (1/0);  total (12/27)
Fisher's exact test for biased strand distribution p-value = 3.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACATGATGTTGCACCGTACGATTGGTGCAAAACATGTTATCGATGATGAGATCGGTTTGAAACTGACTGAGTTGGTTTTCCGCGTCGCTGATGGGCGGGTTGCGCAGTAAAAGCTGCGGATAGTGAGTTTTAATCGCCCGATAGATGACAGGAAGGACCAGTGCTCCGACCGAGGGAGTGGTCGCAATAGTTATCGTTCGTTGCTTATCATAGCTTCCTTCGATATCCAGCGCGCCAAGAATGGACTCAAGGCCCTGACTGATATACTCATGTAGATGCATCGCAAATGCGGTAGGAGTGACACCCTGGCCTTTGCGAATAAACAATGGGTCAGGGAATATAACGCGCAGTTTCTGAATAGACTGACTGATTGCCGAGGGGGTCAGATTAAGCACTTTCGCTGCATTAACGATCCCTTTATGTACATATACAGCTTCAAAAATAGTCAGAAGGTTAAGATCAATATTTCGTAAGGTTGT  >  NC_000913/660831‑661309                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 acaTGATGTTGCACCGTACGATTGGTGCAAAACATGTTATCGATGATGAGATCGGTTTGAAACTGACTGAGTTGGTTTTCCGCGTCGCTGATGGGCGGGTTGCGCAGTAAAAGCTGCGGATAGTGAGTTTTAATCGCCCGATAGATGACAGGAAGGACCAGTGCTCCGACCGAGGGAGTGGTCGCAATAGTTATCGTTCGTTGCTTATCATAGCTTCCTTCGATATCCAGCGCGCCAAAAATGGACTCAAg                                                                                                                                                                                                                                      <  1:432414/251‑1 (MQ=255)                gTACGATTGGTGCAAAACATGTTATCGATGATGAGATCGGTTTGAAACTGACTGAGTTGGTTTTCCGCGTCGCTGATGGGCGGGTTGCGCAGTAAAAGCTGCGGATAGTGAGTTTTAATCGCCCGATAGATGACAGGAAGGACCAGTGCTCCGACCGAGGGAGTGGTCGCAATAGTTATCGTTCGTTGCTTATCATAGCTTCCTTCGATATCCAGCGCGCCAAAAATGGACTCAAGGCCCTGACTGAtata                                                                                                                                                                                                                       <  2:392595/251‑1 (MQ=255)                  atgATTGGTGCAAAACATGTTATCGATGATGAGATCGGTTTGAAACTGACTGAGTTGGTTTTCCGCGTCGCTGATGGGCGGGTTGCGCAGTAAAAGCTGCGGATAGTGAGTTTTAATCGCCCGATAGATGACAGGAAGGACCAGTGCTCCGACCGAGGGAGTGGTCGCAATAGTTATCGTTCGTTGCTTATCATAGCTTCCTTCGATATCCAGCGCGCCAAAAATGGACTCAAGGCCCTGACTGATATAc                                                                                                                                                                                                                      <  2:334117/248‑1 (MQ=255)                        ggTGCAAAACATGTTATCGATGATGAGATCGGTTTGAAACTGACTGAGTTGGTTTTCCGCGTCGCTGATGGGCGGGTTGCGCAGTAAAAGCTGCGGATAGTGAGTTTTAATCGCCCGATAGATGACAGGAAGAACCAGTGCTCCGACCGAGGGAGTGGTCGCAATAGTTATCGTTCGTTGCTTATCATAGCTTCCTTCGATATCCAGCGCGCCAACAATGGACTCAAGGCCCTGACTGATATACTCATGt                                                                                                                                                                                                                >  1:336735/1‑250 (MQ=255)                           gCAAAACATGTTATCGATGATGAGATCGGTTTGAAACTGACTGAGTTGGTTTTCCGCGTCGCTGATGGGCGGGTTGCGCAGTAAAAGCTGCGGATAGTGAGTTTTAATCGCCCGATAGATGACAGGAAGGACCAGTGCTCCGACCGAGGGAGTGGTCGCAATAGTTATCGTTCGTTGCTTATCATAGCTTCCTTCGATATCCAGCGCGCCAAAAATGGACTCAAGGCCCTGACTGATATACTCATGTAGa                                                                                                                                                                                                             <  1:138781/250‑1 (MQ=255)                              aaaCATGTTATCGATGATGAGATCGGTTTGAAACTGACTGAGTTGGTTTTCCGCGTCGCTGATGGGCGGGTTGCGCAGTAAAAGCTGCGGATAGTGAGTTTTAATCGCCCGATAGATGACAGGAAGGACCAGTGCTCCGACCGAGGGAGTGGTCGCAATAGTTATCGTTCGTTGCTTATCATAGCTTCCTTCGATAGCCAGCGCGCCAAAAAGGGACTCAAGGCCCTGACTGATATACTCATGTAGATGCa                                                                                                                                                                                                         <  1:351630/251‑1 (MQ=255)                                 cATGTTATCGATGATGAGATCGGTTTGAAACTGACTGAGTTGGTTTTCCGCGTCGCTGATGGGCGGGTTGCGCAGTAAAAGCTGCGGATAGTGAGTTTTAATCGCCCGATAGATGACAGGAAGGACCAGTGCTCCGACCGAGGGAGTGGTCGCAATAGTTATCGTTCGTTGCTTATCATAGCTTCCTTCGATATCCAGCGCGCCAAAAATGGACTCAAGGCCCTGACTGATATACTCATGTAGATGCATCg                                                                                                                                                                                                      >  1:164944/1‑251 (MQ=255)                                           atgatgAGATCGGTTTGAAACGGACTGAGTTGGTTTTCCGCGTCGCTGATGGGCGGGTTGCGCAGTAAAAGCTGCGGATAGTGAGTTTTAATCGCCCGATAGATGACAGGAAGGACCAGTGCTCCGACCGAGGGAGTGGTCGCAATAGTTATCGTTCGTTGCTTATCATAGCTTCCTTCGATATCCAGCGCGCCAAAAATGGACTCAAGGCCCTGACTGATATACTCATGTAGATGCATCGCAAATGCGGt                                                                                                                                                                                            <  1:250758/251‑1 (MQ=255)                                            tgatgaGATCGGTTTGAAACTGACTGAGTTGGTTTTCCGCGTCGCTGATGGGCGGGTTGCGCAGTAAAAGCTGCGGATAGTGAGTTTTAATCGCCCGATAGATGACAGGAAGGACCAGTGCTCCGACCGAGGGAGTGGTCGCAATAGTTATCGTTCGTTGCTTATCATAGCTTCCTTCGATATCCAGCGCGCCAAAAATGGACTCAAGGCCCTGACTGATATACTCATGTAGATGCATCGCAAATGCGGt                                                                                                                                                                                            >  2:457550/1‑250 (MQ=255)                                                    tCGGTTTGAAACTGACTGAGTTGGTTTTCCGCGTCGCTGATGGGCGGGTTGCGCAGTAAAAGCTGCGGATAGTGAGTTTTAATCGCCCGATAGATGACAGGAAGGACCAGTGCTCCGACCGAGGGAGTGGTCGCAATAGTTATCGTTCGTTGCTTATCATAGCTTCCTTCGATATCCAGCGCGCCAAAAATGGACTCAAGGCCCTGACTGATATACTCATGTAGATGCATCGCAAATGCGGTAGGAGGGac                                                                                                                                                                                   >  2:344503/1‑251 (MQ=255)                                                      ggTTTGAAACTGACTGAGTTGGTTTTCCGCGTCGCTGATGGGCGGGTTGCGCAGTAAAAGCTGCGGATAGTGAGTGTTAATCGCCCGATAGATGACAGGAAGGACCAGTGCTCCGACCGAGGGAGTGGTCGCAATAGTTATCGTTCGTTGCTTATCATAGCTTCCTTCGATATCCAGCGCGCCAAAAATGGACTCAAGGCCCTGACTGATATACTCATGTAGATGCATCGCAAATGCGGTAGGAGTGacac                                                                                                                                                                                 <  2:97765/251‑1 (MQ=255)                                                          tGAAACTTACTGAGTTGGTTTTCCGCGTCGCTGATGGGCGGGTTGCGCAGTAAAAGCTGCGGATAGTGAGTTTTAATCGCCCGATAGATGACAGGAAGGACCAGTGCTCCGACCGAGGGAGTGGTCGCAATAGTTATCGTTCGTTTCTTATCATAGCTTCCTTCGATATCCAGCGCGCCAAAAATGGCCTCAAGGCCCTGACTGATATACTCATGTAGATGCCTCGCAAATGCGGTAGGAGTGACCCCCt                                                                                                                                                                              <  2:159720/250‑1 (MQ=255)                                                                          ggTTTTCCGCGTCGCTGAGGGGCGGGTTGCGCAGGAAAAGCTGCGGATAGTGAGTTTTAATCGCCCGATAGATGACAGGAAGGGCCAGTGCTCCGACCGAGGGAGTGTTCGCAATAGTTATCGTTCGTTGCTTATCATAGCTTCCTTCGATATCCAGCGCGCCAAAAATGGACTCAAGGCCCTGACTGATATACGCATGTAGATGCATCGCAAATGCGGTAGGAGTGACACCCTGGCCTTTGCGAATAAAc                                                                                                                                                             <  2:131049/251‑1 (MQ=255)                                                                              ttCCGCGTCGCTGATGGGCGGGTTGCGCAGTAAAAGCTGCGGATAGTGAGTTTTAATCGCCCGATAGATGACAGGAAGGACCAGTGCTCCGACCGAGGGAGTGGTCGCAATAGTTATCGTTCGTTGCTTATCATAGCTTCCTTCGATATCCAGCGCGCCAAAAATGGACTCAAGGCCCTGACTGATATACTCATGTAGATGCATCGCAAATTCGGTAGGAGTGACACCCTGGCCGTTGCGAATAAACAAt                                                                                                                                                          <  1:394011/250‑1 (MQ=255)                                                                              ttCCGCGTCGCTGATGGGCGGGTTGCGCAGTAAAAGCTGCGGATAGTGAGTTTTAATCGCCCGATAGATGACAGGAAGGACCAGTGCTCCGACCGAGGGAGTGGTCGCAATAGTTATCGTTCGTTGCTTATCATAGCTTCCTTCGATATCCAGCGCGCCAAAAATGGACTCAAGGCCCTGACTGATATACTCATGTAGATGCATCGCAAATGCGGTAGGAGTGACACCCTGGCCTTTGCGAATAAACAAt                                                                                                                                                          <  1:275035/250‑1 (MQ=255)                                                                                     tCGCTGATGGGCGGGGTGCGCAGTAAAAGCTGCGGATAGTGAGTATTAATCGCCAGATAGATGACAGGAAGGACCAGTGCTCCGACCGAGGGAGTGGTCGCAATAGTTATCGTTCGTTGCTTATCATAGCTTCCTTCGATATCAAGCGCGCCAAAAATGGACTCAAGGCCCTGACTGATATACTCATGTAGATGCATCGCAAATGCGGTAGGAGTGACACCCTGGCCTTTGCGAATAAACAATGGGTCAgg                                                                                                                                                  >  1:24189/1‑251 (MQ=255)                                                                                           aTGGGCGGGTTGCGCAGTAAAAGCTGCGGATAGTGAGTTTTAATCGCCCGATAGATGACAGGAATGACCAGTGCTCCGACCGAGGGAGTGGTCGCAATAGTTATCGTTCGTTGCTTATCATAGCTTCCTTCGATATCCAGCGCGCCAAAAATGGACTCAAGGCCCTGACTGATATACTCATGTAGATGCATCGCAAATGCGGTAGGAGTGACTCCCTGGCCTTTGCGAATAAACAATGGGTCAGGGAatat                                                                                                                                            <  1:242054/251‑1 (MQ=255)                                                                                                  ggTTGCGCAGTAAAAGCTGCGGATAGTGAGTTTTAATCGCCCGATAGATGACAGGAAGGACCAGTGCTCCGACCGAGGGAGTGGTCGCAATAGTTATCGTTCGTTGCTTATCATAGCTTCCTTCGATATCCAGCGCGCCAAAAATGGACTCAAGGCCCTGACTGATATACTCATGTAGATGCATCGCAAATGCGGTAGGAGTGACACCCTGGCCTTTGCGAATAAACAATGGGTCAGGGAATATAACgcgc                                                                                                                                     >  2:315381/1‑251 (MQ=255)                                                                                                   gTTGCGCAGTAAAAGCGGCTGATAGTGAGTTTTAAGCGCCCGATAGATGACAGGAAGGACCAGGGCTCCGACCGAGGGAGTGGTCGCAATAGTTATCGTTCGTTGCTTATCATAGCTTCCTTCGATATCCAGCGCGCCAAAAATGGACTCAAGGCCCTGACTGATATACTCATGTAGATGCATCGCAAATGCGGTAGGAGTGACACCCTGGCCTTGGCGAATAAACAATGGGTCAGGGAATATAACGCGCa                                                                                                                                    <  2:242606/251‑1 (MQ=255)                                                                                                                    gCGGATAGTGAGTTTTAATCGCCCGATAGATGACAGGAAGGACCAGTGCTCCGACCGAGGGAGTGGTCGCAATAGTTATCGTTCGTTGCTTATCATAGCTTCCTTCTATATCCAGCGCGCCAAAAATGGACTCAAGGCCCTGACTGATATACTCATGTAGATGCATCGCAAATGCGGTAGGAGTGACACCCTGGCCTTTGCGAATAAACAATGGGTCAGGGAATATAACGCGCAGTTTCTGAATAgact                                                                                                                     <  1:275922/249‑1 (MQ=255)                                                                                                                      ggATAGTGAGTTTTAATCGCCCGATAGATGACAGGAAGGACCAGTGCTCCGACCGAGGGAGTGGTCGCAATAGTTATCGTTCGTTGCTTATCATAGCTTCCTTCGCTATCCAGCGCGCCAAAACTGGACTCAAGGCCCTGAATGATATACTCATGTAGATGCATCGCAAATGCGGTAGGAGTGACACCCTGGCCTTTGCGAATAAACAATGGGACAGGGAATATAACGCGCAGTTTCTGAATAgactgact                                                                                                                 >  2:293898/1‑251 (MQ=255)                                                                                                                        aTAGTGAGTTTTAATCGCCCGATCGATGACAGGAAGGACCAGTGCTCCGACCGAGGGAGTGGTCGCAATAGTTATCGTTCGTTGCTTATCATAGCTTCCTTCGATATCCAGCGCGCCAAAAATGGACTCAAGGCCCTGACTGATATACTCATGTAGATGCATCGCAAATGCGGTAGGAGTGACACCCTGGCCTTTGCGAATAAACAATGGGTCAGGGAATATAACGCGCAGTTTCTGAATAGactgactg                                                                                                                >  2:46405/1‑250 (MQ=255)                                                                                                                                     aTCGCCCGATAGATGACAGGAAGGACCAGTGCTCCGACCGAGGGAGTGGTCGCAATAGTTATCGTTCGTTGCTTATCATAGCTTCCTTCGATATCCAGCGCGCCAAAAATGGACTCAAGGCCCTGACTGATATACTCATGTAGATGCATCGCAAATGCGGTAGGAGTGACACCCTGGCCTTTGCGAATAAACAATGGGTCAGGGAATATAACGCGCAGTTTCTGAATAGACTGACTGATTGCCGAGGGGGt                                                                                                  <  1:381040/251‑1 (MQ=255)                                                                                                                                               aGATGACAGGAAGGACCAGTGCTCCGACCGAGGGAGTGGTCGCAATAGTTATCGTTCGTTGCTTATCATAGCTTCCTTCGATATCCAGCGCGCCAAAAATGGACTCAAGGCCCTGACTGATATACTCATGTAGATGCATCGCAAATGCGGTAGGAGTGACACCCTGGCCTTTGCGAATAAACAATGGGTCAGGGAATATAACGCGCAGTTTCTGAATAGACTGACTGATTGCCGAGGGGGTCAGATTAAGc                                                                                        <  2:37609/251‑1 (MQ=255)                                                                                                                                                    aCAGGAAGGACCAGTGCTCCGACCGAGGGAGTGGTCGCAATAGTTATCGTTCGTTGCTTATCATAGCTTCCTTCGATATCCAGCGCGCCAAAAATGGACTCAAGGCCCTGACTGATATACTCATGTAGATGCATCGCAAATGCGGTAGGAGTGACACCCTGGCCTTTGCGAATAAACAATGGGTCAGGGAATATAACGCGCAGTTTCTGAATAGACTGACTGATTGCCGAGGGGGTCAGATTAAGCACttt                                                                                   <  1:382559/251‑1 (MQ=255)                                                                                                                                                            gaCCAGTGCTCCGACCGAGGGAGTGGTCGCAATAGTTATCGTTCGTTGCTTATCATAGCTTCCTTCGATATCCAGCGCGCCAAAAATGGACTCAAGGCCCTGACTGATATACTCATGTAGATGCATCGCAAATGCGGTAGGAGTGACACCCTGGCCTTTGCGAATAAACAATGGGTCAGGGAATATAACGCGCAGTTTCTGAATAGACTGACTGATTGCCGAGGGGGTCAGATTAAGCACTTTCGCTGCAt                                                                           <  1:145285/251‑1 (MQ=255)                                                                                                                                                                         accgaGGGAGTGGTCGCAATAGTTATCGTTCGTTGCTTATCATAGCTTCCTTCGATATCCAGCGCGCCAAAAATGGACTCAAGGCCCTGACTGATATACTCATGTAGATGCATCGCAAATGCGGTAGGAGTGACACCCTGGCCTTTGCGAATAAACAATGGGTCAGGGAATATAACGCGCAGTTTCTGAATAGACTGACTGATTGCCGAGGGGGTCAGATTAAGCACTTTCGCTGCATTAACGATCCCttt                                                              <  2:384876/251‑1 (MQ=255)                                                                                                                                                                         accgaGGGAGTGGTCGCAATAGTTATCGTTCGTTGCTTATCATAGCTTCCTTCGATATCCAGCGCGCCAAAAATGGACTCAAGGCCCTGACTGATATACTCATGTAGATGCATCGCAAATGCGGTAGGAGTGACACCCTGGCCTTTGCGAATAAACAATGGGTCAGGGAATATAACGCGCAGTTTCTGAATAGACTGACTGATTGCCGAGGGGGTCAGATTAAGCACTTTCGCTGCATTAACGATCCCttt                                                              <  1:432836/251‑1 (MQ=255)                                                                                                                                                                          ccgaGGGAGTGGTCGCAATAGTTATCGTTTGTTGCTTAGCATAGCTTCCTTCGATATCCAGCGCGCCAAAAATGGACTCAAGGCCCTGACGGATATACTCATGTAGATGCATCGCAAATGCGGTAGGAGTGACACCCTGGCCTTTGCGAATAAACAATGGGTCAGGGAATATAACGCGCAGTTTCTGAATATACTGACTGATTGCCGAGGGGGTCAGATTAAGCACTATCTCTTCATTAACGATCCCTTTa                                                             <  2:24189/251‑1 (MQ=255)                                                                                                                                                                                  gTGGTCGCAATAGTTATCGTTCGTTGCTTATCATAGCTTCCTTCGATATCCAGCGCGCCAAAAATGGACTCAAGGCCCTGACTGATATACTCATGTAGATGCATCGCAAATGCGGTAGGAGTGACACCCTGGCCTTTGCGAATAAACAATGGGTCAGGGAATATAACGCGCAGTTTCTGAATAGACTGACTGATTGCCGAGGGGGTCAGATTAAGCACTTTCGCTGCATTAACGATCCCTTTATGTACat                                                      >  1:113141/1‑250 (MQ=255)                                                                                                                                                                                        gCAATAGTTATCGTTCGTTGCTTATCATAGCTTCCTGCGATATCCAGCGCGCCAAAAATGGACTCAAGGCCCTGACTGATATACTCATGTAGATGCATCGCAAATGCGGTAGGAGTGACACCCTGGCCTTTGCGAATAAACAATGGGTCAGGGAATATAACGCGCAGTTTCTGAATAGACTGACTGATTGCCGAGGGGGTCAGATTAAGCACTTGCGCTGCATTAACGATCCCTTTATGTACATAGACAGc                                               <  2:435116/251‑1 (MQ=255)                                                                                                                                                                                              gTTATCGTTCGTTGCTTATCATAGCTTCCTTCGATATCCAGCGCGCCAAAAATGGACTCAAGGCCCTGACTGATATACTCATGTAGATGCATCGCAAATGCGGTAGGAGTGACACCCTGGCCTTTGCGAATAAACAATGGGTCAGGGAATATAACGCGCAGTTTCTGAATAGACTGACTGATTGCCGAGGGGGTCAGATTAAGCACTTTCGCTGCATTAACGATCCCTTTATGTACATATACAGCTTCaa                                          >  2:129997/1‑250 (MQ=255)                                                                                                                                                                                              gTTATCGTTCGTTGCTTAGCATAGCTTCCTTCGATATCCAGCGCGCCAAAAATGGACTCAAGGCCCTGACTTATATACTCATTTAGATGCATCGCAAATGCGGTAGGAGTGACACCCTGGCCTTTGCGAATAAACAATGGGTCAGGGAATATAACGCGCAGTTTCTGAATAGACTGACTGATTGCCGAGGGGGTCAGATTAAGCACTTTCGCTGCATGAACGATCCCTTTATGTACATATCCAGCTTCaaa                                         <  2:318457/251‑1 (MQ=255)                                                                                                                                                                                                       cgttGCTTATCATAGCTTCCTTCGATATCCAGCGCGCCAAAAATGGACTCAAGGCCCTGACTGATATACTCATGTAGATGCATCGCAAATGCGGTAGGAGTGACACCCTGGCCTTTGCGAATAAACAATGGGTCAGGGAATATAACGCGCAGTTTCTGAATAGACTGACTGATTGCCGAGGGGGTCAGATTAAGCACTTTCGCTGCATTAACGATCCCTTTATGTACATATACAGCTTCAAAAATAGTCAg                                <  1:217209/251‑1 (MQ=255)                                                                                                                                                                                                            cgtATCATAGCTTCCTTCGAGAGCCAGCGCGGCAAAAATGCACTCCAGGCCCTTGCGGATATACTCATGTAGATGCATCGCAAATGCGGTAGTAGTGACACCCAGTCCTTTGCGAATAAACAATGGGTCAGGGCATATATCGCGTAGTTTCTGAATAGACGGACTGATTGCCGAGGGGGTCAGGTTCAGCACTTTCGCGGCATTAACGATCCCGGTATGGATATATCCAGCGTCAAAGATAGTCAGCAGGt                           <  2:343165/249‑1 (MQ=255)                                                                                                                                                                                                               aTCATAGCTTCCTTCGATATCCAGCGCGCCAAAAATGGACTCAAGGCCCTGACTGATATACTCATGTAGATGCATCGCAAATGCGGTAGGAGTGACACCCTGGCCTTTGCGAATAAACAATGGGTCAGGGAATATAACGCGCAGTTTCTGAATAGACTGACTGATTGCCGAGGGGGTCAGATTAAGCACTTTCGCTGCATTAACGATCCCTTTATGTACATATACAGCTTCAAAAATAGTCAGAAGGTTa                         >  2:380485/1‑250 (MQ=255)                                                                                                                                                                                                                    aGCTTCCTTCGATATCCAGCGCGCCAAAAATGGACTCAAGGCCCTGACTGATATACTCATGTAGATGCATCGCAAATGCGGTAGGAGTGACACCCTGGCCTTTGCGAATAAACAATGGGTCAGGGAATATAACGCGCAGTTTCTGAATAGACTGACTGATTGCCGAGGGGGTCAGATTAAGCACTTTCGCTGCATTAACGATCCCTTTATGTACATATACAGCTTCAAAAATAGTCAGAAGGTTAAGATCa                   <  2:144736/251‑1 (MQ=255)                                                                                                                                                                                                                              gATATCCAGCGCGCCAAAAATGGACTCAAGGCCCTGACTGATATACTCATGTAGATGCATCGCAAATGCGGTAGGAGTGACACCCTGGCCTTTGCGAATAAACAATGGGTCAGGGAATATAACGCGCAGTTTCTGAATAGACTGACTGATTGCCGAGGGGGTCAGATTAAGCACTTTCGCTGCATTAACGATCCCTTTATGTACATATACAGCTTCAAAAATAGTCAGAAGGTTAAGATCAATATTTCGt          >  2:160886/1‑250 (MQ=255)                                                                                                                                                                                                                                     aGCGCGCCAAAAATGGACTCAAGGCCCTGACTGATATACTCATGTAGATGCATCGCAAATGCGGTAGGAGTGACACCCTGGCCTTTGCGAATAAACAATGGGTCAGGGAATATAACGCGCAGTTTCTGAATAGACTGACTGATTGCCGAGGGGGTCAGATTAAGCACTTTCGCTGCATTAACGATCCCTTTATGTACATATACAGCTTCAAAAATAGTCAGAAGGTTAAGATCAATATTTCGTAAGGTTGt  <  1:344503/251‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 ACATGATGTTGCACCGTACGATTGGTGCAAAACATGTTATCGATGATGAGATCGGTTTGAAACTGACTGAGTTGGTTTTCCGCGTCGCTGATGGGCGGGTTGCGCAGTAAAAGCTGCGGATAGTGAGTTTTAATCGCCCGATAGATGACAGGAAGGACCAGTGCTCCGACCGAGGGAGTGGTCGCAATAGTTATCGTTCGTTGCTTATCATAGCTTCCTTCGATATCCAGCGCGCCAAGAATGGACTCAAGGCCCTGACTGATATACTCATGTAGATGCATCGCAAATGCGGTAGGAGTGACACCCTGGCCTTTGCGAATAAACAATGGGTCAGGGAATATAACGCGCAGTTTCTGAATAGACTGACTGATTGCCGAGGGGGTCAGATTAAGCACTTTCGCTGCATTAACGATCCCTTTATGTACATATACAGCTTCAAAAATAGTCAGAAGGTTAAGATCAATATTTCGTAAGGTTGT  >  NC_000913/660831‑661309

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