BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,554,088	+C	coding (2/624 nt)	yhgKJ ←	RNA phosphage cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,554,088	1	.	C	99.5%	600.2 / ‑5.8	191	S133W (TCG→TGG)	yhgKJ	RNA phosphage cyclase
Reads supporting (aligned to +/- strand): ref base . (0/0); major base C (89/101); minor base A (0/1); total (89/102)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 223 total aligned reads displayed.

CGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCG‑AGGCGGCGCCGACGGGTTATCGGTCGCCG‑CGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTT  >  NC_000913/3553996‑3554174

cGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCg                                                                                          <  2:4280224/93‑1 (MQ=255)

  gCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgca                                                                                        <  1:3239027/93‑3 (MQ=255)

  gCAGAGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgca                                                                                        <  2:233609/93‑3 (MQ=255)

   cAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgcag                                                                                       >  1:1159530/1‑90 (MQ=255)

   cAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgcag                                                                                       <  2:2054435/93‑4 (MQ=255)

      ggTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAggcg                                                                                    <  1:5714153/93‑1 (MQ=255)

      ggTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAggcg                                                                                    >  1:1875681/1‑93 (MQ=255)

        tAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAggcggc                                                                                  <  1:3789498/93‑1 (MQ=255)

        tAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAggcggc                                                                                  <  2:1483341/93‑1 (MQ=255)

        tAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAggcggc                                                                                  <  1:3128961/93‑1 (MQ=255)

        tAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAggcggc                                                                                  >  2:1991612/1‑93 (MQ=255)

              ccGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAc                                                                            <  1:277085/93‑1 (MQ=255)

               cGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACg                                                                           <  1:4889006/93‑1 (MQ=255)

               cGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACg                                                                           <  1:1940719/93‑1 (MQ=255)

                   gCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGtt                                                                       >  1:4400357/1‑93 (MQ=255)

                    cGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTa                                                                      <  1:2508709/93‑1 (MQ=255)

                    cGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTa                                                                      <  2:6117520/93‑1 (MQ=255)

                     gTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTAt                                                                     <  1:2842776/93‑1 (MQ=255)

                     gTAACAGCGTGGCTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTAt                                                                     <  2:5829678/93‑1 (MQ=255)

                      tAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATc                                                                    <  2:168690/93‑1 (MQ=255)

                      tAACAGCGTGGTTAGCTGATGAGTTCCTGTATTCGAACGAAGCGGCTCAAGCACCCGGCAGATAACATCAGAAAAAGGCGCCGACGGGTTATc                                                                    <  2:6307889/93‑1 (MQ=255)

                        aCAGGGTTGTTTGCTGATGAAATCCTATTTGCGGCAGCAGCCGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCgg                                                                  <  1:6155185/93‑1 (MQ=255)

                        aCAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCgg                                                                  <  1:2294436/93‑1 (MQ=255)

                        aCAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCgg                                                                  <  1:3625081/93‑1 (MQ=255)

                        aCAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATACGCAGGCGGCGCCGACGGGGTATCgg                                                                  >  2:3248488/1‑93 (MQ=255)

                         cAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt                                                                 <  1:2855758/93‑1 (MQ=255)

                          aGCGTGGTTTGCTGATGAATTCCTAGGTTCGCGAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtg                                                                <  2:2601717/93‑2 (MQ=255)

                           gCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgc                                                               <  1:3365444/93‑3 (MQ=255)

                           gCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgc                                                               >  2:3418629/1‑91 (MQ=255)

                             gTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgccg                                                             >  1:3490201/1‑89 (MQ=255)

                              tGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCCCCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑c                                                          <  1:1806089/93‑1 (MQ=255)

                               ggTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑cc                                                         >  1:4960523/1‑92 (MQ=255)

                                 tttGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑ccgc                                                       >  2:6240746/1‑90 (MQ=255)

                                    gCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCa                                                     >  1:5543907/1‑93 (MQ=255)

                                    gCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCa                                                     <  1:5841331/93‑1 (MQ=255)

                                     cTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCAc                                                    >  2:4184812/1‑93 (MQ=255)

                                     cTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCAc                                                    >  1:5400190/1‑93 (MQ=255)

                                      tgatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACt                                                   >  2:726003/1‑93 (MQ=255)

                                      tgatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACt                                                   >  1:5732748/1‑93 (MQ=255)

                                       gatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACtt                                                  >  2:4695166/1‑93 (MQ=255)

                                        atgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTc                                                 <  1:3144140/93‑1 (MQ=255)

                                          gaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCaa                                               >  2:1313805/1‑93 (MQ=255)

                                            aTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGGCGC‑CGCTCACTTCAaca                                             >  2:3726791/1‑93 (MQ=255)

                                            aTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAaca                                             <  1:3769593/93‑1 (MQ=255)

                                            aTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAaca                                             <  2:2300845/93‑1 (MQ=255)

                                            aTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGCGTTATCGGT‑GCCGCCGCTCACTTCAaca                                             >  1:2436082/1‑93 (MQ=255)

                                             ttCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAacac                                            <  2:1635375/93‑1 (MQ=255)

                                             ttCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAacac                                            <  2:3133562/93‑1 (MQ=255)

                                              tCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACAcg                                           >  2:2231988/1‑93 (MQ=255)

                                              tCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACAcg                                           <  2:395676/93‑1 (MQ=255)

                                                cTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACAcgcg                                         >  2:3310677/1‑93 (MQ=255)

                                                 tATTTTCGCCAGCAGCGGCTCCGGCACCCGGCGGATAAAATTCGCAGGCGGCCCCGAAGGGTTCTCGGT‑GCCGCCGCTCACTTGAGCAcgcta                                        >  2:4591810/1‑91 (MQ=255)

                                                  aTTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTCCCGC‑CGCTCACTTCAACACGCGaa                                       <  1:2080074/93‑1 (MQ=255)

                                                  aTTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGaa                                       <  2:1123581/93‑1 (MQ=255)

                                                  aTTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGaa                                       <  1:5457611/93‑1 (MQ=255)

                                                   ttttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAg                                      <  2:4376257/93‑1 (MQ=255)

                                                     ttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGt                                    >  1:1304303/1‑93 (MQ=255)

                                                      tCGCCAGCAGGGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCCCTTCAACACGCGAAGGTc                                   <  2:1410112/93‑1 (MQ=255)

                                                        gCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCa                                 >  2:5908859/1‑93 (MQ=255)

                                                        gCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCa                                 <  1:938283/93‑1 (MQ=255)

                                                         ccAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCAt                                >  2:1959279/1‑93 (MQ=255)

                                                         ccAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCAt                                <  2:5457744/93‑1 (MQ=255)

                                                          cagcagCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATc                               >  2:3570187/1‑93 (MQ=255)

                                                          cagcagCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATc                               >  2:872584/1‑93 (MQ=255)

                                                          cagcagCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATc                               >  2:5040354/1‑93 (MQ=255)

                                                           agcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGACGC‑CGCTCACTTCAACACGCGAAGGTCCATCg                              >  2:182925/1‑93 (MQ=255)

                                                           agcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGG‑CGCCGCCGCTCACTTCAACACGCGAACGACCATCg                              >  1:5650423/1‑93 (MQ=255)

                                                            gcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCgg                             >  2:5406360/1‑93 (MQ=255)

                                                              agcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCa                           >  1:5190337/1‑93 (MQ=255)

                                                              agcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGAGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCa                           >  1:2414542/1‑93 (MQ=255)

                                                                 ggCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAc                        <  2:5669513/93‑1 (MQ=255)

                                                                  gCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAcc                       >  1:6052169/1‑93 (MQ=255)

                                                                  gCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAcc                       <  2:20930/93‑1 (MQ=255)

                                                                     ccAGCACCCGGCGGATATAACCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACcaca                    <  2:5554501/93‑1 (MQ=255)

                                                                      cAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAg                   <  1:105556/93‑1 (MQ=255)

                                                                      cAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAg                   <  1:5316136/93‑1 (MQ=255)

                                                                      cAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAg                   >  2:304840/1‑93 (MQ=255)

                                                                      cAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAg                   <  1:4540508/93‑1 (MQ=255)

                                                                      cAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAg                   <  2:2543013/93‑1 (MQ=255)

                                                                       aGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgc                  <  2:4622125/93‑1 (MQ=255)

                                                                       aGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCAACGGCAAACCACAgc                  <  2:5208659/93‑1 (MQ=255)

                                                                         cACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgcgc                <  1:5067627/93‑1 (MQ=255)

                                                                          aCCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGcgcg               >  1:2361986/1‑93 (MQ=255)

                                                                          aCCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGcgcg               >  1:4379970/1‑93 (MQ=255)

                                                                           cccGCCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCgg              <  2:2931073/93‑1 (MQ=255)

                                                                            ccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCggg             <  2:5631872/93‑1 (MQ=255)

                                                                            ccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCggg             >  2:6270485/1‑93 (MQ=255)

                                                                             cggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgc            >  1:775224/1‑93 (MQ=255)

                                                                              ggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgca           >  2:2539358/1‑93 (MQ=255)

                                                                              ggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgca           >  2:763700/1‑93 (MQ=255)

                                                                               gcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcag          <  1:722358/93‑1 (MQ=255)

                                                                                 ggATAGAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcagcc        >  2:5794477/1‑92 (MQ=255)

                                                                                 ggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcagca        <  1:1642989/93‑1 (MQ=255)

                                                                                 ggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcagca        >  1:2973002/1‑93 (MQ=255)

                                                                                 ggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGCCCATCCGCAAACCACAGCGCGGgcagca        >  2:5152509/1‑93 (MQ=255)

                                                                                  gATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAc       >  1:4121911/1‑93 (MQ=255)

                                                                                     aaaaTCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGt    >  2:1069399/1‑93 (MQ=255)

                                                                                     aaaaTCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGt    >  1:1418499/1‑93 (MQ=255)

                                                                                      aaaTCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGtt   <  1:3889667/93‑1 (MQ=255)

                                                                                       aaTCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGttt  <  2:3898301/93‑1 (MQ=255)

CGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCG‑AGGCGGCGCCGACGGGTTATCGGTCGCCG‑CGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTT  >  NC_000913/3553996‑3554174

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