BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,857,508	T→C	E79G (GAA→GGA)	yidK ←	putative cotransporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,857,508	0	T	C	99.4%	497.2 / ‑4.9	153	E79G (GAA→GGA)	yidK	putative cotransporter
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (74/78); total (75/78)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.85e-01

Only 100 of 169 total aligned reads displayed.

AGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTTCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAACGAA  >  NC_000913/3857418‑3857597

aGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTccc                                                                                         >  2:4269217/1‑93 (MQ=255)

 gagaaTCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTcccc                                                                                        <  1:1810793/90‑1 (MQ=255)

 gAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTcccc                                                                                        >  2:1641029/1‑93 (MQ=255)

  aaaaTCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCa                                                                                       >  2:1811320/1‑93 (MQ=255)

   aaaTCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAg                                                                                      >  2:3928676/1‑93 (MQ=255)

    aaTCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGc                                                                                     >  1:295274/1‑93 (MQ=255)

     aTCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGcc                                                                                    >  2:1463259/1‑93 (MQ=255)

       cGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCa                                                                                  <  1:2754071/93‑1 (MQ=255)

         gggATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATc                                                                                <  1:3819575/93‑1 (MQ=255)

         gggATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATc                                                                                >  2:1904614/1‑93 (MQ=255)

         ggaaTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGCAGATCAGCGTCACCGCTTATGGCACTTCCCAGTACATc                                                                                >  2:3823605/4‑93 (MQ=255)

          gtaTGGTGGCAATGCCGCGCTTCAGATAGCGCGTTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCa                                                                               <  2:4611650/91‑1 (MQ=255)

             ttgtggCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCg                                                                            <  2:2697683/91‑1 (MQ=255)

              ggtggCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGa                                                                           <  1:4116963/93‑1 (MQ=255)

                tggCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACa                                                                         <  1:2967854/93‑1 (MQ=255)

                 ggCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACAt                                                                        >  1:4250927/1‑93 (MQ=255)

                 ggCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACAt                                                                        <  1:2078174/93‑1 (MQ=255)

                    aaTGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGcc                                                                     >  1:815280/1‑93 (MQ=255)

                     aTGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCg                                                                    >  1:3171059/1‑93 (MQ=255)

                     aTGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCg                                                                    >  2:1410859/1‑93 (MQ=255)

                     aTGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCg                                                                    >  2:2781411/1‑93 (MQ=255)

                        ccGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCtt                                                                 <  2:4313293/93‑1 (MQ=255)

                         cgcgCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCttt                                                                <  2:2554839/93‑1 (MQ=255)

                         cgcgCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCttt                                                                <  2:2428758/93‑1 (MQ=255)

                         cgcgCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCttt                                                                <  2:1032559/93‑1 (MQ=255)

                         cgcgCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCttt                                                                <  2:734897/93‑1 (MQ=255)

                          gcgcTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCtttt                                                               >  2:3952071/1‑93 (MQ=255)

                           cgcTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTg                                                              <  1:44109/93‑1 (MQ=255)

                            gcTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCTCTCCCCAGCCCATCACCGACATGCCGCTTTTGt                                                             <  1:1096714/93‑1 (MQ=255)

                            gcTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGt                                                             >  1:282279/1‑93 (MQ=255)

                             cTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGtg                                                            >  2:1565213/1‑92 (MQ=255)

                               tCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAgg                                                          <  2:1413946/93‑1 (MQ=255)

                                 aGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGcc                                                        >  1:771345/1‑93 (MQ=255)

                                  gATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCt                                                       <  1:94028/93‑1 (MQ=255)

                                  gATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCt                                                       <  2:520018/93‑1 (MQ=255)

                                    tAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGc                                                     >  2:2645407/1‑93 (MQ=255)

                                     aGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGcc                                                    <  2:2257381/93‑1 (MQ=255)

                                      gcgcGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGccc                                                   <  2:4150916/93‑1 (MQ=255)

                                        gcgGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCgg                                                 >  1:1329646/1‑93 (MQ=255)

                                         cgGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGa                                                <  2:3866332/93‑1 (MQ=255)

                                          ggTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTTCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGaa                                               >  2:3995568/1‑93 (MQ=255)

                                          ggTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGaa                                               >  2:2217139/1‑93 (MQ=255)

                                          ggTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGaa                                               >  1:1148362/1‑93 (MQ=255)

                                          ggTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGaa                                               <  1:252244/93‑1 (MQ=255)

                                           gTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGaaa                                              <  2:2187808/93‑1 (MQ=255)

                                           gTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGaaa                                              <  1:1909732/93‑1 (MQ=255)

                                            tAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAg                                             >  1:841037/1‑93 (MQ=255)

                                             aaaaaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAgg                                            <  1:3465075/93‑1 (MQ=255)

                                              aaaaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGc                                           >  1:113121/1‑93 (MQ=255)

                                              aaaaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGc                                           <  2:4272132/93‑1 (MQ=255)

                                               aaaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGcc                                          <  1:1623895/93‑1 (MQ=255)

                                               aaaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGcc                                          >  1:3740465/1‑93 (MQ=255)

                                               aaaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTACGCCTGCCCGGAAAGGcc                                          >  1:1276120/1‑93 (MQ=255)

                                                aaGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCg                                         >  1:1964479/1‑93 (MQ=255)

                                                   aTTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGAcc                                      >  2:131518/1‑93 (MQ=255)

                                                   aTTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGAcc                                      <  1:1271095/93‑1 (MQ=255)

                                                    ttAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCa                                     <  1:2689234/93‑1 (MQ=255)

                                                      aGCGCGAGGAAGATCAGCATCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAgt                                   <  2:4417705/93‑1 (MQ=255)

                                                          cgAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAgttgtt                               >  2:614613/1‑93 (MQ=255)

                                                           gAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCGTGCCCGGAAAGGCCGACCAGTTGATc                              >  1:2020174/1‑93 (MQ=255)

                                                           gAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTc                              >  2:483845/1‑93 (MQ=255)

                                                           gAGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTc                              <  1:4234082/93‑1 (MQ=255)

                                                            aGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTcc                             <  1:1704373/93‑1 (MQ=255)

                                                            aGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTcc                             <  2:1740621/93‑1 (MQ=255)

                                                            aGGAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTcc                             >  1:3078038/1‑93 (MQ=255)

                                                             ggAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCg                            >  1:4152135/1‑93 (MQ=255)

                                                             ggAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCg                            >  1:3770773/1‑93 (MQ=255)

                                                              gAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGt                           <  2:698726/93‑1 (MQ=255)

                                                              gAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGt                           <  2:412776/93‑1 (MQ=255)

                                                              gAAGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGt                           <  2:2851713/93‑1 (MQ=255)

                                                               aaGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTg                          <  1:2308556/93‑1 (MQ=255)

                                                                aGATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTgg                         >  1:4349700/1‑93 (MQ=255)

                                                                 gATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGa                        >  1:2452063/1‑93 (MQ=255)

                                                                 gATCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGa                        <  1:2283333/93‑1 (MQ=255)

                                                                   tCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGaaa                      >  2:518059/1‑93 (MQ=255)

                                                                   tCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGaaa                      <  2:1431328/93‑1 (MQ=255)

                                                                   tCAGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGaaa                      <  2:112637/93‑1 (MQ=255)

                                                                     aGCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAgg                    <  1:445630/93‑1 (MQ=255)

                                                                      gctgcACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggt                   <  1:661195/89‑1 (MQ=255)

                                                                      gCGTCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggt                   >  1:4339003/1‑93 (MQ=255)

                                                                      gCGTCACCACTGAAGTCGCGCCCCAGCCCATCCCCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggt                   <  1:4635283/93‑1 (MQ=255)

                                                                         tCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggttgg                >  2:87559/1‑93 (MQ=255)

                                                                         tCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggttgg                >  1:2206768/1‑93 (MQ=255)

                                                                         tCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggttgg                >  1:967283/1‑93 (MQ=255)

                                                                         tCACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggttgg                <  2:3149483/93‑1 (MQ=255)

                                                                          cACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAggttggt               >  2:1123217/1‑93 (MQ=255)

                                                                          cACCGCTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCAACCAGTTGTTCCGTGGAAAggttggt               <  2:1778370/93‑1 (MQ=255)

                                                                               cTGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCa          <  2:2055999/93‑1 (MQ=255)

                                                                                tGAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCAt         >  2:310602/1‑93 (MQ=255)

                                                                                 gAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCAtt        >  1:4108480/1‑93 (MQ=255)

                                                                                 gAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCAtt        >  2:659484/1‑93 (MQ=255)

                                                                                 gAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCAtt        >  1:448994/1‑93 (MQ=255)

                                                                                 gAAGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCAtt        >  1:506882/1‑93 (MQ=255)

                                                                                  aaGTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTa       <  2:2549393/93‑1 (MQ=255)

                                                                                    gTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAAc     <  1:3013995/93‑1 (MQ=255)

                                                                                    gTCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAAc     <  1:3862387/93‑1 (MQ=255)

                                                                                     tCACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAACg    >  2:3071562/1‑93 (MQ=255)

                                                                                      cACTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTAGGTTAGCATTAACGa   >  2:3764170/1‑93 (MQ=255)

                                                                                       aCTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAACGaa  >  1:1495801/1‑93 (MQ=255)

                                                                                       aCTCCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAACGaa  <  2:2283201/93‑1 (MQ=255)

AGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTTCCCAGCCCATCACCGACATGCCGCTTTTGTAGGCCTGCCCGGAAAGGCCGACCAGTTGTTCCGTGGAAAGGTTGGTTAGCATTAACGAA  >  NC_000913/3857418‑3857597

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