BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,056,041	A→G	E186E (GAA→GAG)	torT →	TorT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,056,041	0	A	G	99.3%	359.2 / NA	136	E186E (GAA→GAG)	torT	TorT
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (76/59); minor base T (0/1); total (76/60)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 136 total aligned reads displayed.

GCACGGTAAACCACTGAATGTGCTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAAGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTATTGTTGATATTGCGCTTGGTGATAACGATATTGAA  >  NC_000913/1055969‑1056113

gcACGGTAAACCACTGACTGTGCTGTTGATGCCCGGACCCGATAACGCCGCGGGCAGTAAGGATATGGTCGATGGt                                                                       <  1:4604401/76‑1 (MQ=255)

gcACGGTAAACCACTGAATGTGCTGTTGATGCCCGGACCCGCTAACGCCGGGGGCAGTCAGGGGGTGGTCGAGGGt                                                                       >  2:5670634/1‑76 (MQ=255)

gcACGGTAAACCACTGAATGTGCTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGt                                                                       <  2:2450787/76‑1 (MQ=255)

 cACGGTAAACCACTGAATGTGCTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGtt                                                                      >  1:2611148/1‑76 (MQ=255)

  aCGGTAAACCACTGAATGTGCTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGctt                                                                     >  1:1506364/1‑76 (MQ=255)

   cGGTAAACCACTGAATGTGCTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGtttt                                                                    <  1:2211804/76‑1 (MQ=255)

    ggTAAACCACTGAATGTGCTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTc                                                                   >  1:3903289/1‑76 (MQ=255)

    ggTAAACCACTGAATGTGCTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGGTTTc                                                                   >  1:3845988/1‑76 (MQ=255)

     gTAAACCACTGAATGTGCTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTcg                                                                  >  1:4020489/1‑76 (MQ=255)

     gTAAACCACTGAATGTGCTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTcg                                                                  <  2:5015500/76‑1 (MQ=255)

     gTAAACCACTGAATGTGCTGTTGAGGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTcg                                                                  <  2:5274409/76‑1 (MQ=255)

      tagaCCACGGAATGTGCTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTcgc                                                                 <  2:5550195/73‑1 (MQ=255)

      tAAACCACTGAATGTGCTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTcgc                                                                 >  1:151606/1‑76 (MQ=255)

      tAAACCACTGAATGTGCTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTcgc                                                                 >  1:5601423/1‑76 (MQ=255)

       aaaCCACTGAATGTGCTGTTGATGCCCGGACCCGACAACTCCGGGGGCCGTCAGGAGCTGGCCGAGGGCTATcgcg                                                                >  2:641165/1‑76 (MQ=255)

        aaCCACTGAATGTGCTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCgcgc                                                               >  2:4713725/1‑76 (MQ=255)

         aCCACTGAATGTGCTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCa                                                              <  2:1569767/76‑1 (MQ=255)

          ccACTGAATGTGCTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAg                                                             <  2:4544870/76‑1 (MQ=255)

          ccACTGAATGTGCTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAg                                                             >  1:3646012/1‑76 (MQ=255)

          ccACTGAATGTGCTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAg                                                             <  1:3017339/76‑1 (MQ=255)

           cACTGAATGTGCTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGGTTTCGCGCAGc                                                            >  1:352387/1‑76 (MQ=255)

            aCTGAATGTGCTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGcc                                                           >  2:4659861/1‑76 (MQ=255)

             cTGAATGTGCTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCa                                                          <  1:3464609/76‑1 (MQ=255)

              tGAATGTGCTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCCGCCAt                                                         >  2:1347547/1‑76 (MQ=255)

              tGAATGTGCTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCAt                                                         >  2:2614292/1‑76 (MQ=255)

              tGAATGTGCTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCAt                                                         >  1:2822365/1‑76 (MQ=255)

               gAATGTGCTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCAtt                                                        <  1:5821471/76‑1 (MQ=255)

                aaTGTGCTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTg                                                       >  1:1808416/1‑76 (MQ=255)

                aaTGTGCTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTg                                                       >  1:5306414/1‑76 (MQ=255)

                aaTGTGCTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTg                                                       <  1:1174820/76‑1 (MQ=255)

                aaTGTGCTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTg                                                       >  2:4605772/1‑76 (MQ=255)

                 aTGTGCTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGCGGTGGTCGAGGGTCTTCGCGCAGCCATTGc                                                      >  2:4902718/1‑76 (MQ=255)

                 aTGTGCTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGc                                                      >  1:1024212/1‑76 (MQ=255)

                  tgtgCTGTTGATGCCCGGACCCGATAACGCGGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGcc                                                     <  2:5004137/76‑1 (MQ=255)

                  tgtgCTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGcc                                                     <  1:5075852/76‑1 (MQ=255)

                   gtgCTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGGCGAGGGTTTTCGCGCAGCCATTGCCg                                                    >  2:5353206/1‑76 (MQ=255)

                     gCTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGCTGGTCGAGGGATTTCGCGCAGCCATTGCCGGa                                                  >  1:5685314/1‑76 (MQ=255)

                     gCTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGa                                                  <  1:2713057/76‑1 (MQ=255)

                      cTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGaa                                                 <  1:4713726/76‑1 (MQ=255)

                      cTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGaa                                                 <  1:2825641/76‑1 (MQ=255)

                          tGATGCGCGGACCCGATAACGCCGGGGGCAGTAAGGATATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGccc                                             <  2:1806953/76‑1 (MQ=255)

                          tGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGccc                                             <  2:114919/76‑1 (MQ=255)

                          tGAGGCCCGGGCGCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGccc                                             <  2:557892/76‑1 (MQ=255)

                            aTGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCgg                                           <  1:1359707/76‑1 (MQ=255)

                            aTGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCgg                                           >  2:2306636/1‑76 (MQ=255)

                             tGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCggg                                          >  1:2175729/1‑75 (MQ=255)

                              gCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGGAGCCCGGTg                                         >  1:1075088/1‑76 (MQ=255)

                              gCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTg                                         >  2:2629735/1‑76 (MQ=255)

                              gCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTg                                         <  2:4511384/76‑1 (MQ=255)

                              gCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTg                                         >  1:1187792/1‑76 (MQ=255)

                              gCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGGTTTCGCGCAGACATTGCCGGCAGCCCGGTg                                         >  1:480511/1‑76 (MQ=255)

                                cgggACCCGAGAGCGCGGGGGGCAGTAAGGAGAGGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCg                                       <  2:276360/74‑1 (MQ=255)

                                ccGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCg                                       <  1:2837659/76‑1 (MQ=255)

                                ccGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCg                                       <  1:4285706/76‑1 (MQ=255)

                                ccGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCg                                       >  1:411262/1‑76 (MQ=255)

                                ccGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCg                                       >  1:3470345/1‑76 (MQ=255)

                                 cGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGt                                      >  2:4975816/1‑76 (MQ=255)

                                  tgACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGGGCGTa                                     >  2:4206221/2‑76 (MQ=255)

                                  ggACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTa                                     >  1:2613358/1‑76 (MQ=255)

                                   gACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTAt                                    >  1:4596089/1‑76 (MQ=255)

                                    aCCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTAtt                                   <  1:1917858/76‑1 (MQ=255)

                                     cccGATAACGCCGGGGGCAGTAGGGAGTTGGTCGAGGGATTTCGCGCAGCCGTTGCCGGAAGCACGGCGCGTAttg                                  >  1:105091/1‑76 (MQ=255)

                                     cccGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTAttg                                  >  1:3078396/1‑76 (MQ=255)

                                       cGATAGCGGCGGGGGGAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTAttgtt                                <  2:43900/76‑1 (MQ=255)

                                       cGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTAttgtt                                >  1:4398225/1‑76 (MQ=255)

                                       cGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTAttgtt                                >  1:3527741/1‑76 (MQ=255)

                                       cGATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTAttgtt                                >  1:4743171/1‑76 (MQ=255)

                                        gATAACGCCGGGGGCAGTACGGCGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTAttgttg                               >  2:1053425/1‑76 (MQ=255)

                                        gATAACGCCGGGGGCAGTAAGGCGATGGTCGAGGGGTTTCGCGCAGCCATTGCCGGAAGCCCGGAGCGTAttgttg                               >  2:4879209/1‑76 (MQ=255)

                                        gATAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTAttgttg                               >  2:4538819/1‑76 (MQ=255)

                                         aTAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTATTGTTGa                              <  1:3059701/76‑1 (MQ=255)

                                         aTAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTATTGTTGa                              <  2:3045223/76‑1 (MQ=255)

                                          tAACGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTATTGTTGat                             >  1:326803/1‑76 (MQ=255)

                                           aaCGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTATTGTTGata                            <  1:855864/76‑1 (MQ=255)

                                           aaCGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTATTGTTGata                            <  1:2294976/76‑1 (MQ=255)

                                            aCGCGGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTATTGTTGatat                           <  2:2463356/76‑1 (MQ=255)

                                            aCGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTATTGTTGatat                           <  2:542706/76‑1 (MQ=255)

                                             cGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTATTGTTGATAtt                          <  2:2069857/76‑1 (MQ=255)

                                             cGCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATGGCCGGAAGCCCGGTGCGTATTGTTGATAtt                          >  2:4669423/1‑76 (MQ=255)

                                              gCCGGGGGCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTATTGTTGATATTg                         >  1:2568794/1‑76 (MQ=255)

                                                 gggggCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTATTGTTGATATTgcgc                      >  1:1194602/1‑76 (MQ=255)

                                                   gggCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTATTGTTGATATTGCGCtt                    <  1:1803206/76‑1 (MQ=255)

                                                     gCAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTATTGTTGATATTGCGCTTgg                  >  2:2215805/1‑76 (MQ=255)

                                                      cAGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGCGCGTATTGTTGATATTGCGCTTggg                 >  2:4774740/1‑75 (MQ=255)

                                                       aGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTATTGTTGATATTGCGCTTGGtn                <  2:2998204/76‑2 (MQ=255)

                                                       aGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTATTGTTGATATTGCGCTTGGTg                <  2:5074294/76‑1 (MQ=255)

                                                       aGTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTATTGTTGATATTGCGCTTGGTg                >  2:3677188/1‑76 (MQ=255)

                                                        gTAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTATTGTTGATATTGCGCTTGGTGa               >  1:3314244/1‑76 (MQ=255)

                                                         tAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTATTGTTGATATTGCGCTTGGTGAt              >  2:2369833/1‑76 (MQ=255)

                                                         tAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTATTGTTGATATTGCGCTTGGTGAt              <  2:2602975/76‑1 (MQ=255)

                                                         tAAGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTATTGTTGATATTGCGCTTGCTGAt              >  1:5424650/1‑76 (MQ=255)

                                                          aaGGAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTATTGTTGATATTGCGCTTGGTGATa             >  2:2485374/1‑76 (MQ=255)

                                                          aaGGAGATGGTCGAGAGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTATTGTTGATATTGCGCTTGGTGATa             >  1:2106608/1‑76 (MQ=255)

                                                           aGGAGATGGTCGAGGGTTTACGCGCAGCCATTGCCGGAAGCCCGGTGCGTATTGTTGATATTGCTGTTGGTGATaa            <  1:5258702/76‑1 (MQ=255)

                                                            ggAGATGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGCGCGTCTTGTTGCTATTGCGCCTGCTGATCAc           >  2:5488109/1‑76 (MQ=255)

                                                             gagaTGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTATTGTTGATATTGCGCTTGGTGATAACg          <  2:1280308/76‑1 (MQ=255)

                                                               gaTGGTCGAGGTGTTTCGCGCAGCCATTGCCGGACGCCCGGTGCGTATTGTTGATATTGCGCTTGGTGATAACGat        <  2:3619641/76‑1 (MQ=255)

                                                               gaTGGTCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTATTGTTGATATTGCGCTTGGTGATAACGat        >  1:3081193/1‑76 (MQ=255)

                                                                    tCGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTATTGTTGATATTGCGCTTGGTGATAACGATATTGa   <  2:2200159/76‑1 (MQ=255)

                                                                     cGAGGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTATTGTTGATATTGCGCTTGGTGATAACGATATTGaa  <  2:411262/76‑1 (MQ=255)

GCACGGTAAACCACTGAATGTGCTGTTGATGCCCGGACCCGATAACGCCGGGGGCAGTAAGGAGATGGTCGAAGGTTTTCGCGCAGCCATTGCCGGAAGCCCGGTGCGTATTGTTGATATTGCGCTTGGTGATAACGATATTGAA  >  NC_000913/1055969‑1056113

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