Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,214,780	A→C	L232L (CTT→CTG)	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,780	0	A	C	100.0%	235.4 / NA	60	L232L (CTT→CTG)	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 A (0/0);  new base C (30/30);  total (30/30)

GACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGAT  >  NC_000913/1214633‑1214927                                                                                                                                                    |                                                                                                                                                    gACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGt                                                                                                                                                  <  1:462933/151‑1 (MQ=255)                aaaTTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTg                                                                                                                                   <  2:351855/151‑1 (MQ=255)                  aTTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGaa                                                                                                                                 <  2:5391/151‑1 (MQ=255)                  aTTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGa                                                                                                                                  >  2:748812/1‑150 (MQ=255)                       aCGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCaa                                                                                                                            <  2:217303/151‑1 (MQ=255)                         gatgatTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAgc                                                                                                                          <  2:316694/151‑1 (MQ=255)                         gatgatTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAg                                                                                                                           >  2:504733/1‑150 (MQ=255)                               ttGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCtt                                                                                                                     <  1:31670/150‑1 (MQ=255)                                 gTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACt                                                                                                                  >  1:138327/1‑151 (MQ=255)                                    cAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACtttt                                                                                                               >  1:386797/1‑151 (MQ=255)                                    cAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACtttt                                                                                                               >  1:675466/1‑151 (MQ=255)                                    cAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACttt                                                                                                                <  2:184373/150‑1 (MQ=255)                                       gCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTg                                                                                                              >  1:351734/1‑149 (MQ=255)                                                cATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCtgt                                                                                                   <  1:238206/151‑1 (MQ=255)                                                  ttGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCtgtg                                                                                                  >  2:257422/1‑150 (MQ=255)                                                   tGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTgtgt                                                                                                 >  1:804881/1‑150 (MQ=255)                                                        tttATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGAttt                                                                                           >  2:450354/1‑151 (MQ=255)                                                        tttATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGAttt                                                                                           >  1:46472/1‑151 (MQ=255)                                                           aTTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGTAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTccc                                                                                        <  2:770155/151‑1 (MQ=255)                                                             ttCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTcccc                                                                                       <  1:50155/150‑1 (MQ=255)                                                               cATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGt                                                                                     <  1:423790/150‑1 (MQ=255)                                                                 ttAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCtt                                                                                  >  2:679872/1‑151 (MQ=255)                                                                  tAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCttt                                                                                 >  2:100913/1‑151 (MQ=255)                                                                   aaGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCttt                                                                                 <  1:632903/150‑1 (MQ=255)                                                                   aaGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCttt                                                                                 >  2:628161/1‑150 (MQ=255)                                                                    aGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTAc                                                                               >  2:79902/1‑151 (MQ=255)                                                                      tAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGc                                                                             >  2:425743/1‑151 (MQ=255)                                                                       aaaaaaGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGc                                                                             <  1:401153/150‑1 (MQ=255)                                                                          aaaGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGcc                                                                         >  2:172738/1‑151 (MQ=255)                                                                            agagACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGccc                                                                        >  1:314754/1‑150 (MQ=255)                                                                              agaCTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTCCTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAAc                                                                     <  2:370325/151‑1 (MQ=255)                                                                               gaCTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAAcc                                                                    <  1:831069/151‑1 (MQ=255)                                                                                aCTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCg                                                                   <  2:424800/151‑1 (MQ=255)                                                                                        aGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCt                                                           <  1:524300/151‑1 (MQ=255)                                                                                            tCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATg                                                       >  2:565114/1‑151 (MQ=255)                                                                                              gTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATggg                                                     <  2:763852/151‑1 (MQ=255)                                                                                                tAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGc                                                    <  2:875904/150‑1 (MQ=255)                                                                                                 aaCCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCt                                                   <  2:699715/150‑1 (MQ=255)                                                                                                   ccAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATc                                                <  2:45461/151‑1 (MQ=255)                                                                                                      gggTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTc                                             >  2:159208/1‑151 (MQ=255)                                                                                                       ggTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCa                                            <  1:895879/151‑1 (MQ=255)                                                                                                       ggTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCa                                            <  1:452030/151‑1 (MQ=255)                                                                                                          cATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCAtt                                          <  2:88438/150‑1 (MQ=255)                                                                                                             aGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCAttttt                                       <  1:396128/150‑1 (MQ=255)                                                                                                                tAATAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGc                                   >  1:798659/1‑151 (MQ=255)                                                                                                                 aaTAGATTGATTGAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGc                                   <  2:303268/150‑1 (MQ=255)                                                                                                                      attgattgAAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATg                             <  1:212925/151‑1 (MQ=255)                                                                                                                        tgattgaAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATgg                            <  2:424980/150‑1 (MQ=255)                                                                                                                         gattgaAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCt                          <  1:809062/151‑1 (MQ=255)                                                                                                                          attgaAATCATTTTATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCtt                         >  2:780874/1‑151 (MQ=255)                                                                                                                                      ttATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCA‑TTTTTTGCACCATGGCTTCAAAAGCGGtt              >  1:679648/1‑149 (MQ=255)                                                                                                                                       tATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTa             <  1:689540/150‑1 (MQ=255)                                                                                                                                       tATCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTAt            >  2:209222/1‑151 (MQ=255)                                                                                                                                         tCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATcc          >  1:838142/1‑151 (MQ=255)                                                                                                                                         tCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATcc          >  1:757934/1‑151 (MQ=255)                                                                                                                                         tCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATcc          >  1:73415/1‑151 (MQ=255)                                                                                                                                         tCACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATcc          >  2:257893/1‑151 (MQ=255)                                                                                                                                          cACCGAGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATccg         >  2:492081/1‑151 (MQ=255)                                                                                                                                               aGCTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCcgcg    >  1:417686/1‑151 (MQ=255)                                                                                                                                                 cTCCAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGAt  >  1:494020/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GACTTCACTTTCAGTAAATTCAACGATGATTTGTTCAGGCACCAGAGCATTGGCCTTTATTTCATTAAGTAAAAAAGAGACTGCGTCAGGTTCGTTAACCAGGGTCATAGGTAATAGATTGATTGAAATCATTTTATCACCGAGCTCAAGTGCGTGTGCCATCGTGAATGCAAGCGCCTTACTTTTGAGATCCGCTGTGTAGATTTCCCCGTCTTTACGCTGCCCAACCGCTATGGCTGATGGGCTATCTTCATTTTTTTGCACAATGGCTTCAAAAGCGATTATCCGCCGCGAT  >  NC_000913/1214633‑1214927

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