Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,585,101	(C)6→7	coding (1386/2280 nt)	ydeP ←	putative oxidoreductase YdeP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,585,095	1	.	C	100.0%	202.4 / NA	57	G464G (GGC→GGG)	ydeP	putative oxidoreductase YdeP
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (29/28);  total (29/28)

GGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑G‑‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGG  >  NC_000913/1584959‑1585229                                                                                                                                           |                                                                                                                                        ggACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑G‑‑ccc                                                                                                                                     <  2:9704/140‑1 (MQ=255)   aCCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑G‑‑ccccc                                                                                                                                   <  1:797161/140‑1 (MQ=255)   aCCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑G‑‑ccccc                                                                                                                                   <  1:712337/140‑1 (MQ=255)       gCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑g                                                                                                                                          >  2:897126/1‑131 (MQ=255)       gCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑g                                                                                                                                          <  1:897126/131‑1 (MQ=255)       gCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑G‑‑cccccccat                                                                                                                               >  2:550386/1‑137 (MQ=255)       gCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑G‑‑cccccccat                                                                                                                               >  2:12284/1‑137 (MQ=255)       gCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑G‑‑cccccccat                                                                                                                               >  1:944050/1‑137 (MQ=255)                     tGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGTCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTc                                                                                                                  >  2:814940/1‑140 (MQ=255)                     tGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCt                                                                                                                   >  1:1023042/1‑139 (MQ=255)                     tGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCt                                                                                                                   <  2:1023042/139‑1 (MQ=255)                      gCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGTTCg                                                                                                                 >  2:761541/1‑140 (MQ=255)                       ccGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑G‑‑ccccccc                                                                                                                                 >  1:448219/1‑121 (MQ=255)                       ccGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑G‑‑ccccccc                                                                                                                                 <  2:448219/122‑2 (MQ=255)                         gTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGc                                                                                                              <  1:313317/140‑1 (MQ=255)                          tGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGcc                                                                                                             >  1:359182/1‑140 (MQ=255)                          tGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGcc                                                                                                             >  2:492118/1‑140 (MQ=255)                          tGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGcc                                                                                                             >  1:52603/1‑140 (MQ=255)                             ggTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGc                                                                                                          <  2:836607/140‑1 (MQ=255)                             ggTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGc                                                                                                          <  1:550386/140‑1 (MQ=255)                             ggTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGc                                                                                                          <  1:854093/140‑1 (MQ=255)                             ggTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGc                                                                                                          <  2:103234/140‑1 (MQ=255)                              gTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGcc                                                                                                         <  2:846340/140‑1 (MQ=255)                               tCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGccc                                                                                                        <  1:640562/140‑1 (MQ=255)                                cAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGcccc                                                                                                       <  2:359182/140‑1 (MQ=255)                                 aaCAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCg                                                                                                      <  2:52603/140‑1 (MQ=255)                                    aGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTAc                                                                                                   >  2:19499/1‑140 (MQ=255)                                         agagCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGcc                                                                                                         <  2:455001/129‑1 (MQ=255)                                         agagCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGcc                                                                                                         >  1:455001/1‑129 (MQ=255)                                           agCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑GCCCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATAttg                                                                                             >  2:507245/1‑140 (MQ=255)                                             cGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATAttgctt                                                                                          >  2:888226/1‑140 (MQ=255)                                                gtagggttAGTGGCTACGGGTACCGCTAAATCTAATTCCGTTAACGGTACAGCGCTCGCTTCCCGATCGGGCATTGCCAACGCAAAATTAG‑‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCa                                                                                       >  2:290849/7‑140 (MQ=255)                                                              tACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTgca                                                                         >  2:886505/1‑140 (MQ=255)                                                                cGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTgcagc                                                                       <  1:788414/140‑1 (MQ=255)                                                                    tACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATg                                                                    <  2:831534/139‑1 (MQ=255)                                                                    tACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATg                                                                    >  1:831534/1‑139 (MQ=255)                                                                    tACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGt                                                                   >  2:375284/1‑140 (MQ=255)                                                                       cGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCa                                                                <  1:888226/140‑1 (MQ=255)                                                                        gCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTg                                                                           >  2:337668/1‑128 (MQ=255)                                                                        gCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTg                                                                           <  1:337668/128‑1 (MQ=255)                                                                                     gCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTgg                                                  <  2:944050/140‑1 (MQ=255)                                                                                         tAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTg                                              <  1:375284/140‑1 (MQ=255)                                                                                                aCAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTgg                                                  <  1:819337/129‑1 (MQ=255)                                                                                                aCAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTgg                                                  >  2:819337/1‑129 (MQ=255)                                                                                                    cgctcgctTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTg                                              >  1:833053/1‑129 (MQ=255)                                                                                                    cgctcgctTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTg                                              <  2:833053/129‑1 (MQ=255)                                                                                                    cgctcgctTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAgcg                                   <  1:12284/140‑1 (MQ=255)                                                                                                    cgctcgctTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAgcg                                   >  1:719285/1‑140 (MQ=255)                                                                                                    cgctcgctTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAgcg                                   <  2:719285/140‑1 (MQ=255)                                                                                                    cgctcgctTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAgcg                                   >  2:458983/1‑140 (MQ=255)                                                                                                     gctcgctTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAgcgc                                  >  2:375522/1‑140 (MQ=255)                                                                                                      ctcgctTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAgcgc                                  >  1:686906/1‑139 (MQ=255)                                                                                                      ctcgctTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAgcgc                                  <  2:686906/139‑1 (MQ=255)                                                                                                       tcgctTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTc                                >  1:367728/1‑140 (MQ=255)                                                                                                       tcgctTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTc                                >  1:698047/1‑140 (MQ=255)                                                                                                         gctTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCAc                              >  2:386934/1‑140 (MQ=255)                                                                                                          ctTCCCGATCTGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCAcc                             <  1:386934/140‑1 (MQ=255)                                                                                                                    tGGCATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGcc                   <  1:375522/140‑1 (MQ=255)                                                                                                                       cATTGCCAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGa                >  2:19136/1‑140 (MQ=255)                                                                                                                            ccAGCGCAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAActc           >  1:160811/1‑140 (MQ=255)                                                                                                                                cgcAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCa       >  1:726551/1‑140 (MQ=255)                                                                                                                                cgcAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCa       >  2:801438/1‑140 (MQ=255)                                                                                                                                cgcAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCa       >  1:650584/1‑140 (MQ=255)                                                                                                                                  cAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGa     <  2:650584/140‑1 (MQ=255)                                                                                                                                  cAAAATT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGa     <  1:19136/140‑1 (MQ=255)                                                                                                                                    aaaTT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACg   <  1:801438/140‑1 (MQ=255)                                                                                                                                    aaaTT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACg   >  2:119402/1‑140 (MQ=255)                                                                                                                                     aaTT‑GC‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACgg  <  1:1014897/140‑1 (MQ=255)                                                                                                                                           |                                                                                                                                        GGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATT‑G‑‑CCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGG  >  NC_000913/1584959‑1585229

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