BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,456,531	A→G	M82V (ATG→GTG)	b1394 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,456,531	0	A	G	98.9%	564.5 / NA	176	M82V (ATG→GTG)	b1394	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (78/96); minor base T (1/1); total (79/97)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 184 total aligned reads displayed.

TTACTTACAGGTGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCA  >  NC_000913/1456441‑1456620

  tcttacAGGTGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTc                                                                                       <  2:5634986/92‑1 (MQ=255)

  acttacAGGTGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTc                                                                                       >  1:3311315/1‑93 (MQ=255)

  acttacAGGTGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTc                                                                                       >  1:6163334/1‑93 (MQ=255)

  acttacAGGTGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTc                                                                                       >  1:4535272/1‑93 (MQ=255)

   cttacAGGTGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCa                                                                                      <  2:1015091/93‑1 (MQ=255)

    ttacAGGTGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGTCCCCGCACCGGATTTAGGGGTGTCAg                                                                                     <  1:3043867/93‑1 (MQ=255)

      acAGGTGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGtt                                                                                   <  1:1698701/93‑1 (MQ=255)

      acAGGTGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGtt                                                                                   <  1:4767476/93‑1 (MQ=255)

      acAGGTGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGAATTTAGGGGTGTCAGtt                                                                                   <  1:5774762/93‑1 (MQ=255)

         ggTGCCGGGGGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGGGTCGGTTGaa                                                                                >  2:3806640/1‑93 (MQ=255)

         ggTGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGTCCCGCACCGGATTTAGGGGTGTCAGTTGaa                                                                                >  2:776545/1‑93 (MQ=255)

         ggTGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGaa                                                                                >  2:3672338/1‑93 (MQ=255)

             ccGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGtt                                                                            <  2:4474774/93‑1 (MQ=255)

               gggCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTc                                                                          >  1:4376443/1‑93 (MQ=255)

                  cgcgGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCtat                                                                       <  2:3233437/93‑1 (MQ=255)

                  cgcgGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCtat                                                                       >  2:3589468/1‑93 (MQ=255)

                  cgcgAATTTTGTGCT‑GTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCtata                                                                      <  2:1957567/93‑1 (MQ=255)

                    cgGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATaa                                                                     <  2:822526/93‑1 (MQ=255)

                    cgGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATaa                                                                     >  1:5256581/1‑93 (MQ=255)

                    cgGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATaa                                                                     >  1:3993225/1‑93 (MQ=255)

                     ggATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAAc                                                                    <  1:1615779/93‑1 (MQ=255)

                     ggATTTTGTGCTGGTCAGGATCTGAACGATCGTAACGTTGATCCAACAGTCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAAc                                                                    <  2:3935194/93‑1 (MQ=255)

                       aTTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAAccc                                                                  >  2:5808450/1‑93 (MQ=255)

                       aTTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAAccc                                                                  <  2:6007877/93‑1 (MQ=255)

                        ttttGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCa                                                                 >  2:3759522/1‑93 (MQ=255)

                        ttttGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCa                                                                 >  1:1236415/1‑93 (MQ=255)

                        ttttGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCa                                                                 <  2:4537137/93‑1 (MQ=255)

                        ttttGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCa                                                                 <  1:1588636/93‑1 (MQ=255)

                        ttttGT‑CTGGGTAGGGTCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCAc                                                                <  1:3159009/93‑1 (MQ=255)

                         tttGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCAc                                                                <  1:6264704/93‑1 (MQ=255)

                         tttGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCAc                                                                >  1:3954223/1‑93 (MQ=255)

                           tgtgCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTg                                                              >  1:679018/1‑93 (MQ=255)

                            gtgCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTgg                                                             <  1:5419606/93‑1 (MQ=255)

                              gCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTa                                                           >  1:1361326/1‑93 (MQ=255)

                                 ggTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGt                                                        <  1:5604312/93‑1 (MQ=255)

                                 ggTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGt                                                        <  2:2466357/93‑1 (MQ=255)

                                 ggTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGt                                                        <  1:765308/93‑1 (MQ=255)

                                 ggTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGt                                                        <  1:2282231/93‑1 (MQ=255)

                                 ggTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGt                                                        <  1:2598536/93‑1 (MQ=255)

                                 ggTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGt                                                        <  2:117519/93‑1 (MQ=255)

                                 ggTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGt                                                        <  2:5012156/93‑1 (MQ=255)

                                 ggTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGAGTCAGATGAACGGTTCTATAACCCACTGTTAcac                                                        >  2:1743877/1‑91 (MQ=255)

                                  gTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTc                                                       >  2:4450286/1‑93 (MQ=255)

                                  gTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTc                                                       <  2:4379375/93‑1 (MQ=255)

                                    cAGGATCTTAACGATCGTAACGTTGATCCCACAGGACCCGCAGCGGATTTAGGGTTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGc                                                     >  1:5676787/1‑93 (MQ=255)

                                     aGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGcc                                                    <  1:3978820/93‑1 (MQ=255)

                                     aGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGcc                                                    >  2:5504452/1‑93 (MQ=255)

                                      ggATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCt                                                   <  1:91913/93‑1 (MQ=255)

                                      ggATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCt                                                   <  1:1043540/93‑1 (MQ=255)

                                      ggATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCt                                                   <  2:954878/93‑1 (MQ=255)

                                       gATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTg                                                  <  1:4633943/93‑1 (MQ=255)

                                         tCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGc                                                >  1:40680/1‑93 (MQ=255)

                                         tCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGc                                                >  1:1437304/1‑93 (MQ=255)

                                         tCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAAACCACTGGTACGTCGCCTGGc                                                >  2:4796168/1‑93 (MQ=255)

                                           ttAACGATCGTAATGTTGCTCCCACAGGCCCCGCACCGGATTTAAGGGCGTCAGTTGAATGATTCTATAACCCACTGGTACGTCGCCTGGcca                                              >  2:2104898/1‑91 (MQ=255)

                                            tAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCaaa                                             >  2:11461/1‑93 (MQ=255)

                                             aaCGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCaaaa                                            <  2:5614135/93‑1 (MQ=255)

                                                gATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTg                                         <  2:4163205/93‑1 (MQ=255)

                                                 aTCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGc                                        >  1:3441303/1‑93 (MQ=255)

                                                    gTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCaa                                     >  1:5438288/1‑93 (MQ=255)

                                                     tAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCaca                                    >  1:4422485/1‑91 (MQ=255)

                                                       aCGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAAc                                  <  2:544979/93‑1 (MQ=255)

                                                       aCGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAAc                                  <  1:4735181/93‑1 (MQ=255)

                                                       aCGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAAc                                  <  2:2932014/93‑1 (MQ=255)

                                                         gTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCg                                >  1:962220/1‑93 (MQ=255)

                                                         gTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCg                                >  2:4053398/1‑93 (MQ=255)

                                                          ttGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCgg                               >  1:1049876/1‑93 (MQ=255)

                                                          ttGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCgg                               >  2:2051258/1‑93 (MQ=255)

                                                          ttGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCgg                               >  2:5863416/1‑93 (MQ=255)

                                                          ttGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCgg                               >  1:1206691/1‑93 (MQ=255)

                                                           tGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGt                              <  2:1298957/93‑1 (MQ=255)

                                                              tCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGAt                           >  2:4939164/1‑93 (MQ=255)

                                                                ccACAGGCCCCGCACCGGATTTAGGGGTTTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTAATCt                         <  2:3980043/93‑1 (MQ=255)

                                                                ccACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCt                         <  2:3357321/93‑1 (MQ=255)

                                                                ccACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCt                         <  1:5054120/93‑1 (MQ=255)

                                                                 cacaGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCtg                        <  1:4835837/93‑1 (MQ=255)

                                                                 cacaGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCtg                        >  1:2924584/1‑93 (MQ=255)

                                                                 cacaGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCtg                        <  1:491834/93‑1 (MQ=255)

                                                                  acaGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCtgt                       >  2:4347097/1‑93 (MQ=255)

                                                                   caGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCtgtg                      >  2:3141553/1‑93 (MQ=255)

                                                                    aGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTAGCAAAACTGCCAAAACCGGTGATCTGTGc                     <  1:5797543/93‑1 (MQ=255)

                                                                     ggCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCa                    <  1:4371979/93‑1 (MQ=255)

                                                                     ggCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTGTCTGTAATCCACTGGTACGTCGCCTGGCAGAACGGACAAAACCGGTGATCTGTGCa                    >  2:4171206/1‑93 (MQ=255)

                                                                       ccccGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGt                  >  1:4476191/1‑93 (MQ=255)

                                                                           gCATCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAt              <  2:4122673/93‑1 (MQ=255)

                                                                           gCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAt              <  1:2150652/93‑1 (MQ=255)

                                                                           gCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAt              <  2:1211126/93‑1 (MQ=255)

                                                                           gCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAAACCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAt              >  2:2240879/1‑93 (MQ=255)

                                                                            cACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCTGTGATCTGTGCAGTCAAtg             <  2:6126876/93‑1 (MQ=255)

                                                                            cACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtg             >  2:3639573/1‑93 (MQ=255)

                                                                              ccGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtggc           >  2:5959528/1‑93 (MQ=255)

                                                                               cGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtggcg          >  2:2002937/1‑93 (MQ=255)

                                                                               cGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtggcg          <  2:5396199/93‑1 (MQ=255)

                                                                                ggATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtggcgt         <  1:3115283/93‑1 (MQ=255)

                                                                                ggATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtggcgt         >  1:2894189/1‑93 (MQ=255)

                                                                                ggATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtggcgt         <  1:3775277/93‑1 (MQ=255)

                                                                                ggATTTAGGGGTGTCAGTTGAACGTTACTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtggcgt         <  2:2705456/93‑1 (MQ=255)

                                                                                   tttAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtggcgtggc      <  2:177159/93‑1 (MQ=255)

                                                                                       ggggTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCggca  >  2:2592651/1‑93 (MQ=255)

                                                                                       ggggTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCggca  <  1:3406867/93‑1 (MQ=255)

TTACTTACAGGTGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCA  >  NC_000913/1456441‑1456620

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