Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,214,511	G→C	8.9%	P322R (CCT→CGT)	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,511	0	G	C	8.9%	137.9 / 5.1	45	P322R (CCT→CGT)	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 (23/18);  new base C (2/2);  total (25/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 = 9.62e-01

TTTCACGATGATTTACTGAAATCATGTGAAAGAATGTGCTGAAAATTATTTTTTCTCTGGCCACGCTATGGAAGGGATACCATTCAATTTAGCTTTAGCAAACAGATCTCCCTGAAACATCTCAATTCCTGCGGATTCAAGCCACATCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCC  >  NC_000913/1214219‑1214810                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                            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)                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                            TTTCACGATGATTTACTGAAATCATGTGAAAGAATGTGCTGAAAATTATTTTTTCTCTGGCCACGCTATGGAAGGGATACCATTCAATTTAGCTTTAGCAAACAGATCTCCCTGAAACATCTCAATTCCTGCGGATTCAAGCCACATCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCC  >  NC_000913/1214219‑1214810

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: ‑