BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,132,238	A→G	R479R (CGT→CGC)	b2060 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,132,238	0	A	G	99.7%	1069.6 / ‑5.1	311	R479R (CGT→CGC)	b2060	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (138/172); total (139/172)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

Only 100 of 336 total aligned reads displayed.

GCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCACGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCAGTACTTGCGGGCTTTCAATGCCGCGATT  >  NC_000913/2132149‑2132327

gCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCgcgc                                                                                        >  2:1500830/1‑93 (MQ=255)

gCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCgcgc                                                                                        <  1:1487976/93‑1 (MQ=255)

gCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCgcgc                                                                                        >  1:7133385/1‑93 (MQ=255)

gCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCgcgc                                                                                        >  1:7045490/1‑93 (MQ=255)

gCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCgcgc                                                                                        <  1:4849583/93‑1 (MQ=255)

    cAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCtt                                                                                    <  2:713641/93‑1 (MQ=255)

    cAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCtt                                                                                    <  1:7887397/93‑1 (MQ=255)

     aaTCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCACGCGCttt                                                                                   >  1:884533/1‑93 (MQ=255)

      aTCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTc                                                                                  >  1:1288542/1‑93 (MQ=255)

       tCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCt                                                                                 <  1:8238227/93‑1 (MQ=255)

         gCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGc                                                                               >  1:8714313/1‑93 (MQ=255)

          ccAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGcc                                                                              <  2:1872186/93‑1 (MQ=255)

           cAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCa                                                                             <  2:7563503/93‑1 (MQ=255)

           cAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCa                                                                             >  1:1543363/1‑93 (MQ=255)

            aGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCAt                                                                            <  1:5299985/93‑1 (MQ=255)

            aGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCAt                                                                            <  2:695506/93‑1 (MQ=255)

            aGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCAt                                                                            <  2:2812791/93‑1 (MQ=255)

             gATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCAtt                                                                           >  1:2647386/1‑93 (MQ=255)

               tCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTcc                                                                         <  1:5451289/93‑1 (MQ=255)

                cGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCg                                                                        <  1:7293591/93‑1 (MQ=255)

                  gTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGAc                                                                      >  2:6250707/1‑93 (MQ=255)

                     ggATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGc                                                                   >  1:3818593/1‑93 (MQ=255)

                     ggATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGc                                                                   >  1:4162477/1‑93 (MQ=255)

                      gATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCg                                                                  >  1:3112462/1‑93 (MQ=255)

                         tCCCCACCGCCAGTATCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGa                                                               <  1:1503230/93‑1 (MQ=255)

                         tCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGa                                                               >  1:3681013/1‑93 (MQ=255)

                         tCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGa                                                               <  2:1761500/93‑1 (MQ=255)

                         tCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGa                                                               <  2:1099927/93‑1 (MQ=255)

                          ccccACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGAt                                                              <  2:5778933/93‑1 (MQ=255)

                          ccccACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGAt                                                              <  2:2467755/93‑1 (MQ=255)

                            ccACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGc                                                            >  2:1241622/1‑93 (MQ=255)

                            ccACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGc                                                            <  2:2535589/93‑1 (MQ=255)

                             cACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCt                                                           >  1:2119449/1‑93 (MQ=255)

                               ccGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTgg                                                         <  2:5688207/93‑1 (MQ=255)

                                cGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGc                                                        >  1:3352629/1‑93 (MQ=255)

                                    aGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAg                                                    <  2:4165271/93‑1 (MQ=255)

                                    aGTAGCTGGCTCTGTTTAGAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAg                                                    >  2:8624533/1‑93 (MQ=255)

                                     gTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGa                                                   >  1:5875887/1‑93 (MQ=255)

                                     gTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGa                                                   >  2:3720867/1‑93 (MQ=255)

                                      tAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGAc                                                  <  1:2378057/93‑1 (MQ=255)

                                      tAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGAc                                                  <  2:5429792/93‑1 (MQ=255)

                                       aGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACg                                                 >  2:8050334/1‑93 (MQ=255)

                                       aGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACg                                                 <  1:8016942/93‑1 (MQ=255)

                                         cTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCt                                               >  1:792847/1‑93 (MQ=255)

                                          tGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTa                                              >  1:3241268/1‑93 (MQ=255)

                                           ggCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTaa                                             <  1:8561079/93‑1 (MQ=255)

                                           ggCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTaa                                             <  1:2000047/93‑1 (MQ=255)

                                            gCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAAt                                            >  2:3120115/1‑93 (MQ=255)

                                            gCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAAt                                            >  1:5545504/1‑93 (MQ=255)

                                            gCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAAt                                            >  1:5848004/1‑93 (MQ=255)

                                             ctctGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATa                                           >  2:5212910/1‑93 (MQ=255)

                                             ctctGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATa                                           <  2:4363644/93‑1 (MQ=255)

                                              tctGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATAc                                          <  1:395747/93‑1 (MQ=255)

                                               ctGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATAcc                                         <  2:1371977/93‑1 (MQ=255)

                                                tGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCg                                        <  2:5400726/93‑1 (MQ=255)

                                                 gTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCgt                                       >  1:608151/1‑93 (MQ=255)

                                                  tttATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCgtg                                      <  1:6000688/93‑1 (MQ=255)

                                                  tttATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCgtg                                      >  2:3757661/1‑93 (MQ=255)

                                                  tttATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCgtg                                      >  1:180850/1‑93 (MQ=255)

                                                    tataGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGtt                                    <  2:663327/93‑1 (MQ=255)

                                                    tataGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCACGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGtt                                    <  1:4994583/93‑1 (MQ=255)

                                                      taGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGttct                                  <  1:1001082/93‑1 (MQ=255)

                                                       aGCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGttctt                                 <  1:4393094/93‑1 (MQ=255)

                                                        gCGTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGttcttc                                >  1:3878488/1‑93 (MQ=255)

                                                          gTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCa                              >  1:2921548/1‑93 (MQ=255)

                                                          gTTTAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCa                              >  2:7719244/1‑93 (MQ=255)

                                                           tttAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCAg                             <  2:1474897/93‑1 (MQ=255)

                                                            ttAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCAGt                            <  2:3818593/93‑1 (MQ=255)

                                                            ttAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCAGt                            >  2:4020311/1‑93 (MQ=255)

                                                            ttAATCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCAGt                            <  2:7388871/93‑1 (MQ=255)

                                                              aaTCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCAGTAc                          >  2:5364554/1‑93 (MQ=255)

                                                               aTCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCAGTACt                         <  2:2241360/93‑1 (MQ=255)

                                                               aTCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCAGTACt                         <  1:7398060/93‑1 (MQ=255)

                                                               aTCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCAGTACt                         <  2:4995511/93‑1 (MQ=255)

                                                                tCCCTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCAGTACtt                        <  1:1350359/93‑1 (MQ=255)

                                                                 ccttttGATGTTTTTG‑CCCCATCGCGCGCCTTTCGCCCATCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCAGTACTTGc                      <  2:2213989/90‑1 (MQ=255)

                                                                   cTTTGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCAGTACTTGCg                     >  2:1985947/1‑93 (MQ=255)

                                                                    tttGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCAGTACTTGCgg                    >  2:92112/1‑93 (MQ=255)

                                                                    tttGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCAGTACTTGCgg                    <  2:7475745/93‑1 (MQ=255)

                                                                    tttGATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCAGTACTTGCgg                    >  1:682035/1‑93 (MQ=255)

                                                                       gATGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCAGTACTTGCGGGCt                 <  2:5125633/93‑1 (MQ=255)

                                                                        aTGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCAGTACTTGCGGGCtt                <  2:8768735/93‑1 (MQ=255)

                                                                        aTGGTTTTGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCAGTACTTGCGGGCtt                <  2:8515806/93‑1 (MQ=255)

                                                                          ggTTTTGCCGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCAGTACTTGCGGGCTTTc              <  1:7575217/93‑1 (MQ=255)

                                                                            ttttGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCAGTACTTGCGGGCTTTCaa            >  2:3008516/1‑93 (MQ=255)

                                                                            ttttGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGAAGCTAATACCGTGTTCTTACAGTACTTGCGGGCTTTCaa            >  2:4126355/1‑93 (MQ=255)

                                                                             tttGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCAGTACTTGCGGGCTTTCAAt           >  1:1703564/1‑93 (MQ=255)

                                                                             tttGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCAGTACTTGCGGGCTTTCAAt           >  1:6285668/1‑93 (MQ=255)

                                                                             tttGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCAGTACTTGCGGGCTTTCAAt           <  1:5116687/93‑1 (MQ=255)

                                                                             tttGACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCAGTACTTGCGGGCTTTCAAt           <  2:7858476/93‑1 (MQ=255)

                                                                                gACGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCAGTACTTGCGGGCTTTCAATGcc        >  1:6520757/1‑93 (MQ=255)

                                                                                  cGCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCAGTACTTGCGGGCTTTCAATGCcgc      <  2:5983066/93‑1 (MQ=255)

                                                                                   gCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCAGTACTTGCGGGCTTTCAATGCcgcg     >  2:4857145/1‑93 (MQ=255)

                                                                                   gCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCAGTACTTGCGGGCTTTCAATGCcgcg     >  1:3805298/1‑93 (MQ=255)

                                                                                   gCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCAGTACTTGCGGGCTTTCAATGCcgcg     >  1:4351378/1‑93 (MQ=255)

                                                                                   gCTATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCAGTACTTGCGGGCTTTCAATGCcgcg     >  2:735501/1‑93 (MQ=255)

                                                                                     tATCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCAGTACTTGCGGGCTTTCAATGCCGCGAt   <  2:335592/93‑1 (MQ=255)

                                                                                      aTCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCAGTACTTGCGGGCTTTCAATGCCGCGAtt  <  2:6873084/93‑1 (MQ=255)

                                                                                      aTCGCGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCAGTACTTGCGGGCTTTCAATGCCGCGAtt  <  1:2264855/93‑1 (MQ=255)

                                                                                      aTCGCGCGCTTTCTGACATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCAGTACTTGCGGGCTTTCAATGCCGCGAtt  <  1:7838797/93‑1 (MQ=255)

GCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCACGCGCTTTCTGCCATTCCGACAGCGGGATGCTGGCATAGACGCTAATACCGTGTTCTTCCAGTACTTGCGGGCTTTCAATGCCGCGATT  >  NC_000913/2132149‑2132327

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