BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,594,566	A→G	G515G (GGT→GGC)	ydeK ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,594,566	0	A	G	99.6%	752.3 / ‑5.2	237	G515G (GGT→GGC)	ydeK	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (105/130); total (107/130)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01

Only 100 of 244 total aligned reads displayed.

ATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCA  >  NC_000913/1594476‑1594654

aTTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGcc                                                                                        >  1:6796358/1‑93 (MQ=255)

 ttCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCa                                                                                       >  2:6482881/1‑93 (MQ=255)

   ccAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCAtt                                                                                     >  2:643069/1‑93 (MQ=255)

   ccAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCAtt                                                                                     <  1:1251708/93‑1 (MQ=255)

    cAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTa                                                                                    >  1:722752/1‑93 (MQ=255)

     agagACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTAc                                                                                   >  1:1857168/1‑93 (MQ=255)

      gagaCTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACt                                                                                  >  1:1106705/1‑93 (MQ=255)

          cTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACa                                                                              <  1:1719002/93‑1 (MQ=255)

          cTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACa                                                                              >  2:883326/1‑93 (MQ=255)

          cTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACa                                                                              >  1:2754632/1‑93 (MQ=255)

          cTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACa                                                                              >  2:1996322/1‑93 (MQ=255)

           tGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAt                                                                             >  1:9223662/1‑93 (MQ=255)

           tGTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAt                                                                             >  1:4557041/1‑93 (MQ=255)

            gTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAtt                                                                            >  2:8247321/1‑93 (MQ=255)

            gTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAtt                                                                            >  2:782220/1‑93 (MQ=255)

            gTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAtt                                                                            <  1:6289524/93‑1 (MQ=255)

            gTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAtt                                                                            >  2:5006025/1‑93 (MQ=255)

            gTCCGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAtt                                                                            >  2:7168620/1‑93 (MQ=255)

              ccGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTa                                                                          <  1:4035488/93‑1 (MQ=255)

              ccGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTa                                                                          <  2:5683701/93‑1 (MQ=255)

               cGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGCTAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAtttcc                                                                         >  1:9071348/1‑91 (MQ=255)

               cGTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTAc                                                                         <  2:8673850/93‑1 (MQ=255)

                gTTGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACg                                                                        >  1:8521199/1‑93 (MQ=255)

                 ttGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGa                                                                       >  1:4489336/1‑93 (MQ=255)

                 ttGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGa                                                                       <  2:1715727/93‑1 (MQ=255)

                  tGAAATATTAACTATCCCCTTACCGGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGa                                                                       <  1:5096257/93‑1 (MQ=255)

                  tGAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGAt                                                                      >  1:9325485/1‑93 (MQ=255)

                   gAAATATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATa                                                                     >  2:2726253/1‑93 (MQ=255)

                    aaaTATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATAc                                                                    <  2:9204197/93‑1 (MQ=255)

                    aaaTATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATAc                                                                    <  1:1903368/93‑1 (MQ=255)

                    aaaTATTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATAc                                                                    <  2:9920/93‑1 (MQ=255)

                       tatTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGc                                                                 <  1:2353961/93‑1 (MQ=255)

                        atTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCt                                                                <  2:1353281/93‑1 (MQ=255)

                        atTAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGCTGGGGCCATTACTGACATTTACGATACAGCt                                                                >  1:96625/1‑93 (MQ=255)

                         ttAACTATCCCCTTACC‑GGGGGCCTTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTa                                                               <  2:3506518/93‑1 (MQ=255)

                          tAACTATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTAc                                                              <  1:1365963/93‑1 (MQ=255)

                              tATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCaaa                                                          >  1:7684956/1‑93 (MQ=255)

                              tATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCaaa                                                          >  1:1321642/1‑93 (MQ=255)

                              tATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCaaa                                                          >  2:9023413/1‑93 (MQ=255)

                              tATCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCaaa                                                          >  1:5476436/1‑93 (MQ=255)

                               aTCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCaaaa                                                         <  1:1187646/93‑1 (MQ=255)

                               aTCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCaaaa                                                         >  1:8493932/1‑93 (MQ=255)

                                tCCCCTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAAt                                                        >  2:3367325/1‑93 (MQ=255)

                                 ccccTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAAtt                                                       <  1:7897416/93‑1 (MQ=255)

                                    cTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTAcc                                                    <  2:7086861/93‑1 (MQ=255)

                                    cTTACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTAcc                                                    >  2:8355774/1‑93 (MQ=255)

                                     ttACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCg                                                   >  2:4538117/1‑93 (MQ=255)

                                     ttACC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCg                                                   >  2:36120/1‑93 (MQ=255)

                                       aCC‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGAc                                                 >  2:2258755/1‑93 (MQ=255)

                                        cc‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACa                                                <  1:5897583/93‑1 (MQ=255)

                                         c‑GGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACAt                                               <  2:2196681/93‑1 (MQ=255)

                                           gggggCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACAtt                                              >  2:4090610/4‑93 (MQ=255)

                                            ggggCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTg                                             >  1:3680472/3‑93 (MQ=255)

                                            ggggCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTg                                             >  1:4967753/3‑93 (MQ=255)

                                            ggggCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTg                                             >  1:7463353/3‑93 (MQ=255)

                                             gggCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGa                                            <  1:9894824/92‑1 (MQ=255)

                                              ggCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGAt                                           >  2:2455260/1‑93 (MQ=255)

                                               gCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATa                                          <  1:43037/93‑1 (MQ=255)

                                                cATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATAc                                         >  2:5757268/1‑93 (MQ=255)

                                                cATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATAc                                         <  1:8850835/93‑1 (MQ=255)

                                                  tttCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATAccc                                       >  2:3777583/1‑93 (MQ=255)

                                                  tttCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATAccc                                       >  2:4009693/1‑93 (MQ=255)

                                                   ttCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCca                                      >  1:1547922/1‑93 (MQ=255)

                                                      ctccGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAg                                   >  2:673677/1‑93 (MQ=255)

                                                      ctccGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAg                                   <  2:1795065/93‑1 (MQ=255)

                                                      ctccGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAg                                   >  1:3280163/1‑93 (MQ=255)

                                                       tccGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGt                                  >  1:5264604/1‑93 (MQ=255)

                                                        ccGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTa                                 >  2:192751/1‑93 (MQ=255)

                                                        ccGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTa                                 <  2:5409024/93‑1 (MQ=255)

                                                          gATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATc                               <  2:1466889/93‑1 (MQ=255)

                                                           aTAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCg                              >  1:5721076/1‑93 (MQ=255)

                                                               aTCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACg                          >  1:9543311/1‑93 (MQ=255)

                                                                tCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACga                         <  1:6630444/93‑1 (MQ=255)

                                                                tCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACga                         >  1:4710742/1‑93 (MQ=255)

                                                                tCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACga                         <  1:3918350/93‑1 (MQ=255)

                                                                tCCATAGCCCGTTGAATTAAGGGAGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACCGTATCGAACga                         >  2:7556559/1‑93 (MQ=255)

                                                                 ccATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACgat                        >  2:4377337/1‑93 (MQ=255)

                                                                  cATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACgatg                       <  2:8373112/93‑1 (MQ=255)

                                                                   aTAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACgatga                      <  2:412100/93‑1 (MQ=255)

                                                                   aTAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACgatga                      <  2:9409247/93‑1 (MQ=255)

                                                                   aTAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACgatga                      <  1:7746927/93‑1 (MQ=255)

                                                                   aTAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACgatga                      <  2:9845675/93‑1 (MQ=255)

                                                                    tAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACgatgat                     >  1:3104529/1‑93 (MQ=255)

                                                                       cccGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTc                  <  1:115346/93‑1 (MQ=255)

                                                                       cccGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTc                  <  1:5996492/93‑1 (MQ=255)

                                                                         cGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGc                >  2:687649/1‑93 (MQ=255)

                                                                         cGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGc                <  2:6131045/93‑1 (MQ=255)

                                                                             gAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcatcc            >  2:2723339/1‑93 (MQ=255)

                                                                              aaTTAATGGTGGCTCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAAAGATGATGTCGcatcca           >  2:4761927/1‑93 (MQ=255)

                                                                                ttAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcatccatc         >  2:6654745/1‑93 (MQ=255)

                                                                                 tAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCa        <  2:5949823/93‑1 (MQ=255)

                                                                                 tAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCa        <  2:3250029/93‑1 (MQ=255)

                                                                                 tAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCa        <  2:1219281/93‑1 (MQ=255)

                                                                                 tAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGACGATGTCGCATCCATCa        <  2:9555798/93‑1 (MQ=255)

                                                                                   aGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAgg      <  2:4145076/93‑1 (MQ=255)

                                                                                   aGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAgg      >  1:3810925/1‑93 (MQ=255)

                                                                                    gggTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGt     <  1:2179073/93‑1 (MQ=255)

                                                                                      gTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTc   <  1:9737298/93‑1 (MQ=255)

                                                                                      gTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTc   <  1:8612642/93‑1 (MQ=255)

                                                                                      gTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTc   <  1:8148695/93‑1 (MQ=255)

ATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACC‑GGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCA  >  NC_000913/1594476‑1594654

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: `‑`