BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,544,971	T→C	G259G (GGT→GGC)	gntT →	gluconate permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,544,971	0	T	C	99.7%	1051.2 / ‑5.0	332	G259G (GGT→GGC)	gntT	gluconate permease
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (159/172); total (159/173)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.05e-01

Only 100 of 348 total aligned reads displayed.

GGTGCCGGTAGTACTGATGGCGATGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGTGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATGGATCAGATTAACGA  >  NC_000913/3544881‑3545060

  ggCCGGTAGTACTGATGGCGATGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGAcc                                                                                       <  1:7739720/92‑1 (MQ=255)

   gCCGGTAGTACTGATGGCGATGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGAccc                                                                                      <  2:441230/93‑1 (MQ=255)

   ccgggTAGTACTGATTGCGATGCTTGCGATTGCCGAAATGATCCTGCCAAACGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGAccc                                                                                      <  2:2111621/90‑1 (MQ=255)

       gcagtaCGGATGGCGATGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTa                                                                                  <  1:7586971/91‑1 (MQ=255)

            tcTGATGGCGATGCGTGCGATTGCCAAAATGATCCTCCCGAAAGGTCACCCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGc                                                                             <  1:3397436/92‑1 (MQ=255)

             cTGATGGCGATGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCa                                                                            >  2:460997/1‑93 (MQ=255)

             cTGATGGCGATGCGTGCGATTGCCGAAATGATCCCGCCGAAAGGCCACGCTTTCCTGCCGGTAGCGGAGGTCCTCGGCGACCCGGTAATGGCa                                                                            >  2:3185283/1‑93 (MQ=255)

              tGATGGCGATGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCaa                                                                           >  2:3616987/1‑93 (MQ=255)

              tGATGGCGATGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACCCTTTCATGCCGGTAGCGGAGGTCCTCGGCGACCCGGTTATGGCaa                                                                           >  2:7734905/1‑93 (MQ=255)

              tGATGGCGATGCGTGCGATTGCAGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCaa                                                                           <  2:3884404/93‑1 (MQ=255)

              cgATGGCGATGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCaa                                                                           <  2:1088154/92‑1 (MQ=255)

                aTGGCGATGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACg                                                                         >  1:647279/1‑93 (MQ=255)

                aTGGCGATGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACg                                                                         <  2:4822932/93‑1 (MQ=255)

                 tGGCGATGCGTGCGATTGCCGGAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGc                                                                        >  2:9605994/1‑93 (MQ=255)

                 tGGCGATGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGc                                                                        <  1:7112993/93‑1 (MQ=255)

                 tGGCGATGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGc                                                                        >  2:1321000/1‑93 (MQ=255)

                     gATGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGAt                                                                    <  1:4446730/93‑1 (MQ=255)

                     gATGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGAt                                                                    <  1:3653904/93‑1 (MQ=255)

                      aTGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCCCGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACACGGTAATGGCAACGCTGGtt                                                                   >  2:6630456/1‑93 (MQ=255)

                      aTGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGAtt                                                                   >  2:5810840/1‑93 (MQ=255)

                        gcgtgcgATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGc                                                                 >  1:1142235/1‑93 (MQ=255)

                        gcgtgcgATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGc                                                                 <  1:3977219/93‑1 (MQ=255)

                         cgtgcgATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGcc                                                                <  2:287491/93‑1 (MQ=255)

                           ggcgATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGt                                                              <  2:9358101/92‑1 (MQ=255)

                              gATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCt                                                           >  2:634322/1‑93 (MQ=255)

                              gATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGACGTGCt                                                           >  1:2314405/1‑93 (MQ=255)

                                ttGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGa                                                         <  2:8186845/93‑1 (MQ=255)

                                ttGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGa                                                         >  1:9019786/1‑93 (MQ=255)

                                ttGCCGAAATGATCCTGCCGAAAGGACACGCTTTCCTGAAGCTAGCGGAGATCCTCAGCGACCCGGTAATGGCAACGCTCATTGCCGAGCTGa                                                         >  2:8637404/1‑93 (MQ=255)

                                 tGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGAt                                                        >  2:1602129/1‑93 (MQ=255)

                                  gCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGAtt                                                       >  1:1624771/1‑93 (MQ=255)

                                     gAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGTCAACGCTGATTGCCGCTCTGATTGCg                                                    >  2:8992812/1‑93 (MQ=255)

                                      aaaTGATCCTGCCGCAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGa                                                   <  1:8936105/93‑1 (MQ=255)

                                      aaaTGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGa                                                   <  1:9234448/93‑1 (MQ=255)

                                      aaaTGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGa                                                   <  2:3229647/93‑1 (MQ=255)

                                       aaTGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGAt                                                  <  2:9850699/93‑1 (MQ=255)

                                       aaTGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGAt                                                  >  1:5718403/1‑93 (MQ=255)

                                         tGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGt                                                >  1:916617/1‑93 (MQ=255)

                                          gATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGtt                                               <  1:5178499/93‑1 (MQ=255)

                                             ccTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCAc                                            <  2:4189885/93‑1 (MQ=255)

                                             ccTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCAc                                            >  2:7425562/1‑93 (MQ=255)

                                               tGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCt                                          >  1:5934754/1‑93 (MQ=255)

                                                 ccGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCttt                                        <  1:1563949/93‑1 (MQ=255)

                                                  caaaaGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTg                                       <  2:5455816/91‑1 (MQ=255)

                                                  cGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTg                                       <  1:13083/93‑1 (MQ=255)

                                                   gAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTgg                                      >  1:1865385/1‑93 (MQ=255)

                                                    aaaGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGt                                     >  2:3488985/1‑93 (MQ=255)

                                                     aaGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTc                                    <  1:5341445/93‑1 (MQ=255)

                                                     aaGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTc                                    <  1:953240/93‑1 (MQ=255)

                                                     aaGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTc                                    >  2:2252808/1‑93 (MQ=255)

                                                     aaGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTc                                    <  1:9858297/93‑1 (MQ=255)

                                                     aaGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTc                                    <  1:5087519/93‑1 (MQ=255)

                                                       ggTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTa                                  >  2:5359102/1‑93 (MQ=255)

                                                        gTCGCGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTaa                                 <  1:3028666/93‑1 (MQ=255)

                                                        gTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTaa                                 >  2:1878933/1‑93 (MQ=255)

                                                        gTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTaa                                 >  2:8444453/1‑93 (MQ=255)

                                                         tCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTaaa                                >  1:2757852/1‑93 (MQ=255)

                                                           aCGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAAcc                              >  1:7182916/1‑93 (MQ=255)

                                                           aCGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAAcc                              >  1:6917607/1‑93 (MQ=255)

                                                           aCGCTTTCCTGCCGGTAGCGGAGGTCCTCGGCGACCCGGTAATGGCAACGCTGCTTGCCCTGCTGATTGCGATGTTCACCTTTGGTCTACAcc                              >  1:4686432/1‑93 (MQ=255)

                                                            cGCTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCg                             >  2:6971087/1‑93 (MQ=255)

                                                             gCTGTCCTGCGGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGt                            <  1:4198312/93‑1 (MQ=255)

                                                             ccTTTCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGCTTGCGATGTTCACCTTTGGTCTAAACCGt                            <  2:6216251/92‑1 (MQ=255)

                                                               tttCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTgt                          >  2:8338840/1‑92 (MQ=255)

                                                               tttCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTgg                          <  1:2228928/93‑1 (MQ=255)

                                                               tttCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTgg                          >  1:6209487/1‑93 (MQ=255)

                                                               tttCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTgg                          <  1:8358647/93‑1 (MQ=255)

                                                                tgccTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGt                         <  1:9737058/91‑1 (MQ=255)

                                                                 tgcTGGCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGCCAACGCGGAGTGCCGTGTTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTc                        <  1:8631760/91‑1 (MQ=255)

                                                                 tCCTGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTc                        <  1:7435727/93‑1 (MQ=255)

                                                                 tCCTGCCGGTAGCGGAGCTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAGCCGTGGTc                        >  1:8789153/1‑93 (MQ=255)

                                                                    tGCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGtt                     <  2:5724054/93‑1 (MQ=255)

                                                                    cgCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGtt                     <  2:3003067/92‑1 (MQ=255)

                                                                     gCCGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTc                    <  1:7792481/93‑1 (MQ=255)

                                                                       cGGTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCaa                  <  1:1074083/93‑1 (MQ=255)

                                                                        ggTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAAt                 <  1:8468042/93‑1 (MQ=255)

                                                                         gTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATg                <  2:7978357/93‑1 (MQ=255)

                                                                         gTAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATg                >  2:575473/1‑93 (MQ=255)

                                                                          tAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATgg               <  2:5178484/93‑1 (MQ=255)

                                                                          tAGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATgg               <  2:4009647/93‑1 (MQ=255)

                                                                           aGCGTAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATGGa              <  2:6496637/93‑1 (MQ=255)

                                                                           aGCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATGGa              <  1:2483317/93‑1 (MQ=255)

                                                                            gCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATGGAt             <  1:7372091/93‑1 (MQ=255)

                                                                            gCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATGGAt             >  2:4120855/1‑93 (MQ=255)

                                                                            gCGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATGGAt             >  1:7014440/1‑93 (MQ=255)

                                                                             cGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATGGATc            <  2:3710358/93‑1 (MQ=255)

                                                                             cGGAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATGGATc            >  1:5844387/1‑93 (MQ=255)

                                                                               gAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATGGATCAg          >  2:3166568/1‑93 (MQ=255)

                                                                               gAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATGGATCAg          <  1:5560741/93‑1 (MQ=255)

                                                                               gAGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATGGATCAg          >  2:2040961/1‑93 (MQ=255)

                                                                                aGTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATGGATCAGa         <  2:9751052/93‑1 (MQ=255)

                                                                                 gTTCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATGGATCAGAt        >  2:9497453/1‑93 (MQ=255)

                                                                                  ttCCTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATGGATCAGAtt       <  1:7272689/93‑1 (MQ=255)

                                                                                    ccTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATGGATCAGATTaa     >  1:4653819/1‑93 (MQ=255)

                                                                                    ccTCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATGGATCAGATTaa     >  1:9413521/1‑93 (MQ=255)

                                                                                      tCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATGGATCAGATTAACg   <  1:3777288/93‑1 (MQ=255)

                                                                                      tCGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATGGATCAGATTAACg   >  1:8069140/1‑93 (MQ=255)

                                                                                      ccGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATGGATCAGATTAACg   <  2:2400259/92‑1 (MQ=255)

                                                                                       cGGCGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATGGATCAGATTAACGa  <  1:4380796/93‑1 (MQ=255)

                                                                                       cGGCGACCCGGTAATGGCAACCCTGATTGCCGTTCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATGGATCAGATTAACGa  <  1:3350789/93‑1 (MQ=255)

GGTGCCGGTAGTACTGATGGCGATGCGTGCGATTGCCGAAATGATCCTGCCGAAAGGTCACGCTTTCCTGCCGGTAGCGGAGTTCCTCGGTGACCCGGTAATGGCAACGCTGATTGCCGTGCTGATTGCGATGTTCACCTTTGGTCTAAACCGTGGTCGTTCAATGGATCAGATTAACGA  >  NC_000913/3544881‑3545060

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