BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,520,082	A→G	H527H (CAT→CAC)	yjhG ←	putative dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,520,082	0	A	G	99.7%	1263.2 / ‑5.1	389	H527H (CAT→CAC)	yjhG	putative dehydratase
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (193/195); total (193/196)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.22e-01

Only 100 of 398 total aligned reads displayed.

GGGCCTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACATGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCAATAATGACCAG  >  NC_000913/4519993‑4520174

gggCCTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGc                                                                                           <  2:8334817/93‑1 (MQ=255)

    cTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTAc                                                                                       <  1:328522/93‑1 (MQ=255)

     tCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTAcc                                                                                      <  2:4919650/93‑1 (MQ=255)

     tCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTAcc                                                                                      >  1:6045804/1‑93 (MQ=255)

     tCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTAcc                                                                                      >  2:2132784/1‑93 (MQ=255)

     tCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCAGTGATTAACGAAACGTGCTTAcc                                                                                      >  2:6710028/1‑93 (MQ=255)

        ggCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTa                                                                                   <  2:9259381/93‑1 (MQ=255)

        ggCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTa                                                                                   <  2:7779996/93‑1 (MQ=255)

        ggCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCCAGAAACGTGCATCGGTGATTAACGAAACGTGATTACCGTa                                                                                   >  1:6813961/1‑93 (MQ=255)

         gCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTAt                                                                                  <  1:9105234/93‑1 (MQ=255)

          ccAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATg                                                                                 <  1:1886180/93‑1 (MQ=255)

           cAGCGCTGCTGGCACCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGa                                                                                <  2:5400926/93‑1 (MQ=255)

             gcgcTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACa                                                                              >  2:7328523/1‑93 (MQ=255)

             gcgcTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACa                                                                              <  1:5838012/93‑1 (MQ=255)

              cgcTTCTGGACCCACAAGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAg                                                                             <  1:8971130/93‑1 (MQ=255)

                cTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGAt                                                                           >  1:7294078/1‑93 (MQ=255)

                 ttCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATg                                                                          >  1:10041458/1‑93 (MQ=255)

                  tCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGc                                                                         <  2:1415948/93‑1 (MQ=255)

                      gCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCa                                                                     <  2:8198043/93‑1 (MQ=255)

                      gCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCa                                                                     >  2:2988888/1‑93 (MQ=255)

                       ccccACATGGCCGATGCACGCGCCAGTAGAAACGCCTGAGAAACGTGCATCGGGGATTAACGAAAGGTGCTTACCGTATGAGAGATGCTTCAg                                                                    >  2:407829/1‑93 (MQ=255)

                       ccccACATGGCAGGTGCAGACGCCAGTAGAAACGCTCGAGACACGTGCATCGGTGATTAACGAAACGGGCTTACCGTATGTCAGATGCTTCAg                                                                    <  2:4009857/93‑1 (MQ=255)

                         ccACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAggg                                                                  <  2:5230302/93‑1 (MQ=255)

                         ccACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAggg                                                                  <  2:1035871/93‑1 (MQ=255)

                         ccACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAggg                                                                  >  1:720314/1‑93 (MQ=255)

                         ccACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAggg                                                                  >  1:7322355/1‑93 (MQ=255)

                          cacaTGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGc                                                                 <  2:7161847/93‑1 (MQ=255)

                           acaTGGCGGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCa                                                                <  1:374373/93‑1 (MQ=255)

                           acaTGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCa                                                                >  2:418005/1‑93 (MQ=255)

                            caTGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCAc                                                               <  2:3386971/93‑1 (MQ=255)

                              tGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTg                                                             >  2:7898354/1‑93 (MQ=255)

                              tGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTg                                                             >  1:2127552/1‑93 (MQ=255)

                              tGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTg                                                             >  2:946128/1‑93 (MQ=255)

                              tGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTg                                                             >  1:7859101/1‑93 (MQ=255)

                                gCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGt                                                           <  1:9390730/93‑1 (MQ=255)

                                gCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGt                                                           <  2:7081886/93‑1 (MQ=255)

                                 ccGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTa                                                          <  1:8676234/93‑1 (MQ=255)

                                  cGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTaa                                                         <  1:9897060/93‑1 (MQ=255)

                                   gATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAAc                                                        >  1:2330954/1‑93 (MQ=255)

                                   gATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAAc                                                        >  1:50296/1‑93 (MQ=255)

                                    aTGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAAcc                                                       <  2:8813736/93‑1 (MQ=255)

                                    aTGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAAcc                                                       <  2:850535/93‑1 (MQ=255)

                                     tGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGGATGACAGATGCTTCAGGGCACTGGTAAccg                                                      >  1:2590689/1‑92 (MQ=255)

                                      gCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTg                                                     >  1:4723053/1‑93 (MQ=255)

                                           cgcCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAgg                                                <  2:1985512/93‑1 (MQ=255)

                                            gcCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGt                                               >  1:6619203/1‑93 (MQ=255)

                                            gcCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGt                                               >  2:6109472/1‑93 (MQ=255)

                                            gcCAGTAGAAACGCACGAGAAACGCGCAACGGGGATTAAAGAAACGTGATTACCGAATGACAGATGCTTCAGGGCACTGGTAACCTGGGAGGt                                               >  2:4716740/1‑93 (MQ=255)

                                             ccAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGtt                                              >  2:7379169/1‑93 (MQ=255)

                                             ccAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGtt                                              <  2:179634/93‑1 (MQ=255)

                                              cAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTCGTAACCTGGTATGttt                                             >  1:5310017/1‑93 (MQ=255)

                                               aGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTttc                                            >  1:9580024/1‑93 (MQ=255)

                                                gTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTttct                                           >  1:1731761/1‑93 (MQ=255)

                                                 tAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTttctt                                          >  1:2879271/1‑93 (MQ=255)

                                                 tAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTttctt                                          >  2:354078/1‑93 (MQ=255)

                                                 tAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTttctt                                          <  1:3092195/93‑1 (MQ=255)

                                                  aGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTttcttc                                         <  1:1470635/93‑1 (MQ=255)

                                                  aGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTttcttc                                         <  1:5571747/93‑1 (MQ=255)

                                                   gAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTcc                                        >  2:1329663/1‑93 (MQ=255)

                                                     aaCGCCCGAGAAACGTGCGTCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCAt                                      >  2:9834171/1‑93 (MQ=255)

                                                     aaCGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCAt                                      >  1:3453276/1‑93 (MQ=255)

                                                      aCGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATc                                     >  2:6295773/1‑93 (MQ=255)

                                                        gCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATccc                                   >  2:3707089/1‑93 (MQ=255)

                                                        gCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATccc                                   <  2:281755/93‑1 (MQ=255)

                                                        gCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATccc                                   >  2:587685/1‑93 (MQ=255)

                                                        gCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATccc                                   <  2:2676106/93‑1 (MQ=255)

                                                        gCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATccc                                   <  1:8030940/93‑1 (MQ=255)

                                                         cccGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCt                                  <  1:4390347/93‑1 (MQ=255)

                                                          ccGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTg                                 >  1:9065090/1‑93 (MQ=255)

                                                             agaAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTAc                              <  1:4568365/93‑1 (MQ=255)

                                                             agaAACGTGCATCGGTGATTAACGAAAAGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTAc                              >  1:3926604/1‑93 (MQ=255)

                                                              gaAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTAcc                             <  2:6229899/93‑1 (MQ=255)

                                                              gaAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTAcc                             >  1:2194548/1‑93 (MQ=255)

                                                                aaCGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTg                           >  1:829887/1‑93 (MQ=255)

                                                                 aCGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGa                          <  2:6781091/93‑1 (MQ=255)

                                                                   gTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAg                        <  2:6977542/93‑1 (MQ=255)

                                                                      gcgcGGTGATGAAAGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTc                     <  1:12787/90‑1 (MQ=255)

                                                                        tctgTGATTAACGAAATGTGCTTACCGTATGACAGATGCTTCAGGGCACAGGTAACCTGGTAGGTTTCTTCCATCCCTGTTCCTGAAGGTCCa                   >  2:8965433/4‑93 (MQ=255)

                                                                        tCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCa                   >  2:6588288/1‑93 (MQ=255)

                                                                         cGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCaa                  <  1:4097088/93‑1 (MQ=255)

                                                                         cGGTGATTAACAAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCaa                  <  2:7424305/93‑1 (MQ=255)

                                                                           gTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAAcc                >  1:8604514/1‑92 (MQ=255)

                                                                           gTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACg                >  2:770703/1‑93 (MQ=255)

                                                                           gTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACg                >  2:5083174/1‑93 (MQ=255)

                                                                            tGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGc               >  2:3179166/1‑93 (MQ=255)

                                                                            tGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGc               >  2:9726586/1‑93 (MQ=255)

                                                                            tGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGc               >  2:8718433/1‑93 (MQ=255)

                                                                            tGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGc               >  2:6465131/1‑93 (MQ=255)

                                                                               ttAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCaa            >  2:3610271/1‑93 (MQ=255)

                                                                               ttAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCaa            <  2:2386301/93‑1 (MQ=255)

                                                                               ttAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCaa            >  2:9127924/1‑93 (MQ=255)

                                                                                   cGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCaataat        <  2:1525381/93‑1 (MQ=255)

                                                                                    gAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCAATAATg       <  1:6450764/93‑1 (MQ=255)

                                                                                    gAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCAATAATg       >  1:7523917/1‑93 (MQ=255)

                                                                                     aaaCGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCAATAATGa      <  2:382016/93‑1 (MQ=255)

                                                                                     aaaCGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCAATAATGa      <  2:4386968/93‑1 (MQ=255)

                                                                                     aaaCGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCAATAATGa      <  1:567961/93‑1 (MQ=255)

                                                                                     aaaCGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCAATAATGa      <  2:5704988/93‑1 (MQ=255)

                                                                                      aaCGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCAATAATGAc     >  1:5274113/1‑93 (MQ=255)

                                                                                      aaCGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGGAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCAACAATGAc     >  1:7670044/1‑93 (MQ=255)

GGGCCTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACATGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCAATAATGACCAG  >  NC_000913/4519993‑4520174

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