Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,636,493	T→C	L163L (TTG→CTG)	creB →	DNA‑binding transcriptional regulator CreB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,636,493	0	T	C	100.0%	136.2 / NA	45	L163L (TTG→CTG)	creB	DNA‑binding transcriptional regulator CreB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (28/17);  total (28/17)

CTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAGCGCTGACTCGGTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTA  >  NC_000913/4636247‑4636734                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  cTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGGATGAGTTTTTATTGCTGa                                                                                                                                                                                                                                                >  2:435350/1‑250 (MQ=255)                    agGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCCCCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTcgcca                                                                                                                                                                                                                            <  1:372913/251‑5 (MQ=255)                        cGATCGCCTGCTTGGGCTGGAAATGGTTGCTGACGACTATGTGGCTAAACGTTTTTCACCCCGGGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGcc                                                                                                                                                                                                                        <  2:443227/251‑1 (MQ=255)                              ccggcttggggtgagaatgggTGCTGACGGCTGTGTGGCTAAACCGTTTTCCCCCCGCGAAGGGTGCGCCGGGGTGCGCACCTTACGGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCCCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCt                                                                                                                                                                                                                  <  2:78662/232‑1 (MQ=255)                                           caaaTTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGc                                                                                                                                                                                                      >  1:263835/2‑250 (MQ=255)                                                gggTGCGGACGACGATGGGGCAAAACCTTTTTCCCCCCGCGAAGGGTGCGCCAGGGTGCGCACTTTCCTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCCCCCGTCACCCGTATTGGTCATTTTGAATGGAAGGACCCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACtaa                                                                                                                                                                                                <  1:90796/250‑3 (MQ=255)                                                    gCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATgg                                                                                                                                                                                             >  1:262699/1‑250 (MQ=255)                                                          gcctatgTGGCTAAACCGTTTCCCCCCCGCGAAGGGGGCGCCGGGGTGCGCACCTTACTGCGGCGGGTGAAGAAGTTCCCGTCGCCGTCCCCCGTCACCCGTATTGGTCATTTTGAATTGAATGAACCCGGGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAg                                                                                                                                                                                                       <  2:8859/228‑1 (MQ=255)                                                          gacTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAg                                                                                                                                                                                                       >  1:8859/1‑234 (MQ=255)                                                             aatgtgacTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCGTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCCGCTGGTTTGACACGCCATTAACGCTGACGCGGAATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCTAGCCGCGTCTGGTCCCGCCAGCACCTAATGGCTAGCGtgtg                                                                                                                                                                                   >  2:352122/8‑248 (MQ=255)                                                                                   cccGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGCAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGCGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGACGCGTCTGGTCCCGCCAGCCAATAATGGTTAGCGGGTGGGTAGAGGCGCAGGACCCCTAc                                                                                                                                                             >  2:259117/1‑251 (MQ=255)                                                                                   cccGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGaa                                                                                                                                                                                                                                               >  2:337774/1‑169 (MQ=255)                                                                                   cccGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCCCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGaa                                                                                                                                                                                                                                               <  1:337774/169‑1 (MQ=255)                                                                                             gtgCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCg                                                                                                                                                   >  1:272823/1‑251 (MQ=255)                                                                                              ggcgcCGGGGTGCGCCCCTTAGTGGGGCGGGTGAAGAAGTTCTGGCCGCGGGCCCCCGTCATCGGTTTGGGTCATTTTGAATTGAATGACCCCGCGGCGAAGATAGGCGGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGt                                                                                                                                                  <  1:482565/250‑1 (MQ=255)                                                                                                                  aCTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGa                                                                                                                                                                             >  2:337048/1‑204 (MQ=255)                                                                                                                                  aagTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCg                                                                                                                                                        >  2:482331/1‑209 (MQ=255)                                                                                                                                  aagTTCTCGCCGCCGTCCCCCGTCACCCGTATTGGTCATTTTGAATTGAAGGACCCCGCGGCGCAGATCAGCTGGTTGGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCg                                                                                                                                                        <  1:482331/209‑1 (MQ=255)                                                                                                                                    gTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCAAAGCTgcgcg                                                                                                            >  1:79054/1‑251 (MQ=255)                                                                                                                                      tctcGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGCTCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTATGGTCCTGCCAGCAACTAATGGATAGCGTGTGGGAGGATGCGCAGGCCACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGcc                                                                                                          >  2:75592/1‑251 (MQ=255)                                                                                                                                                gTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTACTGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCACGCTGCGCGCCATCACccccc                                                                                                >  2:180867/1‑250 (MQ=255)                                                                                                                                                gTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGCTACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCg                                                                                                >  2:308400/1‑251 (MQ=255)                                                                                                                                                   tcCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGAGCAGATCAGCTGGTTTGACACGCAATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGAACCGTTGATACCAACATTAAAAAGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATc                                                                                             >  1:155630/1‑251 (MQ=255)                                                                                                                                                       gTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTc                                                                                         >  2:116036/1‑251 (MQ=255)                                                                                                                                                         cATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGCTCTTTCAc                                                                                       >  2:366035/1‑251 (MQ=255)                                                                                                                                                                        ttttGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATc                                                                        >  1:70199/1‑251 (MQ=255)                                                                                                                                                                                       ccccggggtgaagatctgttGGTTTGACACGCCATAAACGCTGACGCGGTATGATTTTTTATTGCTGAAGACGTTACCCAAGTCGCCAGGCCGCGTCTGGGCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGAtata                                                         <  2:70199/232‑1 (MQ=255)                                                                                                                                                                                            cggcgCAGGACAGCTGGTTTGACACGCCATTAAAGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGACCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGAGCACTCAACCCCGCTCTTTAACCGATTAATACTCCTCGCGGCATGGGAGATAGccat                                                    >  2:474038/1‑249 (MQ=255)                                                                                                                                                                                            cggcgCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGtt                                                                                                                                                 <  1:411866/158‑2 (MQ=255)                                                                                                                                                                                            cggcgCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGtt                                                                                                                                                 >  2:411866/1‑157 (MQ=255)                                                                                                                                                                                                  ataTCAGTTGGTTTGACACGCCATTAACGTTGACGCGGTATGAGTTTTTAGTGCTGAAGACGTTACTCAAGTCGCCAGGCTGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGc                                              <  1:363699/249‑1 (MQ=255)                                                                                                                                                                                                      cggttgttttGACCCGCCATTACCGTTGACGGGGTAGGATTTTTTATTGCTGAAGCCTTTACCCAAGCCGCCAGGCCGCGTTTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGAGCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGt                                          <  2:452968/244‑1 (MQ=255)                                                                                                                                                                                                      cAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCCTCTCGGCCTGGGATATAGCCTGAGGGGCCTGt                                          >  2:228925/1‑251 (MQ=255)                                                                                                                                                                                                              ttGACACGCCATTAACGTTGTCGCGGTATGAGTTTTAATTGCGGAAGACGTTCCTCAGGTCGCCAGTCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGt                                   <  1:36137/250‑1 (MQ=255)                                                                                                                                                                                                                gcaacGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCAGAAGCTGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATc                                <  1:176256/248‑1 (MQ=255)                                                                                                                                                                                                                aacactccattaacgcggatgcGGTATGAGTTTTTTGTGCGGAAGCCGTTACTCAAGGCGCCAGGCCGCGTCTGTTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACTCCTACGATCGCACTGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTAt                                 <  1:301728/229‑1 (MQ=255)                                                                                                                                                                                                                                tGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCACGATCTTTCACCGATTAATACTCATAGCGGCATGGGATATAGCCTGCGGGCACGGTAATGCGTATCGGACTGCGGTtgctgc                >  2:267636/1‑251 (MQ=255)                                                                                                                                                                                                                                  acgcGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAg                                                                                                                   >  1:221834/4‑150 (MQ=255)                                                                                                                                                                                                                                  acgcGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCGGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAg                                                                                                                   <  2:221834/147‑1 (MQ=255)                                                                                                                                                                                                                                         aTGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGCTATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTAtt       >  2:480911/1‑251 (MQ=255)                                                                                                                                                                                                                                           gAGTTTTTATTGCTGAAGACGGTACCCAAGTCGCCAGTCGGGGTCGTGCCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATTCCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTAtttt     <  1:30569/251‑1 (MQ=255)                                                                                                                                                                                                                                             gTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGCAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGAATCGGCATGCGGTTGCTGCTGGGCTAtttttt   >  1:452621/1‑251 (MQ=255)                                                                                                                                                                                                                                             gTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCgg                              >  1:40646/1‑224 (MQ=255)                                                                                                                                                                                                                                             gTTTTTATTGCTGAAGACGTTACTAAAGTAGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTAtttttt   >  1:394760/1‑251 (MQ=255)                                                                                                                                                                                                                                              tttttATTGCTGAAGGCGTTCCTCAAGGAGCCCTGCCGCGACTGGTCCAGCCAGCAACTAATGGCTAGCGAGTGGGAAGATGCGCAGGACACCTACGATCGCAACGTTGATAACCACATTAAAACGGTACGTGCCAAGCTGAGCGGCATCAACCCCCATCTTTCACCGATTAATACTCATCGCGGCATGGGATAAAGCCTGGAGGGACTGAAATTCGGACCGGAATGCGGTGGCTGCATGGCTAtttttgg  >  2:325083/1‑249 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  CTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAGCGCTGACTCGGTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTA  >  NC_000913/4636247‑4636734

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