BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,039,349	T→C	P77P (CCA→CCG)	mobA ←	MobA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,039,349	0	T	C	99.0%	973.7 / ‑2.8	287	P77P (CCA→CCG)	mobA	MobA
Reads supporting (aligned to +/- strand): ref base T (2/1); new base C (149/135); total (151/136)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

Only 100 of 297 total aligned reads displayed.

AGGGGGAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCTGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGTGAGAGAG  >  NC_000913/4039259‑4039441

aGGGGGAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCgg                                                                                            <  2:7036676/93‑1 (MQ=255)

 gggggAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCggg                                                                                           >  1:6658266/1‑93 (MQ=255)

 gggggAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCggg                                                                                           <  1:3750909/93‑1 (MQ=255)

    ggAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTaa                                                                                        <  1:7461401/93‑1 (MQ=255)

    ggAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTaa                                                                                        >  1:2325045/1‑93 (MQ=255)

     gAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAAt                                                                                       <  1:4935349/93‑1 (MQ=255)

     gAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAAt                                                                                       >  2:3467242/1‑93 (MQ=255)

     gAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAAt                                                                                       <  1:2850432/93‑1 (MQ=255)

     gAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAAt                                                                                       <  2:5838260/93‑1 (MQ=255)

       aTGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATcc                                                                                     <  2:8634840/93‑1 (MQ=255)

       aTGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATcc                                                                                     <  1:682459/93‑1 (MQ=255)

         gTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGc                                                                                   <  1:264549/93‑1 (MQ=255)

            aGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAg                                                                                >  2:7149084/1‑93 (MQ=255)

             ggCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGt                                                                               >  1:5355558/1‑93 (MQ=255)

              gCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTg                                                                              >  2:2633919/1‑93 (MQ=255)

               cGTATCGCGCGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGa                                                                             <  1:8729962/93‑1 (MQ=255)

               cGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGa                                                                             <  2:781434/93‑1 (MQ=255)

               cGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGa                                                                             >  1:8348982/1‑93 (MQ=255)

               cGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGa                                                                             <  2:6599763/93‑1 (MQ=255)

               cGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGa                                                                             <  1:3968036/93‑1 (MQ=255)

                 tATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCTCCAGTGAAt                                                                           >  2:7885109/1‑93 (MQ=255)

                   tCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCt                                                                         <  1:4335667/93‑1 (MQ=255)

                   tCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCt                                                                         <  2:8210490/93‑1 (MQ=255)

                    cGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCtt                                                                        <  1:7745590/93‑1 (MQ=255)

                         gggCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTACTCCGACCGTGAATCTTCAATc                                                                   >  2:2958365/1‑93 (MQ=255)

                          ggCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCa                                                                  >  2:1281015/1‑93 (MQ=255)

                          ggCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCa                                                                  <  1:5177878/93‑1 (MQ=255)

                          ggCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCa                                                                  >  1:3814012/1‑93 (MQ=255)

                           gCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCTGGGTAATCCGCCAGTGAATCTTCAATCAc                                                                 <  1:1326492/93‑1 (MQ=255)

                           gCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCAc                                                                 >  2:2391055/1‑93 (MQ=255)

                            caaaacaaaaACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACt                                                                >  1:1103042/1‑93 (MQ=255)

                            caaaacaaaaACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACt                                                                <  1:7853640/93‑1 (MQ=255)

                             aaaacaaaaACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACtt                                                               <  1:6482777/93‑1 (MQ=255)

                              aaacaaaaACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACttt                                                              >  1:7766206/1‑93 (MQ=255)

                               aacaaaaACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTc                                                             >  2:3708594/1‑93 (MQ=255)

                               aacaaaaACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTc                                                             >  2:3541560/1‑93 (MQ=255)

                               aacaaaaACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTc                                                             <  2:3441342/93‑1 (MQ=255)

                               aacaaaaACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTc                                                             <  1:4549884/93‑1 (MQ=255)

                                acaaaaACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCa                                                            >  1:5877426/1‑93 (MQ=255)

                                 caaaaACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAg                                                           >  1:7057383/1‑93 (MQ=255)

                                 caaaaACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAg                                                           >  1:5812649/1‑93 (MQ=255)

                                   aaaaCCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCCGAGGGCCCGGGTAATCCGCCAGTGCATCTTCAATCCCTTTCAgcc                                                         >  2:4314356/1‑91 (MQ=255)

                                   aaaaCCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGAc                                                         >  2:4664861/1‑93 (MQ=255)

                                     aaCCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCg                                                       >  2:6345071/1‑93 (MQ=255)

                                     aaCCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCg                                                       <  1:6921858/93‑1 (MQ=255)

                                      aCCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGc                                                      <  2:7283322/93‑1 (MQ=255)

                                      aCCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGc                                                      <  1:6246786/93‑1 (MQ=255)

                                        cACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCtt                                                    >  2:1443404/1‑93 (MQ=255)

                                        cACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCtt                                                    <  2:5099433/93‑1 (MQ=255)

                                           tCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGcc                                                 >  1:8595453/1‑93 (MQ=255)

                                           ccACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGcc                                                 <  1:6581283/92‑1 (MQ=255)

                                           ccACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGcc                                                 <  2:2092672/92‑1 (MQ=255)

                                            cACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCt                                                >  1:8165343/1‑93 (MQ=255)

                                            cACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCt                                                >  1:2521956/1‑93 (MQ=255)

                                             aCCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTg                                               >  1:8635700/1‑93 (MQ=255)

                                              ccAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGa                                              >  2:2789230/1‑93 (MQ=255)

                                                aGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATa                                            >  1:4832402/1‑93 (MQ=255)

                                                 gCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAg                                           <  2:2009661/93‑1 (MQ=255)

                                                 gCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAg                                           >  2:3961820/1‑93 (MQ=255)

                                                     ccTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGAttt                                       >  1:5493511/1‑93 (MQ=255)

                                                      ctgctgCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTc                                      >  2:6966333/1‑93 (MQ=255)

                                                      ctgctgCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTc                                      >  1:6434502/1‑93 (MQ=255)

                                                      ctgctgCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTc                                      <  1:7174667/93‑1 (MQ=255)

                                                        gctgcATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCt                                    <  2:3032950/93‑1 (MQ=255)

                                                           gcATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCtgat                                 >  2:1658567/1‑93 (MQ=255)

                                                           gcATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCtgat                                 >  2:3599616/1‑93 (MQ=255)

                                                             aTTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCtgatga                               <  2:1510950/93‑1 (MQ=255)

                                                             aTTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTACCtgatga                               >  2:8380482/1‑93 (MQ=255)

                                                               tACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACg                             >  1:6335753/1‑93 (MQ=255)

                                                                 cTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGAt                           <  1:670030/93‑1 (MQ=255)

                                                                  tGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGAtt                          <  2:4942547/93‑1 (MQ=255)

                                                                   gAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTa                         >  1:3172266/1‑93 (MQ=255)

                                                                    aaaGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAg                        >  2:8026406/1‑93 (MQ=255)

                                                                    aaaGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAg                        <  2:7496775/93‑1 (MQ=255)

                                                                     aaGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGc                       >  1:4648564/1‑93 (MQ=255)

                                                                      aGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCa                      <  1:4088627/93‑1 (MQ=255)

                                                                       gCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCAt                     <  2:770257/93‑1 (MQ=255)

                                                                       gCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCAt                     <  2:557636/93‑1 (MQ=255)

                                                                       gCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCAt                     >  2:5780051/1‑93 (MQ=255)

                                                                        cATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCAtt                    <  1:7500826/93‑1 (MQ=255)

                                                                          ttCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTaa                  >  2:4525059/1‑93 (MQ=255)

                                                                          ttCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTaa                  >  2:8229261/1‑93 (MQ=255)

                                                                          ttCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTaa                  >  2:3812583/1‑93 (MQ=255)

                                                                            ccTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAAcc                >  2:5283520/1‑93 (MQ=255)

                                                                             cTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCa               >  1:3057080/1‑93 (MQ=255)

                                                                             cTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCa               <  2:4474465/93‑1 (MQ=255)

                                                                              tgtgagagGGCCCGGGTAAACCGCCAGTGAATCTTCAATCCCTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCac              <  2:2224180/89‑1 (MQ=255)

                                                                              tGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCac              >  2:6754833/1‑93 (MQ=255)

                                                                              tGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCac              <  1:806620/93‑1 (MQ=255)

                                                                               gCCAGAGGGCCCGGGTAATACGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCacg             >  1:5946344/1‑93 (MQ=255)

                                                                                 cAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCacgac           <  2:181864/93‑1 (MQ=255)

                                                                                 cAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCacgac           >  2:1022578/1‑93 (MQ=255)

                                                                                 cAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCacgac           >  2:4698741/1‑93 (MQ=255)

                                                                                  agagGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCacgacg          >  1:1792776/1‑93 (MQ=255)

                                                                                  agagGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCacgacg          >  1:1693084/1‑93 (MQ=255)

                                                                                   gagGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGt         >  2:6175482/1‑93 (MQ=255)

                                                                                   gagGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGt         >  1:6515110/1‑93 (MQ=255)

                                                                                    agGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGTg        <  2:1547409/93‑1 (MQ=255)

                                                                                     gggCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGTga       <  1:1502206/93‑1 (MQ=255)

                                                                                          tGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGTGagagag  <  1:1462529/93‑1 (MQ=255)

AGGGGGAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCTGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGTGAGAGAG  >  NC_000913/4039259‑4039441

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