BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,465,851	T→C	Y547C (TAC→TGC)	yheB ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,465,851	0	T	C	82.1%	595.0 / 125.7	251	Y547C (TAC→TGC)	yheB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (21/24); new base C (107/99); total (128/123)
Fisher's exact test for biased strand distribution p-value = 6.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

Only 100 of 254 total aligned reads displayed.

CTACGCGCTCTTTGGTTCCGGTCCAGTCAACAAAAACACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGTAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTGCGGCGCATTATTCAGCTCTTCGTTGCTATA  >  NC_000913/3465759‑3465942

cTACGCGCTCTTTGGTTCCGGTCCAGTCAACAAAAACACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGt                                                                                             <  2:6298875/93‑1 (MQ=255)

 tACGCGCTCTTTGGTTCCGGTCCAGTCAACAAAAACACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGTa                                                                                            <  2:5015705/93‑1 (MQ=255)

  aCGCGCTCTTTGGTTCCGGTCCAGTCAACAAAAACACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAg                                                                                           >  1:3835769/1‑93 (MQ=255)

  aCGCGCTCTTTGGTTCCGGTCCAGTCAACAAAAACACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAg                                                                                           <  2:2386204/93‑1 (MQ=255)

   cgcgCTCTTTGGTTCCGGTCCAGTCAACAAAAACACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGt                                                                                          <  1:4979804/93‑1 (MQ=255)

   cgcgCTCTTTGGTTCCGGTCCAGTCAACAAAAACACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGt                                                                                          <  1:1059579/93‑1 (MQ=255)

    gcgcTCTTTGGTTCCGGTCCAGTCAACAAAAACACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGtt                                                                                         >  1:2258454/1‑93 (MQ=255)

      gcTCTTTGGTTCCGGTCCAGTCAACAAAAACACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTca                                                                                       <  2:4413858/93‑1 (MQ=255)

      gcTCTTTGGTTCCGGTCCAGTCAACAAAAACACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTca                                                                                       <  2:4908604/93‑1 (MQ=255)

       ctctTTGGTTCCGGCCCAGTCAACAAAAACACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGTAGTTcac                                                                                      <  1:4581231/93‑1 (MQ=255)

        tttttGGTCCCGGTCCAGTCAACAAAAACACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTcaca                                                                                     <  1:29426/91‑1 (MQ=255)

        tctTTGGTTCCGGTCCAGTCAACAAAAACACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGTAGTTcaca                                                                                     >  1:6107741/1‑93 (MQ=255)

          tttGGTTCCGGTCCAGTCAACAAAAACACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACAcc                                                                                   <  2:1034640/93‑1 (MQ=255)

              gTTCCGGTCCAGTCAACAAAAACACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTa                                                                               >  1:4992750/1‑93 (MQ=255)

                tCCGGTCCAGTCAACAAAAACACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTaaa                                                                             >  1:4231972/1‑93 (MQ=255)

                 ccGGTCCAGTCAACAAAAACACGCCACGGGTTGCTGTCAGAAGGAGCAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAg                                                                            >  1:1000707/1‑93 (MQ=255)

                 ccGGTCCAGTCAACAAAAACACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAg                                                                            >  2:4019844/1‑93 (MQ=255)

                  cGGTCCAGTCAACAAAAACACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGc                                                                           <  1:788683/93‑1 (MQ=255)

                   ggTCCAGTCAACAAAAACACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCg                                                                          <  2:2747226/93‑1 (MQ=255)

                     tCCAGTCAACAAAAACACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGTAGTTCACACCGTTAAAGCGaa                                                                        <  2:6613540/93‑1 (MQ=255)

                       cAGTCAACAAAAACACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATc                                                                      <  1:2465080/93‑1 (MQ=255)

                        aGTCAACAAAAACACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCa                                                                     >  2:4223952/1‑93 (MQ=255)

                        aGTCAACAAAAACACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCa                                                                     <  2:4061558/93‑1 (MQ=255)

                         gTCAACAAAAACACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAg                                                                    >  1:3390076/1‑93 (MQ=255)

                          tCAACAAAAACACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGc                                                                   >  2:5646662/1‑93 (MQ=255)

                          tCAACAAAAACACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGc                                                                   >  2:534135/1‑93 (MQ=255)

                          tCAACAAAAACACGCCACGGGTTGCCGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGCTAAAGCGAATCAGc                                                                   >  2:1658061/1‑93 (MQ=255)

                           caacaaAAACACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCg                                                                  <  2:5136013/93‑1 (MQ=255)

                            aacaaAAACACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGt                                                                 >  2:4338029/1‑93 (MQ=255)

                            aacaaAAACACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGt                                                                 >  2:847710/1‑93 (MQ=255)

                             acaaAAACACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTa                                                                <  1:1668737/93‑1 (MQ=255)

                              caaAAACACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTAGTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTAt                                                               <  2:3557717/93‑1 (MQ=255)

                               aaaaaCACGCCACGGGTTGCTGTCAGACGGAGAAACATTCTGCACTTAGCTATGAGAAAAGCAGTTCACACAGTTAAAGCGAATCAGCGTATc                                                              <  1:6645054/93‑1 (MQ=255)

                               aaaaaCACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATc                                                              <  1:3705051/93‑1 (MQ=255)

                               aaaaaCACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAACCAGCGTATc                                                              <  2:6045265/93‑1 (MQ=255)

                                aaaaCACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGTAGTTCACACCGTTAAAGCGAATCAGCGTATCg                                                             >  1:6029305/1‑93 (MQ=255)

                                 aaaCACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGc                                                            <  2:1542431/93‑1 (MQ=255)

                                 aaaCACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGc                                                            <  1:5245937/93‑1 (MQ=255)

                                 aaaCACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGc                                                            >  1:6187215/1‑93 (MQ=255)

                                  aaCACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGTAGTTCACACCGTTAAAGCGAATCAGCGTATCGCt                                                           <  1:6546572/93‑1 (MQ=255)

                                   acacGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGTAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTg                                                          <  1:1691128/93‑1 (MQ=255)

                                     acGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGc                                                        >  2:204539/1‑93 (MQ=255)

                                     acGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGc                                                        <  1:3067412/93‑1 (MQ=255)

                                      cGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGTAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCa                                                       >  2:3899658/1‑93 (MQ=255)

                                        ccACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGTAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATa                                                     <  1:3321386/93‑1 (MQ=255)

                                        ccACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATa                                                     >  2:74734/1‑93 (MQ=255)

                                         cACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATaa                                                    <  1:2140920/93‑1 (MQ=255)

                                           cGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAg                                                  <  2:257743/93‑1 (MQ=255)

                                             ggTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCg                                                >  2:3933352/1‑93 (MQ=255)

                                                 gCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTc                                            <  2:5555678/93‑1 (MQ=255)

                                                   tGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTTGCATAAAGCGGTTCTg                                          >  1:589527/1‑93 (MQ=255)

                                                    gTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTgg                                         >  2:6017751/1‑93 (MQ=255)

                                                     tCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGa                                        <  1:5988867/93‑1 (MQ=255)

                                                      cAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGAt                                       <  1:2978160/93‑1 (MQ=255)

                                                      cAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGAt                                       <  2:705264/93‑1 (MQ=255)

                                                       aGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGAtt                                      >  1:1033512/1‑93 (MQ=255)

                                                        gAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTa                                     <  2:1622990/93‑1 (MQ=255)

                                                         aaGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTaa                                    >  1:1247586/1‑93 (MQ=255)

                                                          aGGAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTaaa                                   <  2:361220/93‑1 (MQ=255)

                                                           ggAGAGACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAAc                                  >  2:2558232/1‑93 (MQ=255)

                                                           ggAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAAc                                  <  1:3233697/93‑1 (MQ=255)

                                                           ggAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAAc                                  <  1:3019858/93‑1 (MQ=255)

                                                           ggAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAAc                                  <  2:6135650/93‑1 (MQ=255)

                                                           ggAGAAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAAc                                  >  1:3988561/1‑93 (MQ=255)

                                                            gagaAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACt                                 <  1:4255744/93‑1 (MQ=255)

                                                             agaAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTg                                >  2:1279345/1‑93 (MQ=255)

                                                              gaAACTTTCTGCACTTTACTCTGAGAAATGTAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTgc                               >  1:5499098/1‑93 (MQ=255)

                                                              gaAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTgc                               >  2:5163546/1‑93 (MQ=255)

                                                              gaAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTgc                               >  1:1883077/1‑93 (MQ=255)

                                                              gaAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTgc                               <  1:942801/93‑1 (MQ=255)

                                                              gaAACTTTCTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTgc                               >  1:5548295/1‑93 (MQ=255)

                                                                 aCTTTCTGCACTTTACTCTGAGAAATGTAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTgcggc                            >  1:1400439/1‑93 (MQ=255)

                                                                   tttCTGCACTTTACTCTGAGAAATGTAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTGCGgcgc                          <  1:5076262/93‑1 (MQ=255)

                                                                      cTGCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTGCGGCGCatt                       <  2:3117712/93‑1 (MQ=255)

                                                                        gCACTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTGCGGCGCattat                     <  1:2769304/93‑1 (MQ=255)

                                                                          aCTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTGCGGCGCATTATTc                   <  1:4854565/93‑1 (MQ=255)

                                                                          aCTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTGCGGCGCATTATTc                   >  1:252385/1‑93 (MQ=255)

                                                                          aCTTTACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTGCGGCGCATTATTc                   >  2:1168068/1‑93 (MQ=255)

                                                                           cTTTACTCTGAGAAATGTAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTGCGGCGCATTATTCa                  <  1:5708554/93‑1 (MQ=255)

                                                                            tttACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTGCGGCGCATTATTCAg                 >  2:1368959/1‑93 (MQ=255)

                                                                            tttACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTGCGGCGCATTATTCAg                 >  2:6557544/1‑93 (MQ=255)

                                                                             ttACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTGCGGCGCATTATTCAGc                <  2:1783211/93‑1 (MQ=255)

                                                                              tACTCTGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTGCGGCGCATTATTCAGct               <  2:5468028/93‑1 (MQ=255)

                                                                               aCTCTGAGAAATGTAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTGCGGCGCATTATTCAGctc              >  1:2722064/1‑93 (MQ=255)

                                                                                ctctGAGACATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTGCGGCGCATTATTCAGctct             >  1:6452773/1‑93 (MQ=255)

                                                                                ctctGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTGCGGCGCATTATTCAGctct             <  1:2354050/93‑1 (MQ=255)

                                                                                ctctGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTGCGGCGCATTATTCAGctct             >  1:4912161/1‑93 (MQ=255)

                                                                                 tctGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTGCGGCGCATTATTCAGCTCtt            >  1:1465859/1‑93 (MQ=255)

                                                                                 tctGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTGCGGCGCATTATTCAGCTCtt            <  1:5178809/93‑1 (MQ=255)

                                                                                  ctGAGAAATGTAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTGCGGCGCATTATTCAGCTCTTc           <  2:4320872/93‑1 (MQ=255)

                                                                                   tGAGAAATGTAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTGCGGCGCATTATTCAGCTCTTCg          >  2:329924/1‑93 (MQ=255)

                                                                                   tGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTGCGGCGCATTATTCAGCTCTTCg          >  2:579073/1‑93 (MQ=255)

                                                                                   tGAGAAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTGCGGCGCATTATTCAGCTCTTCg          >  1:5660402/1‑93 (MQ=255)

                                                                                    gagaAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTGCGGCGCATTATTCAGCTCTTCGt         >  2:52380/1‑93 (MQ=255)

                                                                                    gagaAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTGCGGCGCATTATTCAGCTCTTCGt         <  1:4594338/93‑1 (MQ=255)

                                                                                      gaAATGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTGCGGCGCATTATTCAGCTCTTCGTTg       <  2:5820144/93‑1 (MQ=255)

                                                                                        aaTGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTGCGGCGCATTATTCAGCTCTTCGTTGCt     <  2:4552288/93‑1 (MQ=255)

                                                                                        aaTGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTGCGGCGCATTATTCAGCTCTTCGTTGCt     >  1:3782638/1‑93 (MQ=255)

                                                                                        aaTGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTGCGGCGCATTATTCAGCTCTTCGTTGCt     >  1:3597431/1‑93 (MQ=255)

                                                                                         aTGCAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTGCGGCGCATTATTCAGCTCTTCGTTGCta    <  1:702034/93‑1 (MQ=255)

CTACGCGCTCTTTGGTTCCGGTCCAGTCAACAAAAACACGCCACGGGTTGCTGTCAGAAGGAGAAACTTTCTGCACTTTACTCTGAGAAATGTAGTTCACACCGTTAAAGCGAATCAGCGTATCGCTGGCATAAAGCGTTTCTGGATTAAACTGCGGCGCATTATTCAGCTCTTCGTTGCTATA  >  NC_000913/3465759‑3465942

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