Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,015,840	A→G	G110G (GGT→GGC)	ysgA ←	putative dienelactone hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,015,840	0	A	G	100.0%	268.8 / NA	71	G110G (GGT→GGC)	ysgA	putative dienelactone hydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (40/31);  total (40/31)

TACAGCCAGGTGATACGTCCACCCCAGCAGAATCCGGTGATCATTAAACGATGAACATCGCCGCCGTTGCGGGAAGCCCAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGACCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACG  >  NC_000913/4015703‑4015969                                                                                                                                          |                                                                                                                                  tACAGCCAGGTGATACGTCCACCCCAGCAGAATCCGGTGATCATTAAACGATGAACATCGCCGCCGTTGCGGGAAGCCCAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGcc                                                                                                                                 <  1:608076/140‑1 (MQ=255)       cAGGTGATACGTCCACCCCAGCAGAATCCGGTGTTCATTAAACGATGAACATCGCCGCCGTGGCGGGAGGCCCAACTGGCGATATGATCGAGATCGGCGAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAg                                                                                                                           <  2:195478/140‑1 (MQ=255)       cAGGTGATACGTCCACCCCAGCAGAATCCGGTGATCATTAAACGATGAACATCGCCGCCGTTGCGGGAAGCCCAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAg                                                                                                                           <  1:1118532/140‑1 (MQ=255)         ggTGATACGTCCACCCCAGCAGAATCCGGTGATCATTAAACGATGAACATCGCCGCCGTTGCGGGAAGCCCAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCa                                                                                                                         >  2:347541/1‑140 (MQ=255)         ggTGATACGTCCACCCCAGCAGAATCCGGTGATCATTAAACGATGAACATCGCCGCCGTTGCGGGAAGCCCAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCa                                                                                                                         >  1:32822/1‑140 (MQ=255)           tGATACGTCCACCCCAGCAGAATCCGGTGATCATTAAACGATGAACATCGCCGCCGTTGCGGGAAGCCCAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAAc                                                                                                                       >  1:573923/1‑140 (MQ=255)              tACGTCCACCCCAGCAGAATCCGGTGATCATTAAACGATGAACATCGCCGCCGTTGCGGGAAGCCCAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTg                                                                                                                    >  2:834536/1‑140 (MQ=255)               aCGTCCACCCCAGCAGAATCCGGTGATCATTAAACGATGAACATCGCCGCCGTTGCGGGAAGCCCAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTgg                                                                                                                   >  1:79455/1‑140 (MQ=255)                cGTCCACCCCAGCAGAATCCGGTGATCATTAAACGATGAACATCGCCGCCGTTGCGGGAAGCCCAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAg                                                                                                                           >  1:809414/1‑131 (MQ=255)                cGTCCACCCCAGCAGAATCCGGTGATCATTAAACGATGAACATCGCCGCCGTTGCGGGAAGCCCAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAg                                                                                                                           <  2:809414/131‑1 (MQ=255)                     cccccAGCAGAATCCGGTGATCATTAAACGATGAACATCGCCGCCGTTGCGGGAAGCCCAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGatat                                                                                                             <  2:79455/139‑1 (MQ=255)                     aCCCCAGCAGAATCCGGTGATCATTAAACGATGAACATCGCCGCCGTTGCGGGAAGCCCAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTTGGGatat                                                                                                             >  2:908987/1‑140 (MQ=255)                     aCCCCAGCAGAATCCGGTGATCATTAAACGATGAACATCGCCGCCGTTGCGGGAAGCCCAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGatat                                                                                                             >  1:546093/1‑140 (MQ=255)                      ccccAGCAGAATCCGGTGATCATTAAACGATGAACATCGCCGCCGTTGCGGGAAGCCCAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATc                                                                                                            >  1:1081943/1‑140 (MQ=255)                            cagAATCCGGTGATCATTAAACGATGAACATCGCCGCCGTTGCGGGAAGCCCAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCaaa                                                                                                      >  2:919707/1‑140 (MQ=255)                                aTCCGGTGATCATTAAACGATGAACATCGCCGCCGTTGCGGGAAGCCCAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCa                                                                                                  <  2:573923/140‑1 (MQ=255)                                 tCCGGTGATCATTAAACGATGAACATCGCCGCCGTTGCGGGAAGCCCAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCAt                                                                                                 >  1:292587/1‑140 (MQ=255)                                  ccGGTGATCATTAAACGATGAACATCGCCGCCGTTGCGGGAAGCCCAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCAtt                                                                                                <  2:14236/140‑1 (MQ=255)                                         tCATTAAACGATGAACATCGCCGCCGTTGCGGGAAGCCCAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCg                                                                                         >  1:1055832/1‑140 (MQ=255)                                               aaCGATGAACATCGCCGCCGTTGCGGGAAGCCCAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTcg                                                                                   <  2:32822/140‑1 (MQ=255)                                               aaCGATGAACATCGCCGCCGTTGCGGGAAGCCCAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTcg                                                                                   <  2:1081943/140‑1 (MQ=255)                                               aaCGATGAACATCGCCGCCGTTGCGGGAAGCCCAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTcg                                                                                   <  1:205448/140‑1 (MQ=255)                                               aaCGATGAACATCGCCGCCGTTGCGGGAAGCCCAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTcg                                                                                   <  1:389866/140‑1 (MQ=255)                                                  gATGAACATCGCCGCCGTTGCGGGAAGCCCAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCgg                                                                                <  2:1055832/140‑1 (MQ=255)                                                    tGAACATCGCCGCCGTTGCGGGAAGCCCAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCgg                                                                                >  2:336674/1‑138 (MQ=255)                                                    tGAACATCGCCGCCGTTGCGGGAAGCCCAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCgg                                                                                <  1:336674/138‑1 (MQ=255)                                                    tGAACATCGCCGCCGTTGCGGGAAGCCCAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGaa                                                                              >  2:1143516/1‑140 (MQ=255)                                                    tGAACATCGCCGCCGTTGCGGGAAGCCCAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGaa                                                                              >  2:1133510/1‑140 (MQ=255)                                                           cgccgcCGTTGCGGGAAGCCCAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGt                                                                       >  1:642673/1‑140 (MQ=255)                                                               gccgTTGCGGGAAGCCCAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAg                                                                   <  1:1033751/140‑1 (MQ=255)                                                               gccgTTGCGGGAAGCCCAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAg                                                                   <  1:834536/140‑1 (MQ=255)                                                                     gCGGGAAGCCCAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGa                                                             >  1:1118412/1‑140 (MQ=255)                                                                      cGGGAAGCCCAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGAt                                                            <  2:365557/140‑1 (MQ=255)                                                                         gAAGCCCAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGc                                                         <  1:443163/140‑1 (MQ=255)                                                                           aGCCCAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCa                                                       <  2:292587/140‑1 (MQ=255)                                                                             cccAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTTGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGa                                                     <  1:908987/140‑1 (MQ=255)                                                                             cccAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGa                                                     >  2:545739/1‑140 (MQ=255)                                                                             cccAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGa                                                     <  2:546093/140‑1 (MQ=255)                                                                                     gCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGa                                                             <  1:755148/124‑1 (MQ=255)                                                                                     gCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGa                                                             >  2:755148/1‑124 (MQ=255)                                                                                     gCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTc                                             >  1:956972/1‑140 (MQ=255)                                                                                       gACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCa                                           >  2:1169969/1‑140 (MQ=255)                                                                                       gACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCa                                           >  2:700151/1‑140 (MQ=255)                                                                                       gACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCa                                           >  1:628672/1‑140 (MQ=255)                                                                                        aCATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAg                                          >  2:163993/1‑140 (MQ=255)                                                                                          aTGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAgcg                                        >  1:211907/1‑140 (MQ=255)                                                                                           tGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAgcgc                                       >  1:641700/1‑140 (MQ=255)                                                                                               cGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGa                                   <  2:211907/140‑1 (MQ=255)                                                                                                 agaTCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACg                                 <  1:347541/140‑1 (MQ=255)                                                                                                 agaTCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGCAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACg                                 <  2:1070971/140‑1 (MQ=255)                                                                                                   aTCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGc                               >  1:565780/1‑140 (MQ=255)                                                                                                      ggCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGAc                            >  2:645270/1‑140 (MQ=255)                                                                                                       gCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGa                                   >  1:10837/1‑132 (MQ=255)                                                                                                       gCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGa                                   <  2:10837/132‑1 (MQ=255)                                                                                                       gCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACa                           >  2:807866/1‑140 (MQ=255)                                                                                                       gCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACa                           >  1:36163/1‑140 (MQ=255)                                                                                                        ccAGCACCTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACaa                          >  1:775901/1‑140 (MQ=255)                                                                                                              ccTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTcg                    <  2:1118412/140‑1 (MQ=255)                                                                                                               cTGCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTcgc                   >  1:754917/1‑140 (MQ=255)                                                                                                                 gCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCgg                 <  2:956972/140‑1 (MQ=255)                                                                                                                 gCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCgg                 >  2:860177/1‑140 (MQ=255)                                                                                                                 gCGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCgg                 <  2:36163/140‑1 (MQ=255)                                                                                                                  cGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGa                >  2:731675/1‑140 (MQ=255)                                                                                                                  cGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGa                >  1:807250/1‑140 (MQ=255)                                                                                                                  cGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGa                >  1:650437/1‑140 (MQ=255)                                                                                                                  cGAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGa                >  1:49180/1‑140 (MQ=255)                                                                                                                   gAGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGat               <  1:545739/140‑1 (MQ=255)                                                                                                                    aGTCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGata              <  1:1143516/140‑1 (MQ=255)                                                                                                                      tCAGGCACTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATg            >  2:19640/1‑140 (MQ=255)                                                                                                                            aCTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATg      <  2:565780/140‑1 (MQ=255)                                                                                                                            aCTTTTGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATg      >  2:1049635/1‑140 (MQ=255)                                                                                                                                ttGCTACCAGGCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACg  >  1:341204/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                  TACAGCCAGGTGATACGTCCACCCCAGCAGAATCCGGTGATCATTAAACGATGAACATCGCCGCCGTTGCGGGAAGCCCAACTGGCGACATGATCGAGATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGACCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACG  >  NC_000913/4015703‑4015969

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