BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	541,457	(G)_5→6	coding (346/1146 nt)	ybbZ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	541,452	1	.	G	98.7%	705.4 / 1.3	227	coding (341/1146 nt)	ybbZ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (1/2); new base G (100/124); total (101/126)
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.99e-01

Only 100 of 268 total aligned reads displayed.

GAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTAC‑GGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCGACGGCGGTAT  >  NC_000913/541360‑541539

gAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTAc                                                                                          >  2:7550625/1‑93 (MQ=255)

 aaaaCGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTAC‑g                                                                                        <  2:519021/93‑1 (MQ=255)

   aaCGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTAC‑ggg                                                                                      >  2:1595253/1‑93 (MQ=255)

   aaCGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTAC‑ggg                                                                                      <  1:7573073/93‑1 (MQ=255)

   aaCGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTAC‑ggg                                                                                      <  2:5254206/93‑1 (MQ=255)

   aaCGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTAC‑ggg                                                                                      >  1:4520/1‑93 (MQ=255)

     cggctgTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTAC‑ggggg                                                                                    <  2:320913/93‑1 (MQ=255)

       gcggTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTAC‑gggggga                                                                                  >  1:3622359/1‑91 (MQ=255)

       gcggTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTAC‑gggggga                                                                                  <  2:3291047/93‑3 (MQ=255)

        cggTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTAC‑ggggggag                                                                                 >  2:7197987/1‑90 (MQ=255)

         ggTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTAC‑ggggggagt                                                                                >  1:7347003/1‑89 (MQ=255)

         ggTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTAC‑ggggggagt                                                                                >  2:5662318/1‑89 (MQ=255)

            gATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTaa                                                                              >  2:5014574/1‑93 (MQ=255)

            gATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTaa                                                                              >  2:1348575/1‑93 (MQ=255)

                  gaTGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTc                                                                        <  1:4315351/93‑1 (MQ=255)

                  gaTGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTc                                                                        <  1:55217/93‑1 (MQ=255)

                  gaTGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTc                                                                        <  2:6349492/93‑1 (MQ=255)

                   aTGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCa                                                                       <  1:7127197/93‑1 (MQ=255)

                   aTGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCa                                                                       <  1:6948418/93‑1 (MQ=255)

                   aTGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCa                                                                       <  1:1168222/93‑1 (MQ=255)

                    tGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAGAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCAt                                                                      >  1:7296429/1‑93 (MQ=255)

                    tGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCAt                                                                      >  1:3308842/1‑93 (MQ=255)

                      gcggcAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATgc                                                                    <  1:584664/93‑1 (MQ=255)

                        ggcAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATgcgc                                                                  >  1:222598/1‑93 (MQ=255)

                            gcaAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGa                                                              <  1:5161483/93‑1 (MQ=255)

                            gcaAGTGGCCTGATGATGGTCGCTCCTGAAAAGCGTAATCCGTGGCTAGCCTCCAGTTTTGGTACGGGGGGCGTTAATTCGTCAATCGGTGGa                                                              <  2:7011038/93‑1 (MQ=255)

                              aaGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATa                                                            >  1:1010512/1‑93 (MQ=255)

                               aGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATaa                                                           >  2:6283352/1‑93 (MQ=255)

                               aGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATaa                                                           >  2:6919933/1‑93 (MQ=255)

                                gTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAAc                                                          <  1:2252610/93‑1 (MQ=255)

                                 tGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACg                                                         <  2:1918967/93‑1 (MQ=255)

                                 tGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACg                                                         <  2:7229384/93‑1 (MQ=255)

                                  ggCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGa                                                        <  1:3991663/93‑1 (MQ=255)

                                   gCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGAc                                                       >  1:2351117/1‑93 (MQ=255)

                                    gcTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACa                                                      <  1:3837201/92‑1 (MQ=255)

                                    ccTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACa                                                      <  1:1410320/93‑1 (MQ=255)

                                     cTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACAt                                                     <  1:5578572/93‑1 (MQ=255)

                                     cTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACAt                                                     >  1:5048242/1‑93 (MQ=255)

                                      tGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACAtt                                                    >  1:6416113/1‑93 (MQ=255)

                                      tGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACAtt                                                    >  1:831948/1‑93 (MQ=255)

                                      tGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACAtt                                                    <  2:389139/93‑1 (MQ=255)

                                      tGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACAtt                                                    <  2:5566097/93‑1 (MQ=255)

                                      tGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACAtt                                                    >  1:7537518/1‑93 (MQ=255)

                                       gATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTc                                                   <  2:7280780/93‑1 (MQ=255)

                                       gATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTc                                                   <  2:2035061/93‑1 (MQ=255)

                                         tGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGGGTTAATTCGTCATGCGCTGGATAACGACATTCGc                                                 >  2:603295/1‑93 (MQ=255)

                                          gCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGcc                                                >  1:2826255/1‑93 (MQ=255)

                                            tGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCAGGCGCTGGAGAACGACATGCGCCaa                                              >  1:6013127/1‑92 (MQ=255)

                                            tGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCGCTAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCat                                              >  1:1375218/1‑93 (MQ=255)

                                             ggTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCata                                             <  2:4269774/93‑1 (MQ=255)

                                             ggTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCata                                             <  2:5935620/93‑1 (MQ=255)

                                              gTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCGAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCatat                                            <  1:5731345/93‑1 (MQ=255)

                                              gTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCatat                                            <  2:4131829/93‑1 (MQ=255)

                                              gTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCatat                                            >  1:3392329/1‑93 (MQ=255)

                                               tCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCAtatt                                           >  2:5345089/1‑93 (MQ=255)

                                                cgcgCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCAtatta                                          <  2:2612874/93‑1 (MQ=255)

                                                   gcCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATtc                                       <  2:838953/93‑1 (MQ=255)

                                                   gcCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATtc                                       <  1:6621083/93‑1 (MQ=255)

                                                     cTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGTGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATtctc                                     >  1:3801436/1‑93 (MQ=255)

                                                     cTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATtctc                                     <  1:206079/93‑1 (MQ=255)

                                                     cTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATtctc                                     >  1:479782/1‑93 (MQ=255)

                                                      tGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCg                                    >  2:2640479/1‑93 (MQ=255)

                                                      tGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCg                                    >  2:5827662/1‑93 (MQ=255)

                                                       gAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCgg                                   <  2:3076382/93‑1 (MQ=255)

                                                         aaaGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCa                                 >  2:1608542/1‑93 (MQ=255)

                                                         aaaGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCa                                 >  1:4728437/1‑93 (MQ=255)

                                                          aaGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCAt                                <  2:1891734/93‑1 (MQ=255)

                                                           aGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCAtt                               <  2:2498195/93‑1 (MQ=255)

                                                              gTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTggc                            <  2:1635829/93‑1 (MQ=255)

                                                              gTAATCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTggc                            >  1:214245/1‑93 (MQ=255)

                                                                  tCCGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCa                        <  2:6430753/93‑1 (MQ=255)

                                                                   ccGTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAg                       >  1:7477821/1‑93 (MQ=255)

                                                                     gTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTg                     >  2:4427018/1‑93 (MQ=255)

                                                                     gTTGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTg                     >  1:529482/1‑93 (MQ=255)

                                                                      ttGCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGc                    <  2:1062883/93‑1 (MQ=255)

                                                                        gCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGa                  >  2:3593995/1‑93 (MQ=255)

                                                                        gCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGa                  <  2:2875806/93‑1 (MQ=255)

                                                                        gCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGa                  >  1:4790845/1‑93 (MQ=255)

                                                                        gCTGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGGGCGa                  >  1:971088/1‑93 (MQ=255)

                                                                          tGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACg                <  2:2269566/93‑1 (MQ=255)

                                                                          tGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACg                >  1:25600/1‑93 (MQ=255)

                                                                          tGGCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATAGTATTCTCGGCATTGGCGGCGGTGGGACg                >  2:672978/1‑93 (MQ=255)

                                                                           ggCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACgg               <  1:5187937/93‑1 (MQ=255)

                                                                           ggCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACgg               <  1:7767225/93‑1 (MQ=255)

                                                                            gCCTCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACggg              >  2:2971914/1‑92 (MQ=255)

                                                                               tCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCgc           >  1:1695407/1‑92 (MQ=255)

                                                                               tCCAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCGa           <  2:2561278/93‑1 (MQ=255)

                                                                                ccAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCGAc          >  2:3492596/1‑93 (MQ=255)

                                                                                ccAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCGAc          >  2:4548311/1‑93 (MQ=255)

                                                                                 cAGTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCGAcg         >  2:7838374/1‑93 (MQ=255)

                                                                                   gTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCGAcggg       >  1:682065/1‑92 (MQ=255)

                                                                                   gTTTTGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCGAcggc       >  1:1978667/1‑93 (MQ=255)

                                                                                      ttGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCGACGGCGGt    >  2:7385838/1‑93 (MQ=255)

                                                                                       tGGTACGGGGGGGGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCGACGGCGGTa   <  1:6877812/93‑1 (MQ=255)

                                                                                       tGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCGACGGCGGTa   <  1:4597367/93‑1 (MQ=255)

                                                                                       tGGTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCGACGGCGGTa   <  1:5854615/93‑1 (MQ=255)

                                                                                        ggTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCGACGGCGGTAt  <  1:2984330/93‑1 (MQ=255)

                                                                                        ggTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCGACGGCGGTAt  <  1:2762326/93‑1 (MQ=255)

                                                                                        ggTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCGACGGCGGTAt  >  2:2337348/1‑93 (MQ=255)

                                                                                        ggTACGGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCGACGGCGGTAt  >  2:7902213/1‑93 (MQ=255)

GAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGTTTTGGTAC‑GGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCGACGGCGGTAT  >  NC_000913/541360‑541539

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