Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,619,913	A→G	E104G (GAA→GGA)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,619,913	0	A	G	100.0%	252.4 / NA	72	E104G (GAA→GGA)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (38/34);  total (38/34)

GCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGG  >  NC_000913/4619786‑4620048                                                                                                                                |                                                                                                                                        gCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGgtg                                                                                                                              >  2:1617386/1‑139 (MQ=255)    aaaGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTcc                                                                                                                           <  2:172728/139‑1 (MQ=255)     aaGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAGGTTACCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTccc                                                                                                                          <  1:1259747/139‑1 (MQ=255)       gCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCgg                                                                                                                        <  1:1617386/139‑1 (MQ=255)        cacaCGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATG‑CATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGtt                                                                                                                      <  2:563215/139‑1 (MQ=255)                 gTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGa                                                                                                              <  1:1535154/139‑1 (MQ=255)                  ttCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAg                                                                                                             >  1:575783/1‑139 (MQ=255)                   tCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGt                                                                                                            <  2:575783/139‑1 (MQ=255)                   tCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGt                                                                                                            <  1:1883536/139‑1 (MQ=255)                    cTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTg                                                                                                           >  1:1699186/1‑139 (MQ=255)                    cTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTg                                                                                                           >  1:1927117/1‑139 (MQ=255)                     tACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTgg                                                                                                          >  2:1197313/1‑139 (MQ=255)                       ccGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTccc                                                                                                                          <  1:1275838/121‑1 (MQ=255)                       ccGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTccc                                                                                                                          >  2:1275838/1‑121 (MQ=255)                          gTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGata                                                                                                     <  2:304927/139‑1 (MQ=255)                           tttCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGa                                                                                                              >  1:1739840/1‑129 (MQ=255)                           tttCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGa                                                                                                              <  2:1739840/129‑1 (MQ=255)                            ttCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATAtt                                                                                                   >  2:622041/1‑139 (MQ=255)                               aTTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTtct                                                                                                >  1:92503/1‑139 (MQ=255)                                       gggAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCAc                                                                                        <  1:1197313/139‑1 (MQ=255)                                                  gCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTATAGATACCCCGTCTGGTCACCGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTTGGGATATTTCTCCGATCACGAAAACAGCtt                                                                             >  1:979956/1‑139 (MQ=255)                                                      cGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCg                                                                         >  2:58997/1‑139 (MQ=255)                                                        cTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCa                                                                       >  2:860795/1‑139 (MQ=255)                                                              ttATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAgct                                                                 <  2:92503/139‑1 (MQ=255)                                                                aTCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTc                                                                                                                            >  1:1307071/1‑78 (MQ=255)                                                                aTCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTc                                                                                                                            <  2:1307071/78‑1 (MQ=255)                                                                aTCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTgg                                                                                                          <  2:1172946/96‑1 (MQ=255)                                                                aTCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTgg                                                                                                          >  1:1172946/1‑96 (MQ=255)                                                                aTCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAagag                                                                   <  2:589561/135‑1 (MQ=255)                                                                aTCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAagag                                                                   >  1:589561/1‑135 (MQ=255)                                                                                gggCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTgg                                                                                                          <  2:263987/80‑1 (MQ=255)                                                                                gggCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTgg                                                                                                          >  1:263987/1‑80 (MQ=255)                                                                                gggCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGa                                               <  1:58997/139‑1 (MQ=255)                                                                                 ggCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGaa                                              >  2:599722/1‑139 (MQ=255)                                                                                   cgcgcgAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCg                                                                                             <  2:1949145/90‑1 (MQ=255)                                                                                   cgcACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGtt                                                                                                                      >  1:1584310/1‑65 (MQ=255)                                                                                   cgcACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGtt                                                                                                                      <  2:1584310/65‑1 (MQ=255)                                                                                   cgcACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGCTATTTCTCCg                                                                                             >  1:1949145/1‑90 (MQ=255)                                                                                     cACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAAt                                                                                                                                      <  2:363984/47‑1 (MQ=255)                                                                                     cACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAAt                                                                                                                                      >  1:363984/1‑47 (MQ=255)                                                                                     cACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAAcgcg                                          >  2:138861/1‑139 (MQ=255)                                                                                            tGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCt                                   <  1:599722/139‑1 (MQ=255)                                                                                              tCATCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGc                                 <  2:1927117/139‑1 (MQ=255)                                                                                                 tCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGa                                               <  2:982571/122‑1 (MQ=255)                                                                                                 tCCGGTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGa                                               >  1:982571/1‑122 (MQ=255)                                                                                                    ggTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAgctgc                                                               <  2:1865437/103‑1 (MQ=255)                                                                                                    ggTAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAgctgc                                                               >  1:1865437/1‑103 (MQ=255)                                                                                                    ggTAAAGATACCCCGTCTGGTCACTGGGGAATTGACGGTGTCCCGGTTCAGTTTGAGTGGGCATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGtt                           >  1:1330129/1‑139 (MQ=255)                                                                                                      tAAAGATACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTAc                         >  2:605554/1‑139 (MQ=255)                                                                                                      gtatctttacccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTAc                         >  2:617535/10‑139 (MQ=255)                                                                                                            tACCCCGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGt                   <  1:860795/139‑1 (MQ=255)                                                                                                              ccccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTaa                 >  2:1784978/1‑139 (MQ=255)                                                                                                               cccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAc                <  2:1699186/139‑1 (MQ=255)                                                                                                               cccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAc                >  2:1935964/1‑139 (MQ=255)                                                                                                                ccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAgct                                                                 >  1:479901/1‑89 (MQ=255)                                                                                                                ccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAgct                                                                 <  2:479901/89‑1 (MQ=255)                                                                                                                ccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATa                                                         <  2:1069391/97‑1 (MQ=255)                                                                                                                ccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATa                                                         >  1:1069391/1‑97 (MQ=255)                                                                                                                ccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGtt                           >  2:1193889/1‑127 (MQ=255)                                                                                                                ccGTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGtt                           <  1:1193889/127‑1 (MQ=255)                                                                                                                  gTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAagag                                                                   >  1:1282392/1‑85 (MQ=255)                                                                                                                  gTCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAagag                                                                   <  2:1282392/85‑1 (MQ=255)                                                                                                                   tCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGa                                               >  2:785520/1‑104 (MQ=255)                                                                                                                   tCTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGa                                               <  1:785520/104‑1 (MQ=255)                                                                                                                    cTGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCa           <  1:605554/139‑1 (MQ=255)                                                                                                                     tGGTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCAc          >  1:1689638/1‑139 (MQ=255)                                                                                                                      ggTCACTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCAct         >  2:1939028/1‑139 (MQ=255)                                                                                                                           cTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGaaaa                                                                                   >  2:885869/1‑60 (MQ=255)                                                                                                                           cTGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGaaaa                                                                                   <  1:885869/60‑1 (MQ=255)                                                                                                                            tGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAAGTGGTCGAACGCGCTGATCTGCCGGGTTACCTCGGTAACTGCCACTCTTcc    <  1:472470/138‑1 (MQ=255)                                                                                                                            tGGGGAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCg   >  2:1882898/1‑139 (MQ=255)                                                                                                                            tGGGGAATTGCCGGTGTCCCGGTTCTGGTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTcc    >  2:472470/1‑138 (MQ=255)                                                                                                                             ggggAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCgg  <  1:1939028/139‑1 (MQ=255)                                                                                                                                |                                                                                                                                        GCGAAAGCACACGAAGGTTCTACCGGTTTCATTCCGGCGGGAATGGACGGCAACGCTGAAGTTATCGGCGCGTACGCATGGGCGCACGAAATGTCATCCGGTAAAGATACCCCGTCTGGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGG  >  NC_000913/4619786‑4620048

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