BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,725,518	(G)_6→5	coding (836/1098 nt)	nemA →	N‑ethylmaleimide reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,725,513	0	G	.	100.0%	421.9 / NA	93	coding (831/1098 nt)	nemA	N‑ethylmaleimide reductase
Reads supporting (aligned to +/- strand): ref base G (0/0); new base . (49/44); total (49/44)

Only 100 of 104 total aligned reads displayed.

GGCCCGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATACACAGTA  >  NC_000913/1725421‑1725600

ggCCCGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCg                                                                                         <  2:2738564/93‑1 (MQ=255)

  ccggAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCggg                                                                                       <  2:2527338/90‑1 (MQ=255)

  cccGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCggg                                                                                       <  2:2811570/93‑1 (MQ=255)

   ccGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCgggg                                                                                      >  2:1256226/1‑93 (MQ=255)

    cGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCggggg                                                                                     >  1:389367/1‑93 (MQ=255)

    cGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCggggg                                                                                     <  1:4448068/93‑1 (MQ=255)

     gAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCgggggt                                                                                    >  2:2254911/1‑92 (MQ=255)

        tGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCgggggggaa                                                                                 >  1:3840370/1‑90 (MQ=255)

         gaagaaGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAAc                                                                               >  1:4594370/1‑93 (MQ=255)

          aagaagCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAAcc                                                                              <  2:4568481/93‑1 (MQ=255)

          aagaagCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAAcc                                                                              <  2:134018/93‑1 (MQ=255)

           agaagCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCg                                                                             >  1:3437033/1‑93 (MQ=255)

           agaagCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCg                                                                             >  1:4011346/1‑93 (MQ=255)

             aagCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGta                                                                           <  2:1320839/93‑1 (MQ=255)

             aagCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGta                                                                           <  2:2810800/93‑1 (MQ=255)

             aagCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGta                                                                           >  2:4235421/1‑93 (MQ=255)

              agCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGtat                                                                          >  1:4014254/1‑93 (MQ=255)

                ccGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATAc                                                                        <  2:1177535/93‑1 (MQ=255)

                  gATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTg                                                                      >  2:260067/1‑93 (MQ=255)

                  gATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTg                                                                      >  1:1133391/1‑93 (MQ=255)

                   aTGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGa                                                                     >  2:4597365/1‑93 (MQ=255)

                   aTGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGa                                                                     >  2:1181545/1‑93 (MQ=255)

                    tGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGAt                                                                    >  1:89709/1‑93 (MQ=255)

                    tGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGAt                                                                    >  2:3303989/1‑93 (MQ=255)

                    tGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGAt                                                                    <  2:426713/93‑1 (MQ=255)

                    tGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGAt                                                                    >  1:1164340/1‑93 (MQ=255)

                     gCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATg                                                                   >  1:1848707/1‑93 (MQ=255)

                     gCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATg                                                                   <  1:1494993/93‑1 (MQ=255)

                     gCACTGTATCTGATTGAACAACTGGGGAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTCTACTGATg                                                                   >  2:1122067/1‑93 (MQ=255)

                      cACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGc                                                                  <  1:2028503/93‑1 (MQ=255)

                       aCTGTATCTGATTGAACAACTGGGTAAACGTGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCg                                                                 >  1:91930/1‑93 (MQ=255)

                         tGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGtt                                                               <  2:2917007/93‑1 (MQ=255)

                         tGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGtt                                                               >  2:4589525/1‑93 (MQ=255)

                         tGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGtt                                                               >  1:302586/1‑93 (MQ=255)

                           tATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTcc                                                             <  1:4540209/93‑1 (MQ=255)

                           tATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTcc                                                             >  2:1242679/1‑93 (MQ=255)

                           tATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTcc                                                             >  1:1681399/1‑93 (MQ=255)

                           tATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTcc                                                             >  2:1517386/1‑93 (MQ=255)

                           tATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTcc                                                             <  1:2464323/93‑1 (MQ=255)

                           tATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTcc                                                             >  2:4375918/1‑93 (MQ=255)

                            aTCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCcg                                                            <  2:3342629/93‑1 (MQ=255)

                                 aTTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGaaa                                                       <  2:1713735/93‑1 (MQ=255)

                                  ttGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGaaaa                                                      >  1:1288701/1‑93 (MQ=255)

                                  ttGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGaaaa                                                      >  2:739945/1‑93 (MQ=255)

                                   tGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGATGAACCGTATACTGATGCGTTCCGCGaaaaa                                                     <  1:924232/93‑1 (MQ=255)

                                    gAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAg                                                    <  2:991530/93‑1 (MQ=255)

                                     aacaacTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGt                                                   <  1:166149/93‑1 (MQ=255)

                                       caacTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTAc                                                 <  2:2653351/93‑1 (MQ=255)

                                       caacTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTAc                                                 <  2:2357344/93‑1 (MQ=255)

                                       caacTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTAc                                                 <  1:3720524/93‑1 (MQ=255)

                                        aacTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTAcg                                                >  2:1892066/1‑93 (MQ=255)

                                         acTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTAcgc                                               <  1:2590354/93‑1 (MQ=255)

                                          cTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTAcgcg                                              >  1:1448081/1‑93 (MQ=255)

                                           tGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACgcgc                                             >  1:448075/1‑93 (MQ=255)

                                                 aaCGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGttt                                       >  1:748753/1‑93 (MQ=255)

                                                 aaCGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGttt                                       >  2:4673454/1‑93 (MQ=255)

                                                  aCGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTc                                      <  2:4303286/93‑1 (MQ=255)

                                                    gcgGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCa                                    >  2:442458/1‑93 (MQ=255)

                                                    gcgGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCa                                    <  1:2109798/93‑1 (MQ=255)

                                                     cgGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCAc                                   <  1:1706024/93‑1 (MQ=255)

                                                      ggCATTGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACg                                  <  2:375335/93‑1 (MQ=255)

                                                          ttGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTcc                              >  2:207905/1‑93 (MQ=255)

                                                          ttGCTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTcc                              <  2:1942723/93‑1 (MQ=255)

                                                             cTTATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGAt                           >  1:2760716/1‑93 (MQ=255)

                                                              ttATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGAtt                          >  2:440760/1‑93 (MQ=255)

                                                              ttATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGAtt                          >  1:1035728/1‑93 (MQ=255)

                                                               tATCTGCATATGTCAGAACCAGGTTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTa                         <  1:3311308/93‑1 (MQ=255)

                                                               tATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTa                         >  1:2588316/1‑93 (MQ=255)

                                                               tATCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTa                         >  1:2884315/1‑93 (MQ=255)

                                                                aTCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTAt                        <  2:2071948/93‑1 (MQ=255)

                                                                aTCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTAt                        <  2:517456/93‑1 (MQ=255)

                                                                 tCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATc                       <  1:3680234/93‑1 (MQ=255)

                                                                 tCTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATc                       <  2:2851856/93‑1 (MQ=255)

                                                                  cTGCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCg                      <  1:1313776/93‑1 (MQ=255)

                                                                    gCATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGgc                    >  2:4447422/1‑93 (MQ=255)

                                                                     cATATGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGgcg                   >  1:2355499/1‑93 (MQ=255)

                                                                      atatGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGCTTCCACGGTCCGATTATCGgcgc                  >  1:3964975/1‑93 (MQ=255)

                                                                       tatGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCa                 >  1:987145/1‑93 (MQ=255)

                                                                       tatGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCa                 >  1:298324/1‑93 (MQ=255)

                                                                        atGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAg                <  2:524865/93‑1 (MQ=255)

                                                                        atGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAg                >  2:3707091/1‑93 (MQ=255)

                                                                        atGTCAGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAg                <  2:1334722/93‑1 (MQ=255)

                                                                          gTCAGAACCAGATTGGGC‑GGGGGTGACCCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGt              <  1:4037301/93‑1 (MQ=255)

                                                                             aGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCa           <  2:2791968/93‑1 (MQ=255)

                                                                             aGAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCa           >  1:3932916/1‑93 (MQ=255)

                                                                              gAACCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCAt          <  1:642686/93‑1 (MQ=255)

                                                                               aaCCAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATa         <  1:2047572/93‑1 (MQ=255)

                                                                                 ccAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATaca       >  2:1484674/1‑93 (MQ=255)

                                                                                 ccAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATaca       >  2:3172384/1‑93 (MQ=255)

                                                                                 ccAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATaca       <  2:916268/93‑1 (MQ=255)

                                                                                  cAGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATacac      <  2:4167801/93‑1 (MQ=255)

                                                                                   aGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATAcaca     <  1:3456303/93‑1 (MQ=255)

                                                                                   aGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATAcaca     <  1:1788885/93‑1 (MQ=255)

                                                                                   aGATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATAcaca     <  2:4121371/93‑1 (MQ=255)

                                                                                    gATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATACACAg    >  2:1337325/1‑93 (MQ=255)

                                                                                    gATTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATACACAg    >  1:3699227/1‑93 (MQ=255)

                                                                                     aTTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAATTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATACACAGt   <  1:4457827/93‑1 (MQ=255)

                                                                                     aTTGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATACACAGt   <  2:1227256/93‑1 (MQ=255)

                                                                                      ttGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCGCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATACACAGTa  <  1:2515002/93‑1 (MQ=255)

                                                                                      ttGGGC‑GGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATACACAGTa  >  1:4464590/1‑93 (MQ=255)

GGCCCGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGGCGGGGGGTGAACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATACACAGTA  >  NC_000913/1725421‑1725600

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 22 ≤` `ATCG/ATCG` `< 29 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`