Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,214,521	G→C	8.3%	R319G (CGC→GGC)	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,521	0	G	C	8.3%	153.3 / 5.2	48	R319G (CGC→GGC)	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/20);  new base C (2/2);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

CACGATGATTTACTGAAATCATGTGAAAGAATGTGCTGAAAATTATTTTTTCTCTGGCCACGCTATGGAAGGGATACCATTCAATTTAGCTTTAGCAAACAGATCTCCCTGAAACATCTCAATTCCTGCGGATTCAAGCCACATCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTG  >  NC_000913/1214222‑1214819                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                           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)                                                                                                                                                                                                                                                                                                          gCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAATCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTg  <  2:273315/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                           CACGATGATTTACTGAAATCATGTGAAAGAATGTGCTGAAAATTATTTTTTCTCTGGCCACGCTATGGAAGGGATACCATTCAATTTAGCTTTAGCAAACAGATCTCCCTGAAACATCTCAATTCCTGCGGATTCAAGCCACATCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTG  >  NC_000913/1214222‑1214819

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