BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,958,357	G→A	S698S (AGC→AGT)	recC ←	RecC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,958,357	0	G	A	96.7%	602.2 / 3.9	179	S698S (AGC→AGT)	recC	RecC
Reads supporting (aligned to +/- strand): ref base G (3/3); new base A (104/69); total (107/72)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.37e-01

Only 100 of 182 total aligned reads displayed.

ATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGA  >  NC_000913/2958267‑2958446

ataGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTTGGTTTCTGACt                                                                                         >  2:21121/1‑93 (MQ=255)

ataGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACt                                                                                         >  1:1194255/1‑93 (MQ=255)

   gagTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACtcat                                                                                      <  1:2170466/93‑1 (MQ=255)

   gagTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACtcat                                                                                      >  1:4320800/1‑93 (MQ=255)

   gagTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACtcat                                                                                      >  1:4242982/1‑93 (MQ=255)

      tttttGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAg                                                                                   >  2:648757/1‑93 (MQ=255)

            ctgcGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCaaa                                                                             >  2:2534878/1‑93 (MQ=255)

             tgcGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAg                                                                            >  1:421904/1‑93 (MQ=255)

               cgcgGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGcc                                                                          <  1:3785228/93‑1 (MQ=255)

                gcgGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGccc                                                                         <  1:3251615/93‑1 (MQ=255)

                gcgGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGccc                                                                         >  1:2341073/1‑93 (MQ=255)

                 cgGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCa                                                                        <  1:3917130/93‑1 (MQ=255)

                  ggAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCaa                                                                       <  2:1491987/93‑1 (MQ=255)

                  ggAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCaa                                                                       >  1:17858/1‑93 (MQ=255)

                  ggAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCaa                                                                       >  2:2177299/1‑93 (MQ=255)

                   gAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAAt                                                                      <  1:3693546/93‑1 (MQ=255)

                   gAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAAt                                                                      >  1:109910/1‑93 (MQ=255)

                   gAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAAt                                                                      >  2:163819/1‑93 (MQ=255)

                   gAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAAt                                                                      <  2:307570/93‑1 (MQ=255)

                   gAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAAt                                                                      >  2:3392351/1‑93 (MQ=255)

                     aaTTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATgg                                                                    >  2:1244515/1‑93 (MQ=255)

                     aaTTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATgg                                                                    >  2:3104777/1‑93 (MQ=255)

                      aTTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGgc                                                                   >  1:3218940/1‑93 (MQ=255)

                         aaCGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCa                                                                >  1:705563/1‑93 (MQ=255)

                         aaCGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCa                                                                >  1:4689401/1‑93 (MQ=255)

                            gCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGc                                                             >  2:2645642/1‑93 (MQ=255)

                             cTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCt                                                            <  2:4579632/93‑1 (MQ=255)

                             cTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCt                                                            <  2:727937/93‑1 (MQ=255)

                              ttCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTg                                                           >  1:680243/1‑93 (MQ=255)

                              ttCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTg                                                           >  1:3575068/1‑93 (MQ=255)

                               tCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATTGTGCAAGCTGa                                                          >  2:4208909/1‑93 (MQ=255)

                                ccAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGAc                                                         <  1:2254645/93‑1 (MQ=255)

                                 cAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACg                                                        >  2:8836/1‑93 (MQ=255)

                                 cAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACg                                                        <  2:357304/93‑1 (MQ=255)

                                 cAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACg                                                        >  2:4294577/1‑93 (MQ=255)

                                 cAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGGTTCTGACTCATCAGGTCAAAGCCAAATGGCGCAAGCTGACg                                                        >  1:98866/1‑93 (MQ=255)

                                 cAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCCTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACg                                                        >  2:4045400/1‑93 (MQ=255)

                                  aGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGt                                                       <  1:1703222/93‑1 (MQ=255)

                                  aGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGt                                                       <  1:108087/93‑1 (MQ=255)

                                   ggAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTg                                                      >  1:3293074/1‑93 (MQ=255)

                                   ggAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTg                                                      >  1:3164579/1‑93 (MQ=255)

                                   ggAACAGAAAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTg                                                      >  1:2052216/1‑93 (MQ=255)

                                     aaCAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGa                                                    >  1:291743/1‑93 (MQ=255)

                                      aCAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGAt                                                   >  2:3671497/1‑93 (MQ=255)

                                        aGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATaa                                                 >  1:3282839/1‑93 (MQ=255)

                                         gATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATaaa                                                >  1:2717562/1‑93 (MQ=255)

                                          aTAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAAc                                               >  1:3016413/1‑93 (MQ=255)

                                          aTAGCGGGCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAAc                                               <  1:2575453/93‑1 (MQ=255)

                                                gtcgtcATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCcgtc                                         <  1:4046645/93‑1 (MQ=255)

                                                gtcgtcATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCcgtc                                         <  2:1649304/93‑1 (MQ=255)

                                                gtcgtcATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCcgtc                                         >  2:2613961/1‑93 (MQ=255)

                                                 tcgtcATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCcgtcg                                        >  1:3085065/1‑93 (MQ=255)

                                                 tcgtcATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCcgtcg                                        >  1:3085064/1‑93 (MQ=255)

                                                 tcgtcATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCcgtcg                                        >  2:1139741/1‑93 (MQ=255)

                                                   gtcATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGtt                                      >  1:475711/1‑93 (MQ=255)

                                                    tcATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAAc                                               <  1:4106709/83‑1 (MQ=255)

                                                    tcATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAAc                                               >  2:4106709/1‑83 (MQ=255)

                                                    tcATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTc                                     >  1:1606519/1‑93 (MQ=255)

                                                      aTCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGCCGTTCAt                                   >  1:1823208/1‑93 (MQ=255)

                                                       tCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCAtt                                  <  2:2690583/93‑1 (MQ=255)

                                                       tCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCAtc                                  >  2:391081/1‑92 (MQ=255)

                                                        cgcgACGGCTATGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTc                                 <  1:1245762/93‑1 (MQ=255)

                                                         gcgACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTcc                                >  2:2527709/1‑93 (MQ=255)

                                                          cgACGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTccc                               >  2:4317889/1‑93 (MQ=255)

                                                            aCGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCcag                             >  1:2529620/1‑93 (MQ=255)

                                                             cGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCcagc                            >  1:2220293/1‑93 (MQ=255)

                                                             cGGCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCcagc                            <  2:4745471/93‑1 (MQ=255)

                                                              ggCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCcagca                           <  1:592757/93‑1 (MQ=255)

                                                               gCTGCGGTCGCCACACTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCcagcag                          <  1:3035694/93‑1 (MQ=255)

                                                               gCTACGGTCTCCACCCTGCGGTTTCTGACTCATCAGGTCAAACCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCcagcag                          <  1:834864/93‑1 (MQ=255)

                                                               gCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCcagcag                          >  1:3155290/1‑93 (MQ=255)

                                                               gCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCcagcag                          <  1:839219/93‑1 (MQ=255)

                                                               gCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCcagcag                          <  1:3995280/93‑1 (MQ=255)

                                                               gCTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCcagcag                          >  2:401193/1‑93 (MQ=255)

                                                                gtACGGTCGCCGCGCTTCGGTTTCTGACTCATCAGGGCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAgg                         <  1:3611257/92‑1 (MQ=255)

                                                                cTACGGTCGCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAgg                         >  1:1319407/1‑93 (MQ=255)

                                                                   cGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCaa                      <  1:891217/93‑1 (MQ=255)

                                                                      tCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAAcc                   >  1:4459211/1‑93 (MQ=255)

                                                                        gCCACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCAc                 >  2:3220444/1‑93 (MQ=255)

                                                                         ccACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACt                >  2:1291877/1‑93 (MQ=255)

                                                                          cACGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACtt               <  1:955792/93‑1 (MQ=255)

                                                                           aCGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACttt              <  2:1720905/93‑1 (MQ=255)

                                                                           aCGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACttt              >  2:1860264/1‑93 (MQ=255)

                                                                            cGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGACGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAACCACCTTTg             >  2:3861435/1‑93 (MQ=255)

                                                                            cGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTg             >  1:1498319/1‑93 (MQ=255)

                                                                            cGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTg             <  1:4075911/93‑1 (MQ=255)

                                                                            cGCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTg             >  2:3759251/1‑93 (MQ=255)

                                                                             gCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGa            >  1:1294826/1‑93 (MQ=255)

                                                                             gCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGa            >  1:1312055/1‑93 (MQ=255)

                                                                             gCTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGa            >  1:3646476/1‑93 (MQ=255)

                                                                              cTTCGGTTTCTGACTCATCATGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGaa           <  2:377767/93‑1 (MQ=255)

                                                                              cTTCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGaa           <  1:1578783/93‑1 (MQ=255)

                                                                               ttCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAAc          <  2:4512751/93‑1 (MQ=255)

                                                                               ttCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAAc          >  1:2462625/1‑93 (MQ=255)

                                                                                tCGGTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACg         >  1:3496993/1‑93 (MQ=255)

                                                                                  ggTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGa       >  2:1597988/1‑93 (MQ=255)

                                                                                   gTTTCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGaa      >  1:2108730/1‑93 (MQ=255)

                                                                                    tttCTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAAt     >  2:4227954/1‑93 (MQ=255)

                                                                                    tttCTGACGCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTGGAACGGAAt     <  2:2274536/93‑1 (MQ=255)

                                                                                       cTGACTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGa  <  1:1444363/93‑1 (MQ=255)

ATAGAGTTTTTGCTGCGCGGAAATTAACGCTTCCAGGAACAGATAGCGGTCGTCATCGCGACGGCTACGGTCGCCACGCTTCGGTTTCTGGCTCATCAGGTCAAAGCCCAATGGCGCAAGCTGACGTGGATAAACGCCGTCGTTCATTCCCAGCAGGCAAACCACTTTGAACGGAATTGA  >  NC_000913/2958267‑2958446

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: `‑`