BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,079,216	T→C	E23E (GAA→GAG)	sbmC ←	SbmC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,079,216	0	T	C	96.4%	397.8 / 3.1	137	E23E (GAA→GAG)	sbmC	SbmC
Reads supporting (aligned to +/- strand): ref base T (2/3); new base C (72/60); total (74/63)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

Only 100 of 137 total aligned reads displayed.

GTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGTTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTT  >  NC_000913/2079129‑2079306

 tAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCAc                                                                                      <  1:1371879/93‑1 (MQ=255)

 tAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCAc                                                                                      >  1:2654386/1‑93 (MQ=255)

  aaTAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACg                                                                                     >  2:2504461/1‑93 (MQ=255)

   aTAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACgg                                                                                    >  2:1669150/1‑93 (MQ=255)

   aTAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACgg                                                                                    >  1:4054731/1‑93 (MQ=255)

    tAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACggg                                                                                   >  1:930102/1‑93 (MQ=255)

    tAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACggg                                                                                   >  2:1589189/1‑93 (MQ=255)

     aGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGc                                                                                  <  1:279188/93‑1 (MQ=255)

      gACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGcc                                                                                 <  2:453497/93‑1 (MQ=255)

        cAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCaa                                                                               <  2:1710291/93‑1 (MQ=255)

          gCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACg                                                                             >  1:2147071/1‑93 (MQ=255)

           cAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACga                                                                            <  2:2505660/93‑1 (MQ=255)

            aaCCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACgag                                                                           >  1:3950002/1‑93 (MQ=255)

            aaCCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACAGTCTGCTCCCACGGGCCAACgag                                                                           >  1:2933244/1‑93 (MQ=255)

             aCCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGTTCCCACGGGCCAACgaga                                                                          <  1:3858808/93‑1 (MQ=255)

             aCCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACgaga                                                                          <  1:992921/93‑1 (MQ=255)

              cccACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGAt                                                                         >  1:2829239/1‑93 (MQ=255)

                cACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATgg                                                                       >  1:2434684/1‑93 (MQ=255)

                 aCTCCTTTGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGa                                                                      >  2:2773241/1‑93 (MQ=255)

                 aCTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGa                                                                      >  2:3296940/1‑93 (MQ=255)

                  cTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCCACGAGATGGaa                                                                     >  2:2070446/1‑93 (MQ=255)

                  cTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGaa                                                                     >  1:3875488/1‑93 (MQ=255)

                  cTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGaa                                                                     >  2:3284885/1‑93 (MQ=255)

                   tCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGaaa                                                                    >  1:2808160/1‑93 (MQ=255)

                   cccTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGaaa                                                                    <  1:4145206/92‑1 (MQ=255)

                         ggCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCa                                                              >  2:1507686/1‑93 (MQ=255)

                         ggCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCa                                                              <  2:541717/93‑1 (MQ=255)

                         ggCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCa                                                              >  1:3355705/1‑93 (MQ=255)

                           cacaATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAAc                                                            >  2:1201283/1‑93 (MQ=255)

                           cacaATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAAc                                                            <  2:3161519/93‑1 (MQ=255)

                             caATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACgg                                                          <  1:4386538/93‑1 (MQ=255)

                             caATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACgg                                                          <  1:1542593/93‑1 (MQ=255)

                              aaTATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGt                                                         >  2:1451990/1‑93 (MQ=255)

                               atatTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTa                                                        >  2:1577176/1‑93 (MQ=255)

                               atatTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTa                                                        <  1:2322845/93‑1 (MQ=255)

                               atatTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTa                                                        <  2:1699482/93‑1 (MQ=255)

                                tatTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCTTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTAc                                                       <  2:1148128/93‑1 (MQ=255)

                                tatTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTAc                                                       >  1:1012390/1‑93 (MQ=255)

                                       gCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTtctc                                                >  2:3124113/1‑93 (MQ=255)

                                       gCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTtctc                                                <  1:2319359/93‑1 (MQ=255)

                                          aTCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTc                                             >  1:2104370/1‑93 (MQ=255)

                                           tCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTcc                                            >  1:2124350/1‑93 (MQ=255)

                                           tCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTcc                                            <  1:515117/93‑1 (MQ=255)

                                           tCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTcc                                            >  1:230131/1‑93 (MQ=255)

                                              aCCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGc                                         <  2:367236/93‑1 (MQ=255)

                                              aCCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGc                                         >  2:1921872/1‑93 (MQ=255)

                                              aCCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGc                                         >  1:3688976/1‑93 (MQ=255)

                                               cccACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCt                                        >  1:3136117/1‑93 (MQ=255)

                                               cccACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCt                                        >  2:1975665/1‑93 (MQ=255)

                                               cccACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCt                                        >  2:1410989/1‑93 (MQ=255)

                                                ccACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCtt                                       >  2:3235505/1‑93 (MQ=255)

                                                ccACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCtt                                       >  2:4447713/1‑93 (MQ=255)

                                                  acaTCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTaa                                     >  2:2981149/1‑93 (MQ=255)

                                                  acaTCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTaa                                     <  1:1628393/93‑1 (MQ=255)

                                                    atcatcAACTGCTCAAAGCCTTTCTTTACTGTCTGCTCCCACGGTCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAAtc                                   <  1:2941807/93‑1 (MQ=255)

                                                     tcatcaACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAAtct                                  >  2:3540520/1‑93 (MQ=255)

                                                       atcaACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCg                                >  2:2892317/1‑93 (MQ=255)

                                                       atcaACTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCg                                <  2:3679604/93‑1 (MQ=255)

                                                           aCTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGt                            >  1:1666807/1‑93 (MQ=255)

                                                           aCTGCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGt                            <  1:3854035/93‑1 (MQ=255)

                                                           aCTGCTCAAAGCCTTTCTTCACCGTCTGCTCCCACGGGCCAACGCGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGt                            >  2:4584403/1‑93 (MQ=255)

                                                              gCTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCa                         >  2:2403754/1‑93 (MQ=255)

                                                               cTCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCAt                        <  1:691146/93‑1 (MQ=255)

                                                                tCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATg                       <  1:2310669/93‑1 (MQ=255)

                                                                tCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATg                       >  1:2787426/1‑93 (MQ=255)

                                                                tCAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTATTCCTGCTTAATCTCGTCGTTCATg                       >  2:2595051/1‑93 (MQ=255)

                                                                 cAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGa                      <  1:3311177/93‑1 (MQ=255)

                                                                 cAAAGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGa                      <  2:3607441/93‑1 (MQ=255)

                                                                  aaaGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAACCCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGAc                     <  2:309036/93‑1 (MQ=255)

                                                                    aGCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGc                   >  2:2359446/1‑93 (MQ=255)

                                                                     gCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGcc                  >  1:1747611/1‑93 (MQ=255)

                                                                     gCCTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGcc                  <  1:172680/93‑1 (MQ=255)

                                                                      ccTTTCTTTACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCt                 >  1:2903965/1‑93 (MQ=255)

                                                                       ctttctttACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTc                >  2:1149558/1‑93 (MQ=255)

                                                                        tttctttACCGTCTGTTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTcc               <  1:4409715/93‑1 (MQ=255)

                                                                        tttctttACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTcc               >  2:658775/1‑93 (MQ=255)

                                                                        tttctttACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTcc               >  2:2259199/1‑93 (MQ=255)

                                                                        tttctttACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTcc               >  2:2850389/1‑93 (MQ=255)

                                                                        tttctttACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTcc               >  2:2908636/1‑93 (MQ=255)

                                                                         ttctttACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCa              <  1:2256470/93‑1 (MQ=255)

                                                                         ttctttACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCa              <  1:1594955/93‑1 (MQ=255)

                                                                         ttctttACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCa              <  1:153816/93‑1 (MQ=255)

                                                                           ctttACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCAtt            <  1:2949597/93‑1 (MQ=255)

                                                                            tttACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTg           <  1:899854/93‑1 (MQ=255)

                                                                            tttACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTg           >  2:1792645/1‑93 (MQ=255)

                                                                            tttACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTg           >  2:1913690/1‑93 (MQ=255)

                                                                             ttACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGa          <  1:1526017/93‑1 (MQ=255)

                                                                              tACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGAt         <  2:4206158/93‑1 (MQ=255)

                                                                              tACCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGAt         <  2:1979293/93‑1 (MQ=255)

                                                                               aCCGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATa        <  1:2549101/93‑1 (MQ=255)

                                                                                ccGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATAc       <  1:610223/93‑1 (MQ=255)

                                                                                ccGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATAc       >  1:3463377/1‑93 (MQ=255)

                                                                                ccGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATAc       <  1:2739788/93‑1 (MQ=255)

                                                                                ccGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATAc       <  1:2739787/93‑1 (MQ=255)

                                                                                ccGTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATAc       <  2:2580655/93‑1 (MQ=255)

                                                                                  gTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTg     <  2:3073685/93‑1 (MQ=255)

                                                                                  gTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTg     <  1:1436896/93‑1 (MQ=255)

                                                                                  gTCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTg     >  2:2512880/1‑93 (MQ=255)

                                                                                   tCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGt    <  2:1728392/93‑1 (MQ=255)

                                                                                   tCTGCTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGt    >  1:2407052/1‑93 (MQ=255)

GTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGTTCCCACGGGCCAACGAGATGGAAACCTGCAACGGTACGTTTCTCTTCCTGCTTAATCTCGTAGTTCATGACGCCTCCATTGATACTGTTT  >  NC_000913/2079129‑2079306

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