BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	19,161	T→C	I149I (ATT→ATC)	nhaR →	NhaR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	19,161	0	T	C	86.9%	776.6 / 102.6	312	I149I (ATT→ATC)	nhaR	NhaR
Reads supporting (aligned to +/- strand): ref base T (21/20); new base C (152/119); total (173/139)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.29e-02

Only 100 of 323 total aligned reads displayed.

AGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATCCA  >  NC_000913/19071‑19251

 gAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATctct                                                                                         >  1:1318401/1‑93 (MQ=255)

 gAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATctct                                                                                         >  2:2273060/1‑93 (MQ=255)

   aGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGa                                                                                       >  1:1186095/1‑93 (MQ=255)

    ggCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGAc                                                                                      <  2:221171/93‑1 (MQ=255)

        aGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTc                                                                                  <  2:3443046/93‑1 (MQ=255)

        aGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTc                                                                                  >  2:4315075/1‑93 (MQ=255)

         gCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTcc                                                                                 >  1:7483028/1‑93 (MQ=255)

         gCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTcc                                                                                 <  1:574738/93‑1 (MQ=255)

         gCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTcc                                                                                 >  1:2827452/1‑93 (MQ=255)

          cccATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCg                                                                                <  2:3704410/93‑1 (MQ=255)

           ccATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGa                                                                               <  2:2919233/93‑1 (MQ=255)

           ccATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGa                                                                               >  1:5923404/1‑93 (MQ=255)

           ccATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGa                                                                               <  2:7259352/93‑1 (MQ=255)

            cATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGAt                                                                              <  2:6384173/93‑1 (MQ=255)

             aTTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATa                                                                             >  2:4843499/1‑93 (MQ=255)

             aTTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATa                                                                             >  2:5088307/1‑93 (MQ=255)

              ttCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAg                                                                            >  1:4012143/1‑93 (MQ=255)

               tCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGa                                                                           <  2:4501906/93‑1 (MQ=255)

               tCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAgg                                                                           >  1:225217/1‑92 (MQ=255)

               tCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGa                                                                           <  1:3756152/93‑1 (MQ=255)

               tCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGa                                                                           <  2:1070059/93‑1 (MQ=255)

                cATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGAc                                                                          <  2:1482876/93‑1 (MQ=255)

                cATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGAc                                                                          >  1:5695706/1‑93 (MQ=255)

                 aTCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGAct                                                                         <  2:4546269/93‑1 (MQ=255)

                  tCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGActc                                                                        <  2:5207567/93‑1 (MQ=255)

                  tCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGActc                                                                        <  2:782268/93‑1 (MQ=255)

                   cTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGActct                                                                       <  1:1484382/93‑1 (MQ=255)

                   cTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGActct                                                                       <  1:3027475/93‑1 (MQ=255)

                    ttCGCTGCTTCGAATCCACCCAGGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTa                                                                      <  1:5654403/93‑1 (MQ=255)

                    ttCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTa                                                                      <  2:4146234/93‑1 (MQ=255)

                          gctTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACgcagc                                                                >  2:5542444/1‑93 (MQ=255)

                           ctTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACgcagca                                                               >  2:1357711/1‑93 (MQ=255)

                           ctTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACgcagca                                                               >  2:2549378/1‑93 (MQ=255)

                            ttCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGcagcag                                                              <  1:4683998/93‑1 (MQ=255)

                            ttCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGcagcag                                                              >  1:5333085/1‑93 (MQ=255)

                             tCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAgg                                                             <  1:2375574/93‑1 (MQ=255)

                              cGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGa                                                            <  1:495821/93‑1 (MQ=255)

                               gAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGaa                                                           <  1:4134714/93‑1 (MQ=255)

                               gAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGaa                                                           >  2:2780753/1‑93 (MQ=255)

                                aaTCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAg                                                          >  1:6504062/1‑93 (MQ=255)

                                  tCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGc                                                        <  1:6729083/93‑1 (MQ=255)

                                    cacccacGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACAGTCCGATAGACTCTACGCAGCAGGAAGGCCt                                                      >  1:6333235/1‑93 (MQ=255)

                                      cccacGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGt                                                    >  1:457205/1‑93 (MQ=255)

                                        cacGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTtc                                                  <  1:6587096/93‑1 (MQ=255)

                                        cacGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTtc                                                  >  1:6817216/1‑93 (MQ=255)

                                        cacGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTtc                                                  >  1:2146213/1‑93 (MQ=255)

                                         acGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTtct                                                 <  1:5189944/93‑1 (MQ=255)

                                          cGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTtctc                                                >  1:5110841/1‑93 (MQ=255)

                                           gAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTcc                                               >  2:5811020/1‑93 (MQ=255)

                                            aaaTGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCg                                              >  1:526547/1‑93 (MQ=255)

                                             aaTGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGt                                             >  1:4251047/1‑93 (MQ=255)

                                             aaTGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGt                                             <  1:7098452/93‑1 (MQ=255)

                                             aaTGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGt                                             >  1:5291501/1‑93 (MQ=255)

                                             aaTGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGt                                             <  2:4018239/93‑1 (MQ=255)

                                                  tgctgGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGaa                                        <  2:1018754/93‑1 (MQ=255)

                                                  tgctgGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGaa                                        <  1:5829415/93‑1 (MQ=255)

                                                  tgctgGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGaa                                        <  1:3317693/93‑1 (MQ=255)

                                                   gctgGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAAt                                       >  1:925723/1‑93 (MQ=255)

                                                     tgGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTg                                     <  2:410011/93‑1 (MQ=255)

                                                     tgGAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTg                                     >  1:5520953/1‑93 (MQ=255)

                                                      ggAGCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTgg                                    >  1:5068069/1‑93 (MQ=255)

                                                      ggAGCAACTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTgg                                    <  2:316906/93‑1 (MQ=255)

                                                         gCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGa                                 >  1:6408101/1‑93 (MQ=255)

                                                         gCAATTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGa                                 >  2:3529819/1‑93 (MQ=255)

                                                           aaTTAAGTCAGCATAAACTGGATATGATCTTCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAAt                               >  1:4470042/1‑93 (MQ=255)

                                                           aaTTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGAGTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAAt                               >  1:7303246/1‑93 (MQ=255)

                                                           aaTTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAAt                               >  2:5313676/1‑93 (MQ=255)

                                                            aTTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAAtg                              <  2:1654229/93‑1 (MQ=255)

                                                            aTTAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAAtg                              >  1:293890/1‑93 (MQ=255)

                                                             ttAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAAtgt                             >  2:1562633/1‑93 (MQ=255)

                                                             ttAAGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAAtgt                             >  1:6021959/1‑93 (MQ=255)

                                                               aaGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTgg                           <  2:7036779/93‑1 (MQ=255)

                                                                aGTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGc                          >  2:6809358/1‑93 (MQ=255)

                                                                 gTCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCg                         <  1:3797147/93‑1 (MQ=255)

                                                                  tCAGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGt                        >  1:3356744/1‑93 (MQ=255)

                                                                   cAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTg                       <  2:6117627/93‑1 (MQ=255)

                                                                    aGCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGa                      >  2:213020/1‑93 (MQ=255)

                                                                     gCATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAg                     <  2:7316439/93‑1 (MQ=255)

                                                                      cATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGt                    >  2:1984907/1‑93 (MQ=255)

                                                                      cATAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGt                    >  2:4983881/1‑93 (MQ=255)

                                                                        tAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGttt                  >  1:1056036/1‑93 (MQ=255)

                                                                        tAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGttt                  >  1:720113/1‑93 (MQ=255)

                                                                        tAAACTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGttt                  >  1:2314138/1‑93 (MQ=255)

                                                                          aaCTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATTTGGCGTGAGTTTCt                <  2:6442286/93‑1 (MQ=255)

                                                                          aaCTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCt                <  2:2589857/93‑1 (MQ=255)

                                                                           aCTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTg               >  1:3226544/1‑93 (MQ=255)

                                                                            cTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTgg              <  1:5392103/93‑1 (MQ=255)

                                                                            cTGGATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTgg              <  1:3675307/93‑1 (MQ=255)

                                                                               gATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGgtgt           >  2:6025874/1‑93 (MQ=255)

                                                                               gATATGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGgtgt           <  1:2049171/93‑1 (MQ=255)

                                                                                atatGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTa          <  2:715898/93‑1 (MQ=255)

                                                                                 tatGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTAc         >  2:6619901/1‑93 (MQ=255)

                                                                                  atGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACa        <  1:2604298/93‑1 (MQ=255)

                                                                                  atGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACa        >  2:6327933/1‑93 (MQ=255)

                                                                                   tGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACaa       >  1:3496538/1‑93 (MQ=255)

                                                                                   tGATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACaa       >  2:3831826/1‑93 (MQ=255)

                                                                                    gATCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACaaa      <  2:595721/93‑1 (MQ=255)

                                                                                     aTCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAAt     >  2:201422/1‑93 (MQ=255)

                                                                                      tCATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATc    <  2:542943/93‑1 (MQ=255)

                                                                                       cATCTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATcc   <  1:3017288/93‑1 (MQ=255)

AGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATCCA  >  NC_000913/19071‑19251

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