Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	289,154	A→C	I3M (ATT→ATG)	yagI ←	CP4‑6 prophage; putative DNA‑binding transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	289,154	0	A	C	100.0%	225.0 / NA	60	I3M (ATT→ATG)	yagI	CP4‑6 prophage; putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (28/32);  total (28/32)

AGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACC  >  NC_000913/289009‑289297                                                                                                                                                  |                                                                                                                                                aGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCg                                                                                                                                             <  2:790057/150‑1 (MQ=255)  gCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGc                                                                                                                                           >  2:534479/1‑151 (MQ=255)  gCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGc                                                                                                                                           <  2:80939/151‑1 (MQ=255)    cGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCAt                                                                                                                                         >  1:72682/1‑151 (MQ=255)     gTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATa                                                                                                                                        <  2:848449/151‑1 (MQ=255)             gCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTg                                                                                                                                  <  1:142827/149‑1 (MQ=255)              cAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGc                                                                                                                                >  2:917019/1‑150 (MQ=255)                gggTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATg                                                                                                                             >  2:805432/1‑151 (MQ=255)                 ggTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGt                                                                                                                            <  1:402964/151‑1 (MQ=255)                       tAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACg                                                                                                                       <  2:252187/150‑1 (MQ=255)                       tAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGc                                                                                                                      >  2:529718/1‑151 (MQ=255)                       tAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGc                                                                                                                      >  2:487123/1‑151 (MQ=255)                       tAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGc                                                                                                                      >  1:716519/1‑151 (MQ=255)                          cagcGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTgg                                                                                                                   >  1:917997/1‑151 (MQ=255)                              gAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAAtt                                                                                                               >  2:4103/1‑151 (MQ=255)                              gAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAAtt                                                                                                               >  1:899457/1‑151 (MQ=255)                                    aGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCg                                                                                                         >  1:220014/1‑151 (MQ=255)                                     gggTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCg                                                                                                         <  2:165895/150‑1 (MQ=255)                                        tACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGaaa                                                                                                      >  1:552433/1‑150 (MQ=255)                                            ctTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACg                                                                                                   <  2:104007/149‑1 (MQ=255)                                              tGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTAc                                                                                                <  1:644029/150‑1 (MQ=255)                                                  cAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCttt                                                                                            <  1:88340/150‑1 (MQ=255)                                                      cccATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCgg                                                                                       <  2:22881/151‑1 (MQ=255)                                                               ttGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCCCGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGaaaa                                                                               >  1:319697/1‑150 (MQ=255)                                                               ttGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGaaaa                                                                               <  1:735636/150‑1 (MQ=255)                                                                    gATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCaa                                                                          <  1:768739/150‑1 (MQ=255)                                                                      tGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATg                                                                        >  1:664282/1‑150 (MQ=255)                                                                         cGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAAt                                                                    <  1:404508/151‑1 (MQ=255)                                                                         cGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAAt                                                                    >  1:95623/1‑151 (MQ=255)                                                                         cGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAAt                                                                    <  2:862948/151‑1 (MQ=255)                                                                          ggTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATg                                                                   <  2:661756/151‑1 (MQ=255)                                                                           gTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATgg                                                                  >  1:120349/1‑151 (MQ=255)                                                                             gATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGa                                                                 <  2:207169/150‑1 (MQ=255)                                                                               tCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGcc                                                              >  2:667191/1‑151 (MQ=255)                                                                                cTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCg                                                             >  2:167683/1‑151 (MQ=255)                                                                                  tAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCg                                                             <  1:374660/149‑1 (MQ=255)                                                                                  tAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCg                                                             <  1:513454/149‑1 (MQ=255)                                                                                   aaGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGgtg                                                          >  1:544539/1‑151 (MQ=255)                                                                                       tCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTtct                                                      >  2:765329/1‑151 (MQ=255)                                                                                                 gCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTctccctctccct                                            <  1:833260/151‑1 (MQ=255)                                                                                                    cGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCtgtg                                         <  1:609298/151‑1 (MQ=255)                                                                                                       tGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTggg                                       <  1:406623/150‑1 (MQ=255)                                                                                                        gAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGgag                                     <  2:95132/151‑1 (MQ=255)                                                                                                        gAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGgag                                     <  2:855769/151‑1 (MQ=255)                                                                                                         aaCAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGgaga                                    <  2:587056/151‑1 (MQ=255)                                                                                                         aaCAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGgaga                                    <  2:716501/151‑1 (MQ=255)                                                                                                          aCAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGagag                                   <  1:335951/151‑1 (MQ=255)                                                                                                            aGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAggg                                 <  1:499957/151‑1 (MQ=255)                                                                                                                cGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCgg                             <  2:755046/151‑1 (MQ=255)                                                                                                                cGAGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCgg                             >  2:705341/1‑151 (MQ=255)                                                                                                                  aGGATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCgggg                           <  1:634961/151‑1 (MQ=255)                                                                                                                   ggATCTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGt                          <  2:316625/151‑1 (MQ=255)                                                                                                                       cTGCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAgg                      >  2:226972/1‑151 (MQ=255)                                                                                                                         gCAACGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                       <  2:111042/148‑1 (MQ=255)                                                                                                                             cGCACGTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAcc                 >  2:888036/1‑150 (MQ=255)                                                                                                                                  gTTCAACAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAAcgcg           >  2:197355/1‑151 (MQ=255)                                                                                                                                        cAGACTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCgg     >  1:775443/1‑151 (MQ=255)                                                                                                                                           aCTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGAcc  >  1:886782/1‑151 (MQ=255)                                                                                                                                           aCTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGAcc  >  1:806173/1‑151 (MQ=255)                                                                                                                                           aCTGAATCATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGAcc  >  1:204853/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                AGCCGTGAAGCTGCAGGGTTTTTAGCAGCGAGTGGAGGGTACTCTTGCTCAGCCCCATCAGTTTGCTGATGTCGGTGATCTTAAGCTCGGTGGCCTGCTCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACC  >  NC_000913/289009‑289297

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