Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,016,062	A→G	P36P (CCT→CCC)	ysgA ←	putative carboxymethylenebutenolidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,016,062	0	A	G	100.0%	214.8 / NA	57	P36P (CCT→CCC)	ysgA	putative carboxymethylenebutenolidase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (30/27);  total (30/27)

CCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGC‑ATGGTAAGCAGGCATGTTAT‑CCCCTT‑GTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTT  >  NC_000913/4015914‑4016208                                                                                                                                                        |                                                                                                                                                   ccAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCGCTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGC‑ATGGTAAGCAGGCATGTTAT‑CCCCTT‑GTGAggg                                                                                                                                                  <  1:898969/151‑1 (MQ=255)    gATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGC‑ATGGTAAGCAGGCATGTTAT‑CCCCTT‑GTGAGGGGAt                                                                                                                                               >  1:155943/1‑151 (MQ=255)     aTACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGC‑ATGGTAAGCAGGCATGTTAT‑CCCCTT‑GTGAGGGGAt                                                                                                                                               >  1:726780/1‑150 (MQ=255)     aTACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGC‑ATGGTAAGCAGGCATGTTAT‑CCCCTT‑GTGAGGGGAt                                                                                                                                               >  1:272839/1‑150 (MQ=255)       aCCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCAAGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGC‑ATGGTAAGCAGGCATGTTAT‑CCCCTT‑GTGAGGGGAt                                                                                                                                               >  2:643943/1‑148 (MQ=255)                     aGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGC‑ATGGTAAGCAGGCATGTTAT‑CCCCTT‑GTGAGGGGATAGAGGTGAAGCCCGcc                                                                                                                               >  1:842999/1‑150 (MQ=255)                         ggCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGC‑ATGGTAAGCAGGCATGTTAT‑CCCCTT‑GTGAGGGGATAGAGGTGAAGCCCGCCACAAt                                                                                                                          >  2:108518/1‑151 (MQ=255)                           cGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGC‑ATGGTAAGCAGGCATGTTAT‑CCCCTT‑GTGAGGGGATAGAGGTGAAGCCCGCCACAATcg                                                                                                                        >  2:1007146/1‑151 (MQ=255)                              cAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGC‑ATGGTAAGCAGGCATGTTAT‑CCCCTT‑GTGAGGGGATAGAGGTGAAGCCCGCCACAATCGcgt                                                                                                                     >  2:776638/1‑151 (MQ=255)                              cAAATGTCGCGGATATGTTCATGCACACCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGC‑ATGGTAAGCAGGCATGTTAT‑CCCCTT‑GTGAGGGGATAGAGGTGAAGCCCGCCACAATCGcgt                                                                                                                     >  1:539193/1‑151 (MQ=255)                                  tGTCGCGGATATTTTCAGGCCCGCCAAAAATTTCCTGCACTAAAAGGCCCACTGGCAGTGGGCAATCGTTTTGTTTTGGTCGGGC‑ATGGTAAGCGGGCATGTTAC‑CCCCTT‑GTGAGGGGATAGAGGTGAAGCCCGCCACAATCGcgtcgt                                                                                                                  <  1:758914/150‑1 (MQ=255)                                     cgcgGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGC‑ATGGTAAGCAGGCATGTTAT‑CCCCTT‑GTGAGGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCgg                                                                                                              >  1:493553/1‑151 (MQ=255)                                       cgGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGC‑ATGGTAAGCAGGCATGTTAT‑CCCCTT‑GTGAGGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCgggg                                                                                                            >  1:640050/1‑151 (MQ=255)                                          atatGTTCATGCACGCCAAAAATTTCCTGCACTACAACGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGC‑ATGGTAAGCAGGCATGTTAT‑CCCCTT‑GTGAGGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCt                                                                                                         >  2:452604/1‑151 (MQ=255)                                             tGTTCATGCACGCAAAAATTTTCCTGCACTACAAGGCCCACGGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGC‑TTGTTAAGCAGGCATGTTAC‑CCCCTT‑GTGAGGGGATAGAGGTAAAGCCCGCCACAATCGCGTCGTCCGGGGTCTg                                                                                                        <  1:132728/149‑1 (MQ=255)                                               ttCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGC‑ATGGTAAGCAGGCATGTTAT‑CCCCTT‑GTGAGGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACg                                                                                                    >  2:290239/1‑151 (MQ=255)                                                       cGCAAAAATTTTCCGGCACTCAAAGCCCCAGGGGAAGTGGGCCATCGTTTTGTTTTGGTCTGGC‑TTGTTAAGCGGTCAGTTTCC‑CCCCTT‑GTGGGGGAATGGGGGTAAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGt                                                                                               <  2:549327/148‑1 (MQ=255)                                                          cAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGC‑ATGGTAAGCAGGCATGTTAT‑CCCCTT‑GTGAGGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGc                                                                                         >  1:1012493/1‑151 (MQ=255)                                                           aaaaaTTCCCTGCACTACAAGGACCACTGGCATTGGGCCATCGTTTTGCTTTGGTCTGGC‑ATGATAGGCAGGTATGTTTC‑CCCCTT‑GTGGGGGGATAGAGGTGAAGGCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGc                                                                                         <  1:535056/150‑1 (MQ=255)                                                             aaaTTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGC‑ATGGTAAGCAGGCATGTTAT‑CCCCTT‑GTGAGGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCg                                                                                        >  1:30760/1‑149 (MQ=255)                                                              attttCCTGCACTAAAATGCCCACGGGCAGTGGCCCATCCTTTTGTTTTGGTCTGGC‑ATGTTAAGCAGGCTTTTTTT‑CCCCTG‑GTGAGGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCgag                                                                                      <  2:549391/148‑1 (MQ=255)                                                                   ccTGCACTACAATGACCGCTGGCGGTGGGACATCTTTTTTCTTTGGTCTGGC‑ATGGAAAGCAGGAATGTTCT‑CCCTTT‑GGGAGGGGATAGAGGGGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGa                                                                                <  2:274853/151‑1 (MQ=255)                                                                       cACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGC‑ATGGTAAGCAGGCATGTTAT‑CCCCTT‑GTGAGGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGCAGCGAGAGGCGATgc                                                                             >  2:464321/1‑150 (MQ=255)                                                                         cTAAAATGCCCCCGGGCGGTGGGCCATCGTTTTGCTTTGGTCTGGC‑ATGGTAAGCAGGCATGTTAT‑CCCCTT‑GGGGGGGGAAAGGGGGGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATgcag                                                                           <  1:70453/150‑1 (MQ=255)                                                                           aCAATGACCGCGGGAGGGGGGCCTTGTTTTTGTTTGGGTCTGGC‑AGGGAAAGCAGGAATGTTAC‑CCCTTT‑GGGAGGGGATAGAGGTGAGGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATgca                                                                            <  2:18401/147‑1 (MQ=255)                                                                           aCAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGC‑ATGGTAAGCAGGCATGTTAT‑CCCCTT‑GTGAGGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATgcagc                                                                          >  2:473519/1‑149 (MQ=255)                                                                              atccccACTGGCGGGGGCCCATCGTTTTGTTTTGGTCGGGA‑TTGGTAAGGAGGAAT‑TTTT‑CCCCTTGGGGGGGGGATAGGGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTg                                                                     <  2:871429/147‑1 (MQ=255)                                                                              aTGCCCACGGGCGGTGGGCCATCGTTTTGTTTGGGTCTGGC‑ATGGTAAGCAGGCATTTTAT‑CCCCTT‑GTGAGGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTg                                                                     <  1:473544/151‑1 (MQ=255)                                                                                  ccACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGC‑ATGGTAAGCAGGCATGTTAT‑CCCCTT‑GTGAGGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCaa                                                                  >  2:797950/1‑150 (MQ=255)                                                                                      tGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGC‑ATGGTAAGCAGGCATGTTAT‑CCCCTT‑GTGAGGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATcc                                                              >  2:81489/1‑150 (MQ=255)                                                                                      gggCTGTGGGCCCCGGCTTTGCTTTGG‑CGGGCAATGGTAAGCAGGCAGTTTAT‑CCCCTT‑GGGAGGGGATAGGGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATcc                                                              <  2:25788/149‑1 (MQ=255)                                                                                          aGTGGGCCATCGCTTTGCTTTGGTCTGGC‑ATGGTAAGCAGGCATGTTAT‑CCCCTT‑GTGAGGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGAtt                                                         >  1:618471/1‑151 (MQ=255)                                                                                           gTGGGCCATCGCTTTGCTTTGGTCTGGC‑ATGATAAGCAGGCATGTTAT‑CCCCTT‑GTGAGGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGAttg                                                        >  1:939188/1‑151 (MQ=255)                                                                                              gcccACCGTTTTGTTTTGGTCTGGT‑ATGGAAGGCAGGCATGTTAC‑CCCTTT‑GGGAGGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGAttgttg                                                     <  2:1025679/149‑1 (MQ=255)                                                                                                 cATCGCTTTGCTTTGGTCTGGC‑ATGGTAAGCAGGCATGTTAT‑CCCCTT‑GTGAGGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgttgtgt                                                  >  2:167401/1‑151 (MQ=255)                                                                                                  aTCGCTTTGCTTTGGTCTGGC‑ATGGTAAGCAGGCATGTTAT‑CCCCTT‑GTGAGGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgttgtgtt                                                 >  1:1024782/1‑151 (MQ=255)                                                                                                  aTCGCTTTGCTTTGGTCTGGC‑ATGGTAAGCAGGCATGTTAT‑CCCCTT‑GTGAGGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgttgtgtt                                                 >  1:178267/1‑151 (MQ=255)                                                                                                   tCGCTTTGCTTTGGTCTGGC‑ATGGTAAGCAGGCATGTTAT‑CCCCTT‑GTGAGGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgttgtgttg                                                >  1:445371/1‑151 (MQ=255)                                                                                                   tCGCTTTGCTTTGGTCTGGC‑ATGGTAAGCAGGCATGTTAT‑CCCCTT‑GTGAGGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGCGCAAATCCAGATtgttgtgttg                                                >  2:234629/1‑151 (MQ=255)                                                                                                    cGCTTTGTTTTGGTCTGGC‑ATGGAAAGCAGGCATTTTAT‑CCCCTT‑GTGGGGGGATAGGGGTGAGGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgttgtgtt                                                 <  1:227247/149‑1 (MQ=255)                                                                                                      ctttgctttgGTCTGGC‑ATGGTAAGCAGGCATGTTAT‑CCCCTT‑GTGAGGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGttgttg                                             >  1:573243/1‑151 (MQ=255)                                                                                                       tttgttttgGTGGGGC‑GTGGTAAGCGGGCATGTTAC‑CCCTTT‑GGGAGGGGATAGGGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTgttgtgttgt                                               <  1:712193/148‑1 (MQ=255)                                                                                                       tttgctttgGTCGGGA‑ATGGTAAGCAGGCATGTTAC‑CCCCTT‑GGGGGGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTgttgtt                                              <  1:492932/149‑1 (MQ=255)                                                                                                            tttgTTGTGCC‑ATGGAAAGCAGGCTGGTTAC‑CCCCTG‑GTGAGGGGATAGAGGTGACGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCa                                          <  2:1023836/148‑1 (MQ=255)                                                                                                                gTCTGGC‑ATGGTAAGCAGGCATGTTAT‑CCCCTT‑GTGAGGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTAtt                                   >  1:630044/1‑151 (MQ=255)                                                                                                                   tGGC‑ATGGTAAGCAGGCATGTTAT‑CCCCTT‑GTGAGGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATtctc                                >  2:357186/1‑151 (MQ=255)                                                                                                                    ggC‑ATGGTAAGCAGGCATGTTAC‑CCCTTT‑GTGGGGGGATAGGGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATtctc                                <  1:378740/150‑1 (MQ=255)                                                                                                                       cagggTAAGCAGG‑ATGTTTTCCCCCTG‑GGGGGGGGAGAGGGGTGAACCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCg                              <  2:305803/147‑1 (MQ=255)                                                                                                                          ggTAGGCGGGCATGTTCC‑CCCCTT‑GTGGGGGGATAGGGGGGAGGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATtc                                  <  2:631475/143‑1 (MQ=255)                                                                                                                               gCAGGCATGTTAT‑CCCTTG‑GGGAGGGGATAGGGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTta                      <  2:390732/150‑1 (MQ=255)                                                                                                                                  ggCTTGTTAT‑CCCCTG‑GTGGGGGGATAGGGGGAAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTAccc                         <  2:301584/144‑1 (MQ=255)                                                                                                                                     aTGTTAT‑CCCCTT‑GTGAGGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAAt               >  2:47114/1‑151 (MQ=255)                                                                                                                                     aTGTTAC‑CCCTTT‑GTGAGGGGAGAGGGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaa                   <  2:696619/147‑1 (MQ=255)                                                                                                                                      ttttAC‑CCCTTT‑GTGGGGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAAt               <  1:597289/148‑1 (MQ=255)                                                                                                                                       ttttccccccTT‑GTGGGGGGATAGAGGTGAGGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAAt               <  1:110148/144‑1 (MQ=255)                                                                                                                                        ttAT‑CCCCTG‑GTGGGGGGATAGGGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAAt               <  1:324025/148‑1 (MQ=255)                                                                                                                                             ccccTT‑GTGAGGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTgc          <  1:940555/149‑1 (MQ=255)                                                                                                                                                     tgAGGGGATAGAGGTGAAGCCCGCCACGATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGtt  <  1:959968/150‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                   CCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGC‑ATGGTAAGCAGGCATGTTAT‑CCCCTT‑GTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTT  >  NC_000913/4015914‑4016208

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