Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,104,065	G→A	*351* (TAG→TAA)	mutY →	adenine DNA glycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,104,065	0	G	A	100.0%	254.5 / NA	78	*351* (TAG→TAA)	mutY	adenine DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (36/42);  total (36/42)

CCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAGCGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGAAAACGCATCTATAACGAG  >  NC_000913/3103931‑3104191                                                                                                                                       |                                                                                                                               cccGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAAcgcg                                                                                                                            <  2:1275746/139‑1 (MQ=255) cccGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAAcgcg                                                                                                                            >  1:238987/1‑139 (MQ=255)   cGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTg                                                                                                                          <  1:494595/139‑1 (MQ=255)      gtcgtcATTCACCGGCTGCTTGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGt                                                                                                                       >  1:1339696/1‑139 (MQ=255)              tCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAag                                                                                                               <  2:182202/139‑1 (MQ=255)                aCCGGCTGCTTGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAagag                                                                                                             <  2:1339696/139‑1 (MQ=255)                    gCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTg                                                                                                                          <  2:648342/122‑1 (MQ=255)                    gCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTg                                                                                                                          >  1:648342/1‑122 (MQ=255)                    gCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAg                                                                                                                        >  1:109363/1‑124 (MQ=255)                    gCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAg                                                                                                                        <  2:109363/124‑1 (MQ=255)                    gCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGgatg                                                                                                         >  1:1604894/1‑139 (MQ=255)                    gCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGgatg                                                                                                         >  1:858101/1‑139 (MQ=255)                        cATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTg                                                                                                                          <  2:42711/118‑1 (MQ=255)                        cATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTg                                                                                                                          >  1:42711/1‑118 (MQ=255)                               gAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAg                                                                                                <  1:1181465/137‑1 (MQ=255)                               gAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAACGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAg                                                                                                >  2:1181465/1‑137 (MQ=255)                                      aTGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGAt                                                                                       >  2:895615/1‑139 (MQ=255)                                       tGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAag                                                                                                               >  2:1438765/1‑114 (MQ=255)                                       tGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAag                                                                                                               <  1:1438765/114‑1 (MQ=255)                                       tGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGAtt                                                                                      >  1:1563343/1‑139 (MQ=255)                                       tGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGAtt                                                                                      >  1:931574/1‑139 (MQ=255)                                            tctGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGGCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTg                                                                                 <  1:553190/139‑1 (MQ=255)                                              tGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAagag                                                                                                             >  1:933796/1‑109 (MQ=255)                                              tGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAagag                                                                                                             <  2:933796/109‑1 (MQ=255)                                                      cTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTg                                                                       <  2:858101/139‑1 (MQ=255)                                                      cTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTg                                                                       >  1:994316/1‑139 (MQ=255)                                                      cTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTg                                                                       >  2:619103/1‑139 (MQ=255)                                                      cTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTg                                                                       <  2:931574/139‑1 (MQ=255)                                                        tAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGgat                                                                                                          <  2:1123601/102‑1 (MQ=255)                                                        tAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGgat                                                                                                          >  1:1123601/1‑102 (MQ=255)                                                          gcgcAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTa                                                                               >  1:1132871/1‑127 (MQ=255)                                                          gcgcAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTa                                                                               <  2:1132871/127‑1 (MQ=255)                                                          gcgcAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAAc                                                                   >  1:1593777/1‑139 (MQ=255)                                                           cgcAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAAcgcg                                                                                                                            >  1:1687376/1‑81 (MQ=255)                                                           cgcAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAAcgcg                                                                                                                            <  2:1687376/81‑1 (MQ=255)                                                           cgcAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACg                                                                  <  1:1679349/139‑1 (MQ=255)                                                             cAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTg                                                                <  2:994316/139‑1 (MQ=255)                                                              aaCCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGa                                                               <  2:1071760/139‑1 (MQ=255)                                                                ccgccgTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACg                                                                                         >  1:153845/1‑111 (MQ=255)                                                                ccgccgTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACg                                                                                         <  2:153845/111‑1 (MQ=255)                                                                  gccgTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCa                                                           <  2:1386408/139‑1 (MQ=255)                                                                    cgTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGa                                                         <  2:1563343/139‑1 (MQ=255)                                                                          ttGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCa                                                   <  2:238987/139‑1 (MQ=255)                                                                            ggCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAgg                                                 >  1:1025199/1‑139 (MQ=255)                                                                             gCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTa                                                                               >  2:299899/1‑108 (MQ=255)                                                                             gCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTa                                                                               <  1:299899/108‑1 (MQ=255)                                                                                  gCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTc                                           <  1:619103/139‑1 (MQ=255)                                                                                       tCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACg                                                                                         >  2:196185/1‑88 (MQ=255)                                                                                       tCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACg                                                                                         <  1:196185/88‑1 (MQ=255)                                                                                       tCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTg                                      >  2:248041/1‑139 (MQ=255)                                                                                        cccGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGt                                     >  2:1670770/1‑139 (MQ=255)                                                                                            tGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGgat                                                                                                          <  1:1619343/66‑1 (MQ=255)                                                                                            tGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGgat                                                                                                          >  2:1619343/1‑66 (MQ=255)                                                                                             ggAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTAcccc                                >  2:1574981/1‑139 (MQ=255)                                                                                               aGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCgg                              <  1:248041/139‑1 (MQ=255)                                                                                                   tttGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAg                          >  2:1163538/1‑139 (MQ=255)                                                                                                       ttACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTgg                      <  2:1604894/139‑1 (MQ=255)                                                                                                        tACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTa                                                                               <  1:518593/81‑1 (MQ=255)                                                                                                        tACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTa                                                                               >  2:518593/1‑81 (MQ=255)                                                                                                        tACAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTggg                     <  2:1288739/139‑1 (MQ=255)                                                                                                         aCAGCAGTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGAGAGCTGGGa                    <  1:1163538/139‑1 (MQ=255)                                                                                                             cagTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAag                                                                                                               <  2:837657/44‑1 (MQ=255)                                                                                                             cagTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAag                                                                                                               <  2:511647/44‑1 (MQ=255)                                                                                                             cagTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAag                                                                                                               >  1:837657/1‑44 (MQ=255)                                                                                                             cagTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAag                                                                                                               >  1:511647/1‑44 (MQ=255)                                                                                                              agTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATaaa                                                                                                                >  2:291793/1‑42 (MQ=255)                                                                                                              agTTACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATaaa                                                                                                                <  1:291793/42‑1 (MQ=255)                                                                                                                 tACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGAAAACGCAt            <  1:1574981/139‑1 (MQ=255)                                                                                                                 tACGCACTGGCGCGCCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGAAAACGCAt            <  2:1593777/139‑1 (MQ=255)                                                                                                                          gcgcgcCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAg                                                                                                <  2:717150/46‑1 (MQ=255)                                                                                                                          gcgcgcCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAg                                                                                                >  1:717150/1‑46 (MQ=255)                                                                                                                          gcgcgcCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAgg                                                      <  2:1190289/88‑1 (MQ=255)                                                                                                                          gcgcgcCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAgg                                                      >  1:1190289/1‑88 (MQ=255)                                                                                                                          gcgcgcCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGAAAACGCATCTATAACga   >  2:485700/1‑139 (MQ=255)                                                                                                                           cgcgcCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGAAAACGCATCTATAACgag  <  1:895615/139‑1 (MQ=255)                                                                                                                           cgcgcCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTCTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGAAAACGCATCTATAACgag  <  2:1025199/139‑1 (MQ=255)                                                                                                                             cgcCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTa                                                                               <  2:818144/60‑1 (MQ=255)                                                                                                                             cgcCGGTTTAACGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTa                                                                               >  1:818144/1‑60 (MQ=255)                                                                                                                                       |                                                                                                                               CCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGCACTGGCGCGCCGGTTTAGCGCGTGAGTCGATAAAGAGGATGATTTATGAGCAGAACGATTTTTTGTACTTTCCTGCAACGTGAAGCAGAAGGTCAGGATTTTCAGCTGTACCCCGGCGAGCTGGGAAAACGCATCTATAACGAG  >  NC_000913/3103931‑3104191

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