BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	161,581	C→T	E151K (GAG→AAG)	yadP ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	161,581	0	C	T	99.4%	1228.9 / ‑5.9	335	E151K (GAG→AAG)	yadP	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (1/1); new base T (152/181); total (153/182)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 347 total aligned reads displayed.

ACAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTCCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACG  >  NC_000913/161491‑161670

 aaaTCCTTGTTATTGCGTTAGCGCCCAGCGATTTAGCGGCGTGTAACGTGTGCGTGCACGGGCAAACGAGGAGGCGGAAACGGCGAACTTCGt                                                                                        >  2:6076654/2‑93 (MQ=255)

  aaTCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTc                                                                                       <  2:6607366/93‑1 (MQ=255)

  aaTCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTc                                                                                       >  2:690147/1‑93 (MQ=255)

    tCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCAc                                                                                     >  2:5311267/1‑93 (MQ=255)

    tCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCAc                                                                                     >  2:5536679/1‑93 (MQ=255)

     ccTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCAcc                                                                                    >  1:6124565/1‑93 (MQ=255)

       ttGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGc                                                                                  >  1:6805135/1‑93 (MQ=255)

       ttGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGc                                                                                  >  2:7196410/1‑93 (MQ=255)

       ttGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGCGCAAACGAGGAGGCGGAAAGGGGGAACTCCGTCACCGc                                                                                  >  2:5917379/1‑93 (MQ=255)

         tttATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCAt                                                                                <  2:518971/92‑1 (MQ=255)

          ttATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATa                                                                               >  2:1672718/1‑93 (MQ=255)

          ttATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATa                                                                               >  2:2261963/1‑93 (MQ=255)

           tATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACCAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATAc                                                                              >  1:9767856/1‑93 (MQ=255)

            tttGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACg                                                                             <  1:7310646/92‑1 (MQ=255)

             ttGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGa                                                                            <  1:9312295/93‑1 (MQ=255)

             ttGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGa                                                                            <  2:8499415/93‑1 (MQ=255)

              tGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGAc                                                                           <  2:6335122/93‑1 (MQ=255)

                 gTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAg                                                                        <  1:3405827/93‑1 (MQ=255)

                 gTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAg                                                                        <  2:9675661/93‑1 (MQ=255)

                  ttAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGt                                                                       >  1:1755471/1‑93 (MQ=255)

                  ttAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGt                                                                       >  1:4953163/1‑93 (MQ=255)

                   tAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGtt                                                                      <  2:4038282/93‑1 (MQ=255)

                       gcCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGGCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTaaaa                                                                  <  2:2425999/93‑1 (MQ=255)

                        cccAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAAc                                                                 >  1:7743046/1‑93 (MQ=255)

                         ccAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAAcc                                                                <  2:4277540/93‑1 (MQ=255)

                          cAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCt                                                               <  1:7052167/93‑1 (MQ=255)

                           aGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTg                                                              <  2:6130431/93‑1 (MQ=255)

                           aGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTg                                                              <  1:9990316/93‑1 (MQ=255)

                           aGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTg                                                              >  2:4472091/1‑93 (MQ=255)

                           aGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTg                                                              <  1:6675353/93‑1 (MQ=255)

                           aGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTg                                                              >  1:459437/1‑93 (MQ=255)

                            gCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTgg                                                             >  1:10012092/1‑93 (MQ=255)

                             cGTTTTAGCGGCGTGTAGCGTGTGCGTCCATGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTggc                                                            <  2:6732625/93‑1 (MQ=255)

                               ttttAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTggcgg                                                          >  1:392585/1‑93 (MQ=255)

                               ttttAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTggcgg                                                          <  2:9661044/93‑1 (MQ=255)

                                tttAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTggcggc                                                         <  1:8444261/93‑1 (MQ=255)

                                tttAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTggcggc                                                         <  1:7911729/93‑1 (MQ=255)

                                    gcggcgTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGa                                                     >  2:411742/1‑93 (MQ=255)

                                    gcggcgTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGa                                                     <  1:2812059/93‑1 (MQ=255)

                                     cggcgTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGAt                                                    <  1:8821465/93‑1 (MQ=255)

                                       gcgTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTg                                                  >  2:1054797/1‑93 (MQ=255)

                                         gtgtAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTc                                                >  2:679192/1‑93 (MQ=255)

                                           gtAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCAc                                              <  1:5925846/93‑1 (MQ=255)

                                           gtAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCAc                                              <  2:2172151/93‑1 (MQ=255)

                                           gtAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCAc                                              <  2:1635973/93‑1 (MQ=255)

                                           gtAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCGGGCGGCGGGATTGTCAc                                              >  1:8692977/1‑93 (MQ=255)

                                                gtgtGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCt                                         <  2:1236873/93‑1 (MQ=255)

                                                gtgtGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCt                                         <  1:9576187/93‑1 (MQ=255)

                                                gtgtGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCt                                         <  2:5982733/93‑1 (MQ=255)

                                                 tgtgCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTc                                        >  1:6536086/1‑93 (MQ=255)

                                                 tgtgCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTc                                        >  2:6952944/1‑93 (MQ=255)

                                                 tgtgCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTc                                        >  2:4523443/1‑93 (MQ=255)

                                                   tgCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGc                                      >  1:7482322/1‑93 (MQ=255)

                                                    gCGTCCACGGGCAATCGAGGAGGCGTAACGGGTGAACTTCGTCACCGAATCCGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCt                                     <  2:1877009/93‑1 (MQ=255)

                                                    gCGTCCACGGGCAAACGAGGAGGCTTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCt                                     <  2:3111773/93‑1 (MQ=255)

                                                       gccACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGc                                  <  1:5844893/92‑1 (MQ=255)

                                                         cACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGt                                >  2:857430/1‑93 (MQ=255)

                                                         cACGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGt                                >  1:9052277/1‑93 (MQ=255)

                                                          aCGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTc                               <  2:5956494/93‑1 (MQ=255)

                                                          aCGGGCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTc                               <  2:3297374/93‑1 (MQ=255)

                                                            gggCAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTcgc                             >  1:7224402/1‑93 (MQ=255)

                                                               cAAACGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCa                          <  1:5617667/93‑1 (MQ=255)

                                                                aaaCGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCaa                         <  1:1435332/93‑1 (MQ=255)

                                                                 aaCGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAAt                        >  1:1495081/1‑93 (MQ=255)

                                                                 aaCGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAAt                        >  1:2158099/1‑93 (MQ=255)

                                                                  aCGAGGAGGCGTAAAGGGGTAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATa                       <  2:3542293/93‑1 (MQ=255)

                                                                   cGAGGAGGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATaa                      <  2:8805387/93‑1 (MQ=255)

                                                                    gaggagGCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAg                     >  1:1131525/1‑93 (MQ=255)

                                                                     aggaggCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAgg                    >  1:1337246/1‑93 (MQ=255)

                                                                      ggaggCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAgga                   >  1:319604/1‑92 (MQ=255)

                                                                      ggaggCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAgga                   >  2:2734875/1‑92 (MQ=255)

                                                                      ggagcCGTAAAGGATGAACTTCGTCACCGCATACGACCAGTTAAATCCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGt                   <  1:2512567/93‑1 (MQ=255)

                                                                       gaggCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTa                  >  2:3640366/1‑93 (MQ=255)

                                                                        aggCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTaa                 <  1:2939024/93‑1 (MQ=255)

                                                                        aggCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTaa                 <  1:1607610/93‑1 (MQ=255)

                                                                        aggCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGGaa                 >  2:8883205/1‑93 (MQ=255)

                                                                         ggCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAat                >  2:8167291/1‑93 (MQ=255)

                                                                          gCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAata               <  2:176915/93‑1 (MQ=255)

                                                                          gCGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAata               <  1:8045920/93‑1 (MQ=255)

                                                                           cGTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAatat              <  2:9814052/93‑1 (MQ=255)

                                                                            gTAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATAtg             <  1:3527928/93‑1 (MQ=255)

                                                                             tAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATAtgt            >  2:6176877/1‑93 (MQ=255)

                                                                             tAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATAtgt            >  2:3378904/1‑93 (MQ=255)

                                                                             tAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATAtgt            <  2:4986376/93‑1 (MQ=255)

                                                                             tAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATAtgt            >  2:7347490/1‑93 (MQ=255)

                                                                             tAAAGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATAtgt            >  1:4026685/1‑93 (MQ=255)

                                                                              aaaGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATAtgtg           <  2:9847222/93‑1 (MQ=255)

                                                                              aaaGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATAtgtg           <  1:5238663/93‑1 (MQ=255)

                                                                              aaaGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATAtgtg           <  1:8206783/93‑1 (MQ=255)

                                                                               aaGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTgg          <  2:3730661/93‑1 (MQ=255)

                                                                               aaGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTgg          >  2:6438394/1‑93 (MQ=255)

                                                                               aaGGGTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTgg          >  2:9038077/1‑93 (MQ=255)

                                                                                 gggTGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGAt        >  2:1444551/1‑93 (MQ=255)

                                                                                  ggTGAACTTCGTCACCGCATACGACCAGTTAAACCCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATg       <  2:4912807/93‑1 (MQ=255)

                                                                                    tGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGaa     >  1:7406171/1‑93 (MQ=255)

                                                                                    tGAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGaa     <  1:8629831/93‑1 (MQ=255)

                                                                                     gAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGaaa    <  2:5689013/93‑1 (MQ=255)

                                                                                     gAACTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGGAATATGTGGATGaaa    >  1:4653893/1‑93 (MQ=255)

                                                                                       aCTTCGTCACCGCATACGACCAGTTAAAGCCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACg  <  1:3806265/93‑1 (MQ=255)

                                                                                       aCTTCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACg  <  1:6310473/93‑1 (MQ=255)

ACAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGTAGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTCCGTCACCGCATACGACCAGTTAAAACCTGGCGGCGGGATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACG  >  NC_000913/161491‑161670

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