BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,642,279	A→G	V30A (GTC→GCC)	b1560 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,642,279	0	A	G	99.6%	727.4 / ‑5.3	226	V30A (GTC→GCC)	b1560	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (129/96); total (129/97)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 240 total aligned reads displayed.

GATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACG  >  NC_000913/1642189‑1642368

gATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGgcc                                                                                         <  1:4610327/93‑1 (MQ=255)

gATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGgcc                                                                                         <  1:5617606/93‑1 (MQ=255)

gATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGgcc                                                                                         >  1:6004309/1‑93 (MQ=255)

 aTGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGgcca                                                                                        >  1:1735430/1‑93 (MQ=255)

  tGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGgccac                                                                                       >  1:5373184/1‑93 (MQ=255)

   ggCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGGACGGCTTAAGGgccacc                                                                                      >  1:5919084/1‑93 (MQ=255)

     cAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAg                                                                                    <  2:3884475/93‑1 (MQ=255)

       gATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGcc                                                                                  >  2:3502112/1‑93 (MQ=255)

       gATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGcc                                                                                  >  2:4116689/1‑93 (MQ=255)

       gATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGACTATGAACTATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGcc                                                                                  >  2:2576142/1‑93 (MQ=255)

        aTTACGCATGCTTTTTGGTCCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGccc                                                                                 <  1:1200594/93‑1 (MQ=255)

          tACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCaa                                                                               <  1:297720/93‑1 (MQ=255)

          tACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCaa                                                                               >  2:5567606/1‑93 (MQ=255)

          tACGCATGCTTTTTGGTTCCGGATCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACCGCCCGGCTTAAGGGCCACCCGCCCaa                                                                               >  1:1070235/1‑93 (MQ=255)

            cGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGc                                                                             <  1:3072289/93‑1 (MQ=255)

            cGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCAGGCTTAAGGGCCACCAGCCCAAGc                                                                             <  1:6460110/93‑1 (MQ=255)

             gCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCt                                                                            >  1:2720963/1‑93 (MQ=255)

             gCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCt                                                                            <  1:373777/93‑1 (MQ=255)

               aTGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCa                                                                          <  2:1232156/93‑1 (MQ=255)

               aTGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCa                                                                          <  2:1974406/93‑1 (MQ=255)

               aTGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCa                                                                          <  1:4144263/93‑1 (MQ=255)

               aTGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCa                                                                          >  2:5263144/1‑93 (MQ=255)

               aTGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCa                                                                          <  1:3686659/93‑1 (MQ=255)

                tGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAg                                                                         <  2:5997182/93‑1 (MQ=255)

                   tttttGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCa                                                                      >  2:763161/1‑93 (MQ=255)

                   tttttGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCa                                                                      >  2:680702/1‑93 (MQ=255)

                    ttttGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCAc                                                                     >  1:3515907/1‑93 (MQ=255)

                    ttttGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCAc                                                                     <  2:4348497/93‑1 (MQ=255)

                     tttGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCAcc                                                                    >  1:3552790/1‑93 (MQ=255)

                      ttGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCa                                                                   <  2:4699146/93‑1 (MQ=255)

                          ttccggttccACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAAc                                                               <  1:6212725/93‑1 (MQ=255)

                             cggttccACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAgg                                                            <  1:1799006/93‑1 (MQ=255)

                              ggttccACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGt                                                           <  1:5445777/93‑1 (MQ=255)

                               gttccACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTc                                                          >  1:1492189/1‑93 (MQ=255)

                               gttccACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTc                                                          >  2:6078802/1‑93 (MQ=255)

                                  ccaccaGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTaa                                                       <  2:3079297/93‑1 (MQ=255)

                                  ccaccaGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTaa                                                       <  2:4672974/93‑1 (MQ=255)

                                   caccaGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAg                                                      <  2:330140/93‑1 (MQ=255)

                                   caccaGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAg                                                      <  1:4564329/93‑1 (MQ=255)

                                   caccaGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGACTAAg                                                      >  2:5655943/1‑93 (MQ=255)

                                    accaGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGt                                                     <  1:3375609/93‑1 (MQ=255)

                                    accaGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGt                                                     >  1:1925205/1‑93 (MQ=255)

                                     ccaGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTa                                                    >  1:2842402/1‑93 (MQ=255)

                                     ccaGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTa                                                    >  1:5663579/1‑93 (MQ=255)

                                      caGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTaa                                                   <  2:2076656/93‑1 (MQ=255)

                                       aGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAAt                                                  >  2:4528238/1‑93 (MQ=255)

                                        gCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTGAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATa                                                 >  2:2661548/1‑93 (MQ=255)

                                         cACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATAc                                                <  2:6270564/93‑1 (MQ=255)

                                          aCTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATAcc                                               >  2:5122064/1‑93 (MQ=255)

                                          aCTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATAcc                                               <  2:822053/93‑1 (MQ=255)

                                           cTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATAccc                                              <  2:3027015/93‑1 (MQ=255)

                                            tCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCg                                             >  2:2241749/1‑93 (MQ=255)

                                            tCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCg                                             >  1:4006297/1‑93 (MQ=255)

                                               aGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCAc                                          >  1:6317337/1‑93 (MQ=255)

                                               aGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCAc                                          <  2:19098/93‑1 (MQ=255)

                                               aGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCAc                                          <  2:1461647/93‑1 (MQ=255)

                                                ggATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACg                                         >  2:5938887/1‑93 (MQ=255)

                                                 gATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGt                                        >  2:2611549/1‑93 (MQ=255)

                                                   ttATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTa                                      >  2:3543782/1‑93 (MQ=255)

                                                   ttATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTa                                      >  2:4973809/1‑93 (MQ=255)

                                                   ttATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTa                                      >  1:4400655/1‑93 (MQ=255)

                                                      tGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCt                                   >  2:313228/1‑93 (MQ=255)

                                                      tGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCt                                   <  1:287836/93‑1 (MQ=255)

                                                       gAAATATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTc                                  >  2:1571492/1‑93 (MQ=255)

                                                         aaTATCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCa                                >  1:2380102/1‑93 (MQ=255)

                                                           tatCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGa                              <  1:5706165/93‑1 (MQ=255)

                                                           tatCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGa                              >  2:4497390/1‑93 (MQ=255)

                                                            atCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGAt                             >  2:5035561/1‑93 (MQ=255)

                                                            atCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGAt                             >  2:6678441/1‑93 (MQ=255)

                                                             tCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATc                            <  2:140800/93‑1 (MQ=255)

                                                             tCTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATc                            >  1:7227/1‑93 (MQ=255)

                                                              cTGTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATcc                           >  1:1455105/1‑93 (MQ=255)

                                                                gTATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGt                         <  1:235885/93‑1 (MQ=255)

                                                                 tATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGtt                        >  1:121339/1‑93 (MQ=255)

                                                                 tATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGtt                        <  2:4023216/93‑1 (MQ=255)

                                                                 tATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGtt                        >  2:1361966/1‑93 (MQ=255)

                                                                 tATGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGtt                        <  1:6632008/93‑1 (MQ=255)

                                                                  aTGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTg                       <  1:5147633/93‑1 (MQ=255)

                                                                   tGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGc                      >  2:5421831/1‑93 (MQ=255)

                                                                   tGGATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGc                      >  2:1654415/1‑93 (MQ=255)

                                                                     gATTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTg                    >  2:1861618/1‑93 (MQ=255)

                                                                      aTTCACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTgg                   <  2:6618159/93‑1 (MQ=255)

                                                                         cACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGaaa                >  2:1795778/1‑93 (MQ=255)

                                                                         cACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGaaa                >  1:3853692/1‑93 (MQ=255)

                                                                         cACGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGaaa                >  2:1318997/1‑93 (MQ=255)

                                                                           cGGCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGt              >  1:4104020/1‑93 (MQ=255)

                                                                           aggCCCGGCTTAAGGGCCACCAGCCCAAGATCAGGCACCAGAACAGGTCTAAGTAATAACCGCACGTTACCTCCAGATCCGTTGCTGGAAAgg              >  2:1703575/2‑92 (MQ=255)

                                                                            ggCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTg             <  2:3848462/93‑1 (MQ=255)

                                                                            ggCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTg             <  1:6604121/93‑1 (MQ=255)

                                                                            ggCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTg             >  2:2800983/1‑93 (MQ=255)

                                                                             gCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGc            <  2:2704624/93‑1 (MQ=255)

                                                                             gCCCGGCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGc            >  2:3269533/1‑93 (MQ=255)

                                                                                  gCTTAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGAc       >  2:6218792/1‑93 (MQ=255)

                                                                                  gCTTAAGGACCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGAc       <  2:5288411/93‑1 (MQ=255)

                                                                                    ttAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGc     >  1:844990/1‑93 (MQ=255)

                                                                                    ttAAGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGc     >  1:4835587/1‑93 (MQ=255)

                                                                                       aGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACg  >  1:3979129/1‑93 (MQ=255)

                                                                                       aGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACg  >  1:6080253/1‑93 (MQ=255)

                                                                                       aGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACg  >  2:2221380/1‑93 (MQ=255)

                                                                                       aGGGCCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACg  >  2:3410056/1‑93 (MQ=255)

GATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACG  >  NC_000913/1642189‑1642368

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