BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	673,104	T→C	K301K (AAA→AAG)	leuS ←	leucine tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	673,104	0	T	C	97.8%	407.9 / 0.0	134	K301K (AAA→AAG)	leuS	leucine tRNA synthetase
Reads supporting (aligned to +/- strand): ref base T (2/1); new base C (67/64); total (69/65)
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.21e-01

Only 100 of 152 total aligned reads displayed.

GCCGTACTCCATCAATACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCTTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAGGATTATTTTCCGC  >  NC_000913/673014‑673193

gCCGTACTCCATCAATACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCctt                                                                                         >  2:2130795/1‑93 (MQ=255)

 ccGTACTCCATCAATACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCtttt                                                                                        <  2:4276108/93‑1 (MQ=255)

 ccGTACTCCATCAATACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCcttt                                                                                        <  1:3519941/93‑1 (MQ=255)

  cGTACTCCATCAATACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCctttt                                                                                       <  2:1825336/93‑1 (MQ=255)

   gTACTCCATCAATACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCcttttt                                                                                      >  1:4362354/1‑93 (MQ=255)

   gTACTCCATCAATACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTCTAAAGCCAGTATCGACGTCcttttt                                                                                      >  2:203162/1‑93 (MQ=255)

    tACTCCATCAATACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTtc                                                                                     >  1:541429/1‑93 (MQ=255)

    tACTCCATCAATACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTtc                                                                                     <  1:3785500/93‑1 (MQ=255)

       tCCATCAATACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTcc                                                                                  >  1:4032247/1‑93 (MQ=255)

        ccATCAATACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCa                                                                                 <  1:3268866/93‑1 (MQ=255)

        ccATCAATACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCa                                                                                 <  1:2099853/93‑1 (MQ=255)

         cATCAATACTAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCAt                                                                                <  1:3320805/93‑1 (MQ=255)

          aTCAATACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCAtc                                                                               <  2:166339/93‑1 (MQ=255)

          aTCAATACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCAtc                                                                               <  1:4737670/93‑1 (MQ=255)

          aTCAATACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCAtc                                                                               >  2:3119212/1‑93 (MQ=255)

           tCAATACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCAtcg                                                                              >  1:4274498/1‑93 (MQ=255)

             aaTACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCAtcgtc                                                                            <  2:4455014/93‑1 (MQ=255)

               tACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGc                                                                          <  2:1417822/93‑1 (MQ=255)

                  gAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCAt                                                                       >  2:2995771/1‑93 (MQ=255)

                   aaGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCAtt                                                                      >  1:3479736/1‑93 (MQ=255)

                      tttGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCa                                                                   >  1:703491/1‑93 (MQ=255)

                      gttGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCa                                                                   <  1:4069269/92‑1 (MQ=255)

                        tgctgcCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCGTCGTCGCCATTTCAGc                                                                 <  1:531501/93‑1 (MQ=255)

                          ctgcCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCtt                                                               <  1:878532/93‑1 (MQ=255)

                          ctgcCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCtt                                                               <  2:2540760/93‑1 (MQ=255)

                           tgcCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTc                                                              <  1:3932195/93‑1 (MQ=255)

                           tgcCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTc                                                              >  2:1545858/1‑93 (MQ=255)

                           cgcCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTc                                                              <  1:3932197/92‑1 (MQ=255)

                              ccAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGc                                                           <  1:3010395/93‑1 (MQ=255)

                              ccAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGc                                                           <  1:1317250/93‑1 (MQ=255)

                               cAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCa                                                          >  2:1715042/1‑93 (MQ=255)

                               cAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCa                                                          >  2:3069667/1‑93 (MQ=255)

                                aaaCGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCaa                                                         >  2:506667/1‑93 (MQ=255)

                                aaaCGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCaa                                                         >  2:864343/1‑93 (MQ=255)

                                  aCGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACt                                                       >  1:4685884/1‑93 (MQ=255)

                                     ggAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTg                                                    >  2:809178/1‑93 (MQ=255)

                                     ggAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTg                                                    >  2:4152586/1‑93 (MQ=255)

                                      gAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTgg                                                   <  2:15137/93‑1 (MQ=255)

                                          ttcttcTCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGtt                                               <  2:326322/93‑1 (MQ=255)

                                           tcttcGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTa                                              >  2:4578312/1‑93 (MQ=255)

                                             ttcGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACg                                            <  1:2060565/93‑1 (MQ=255)

                                                  ccGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCAtt                                       >  2:2609676/1‑93 (MQ=255)

                                                  ccGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCAtt                                       <  2:2771937/93‑1 (MQ=255)

                                                  ccGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCAtt                                       >  2:1441167/1‑93 (MQ=255)

                                                   cGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTc                                      <  1:3753403/93‑1 (MQ=255)

                                                   cGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTc                                      >  2:2198628/1‑93 (MQ=255)

                                                   cGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTc                                      <  2:3638614/93‑1 (MQ=255)

                                                     ttAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGt                                    >  1:4456496/1‑93 (MQ=255)

                                                     ttAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGt                                    <  1:2478538/93‑1 (MQ=255)

                                                     ttAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGt                                    >  2:641880/1‑93 (MQ=255)

                                                       aaTGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCa                                  >  2:2319314/1‑93 (MQ=255)

                                                       aaTGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGACAGTTGAGCTTCGGCAGCTTTGGTGTTGGGGCATTCGTCa                                  >  1:2541939/1‑93 (MQ=255)

                                                        aTGTGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCaa                                 <  1:1135888/93‑1 (MQ=255)

                                                        aTGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCaa                                 <  1:1135889/93‑1 (MQ=255)

                                                         tGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAAt                                <  1:460557/93‑1 (MQ=255)

                                                         tGGGTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAAt                                >  1:1567537/1‑93 (MQ=255)

                                                          gggTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATa                               <  2:3539917/93‑1 (MQ=255)

                                                          gggTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATa                               >  2:4061734/1‑93 (MQ=255)

                                                            gTGAACCGCTTTAAAGCCAGTATCGACGCCTTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATaaa                             <  1:4327440/93‑1 (MQ=255)

                                                            gTGAACCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATaaa                             >  1:24244/1‑93 (MQ=255)

                                                             tGAACCGCTTTAAAGCCAGTATCGACTCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAg                            <  2:2836771/93‑1 (MQ=255)

                                                               aaCCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGgc                          <  2:329293/93‑1 (MQ=255)

                                                               aaCCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGgc                          <  1:1353908/93‑1 (MQ=255)

                                                               aaCCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGgc                          <  1:3139973/93‑1 (MQ=255)

                                                               aaCCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGgc                          <  2:2544568/93‑1 (MQ=255)

                                                               aaCCGCTTAAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGgc                          <  2:3969719/93‑1 (MQ=255)

                                                                aCCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGgcc                         >  1:1560390/1‑93 (MQ=255)

                                                                aCCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGgcc                         >  2:2646776/1‑93 (MQ=255)

                                                                aCCGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGgcc                         <  2:601830/93‑1 (MQ=255)

                                                                  cGCTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGgccgc                       >  2:2779555/1‑93 (MQ=255)

                                                                   gCTTTAAAGCCATTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCATCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGgccgcc                      <  2:2002461/93‑1 (MQ=255)

                                                                    cTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCa                     >  2:2415584/1‑93 (MQ=255)

                                                                    cTTTAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCa                     <  1:937136/93‑1 (MQ=255)

                                                                     tttAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAg                    >  2:1893706/1‑93 (MQ=255)

                                                                      ttAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGt                   >  2:927648/1‑93 (MQ=255)

                                                                      ttAAAGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGt                   <  2:1298897/93‑1 (MQ=255)

                                                                       tAAAGCCAGTAGCGACGCCCTTTTACTCCATCGTCGCCATTTCATCTTCGGAAACTTTCGTGTAACGGCATTCGTCAATAAAAGCCGCCAGtt                  <  2:3160850/93‑1 (MQ=255)

                                                                         aaGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCa                <  2:3270504/93‑1 (MQ=255)

                                                                         aaGCCAGTATCGACGCCCTTTGTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCa                <  2:3270503/93‑1 (MQ=255)

                                                                          aGCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAg               <  1:714275/93‑1 (MQ=255)

                                                                          aGCCAGTATCGACACCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAg               <  1:295410/93‑1 (MQ=255)

                                                                           gCCAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCacg              >  1:4286476/1‑91 (MQ=255)

                                                                            ccAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAGGa             >  1:2700995/1‑93 (MQ=255)

                                                                             cAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAGGat            <  2:4104983/93‑1 (MQ=255)

                                                                             cAGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAGGat            >  2:2855613/1‑93 (MQ=255)

                                                                              agcaTCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAGGatt           >  2:4094987/4‑93 (MQ=255)

                                                                              aGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAGGatt           >  1:3539927/1‑93 (MQ=255)

                                                                              aGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAGGatt           <  1:2204386/93‑1 (MQ=255)

                                                                              aGTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGACCGCCAGTTCAGGatt           >  2:1361944/1‑93 (MQ=255)

                                                                               gTATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAGGatta          <  1:2542646/93‑1 (MQ=255)

                                                                                tATCGACGCCTTTTTTCTCCCTCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAGGattat         <  1:2677481/93‑1 (MQ=255)

                                                                                tATCGACGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAGGattat         >  2:3948973/1‑93 (MQ=255)

                                                                                    gACGCCCTTTTTCTCCCTGGTCGCCATTTCAGCTTCGTGAACTTTGGAGTTACGACATTCGTCAATAGAGGCAGCCAGTTCAGGATTATTTTc     >  2:4342488/1‑93 (MQ=255)

                                                                                    gACGCCCTTTTTCTCCATCGTCGCCCTTTCAGCTTCGGCAACTTTGGTGTTACGGCAGTCGTCAATAAAGGCCGCCAGTTCAGGATTATTTTc     >  1:2074802/1‑93 (MQ=255)

                                                                                     aCGCCTTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAGGATTATTTTcc    >  1:3127498/1‑93 (MQ=255)

                                                                                     aCGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAGGATTATTTTcc    >  2:146817/1‑93 (MQ=255)

                                                                                      cGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAGGATTATTTTccg   >  2:4559263/1‑93 (MQ=255)

                                                                                      cGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAGGATTATTTTccg   >  1:2103802/1‑93 (MQ=255)

                                                                                      cGCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAGGATTATTTTccg   <  1:2359137/93‑1 (MQ=255)

                                                                                       gCCCTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAGGATTATTTTccgc  >  2:379480/1‑93 (MQ=255)

GCCGTACTCCATCAATACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAATGGGTGAACCGCTTTAAAGCCAGTATCGACGCCTTTTTTCTCCATCGTCGCCATTTCAGCTTCGGCAACTTTGGTGTTACGGCATTCGTCAATAAAGGCCGCCAGTTCAGGATTATTTTCCGC  >  NC_000913/673014‑673193

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