Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,214,515	G→T	8.7%	Q321K (CAG→AAG)	bluF ←	anti‑repressor for YcgE, blue light‑responsive; FAD‑binding; inactive c‑di‑GMP phosphodiesterase‑like EAL domain protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,214,515	0	G	T	8.7%	142.5 / 4.5	46	Q321K (CAG→AAG)	bluF	anti‑repressor for YcgE, blue light‑responsive; FAD‑binding; inactive c‑di‑GMP phosphodiesterase‑like EAL domain protein
Reads supporting (aligned to +/- strand):  ref base G (24/18);  new base T (2/2);  total (26/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.46e-01

TTTCACGATGATTTACTGAAATCATGTGAAAGAATGTGCTGAAAATTATTTTTTCTCTGGCCACGCTATGGAAGGGATACCATTCAATTTAGCTTTAGCAAACAGATCTCCCTGAAACATCTCAATTCCTGCGGATTCAAGCCACATCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACT  >  NC_000913/1214219‑1214815                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                             tttCACGATGATTTACTGAAATCATGTGAAAGAATGTGCTGAAAATTATTTTTTCTCTGGCCACGCTATGGAAGGGATACCATTCAATTTAGCTTTAGCAAACAGATCTCCCTGAAACATCTAAATTCCTGCGGATTCAAGCCACATCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCAATTTTATGAAAAATGGAAACCAATAACTAACTGATTTTAATTCTGTCAGGCTGGAAg                                                                                                                                                                                                                                                                                                          >  1:147969/1‑301 (MQ=255)    taCGATGATTTACTGAAATCATGTGAAAGAATGTGCTGAAAATTATTTTTTCTCTGGCCACGCTATGGAAGGGATACCATTCAATTTAGCTTTAGCAAACAGATCTCCCTGAAACATCTCAATTCCTGCGGATTCAAGCCACATCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGt                                                                                                                                                                                                                                                                                                       <  2:61000/300‑1 (MQ=255)    cACGATGATTTACTGAAATCATGTGAAAGAATGTGCTGAAAATTATTTTTTCTCTGGCCACGCTATGGAAGGGATACCATTCAATTTAGCTTTAGCAAACAGATCTCCCTGAAACATCTCAATTCCTGCGGATTCAAGCCACATCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGt                                                                                                                                                                                                                                                                                                       >  1:120195/1‑301 (MQ=255)     aCGATGATTTACTGAAATCATGTGAATGAATGTTCTGAAAATTATTTTTTCTCGTGTCACGCTATGGAAGGGATACCATTCAATTTAGCTTTAGCAAACAGATCTCCCTGAAACATCTCAATTCCTGCGGATTCAAGCCACATCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTg                                                                                                                                                                                                                                                                                                      <  1:193527/301‑1 (MQ=255)         tgatTTTCTGAAATCATGTGAAAGAATGTGCTGAAAATTATTTTTTCTCTGGCCACGCTATGGAAGGGATACCATTCAATTTAGCTTTAGCAAACAGATCTCCCTGAAACATCTCAATTCCTGCGGATTCAAGCCACATCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAg                                                                                                                                                                                                                                                                                                  <  2:10923/301‑1 (MQ=255)          ccgacgctcttccgatctgAAAGAATGTGCTGAAAATTATTTTTTTTCTGGCCACGCTATGGAAGGGATACCATTCAATTTAGCTTTAGCAAACAGATCTCCCTGAAACATCTCAATTCCTGCGGATTCAAGCCACATCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCACGCTTCCAGCCTGACAga                                                                                                                                                                                                                                                                                                 <  1:406912/284‑2 (MQ=255)             ttACTGAAATCATGTTAAAGAATGTGCTGAAAATTATTTTTTCTCTGGCCACGCTATGGAAGGGATACCATTCAATTTAGCTTTAGCAAACAGATCTCCCTGAAACATCTCAATTCCTGCGGATTCAAGCCACATCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGt                                                                                                                                                                                                                                                                                              <  2:120319/301‑1 (MQ=255)              tACTGAAATCATGTGAAAGAATGTGCTGAAAATTATTTTTTCTCTGGCCACGCTATGGAAGGGATACCATTCAATTTAGCTTTAGCAAACAGATCTCCCTGAAACATCTCAATTCCTGCGGATTCAAGCCACATCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTa                                                                                                                                                                                                                                                                                             >  1:372908/1‑301 (MQ=255)                            gAAAGAATGTGCTGAAAATTATTTTTTCTCTGGCCACGCTATGGAAGGGATACCATTCAATTTAGCTTTAGCAAACAGATCTCCCTGAAACATCTCAATTCCTGCGGATTCAAGCCACATCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCACGCTTCCAGCCTGACAGAAgatcggaagagcacacg                                                                                                                                                                                                                                                                               >  2:407475/1‑285 (MQ=255)                                        tGAAAATTATTTTTTCTCTGGCCACGCTATGGAAGGGATACCATTCAATTTAGCTTTAGCAAACAGATCTCCCTGAAACATCTCAATTCCTGCGGATTCAAGCCACATCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGAAAGGAGTAAAAAACCAGAAAAACCAACGCCaaaa                                                                                                                                                                                                                                                                   >  1:193645/1‑301 (MQ=255)                                                 ttttttCTCTGGCCACGCTATGGAAGGGATACCATTCAATTTAGCTTTAGCAAACAGATCTCCCTGAAACATCTCAATTCCTGCGGATTCAAGCCACATCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGAACAAtt                                                                                                                                                                                                                                                          >  2:146657/1‑301 (MQ=255)                                                               acttgaTGGAAGTGGTATCATTCAATTTATCTTTAGCAAACAGATCTCCCTGAAGCTTCTCAATGCCTGCGGATTCAAGCCACATCCACTCTTCTGGNNTTGCCACGCCCATAGGNCTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGNCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATAcc                                                                                                                                                                                                                                            <  1:98156/296‑1 (MQ=255)                                                                                           agctttagcAAACAGATCTCTCTGAAACATCTCAATTCCTGCGGATTCAAGCCACATCCACTCGTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCg                                                                                                                                                                                                                <  1:29536/301‑1 (MQ=255)                                                                                            tctttagcAAACAGATCTCCCTGAAACATCTCAATTCCTGCGGATTCAAGCCACATCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTTATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTAGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTAAATCGCtaaag                                                                                                                                                                                                               <  2:395697/300‑5 (MQ=255)                                                                                            tctttagcAAACAGATCTCCCTGAAACATCTCAATTCCTGCGGATTCAAGCCACATCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCgg                                                                                                                                                                                                               <  2:193888/300‑1 (MQ=255)                                                                                                              cccTGAAACATCTCAATTCCTGCGGATTCAAGCCACATCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACATGAATTTCAAGTGATGTACAGCATCTTATGATCGCCTGAATAACTGCCTGCCGTGGCACACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTTACAGGAGTAACAAACCAGCAAAAACTGCGCAAAAATGATCAATTGCTACAATGATACCAACAACCTTTAGCCATTTAATAGCTTAAGCaaaacaataaaaacgag                                                                                                                                                                                             >  1:227122/1‑287 (MQ=255)                                                                                                                           ttcttccgttttaTTCAAGCCACATCCACTCTTCTGGTGTTGCCACGCCCCTAGCACTGACTTTAATTTCAAGTGATGTACAGCATTTTATGTTCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGTTAATCAATTCCTGACTGATTTTAAGTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACt                                                                                                                                                                                <  1:231517/289‑1 (MQ=255)                                                                                                                               tCCTGCGGATTCAAGCCACTTCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCa                                                                                                                                                                            <  2:148095/301‑1 (MQ=255)                                                                                                                                       aTTCAAGCCACATCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACtagatcggaagaa                                                                                                                                                                    >  2:231760/1‑289 (MQ=255)                                                                                                                                        ttCAAGCCACATCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTACAGTAAATTCa                                                                                                                                                                   >  2:417702/1‑301 (MQ=255)                                                                                                                                                                    gCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGAATTGTTACGGCACCAGACCAt                                                                                                                                       >  1:453948/1‑301 (MQ=255)                                                                                                                                                                       ccGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACATCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTgg                                                                                                                                    <  2:269476/300‑1 (MQ=255)                                                                                                                                                                                  aCTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACCCATCAAACCGACAACTGACTTCACTTTCAGTAAATTCAACGATGCTTTGTTCAGGCCCCATAACTTTGGCCTtttttttt                                                                                                                         >  1:324745/1‑299 (MQ=255)                                                                                                                                                                                       cTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTaaa                                                                                                                    >  2:376728/1‑300 (MQ=255)                                                                                                                                                                                                        gTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTCAGTAAAAAAAAGAACGCGt                                                                                                   >  1:559648/1‑301 (MQ=255)                                                                                                                                                                                                                   tATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTa                                                                                        >  1:52523/1‑301 (MQ=255)                                                                                                                                                                                                                    aTGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATAATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAAAGACTGCGTCAGGTGAGTTaa                                                                                       >  1:589892/1‑301 (MQ=255)                                                                                                                                                                                                                      tccgatcTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGGGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAAcc                                                                                     <  2:448107/295‑1 (MQ=255)                                                                                                                                                                                                                            cTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTAACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTAAGGTACGTTAACCCAGGTc                                                                               >  1:248704/1‑301 (MQ=255)                                                                                                                                                                                                                              gAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGGGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAgatcgga                                                                             >  1:447543/1‑295 (MQ=255)                                                                                                                                                                                                                                 taatTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAgg                                                                          >  1:136683/1‑301 (MQ=255)                                                                                                                                                                                                                                 gacgctcttcctatctcccACTTTTATTAACATTTGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAgg                                                                          <  1:593191/285‑1 (MQ=255)                                                                                                                                                                                                                                   atTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTCTTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGATCATAGGTa                                                                        >  1:609808/1‑301 (MQ=255)                                                                                                                                                                                                                                           ccGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAgattg                                                                <  2:373404/301‑1 (MQ=255)                                                                                                                                                                                                                                              tGGCCTACTTTTATTATCATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGTAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAgattgatt                                                             <  2:227365/301‑1 (MQ=255)                                                                                                                                                                                                                                                 cccACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAggagatcggaagagcaca                                                          >  2:593952/1‑285 (MQ=255)                                                                                                                                                                                                                                                  ccACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGaaa                                                         <  2:136808/301‑1 (MQ=255)                                                                                                                                                                                                                                                  ccACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTAATTGaaa                                                         >  2:496168/1‑301 (MQ=255)                                                                                                                                                                                                                                                                   ggTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCATAACCTGCGCAAACATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGCTATGACTTCTTTTTCAGTAAATTCAACGATGTTTTGTTCAGGCACCAGTTCATTGACCTTTATTACATTAAGTAAAAACTATACTACCTCAGTTACGTTAACGAAGGTAATAGATAATAAAATAATTGAAATCAtttaaaaaccaacc                                        >  1:582066/1‑290 (MQ=255)                                                                                                                                                                                                                                                                          aaTTCCTGACTTATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTg                                 <  1:417139/301‑1 (MQ=255)                                                                                                                                                                                                                                                                               ctgactgaTTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATAGAAATCATTTTATCACCGAGCCCAAGTgaatat                            >  1:532101/1‑296 (MQ=255)                                                                                                                                                                                                                                                                                tgactgaTTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGtgtg                           <  2:52591/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                tctggtttgttactccTGTCACGCTTCCAGCCTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTTACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGtgtg                           <  2:74492/286‑1 (MQ=255)                                                                                                                                                                                                                                                                                 ctggtttgttactccTGTCACGCTTCCAGCCTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTTACCAGGGTCATAGGTAATCAATTGATTGAAATCATTTAATCACCGAGCTCAAGTGCGtgtga                          >  1:74402/15‑300 (MQ=255)                                                                                                                                                                                                                                                                                                    aGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGcc       >  2:27545/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                                        tGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGTTTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTAc   >  2:518885/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                                         ggAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACt  >  2:190253/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                                         ggAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTACAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACt  <  1:146531/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                             TTTCACGATGATTTACTGAAATCATGTGAAAGAATGTGCTGAAAATTATTTTTTCTCTGGCCACGCTATGGAAGGGATACCATTCAATTTAGCTTTAGCAAACAGATCTCCCTGAAACATCTCAATTCCTGCGGATTCAAGCCACATCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACT  >  NC_000913/1214219‑1214815

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑