BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,601,332	A→G	Q96R (CAG→CGG)	yjjQ →	putative regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,601,332	0	A	G	81.5%	688.4 / 166.1	335	Q96R (CAG→CGG)	yjjQ	putative regulator
Reads supporting (aligned to +/- strand): ref base A (42/20); new base G (148/125); total (190/145)
Fisher's exact test for biased strand distribution p-value = 6.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01

Only 100 of 344 total aligned reads displayed.

CCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCAGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGACCGTACGTTCCGGCAATA  >  NC_000913/4601241‑4601424

 cACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCggg                                                                                            >  2:5101743/1‑93 (MQ=255)

  aCGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGc                                                                                           >  1:1373642/1‑93 (MQ=255)

  aCGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGc                                                                                           >  2:7457018/1‑93 (MQ=255)

  aCGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGc                                                                                           >  2:6929963/1‑93 (MQ=255)

   cGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCa                                                                                          <  1:2068376/93‑1 (MQ=255)

    gTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAg                                                                                         >  2:5460712/1‑93 (MQ=255)

     ttCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCAGGCAGt                                                                                        >  1:3947839/1‑93 (MQ=255)

      tCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCAGGCAGTa                                                                                       <  1:6119034/93‑1 (MQ=255)

      tCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCAGGCAGTa                                                                                       <  1:5592914/93‑1 (MQ=255)

       cTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAg                                                                                      <  2:650574/93‑1 (MQ=255)

       cTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAg                                                                                      <  1:3621851/93‑1 (MQ=255)

        tGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGa                                                                                     <  1:4288033/93‑1 (MQ=255)

         gTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTGACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGaa                                                                                    >  2:927685/1‑93 (MQ=255)

          tttGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCAGGCAGTAGAAc                                                                                   >  1:3413413/1‑93 (MQ=255)

             gtgAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGGGTCACCCTCCTGGTATTCCCGGGCAGTAGAACTGc                                                                                >  1:3663409/1‑93 (MQ=255)

              tgAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCt                                                                               <  1:5900417/93‑1 (MQ=255)

                aaCATTATTATACTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCa                                                                             >  2:4655691/1‑93 (MQ=255)

                 aCATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCAt                                                                            >  1:5745413/1‑93 (MQ=255)

                   attattattCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCAGGCAGTAGAACTGCTCAtgt                                                                          >  2:2854725/1‑93 (MQ=255)

                    ttattattCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCAtgtg                                                                         >  1:911110/1‑93 (MQ=255)

                    ttattattCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCAGGCAGTAGAACTGCTCAtgtg                                                                         >  2:5837932/1‑93 (MQ=255)

                      attattCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGcc                                                                       >  1:6933856/1‑93 (MQ=255)

                       ttattCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGccc                                                                      >  1:7121090/1‑93 (MQ=255)

                        tattCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCt                                                                     <  2:3240278/93‑1 (MQ=255)

                         attCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTa                                                                    <  1:2172076/93‑1 (MQ=255)

                           tCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACg                                                                  >  1:2767207/1‑93 (MQ=255)

                           tCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCAGGCAGTAGAACTGCTCATGTGCCCTACg                                                                  >  1:2849995/1‑93 (MQ=255)

                            cTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCAGGCAGTAGAACTGCTCATGTGCCCTACgg                                                                 <  1:487114/93‑1 (MQ=255)

                             tttAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCAGGCAGTAGAACTGCTCATGTGCCCTACGGc                                                                >  2:4134167/1‑93 (MQ=255)

                               tAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGa                                                              <  2:2327297/93‑1 (MQ=255)

                               tAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCAGGCAGTAGAACTGCTCATGTGCCCTACGGCGa                                                              <  1:6820100/93‑1 (MQ=255)

                                 aTCTCATAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCAGGCAGTAGAACTGCTCATGTGCCCTACGGCGACg                                                            >  2:4776470/1‑93 (MQ=255)

                                 aTCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACg                                                            <  1:1301031/93‑1 (MQ=255)

                                 aTCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACg                                                            <  2:5245166/93‑1 (MQ=255)

                                 aTCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCAGGCAGTAGAACTGCTCATGTGCCCTACGGCGACg                                                            >  1:3351267/1‑93 (MQ=255)

                                 aTCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCAGGCAGTAGAACTGCTCATGTGCCCTACGGCGACg                                                            >  1:3494063/1‑93 (MQ=255)

                                    tcACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTa                                                         >  1:4621230/1‑93 (MQ=255)

                                    tcACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTa                                                         <  2:7180490/93‑1 (MQ=255)

                                    tcACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCAGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTa                                                         <  1:6564913/93‑1 (MQ=255)

                                     cacaATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTAt                                                        >  2:7307785/1‑93 (MQ=255)

                                     cacaATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTAt                                                        <  2:5279628/93‑1 (MQ=255)

                                      acaATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTAtt                                                       <  2:1145641/93‑1 (MQ=255)

                                      acaATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTAtt                                                       >  1:5805553/1‑93 (MQ=255)

                                      acaATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTAtt                                                       <  2:6795907/93‑1 (MQ=255)

                                         atatCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTc                                                    <  2:2711450/93‑1 (MQ=255)

                                           atCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTg                                                  >  1:1221783/1‑93 (MQ=255)

                                            tCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGa                                                 <  2:4332580/93‑1 (MQ=255)

                                                gAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGtt                                             <  1:2272197/93‑1 (MQ=255)

                                                gAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGtt                                             >  2:6106705/1‑93 (MQ=255)

                                                gAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGtt                                             >  1:1987369/1‑93 (MQ=255)

                                                 aaaTTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTg                                            >  2:3499398/1‑93 (MQ=255)

                                                   aTTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGaa                                          >  2:1877298/1‑93 (MQ=255)

                                                   aTTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCAGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGaa                                          <  2:4960442/93‑1 (MQ=255)

                                                    ttaaCTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCATACGGCGACGTTATTGTCTGATGTTGAAc                                         >  1:2713262/4‑93 (MQ=255)

                                                    ttCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAAc                                         <  2:1775113/93‑1 (MQ=255)

                                                        cTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCAt                                     >  2:3429588/1‑93 (MQ=255)

                                                        cTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCAt                                     >  1:5233838/1‑93 (MQ=255)

                                                        cTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCAt                                     >  1:5048480/1‑93 (MQ=255)

                                                         tGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCAtt                                    >  1:4406178/1‑93 (MQ=255)

                                                         tGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCAtt                                    <  2:3666908/93‑1 (MQ=255)

                                                         tGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCAtt                                    <  1:1614777/93‑1 (MQ=255)

                                                           ggTTTTCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTga                                  <  2:3102593/93‑1 (MQ=255)

                                                           ggTTTTCATGGTGTCACGCTCCTGGTATTCCCAGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTga                                  >  1:7519241/1‑93 (MQ=255)

                                                              tttCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGaa                               <  2:7095917/93‑1 (MQ=255)

                                                              tttCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGaa                               <  1:6676694/93‑1 (MQ=255)

                                                              tttCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGaa                               <  2:6631379/93‑1 (MQ=255)

                                                                tCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATc                             <  1:3074230/93‑1 (MQ=255)

                                                                tCATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATc                             <  1:929126/93‑1 (MQ=255)

                                                                 cATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCt                            <  1:3807194/93‑1 (MQ=255)

                                                                 cATGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCt                            <  1:2432671/93‑1 (MQ=255)

                                                                  aTGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTg                           >  1:500987/1‑93 (MQ=255)

                                                                   tGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTgg                          >  2:3287410/1‑93 (MQ=255)

                                                                   tGGTGTCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTgg                          >  1:4138035/1‑93 (MQ=255)

                                                                     gtgtCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTggnn                        >  2:1504297/1‑91 (MQ=255)

                                                                     gtgtCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTc                        <  2:5588545/93‑1 (MQ=255)

                                                                       gtCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCaa                      <  1:4795714/93‑1 (MQ=255)

                                                                       gtCACGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAACCTGGTCaa                      <  2:4857111/93‑1 (MQ=255)

                                                                          aCGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGAc                   >  1:1124520/1‑93 (MQ=255)

                                                                          aCGCTCCTGGTATTCCCAGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGAc                   <  2:1394712/93‑1 (MQ=255)

                                                                           cGCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGAcc                  >  2:7411964/1‑93 (MQ=255)

                                                                            gCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGACCg                 >  1:1031556/1‑93 (MQ=255)

                                                                            gCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGACCg                 <  1:7217479/93‑1 (MQ=255)

                                                                            gCTCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGACCg                 >  1:3534595/1‑93 (MQ=255)

                                                                             cTCCTGGTATTCCCAGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGACCGt                >  2:2157270/1‑93 (MQ=255)

                                                                              tCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGACCGTa               >  2:2967278/1‑93 (MQ=255)

                                                                              tCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGACCGTa               >  2:84542/1‑93 (MQ=255)

                                                                              tCCTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGACCGTa               <  1:181342/93‑1 (MQ=255)

                                                                               ccTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGACCGTAc              <  2:786665/93‑1 (MQ=255)

                                                                               ccTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGACCGTAc              <  1:4629952/93‑1 (MQ=255)

                                                                               ccTGGTATTCCCAGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGACCGTAc              <  1:3251423/93‑1 (MQ=255)

                                                                                cTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGACCGTACg             >  1:6852709/1‑93 (MQ=255)

                                                                                cTGGTATTCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGACCGTACg             <  1:6433859/93‑1 (MQ=255)

                                                                                      ttCCCAGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGACCGTACGTTCCgg       >  1:5532031/1‑93 (MQ=255)

                                                                                      ttCCCAGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGACCGTACGTTCCgg       <  1:848315/93‑1 (MQ=255)

                                                                                       tCCCGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGACCGTACGTTCCGGc      <  2:4449532/93‑1 (MQ=255)

                                                                                        cccGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGACCGTACGTTCCGGCa     >  2:1601501/1‑93 (MQ=255)

                                                                                        cccGGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGACCGTACGTTCCGGCa     <  2:536073/93‑1 (MQ=255)

                                                                                         ccAGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGACCGTACGTTCCGGCaa    <  2:258530/93‑1 (MQ=255)

                                                                                           aGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGACCGTACGTTCCGGCAATa  >  2:6430346/1‑93 (MQ=255)

                                                                                           aGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGACCGTACGTTCCGGCAATa  <  2:3375093/93‑1 (MQ=255)

CCACGTTCTGTTTGTGAACATTATTATTCTTTAATCTCACAATATCGGGAAATTCACTGGGTTTTCATGGTGTCACGCTCCTGGTATTCCCAGGCAGTAGAACTGCTCATGTGCCCTACGGCGACGTTATTGTCTGATGTTGAACCCATTGAGAATCTGGTCAAGACCGTACGTTCCGGCAATA  >  NC_000913/4601241‑4601424

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