Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,271,600	(C)5→6	coding (171/1146 nt)	garK ←	glycerate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,271,595	1	.	C	100.0%	228.6 / NA	61	A59G (GCT→GGT)	garK	glycerate kinase I
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (32/29);  total (32/29)

TCGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAG‑‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCT  >  NC_000913/3271457‑3271722                                                                                                                                            |                                                                                                                                 tCGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAg                                                                                                                                   <  1:370366/139‑1 (MQ=255)    cgTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAG‑‑ccc                                                                                                                              >  2:1242906/1‑139 (MQ=255)       ttttCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAg                                                                                                                                   >  1:1172103/1‑133 (MQ=255)       ttttCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAg                                                                                                                                   <  2:1172103/133‑1 (MQ=255)                         tCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAatcatc                                                                                                          <  2:60695/139‑1 (MQ=255)                                        gccgccATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAg                                                                                                                                   <  2:824869/100‑1 (MQ=255)                                        gccgccATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAg                                                                                                                                   >  1:824869/1‑100 (MQ=255)                                          cgccATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAg                                                                                                                                   >  1:1104014/1‑98 (MQ=255)                                          cgccATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAg                                                                                                                                   <  2:1104014/98‑1 (MQ=255)                                                       aaCGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCa                                                                            >  1:1173037/1‑139 (MQ=255)                                                         cgcgGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAAc                                                                          >  2:1120782/1‑139 (MQ=255)                                                           cgGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCg                                                                        >  1:1331778/1‑139 (MQ=255)                                                             gTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGa                                                                      >  1:749942/1‑139 (MQ=255)                                                             gTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAG‑‑c                                                                                                                                <  1:86016/80‑1 (MQ=255)                                                             gTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAG‑‑c                                                                                                                                >  2:86016/1‑80 (MQ=255)                                                              ttttGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGaa                                                                     <  1:1242906/139‑1 (MQ=255)                                                               tttGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAAc                                                                    >  1:1572264/1‑139 (MQ=255)                                                                 tGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTc                                                                                             <  2:277640/112‑1 (MQ=255)                                                                 tGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTc                                                                                             >  1:277640/1‑112 (MQ=255)                                                                   ccATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGaa                                                                <  2:1331778/139‑1 (MQ=255)                                                                   ccATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGaa                                                                <  2:1421151/139‑1 (MQ=255)                                                                     aTCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAAc                                                              <  2:499464/139‑1 (MQ=255)                                                                       cGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGt                                                            >  2:989155/1‑139 (MQ=255)                                                                        gCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAatc                                                                                                             <  1:1301137/89‑1 (MQ=255)                                                                        gCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAatc                                                                                                             >  2:1301137/1‑89 (MQ=255)                                                                         ccGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTAc                                                          <  1:1120782/139‑1 (MQ=255)                                                                            gagaTCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACtgt                                                       >  1:1616864/1‑139 (MQ=255)                                                                                  ccccAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAg                                                                                                                                   >  1:1287993/1‑58 (MQ=255)                                                                                  ccccAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAg                                                                                                                                   <  2:1287993/58‑1 (MQ=255)                                                                                  ccccAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCa                                                 >  2:1014204/1‑139 (MQ=255)                                                                                     cAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGa                                              <  1:140690/139‑1 (MQ=255)                                                                                          ggCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCAt                                                   <  1:269552/129‑1 (MQ=255)                                                                                          ggCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCAt                                                   >  2:269552/1‑129 (MQ=255)                                                                                           gCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTc                                                                                                     <  2:79548/78‑1 (MQ=255)                                                                                           gCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTc                                                                                                     >  1:79548/1‑78 (MQ=255)                                                                                           gCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCgg                                                                       <  2:878145/108‑1 (MQ=255)                                                                                           gCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCgg                                                                       >  1:878145/1‑108 (MQ=255)                                                                                           gCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAAtt                                        >  2:1606063/1‑139 (MQ=255)                                                                                           gCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAAtt                                        >  1:996703/1‑139 (MQ=255)                                                                                            cATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAAttt                                       >  2:480280/1‑139 (MQ=255)                                                                                             aTTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCa                                                                                     <  2:1041227/92‑1 (MQ=255)                                                                                             aTTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCa                                                                                     >  1:1041227/1‑92 (MQ=255)                                                                                             aTTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAG‑‑cccccctgg                                                                                                                        >  1:1303911/1‑52 (MQ=255)                                                                                             aTTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAG‑‑cccccctgg                                                                                                                        <  2:1303911/56‑5 (MQ=255)                                                                                                 aCTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGaa                                                                     <  2:1527745/104‑1 (MQ=255)                                                                                                 aCTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGaa                                                                     >  1:1527745/1‑104 (MQ=255)                                                                                                   tttCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAAc                                                                          <  1:404297/97‑1 (MQ=255)                                                                                                   tttCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAAc                                                                          >  2:404297/1‑97 (MQ=255)                                                                                                   tttCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGaa                                                                     <  2:571856/102‑1 (MQ=255)                                                                                                   tttCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGaa                                                                     >  1:571856/1‑102 (MQ=255)                                                                                                    ttCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACATAAACGTACTGTGCATCAGGAAAAATTTCCCGAAAt                               <  1:817232/139‑1 (MQ=255)                                                                                                    ttCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAAt                               <  1:989155/139‑1 (MQ=255)                                                                                                     tctcGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATc                              >  1:1399563/1‑139 (MQ=255)                                                                                                     tctcGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATc                              >  1:973398/1‑139 (MQ=255)                                                                                                      ctcGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATcc                             <  1:572654/139‑1 (MQ=255)                                                                                                         gCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTa                                                           >  2:816282/1‑106 (MQ=255)                                                                                                         gCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTa                                                           <  1:816282/106‑1 (MQ=255)                                                                                                          cccAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCtttt                         >  2:810250/1‑139 (MQ=255)                                                                                                           ccAGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCttttt                        >  2:914504/1‑139 (MQ=255)                                                                                                           ccAGCGGCCCTGTAACCCAGGCGTGACGTTCAG‑‑ccc                                                                                                                              <  2:675857/36‑1 (MQ=255)                                                                                                           ccAGCGGCCCTGTAACCCAGGCGTGACGTTCAG‑‑ccc                                                                                                                              >  1:675857/1‑36 (MQ=255)                                                                                                             aGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACa                                                                   <  1:919616/94‑1 (MQ=255)                                                                                                             aGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACa                                                                   >  2:919616/1‑94 (MQ=255)                                                                                                             aGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAAtt                                        <  2:364477/121‑1 (MQ=255)                                                                                                             aGCGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAAtt                                        >  1:364477/1‑121 (MQ=255)                                                                                                               cGGCCCTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTAt                    >  1:568537/1‑139 (MQ=255)                                                                                                                  cccTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGc                 <  1:284067/139‑1 (MQ=255)                                                                                                                    cTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGaa                                                                     <  2:545739/85‑1 (MQ=255)                                                                                                                    cTGTAACCCAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGaa                                                                     >  1:545739/1‑85 (MQ=255)                                                                                                                          cccAGGCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCaa         <  1:914504/139‑1 (MQ=255)                                                                                                                              ggCGTGACGTTCAGCCCCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCt      <  1:844702/139‑1 (MQ=255)                                                                                                                              ggCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTc     >  2:165418/1‑139 (MQ=255)                                                                                                                               gCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCtt                                                                                                      <  1:1291924/41‑1 (MQ=38)                                                                                                                               gCGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCtt                                                                                                      >  2:1291924/1‑41 (MQ=38)                                                                                                                                cGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCt                                                                                                       <  1:1120700/39‑1 (MQ=255)                                                                                                                                cGTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCt                                                                                                       >  2:1120700/1‑39 (MQ=255)                                                                                                                                 gTGACGTTCAGC‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCt  <  2:973398/139‑1 (MQ=255)                                                                                                                                            |                                                                                                                                 TCGCGTTTTTCCGCAGGTACCAGCTCCAGCCCACTGGCCGCCGCCATTTCAATAAACGCGGTTTTGCCATCGCCGGAGATCCCCCAACTGGCATTCACTTTCTCGCCCAGCGGCCCTGTAACCCAGGCGTGACGTTCAG‑‑CCCCCTGGGTGGCTGCAATCATCGCTTCCACCGTTCCTTCGCCACCGTCGGCAACCGGAACAGAAACGTACTGTGCATCAGGAAAAATTTCCCGAAATCCTTTTTCTATCGCCTGCGCAACCTCGCT  >  NC_000913/3271457‑3271722

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