BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,011,893	G→A	A123V (GCG→GTG)	b2875 ←	putative synthases

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,011,893	0	G	A	99.1%	1135.1 / ‑4.8	336	A123V (GCG→GTG)	b2875	putative synthases
Reads supporting (aligned to +/- strand): ref base G (3/0); new base A (169/161); total (172/164)
Fisher's exact test for biased strand distribution p-value = 2.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 356 total aligned reads displayed.

GCATGGAGAAGCGTGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCGCCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGACGCTGATAAATACTGAC  >  NC_000913/3011803‑3011982

gCATGGAGAAGCGTGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGcacc                                                                                         >  1:5806048/1‑93 (MQ=255)

gCATGGAGAAGCGTGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGcacc                                                                                         <  1:5682655/93‑1 (MQ=255)

gCATGGAGAAGCGTGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGcacc                                                                                         <  1:4436009/93‑1 (MQ=255)

gCATGGAGAAGCGTGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGcacc                                                                                         <  2:5959565/93‑1 (MQ=255)

  aTGGAGAAGCGTGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGcaccgc                                                                                       <  1:183604/93‑1 (MQ=255)

   tGGAGAAGCGTGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCa                                                                                      >  2:6134572/1‑93 (MQ=255)

    ggAGAAGCGTGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCAc                                                                                     >  1:954344/1‑93 (MQ=255)

     gagaAGCGTGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCAct                                                                                    >  2:2284292/1‑93 (MQ=255)

     gagaAGCGTGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCAct                                                                                    >  1:5345496/1‑93 (MQ=255)

      agaAGCGTGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCActc                                                                                   <  2:3367966/93‑1 (MQ=255)

       gaAGCGTGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCActct                                                                                  <  1:9884809/93‑1 (MQ=255)

       gaAGCGTGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCActct                                                                                  <  1:427696/93‑1 (MQ=255)

        aaGCGTGGTTGATGGTGGGAATAGTTCAGGATTGATGCTTTCGCTATAAATATCAGCAACTGCAATATCAAATCCTAATAGCCCCGCACTCTg                                                                                 <  2:8673813/93‑1 (MQ=255)

        aaGCGTGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTg                                                                                 <  1:7280749/93‑1 (MQ=255)

        aaGCGTGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTg                                                                                 >  2:5864143/1‑93 (MQ=255)

        aaGCGTGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTg                                                                                 <  2:2683816/93‑1 (MQ=255)

         aGCGTGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTgg                                                                                >  2:1426743/1‑93 (MQ=255)

             tGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTa                                                                            >  1:9293210/1‑93 (MQ=255)

              ggTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTAt                                                                           >  2:6932824/1‑93 (MQ=255)

              ggTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTAt                                                                           >  2:5724956/1‑93 (MQ=255)

               gTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATc                                                                          <  2:6639158/93‑1 (MQ=255)

                ttGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCCCCGCACTCTGGGCTATCg                                                                         <  1:8265041/93‑1 (MQ=255)

                 tGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGc                                                                        >  1:6089852/1‑93 (MQ=255)

                  gATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGcc                                                                       <  1:6413813/93‑1 (MQ=255)

                  gATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGcc                                                                       <  2:4450584/93‑1 (MQ=255)

                   aTGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGccc                                                                      <  1:400351/93‑1 (MQ=255)

                   aTGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGccc                                                                      >  2:7200809/1‑93 (MQ=255)

                    tggtggGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCg                                                                     <  2:4647941/93‑1 (MQ=255)

                    tggtggGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCg                                                                     >  1:7382309/1‑93 (MQ=255)

                       tggGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGtt                                                                  >  1:6266622/1‑93 (MQ=255)

                       tggGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGtt                                                                  >  1:3322294/1‑93 (MQ=255)

                        gggAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTg                                                                 >  2:1563336/1‑93 (MQ=255)

                           aaTAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACa                                                              <  1:3529322/93‑1 (MQ=255)

                            aTAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCACATCCTAATAGCACCGCACTCTGGGCTATACCCCGGTGGACac                                                             >  1:6774606/1‑92 (MQ=255)

                            aTAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACAt                                                             >  1:5160971/1‑93 (MQ=255)

                            aTAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACAt                                                             <  2:1869567/93‑1 (MQ=255)

                             tAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATg                                                            <  2:1257567/93‑1 (MQ=255)

                              aGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGc                                                           >  2:5983320/1‑93 (MQ=255)

                                ttCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGccc                                                         >  1:4369137/1‑93 (MQ=255)

                                ttCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGccc                                                         <  1:7842202/93‑1 (MQ=255)

                                 tgaGGATTGAGGCGTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGcccc                                                        <  1:2600399/91‑1 (MQ=255)

                                  cAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCcg                                                       <  1:6909847/93‑1 (MQ=255)

                                   aGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCcgc                                                      >  1:3682963/1‑93 (MQ=255)

                                   aGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCcgc                                                      <  2:8512608/93‑1 (MQ=255)

                                     gATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCgcgc                                                    >  2:8855326/1‑93 (MQ=255)

                                      aTTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGcc                                                   >  2:3777758/1‑93 (MQ=255)

                                      aTTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGcc                                                   >  1:2324978/1‑93 (MQ=255)

                                       ttGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCa                                                  >  1:3880190/1‑93 (MQ=255)

                                        tGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCGCCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCaa                                                 >  1:7601155/1‑93 (MQ=255)

                                        tGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCaa                                                 >  1:5077626/1‑93 (MQ=255)

                                         gAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAAt                                                >  2:3663914/1‑93 (MQ=255)

                                         gAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAAt                                                <  2:4245229/93‑1 (MQ=255)

                                          aGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATc                                               >  2:321921/1‑93 (MQ=255)

                                          aGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATc                                               >  2:4520177/1‑93 (MQ=255)

                                            gCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCaa                                             >  1:2210288/1‑93 (MQ=255)

                                             cTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAac                                            >  2:5371323/1‑93 (MQ=255)

                                              tttCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATTCCCCGCGCCAATCAaca                                           >  2:8140857/1‑93 (MQ=255)

                                              tttCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAaca                                           >  2:4745617/1‑93 (MQ=255)

                                              tttCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAaca                                           <  1:4710463/93‑1 (MQ=255)

                                               ttCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAacac                                          <  2:1030652/93‑1 (MQ=255)

                                                tCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACAcc                                         >  2:5231811/1‑93 (MQ=255)

                                                tCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACAcc                                         <  1:6296402/93‑1 (MQ=255)

                                                      tAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCCCCCGGGTGACATGCCAGGCGCCAATCAACACCACACgg                                   >  1:2937681/1‑92 (MQ=255)

                                                        aaTATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTgg                                 <  1:978659/93‑1 (MQ=255)

                                                         atatCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGa                                <  2:1685904/93‑1 (MQ=255)

                                                         atatCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGa                                >  2:1531305/1‑93 (MQ=255)

                                                          tatCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGAc                               >  1:1219763/1‑93 (MQ=255)

                                                            tCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGc                             >  2:7024108/1‑93 (MQ=255)

                                                             cAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCa                            <  2:9058688/93‑1 (MQ=255)

                                                              aGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCAt                           >  2:7471871/1‑93 (MQ=255)

                                                               gCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATa                          <  1:3595800/93‑1 (MQ=255)

                                                                cAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATAc                         <  2:7658237/93‑1 (MQ=255)

                                                                 aaCGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATAcc                        <  2:9393836/93‑1 (MQ=255)

                                                                 aaCGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATAcc                        <  2:7843808/93‑1 (MQ=255)

                                                                 aaCGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATAcc                        >  1:6554171/1‑93 (MQ=255)

                                                                  aCGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCa                       >  2:4153424/1‑93 (MQ=255)

                                                                   cGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCAt                      <  1:808883/93‑1 (MQ=255)

                                                                      cAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGa                   <  1:10008990/93‑1 (MQ=255)

                                                                       aaTATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGAc                  >  2:8292273/1‑93 (MQ=255)

                                                                       aaTATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGAc                  <  1:2026163/93‑1 (MQ=255)

                                                                        atatCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGACg                 >  1:1460749/1‑93 (MQ=255)

                                                                         tatCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGACGc                <  2:928922/93‑1 (MQ=255)

                                                                         tatCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGACGc                >  2:2536483/1‑93 (MQ=255)

                                                                           tCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGACGCTg              >  1:9016756/1‑93 (MQ=255)

                                                                            cAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGACGCTGa             >  2:9557503/1‑93 (MQ=255)

                                                                             aaaTCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGACGCTGAt            >  1:8457333/1‑93 (MQ=255)

                                                                              aaTCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGACGCTGATa           >  2:296881/1‑93 (MQ=255)

                                                                                 ccTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGACGCTGATAAAt        <  1:6991253/93‑1 (MQ=255)

                                                                                  cTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGACGCTGATAAATa       <  1:2907545/93‑1 (MQ=255)

                                                                                  cTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGACGCTGATAAATa       <  2:2909875/93‑1 (MQ=255)

                                                                                  cTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGACGCTGATAAATa       <  1:5484245/93‑1 (MQ=255)

                                                                                   tAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGACGCTGATAAATAc      >  2:5032774/1‑93 (MQ=255)

                                                                                   tAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGACGCTGATAAATAc      >  2:357204/1‑93 (MQ=255)

                                                                                    aaTAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGACGCTGATAAATACt     >  1:1563823/1‑93 (MQ=255)

                                                                                     aTAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGACGCTGATAAATACTg    <  2:7375202/93‑1 (MQ=255)

                                                                                     aTAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGACGCTGATAAATACTg    >  2:8041827/1‑93 (MQ=255)

                                                                                     aTAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGACGCTGATAAATACTg    <  1:8001680/93‑1 (MQ=255)

                                                                                      tAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGACGCTGATAAATACTGa   >  1:9131370/1‑93 (MQ=255)

                                                                                      tAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGACGCTGATAAATACTGa   >  1:4575605/1‑93 (MQ=255)

                                                                                       aGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGACGCTGATAAATACTGAc  <  1:7134641/93‑1 (MQ=255)

GCATGGAGAAGCGTGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCGCCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGACGCTGATAAATACTGAC  >  NC_000913/3011803‑3011982

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