BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,057,443	G→A	100%	L198L (CTG→TTG)	glnA ←	glutamine synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,057,443	0	G	A	59.1%	111.0 / 213.7	194	L198L (CTG→TTG)	glnA	glutamine synthetase
Reads supporting (aligned to +/- strand): ref base G (42/37); new base A (61/53); total (103/91)
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 197 total aligned reads displayed.

TGGTCATGGTATTGAAGCGGGTAGCCACTTCGTTCTGACCAGCAGTCGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAGACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCGGACGGTGACCTTTGTTACCACCTTCG  >  NC_000913/4057344‑4057540

    cATGGTATTGAAGCGGGTAGCCACTTCGTTCTGACCAGCAGTCGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAAacac                                                                                               >  1:517405/1‑100 (MQ=255)

     aTGGTATTGAAGCGGGTAGCCACTTCGTTCTGACCAGCAGTCGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAGAcaca                                                                                              >  1:3477614/1‑100 (MQ=255)

      tGGTATTGAAGCGGGTAGCCACTTCGTTCTAACCAGCAGTCGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAAACACAt                                                                                             >  2:1847382/1‑100 (MQ=255)

       ggTATTGAAGCGGGTAGCCACTTCGTTCTGACCAGCAGTCGCTACTTCGTGATGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAAACACAtt                                                                                            >  1:1942966/1‑100 (MQ=255)

        gTATTGAAGCGGGTAGCCACTTCGTTCTGACCAGCAGTCGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAGACACAttt                                                                                           >  2:1366575/1‑100 (MQ=255)

         tATTGAAGCGGGTAGCCACTTCGTTCTGACCAGCAGTCGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAGACACATTTc                                                                                          >  2:2940283/1‑100 (MQ=255)

         tATTGAAGCGGGTAGCCACTTCGTTCTGACCAGCAGTCGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAGACACATTTc                                                                                          <  1:2564446/100‑1 (MQ=255)

         tATTGAAGCGGGTAGCCACTTCGTTCTGACCAGCAGTCGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAAACACATTTc                                                                                          <  1:2705268/100‑1 (MQ=255)

         tATTGAAGCGGGTAGCCACTTCGTTCTGACCAGCAGTCGCTACTTCGTGGTGATGGGCTTCAAACACCAGACCCATCTGTTCCATCACCAAACACATTTc                                                                                          >  1:1009199/1‑100 (MQ=255)

          aTTGAAGCGGGTAGCCACTTCGTTCTGACCAGCAGTCGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAAACACATTTCa                                                                                         >  1:1860478/1‑100 (MQ=255)

           ttGAAGCGGGTAGCCACTTCGTTCTGACCAGCAGTCGCTACTTCGTGGTGATGGGCTTCAACCACAAGACCGATCTGTTCCATCACCAAAAAAATTTCa                                                                                         >  1:589032/1‑99 (MQ=255)

            tGAAGCGGGTAGCCACTTCGTTCTGACCAGCAGTCGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAAACACATTTCAGa                                                                                       <  2:1186830/100‑1 (MQ=255)

             gAAGCGGGTAGCCACTTCGTTCTGACCAGCAGTCGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAGACACATTTCAGaa                                                                                      >  1:1397813/1‑100 (MQ=255)

             gAAGCGGGTAGCCACTTCGTTCTGACCAGCAGTCGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAAACACATTTCAGaa                                                                                      >  1:230957/1‑100 (MQ=255)

             gAAGCGGGTAGCCACTTCGTTCTGACCAGCAGTCGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAAACACATTTCAGaa                                                                                      >  1:2270964/1‑100 (MQ=255)

              aaGCGGGTAGCCACTTCGTTCTGACCAGCAGTCGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAGACACATTTCAGAAc                                                                                     <  1:379670/100‑1 (MQ=255)

              aaGCGGGTAGCCACTTCGTTCTGACCAGCAGTCGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAAACACATTTCAGAAc                                                                                     <  2:3536736/100‑1 (MQ=255)

               aGCGGGTAGCCACTTCGTTCTGACCAGCAGTCGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAGACACATTTCAGAACg                                                                                    >  1:1681929/1‑100 (MQ=255)

               aGCGGGTAGCCACTTCGTTCTGACCAGCAGTCGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAAACACATTTCAGAACg                                                                                    <  1:3200359/100‑1 (MQ=255)

                   ggTAGCCACTTCGTTCTGACCAGCAGTCGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAGACACATTTCAGCACGAata                                                                                >  2:248038/1‑100 (MQ=255)

                      aGCCACTTCGTTCTGACCAGCAGTCGCTACTTCGTGGTGATGGGCTTCAACCACTAGACCCATCTGTTCCATCACCAAACACATTTCAGAACGAATATcc                                                                             <  2:3234115/100‑1 (MQ=255)

                      aGCCACTTCGTTCTGACCAGCAGTCGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAGACACATTTCAGAACGAATATcc                                                                             >  2:2535245/1‑100 (MQ=255)

                         cACTTCGTTCTGACCAGCAGTCGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAAACACATTTCAGAACGAATATCCTGa                                                                          >  2:1526701/1‑100 (MQ=255)

                         cACTTCGTTCTGACCAGCAGTCGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAAACACATTTCAGAACGAATATCCTGa                                                                          >  2:1059201/1‑100 (MQ=255)

                         atattCGTTCTGACCAGCAGTCGCTACTTCGTGGTGATGGGCTTCAACCATCAGACCCATCTGTTCCATCACCAAACACATTTCAGAACGAATATCCTGa                                                                          <  1:1004273/97‑1 (MQ=255)

                           cTTCGTTCTGACCAGCAGTCGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAAACACATTTCAGAACGAATATCCTGAGc                                                                        <  1:2117188/100‑1 (MQ=255)

                               gTTCTGACCAGCAGTCGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAAACACATTTCAGAACGAATATCCTGAGCagag                                                                    >  2:1280632/1‑100 (MQ=255)

                               gTTCTGACCAGCAGTCGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAAACACATTTCAGAACGAATATCCTGAGCCGAg                                                                    >  2:1657600/1‑100 (MQ=255)

                               gTTCTGACCAGCAGTCGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAAACACATTTCAGAACGAATATCCTGAGCCGAg                                                                    <  2:1994047/100‑1 (MQ=255)

                                 tCTGACCAGCAGTCGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAGACACATTTCAGAACGAATATCCTGAGCCGAGTc                                                                  <  2:618952/100‑1 (MQ=255)

                                    gACCAGCAGTCGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAAACACATTTCAGAACGAATACCCTGAGCCGAGTCTAc                                                               >  2:836674/1‑100 (MQ=255)

                                        agcagTCGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAAACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGt                                                           >  1:2825696/1‑100 (MQ=255)

                                         gtagtccatACTTCGTGGTAATGGGCTTCAACCACCAGATCCATCTGTTCCATCACCAAACACATTTCAGAACGAATATTCTGAGCCGAGTCTACCgg                                                            <  2:3280499/90‑1 (MQ=255)

                                         gcagTCGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAGACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTg                                                          >  2:2426219/1‑100 (MQ=255)

                                              cGCTACTTCGTTGTGATGGGCTTCAACCATCAGACCCATCTGTTCCATCACCAAACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAAcc                                                     <  1:789047/100‑1 (MQ=255)

                                              cGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAGACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAAcc                                                     >  2:1324639/1‑100 (MQ=255)

                                              cGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAGACACATTTCAGAACGAATATCATGAGCCGAGTCTACCGGTGGAAcc                                                     >  2:1555805/1‑100 (MQ=255)

                                              cGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAAACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAAcc                                                     >  1:3921189/1‑100 (MQ=255)

                                               gCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAGACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCg                                                    >  1:1223654/1‑100 (MQ=255)

                                               gCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACAAGACACATTTCAGAAAAAATATCCTGAGCCGAGTCTACCGGTGGAAc                                                      >  1:207884/1‑98 (MQ=255)

                                                cTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAGACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCgg                                                   >  2:2896342/1‑100 (MQ=255)

                                                 tACTTCGTGTTGATGGGCTTCAACCATCAGACCTATCTGTTCCATCACCAGAAACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCggg                                                  <  1:1179395/100‑1 (MQ=255)

                                                  aCTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAGACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGa                                                 <  1:3423981/100‑1 (MQ=255)

                                                  aCTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAAACACATTTCAGAACGAATATCCTAAGCCGAGTCTACCGGTGGAACCGGGa                                                 >  1:412072/1‑100 (MQ=255)

                                                  aCTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGCTCCATCACCAAACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGa                                                 <  2:2121775/100‑1 (MQ=255)

                                                    ttCGTGGTGATGGGCTTCAACCACCAGACCCCTCTGTCCCATCACCAAACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAg                                               <  2:2636165/100‑1 (MQ=255)

                                                    ttCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAGACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAg                                               <  2:1729561/100‑1 (MQ=255)

                                                    ttCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAGACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAg                                               <  1:3296891/100‑1 (MQ=255)

                                                      cGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAAACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAg                                               <  2:589025/98‑1 (MQ=255)

                                                      cGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAAACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTa                                             >  2:1081926/1‑100 (MQ=255)

                                                         ggtgATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAAACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAAcc                                          <  1:2367417/100‑1 (MQ=255)

                                                           tgATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAGACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGc                                        <  2:1397798/100‑1 (MQ=255)

                                                            gATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAGACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGcc                                       <  1:1366501/100‑1 (MQ=255)

                                                            gATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAAACACATTTCAGAACGAATATCTTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGc                                        >  2:1629737/1‑99 (MQ=255)

                                                              tGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAGACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGCCtt                                     >  2:2602095/1‑100 (MQ=255)

                                                                ggCTTCAACCACCAGACCCATCTGTTCCATCACCAGACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTc                                   <  2:203504/100‑1 (MQ=255)

                                                                 gCTTCAACCACCAGACCCATCTGTTCCATCACCAGACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCa                                  >  2:311303/1‑100 (MQ=255)

                                                                  cTTCAACCACCATACCCATCTGTTCCATCACCAAACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCAc                                 <  1:917939/100‑1 (MQ=255)

                                                                  cTTCAACCACCAGACCCATCTGTTCCATCACCAGACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCAc                                 >  2:3899272/1‑100 (MQ=255)

                                                                  cTTCAACCACCAGACCCATCTGTTCCATCACCAGACACATTTCAGAACGAATATCCTAAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCa                                  >  1:966104/1‑99 (MQ=255)

                                                                    tCAACCACCAGACCCATCTGTTCCATCACCAGACACATTTCAGAACGAATATCCTGAGCCGAGACTACCGGTGGAACCGGGAAGTAACCGCCTTTCACTg                               >  1:3378039/1‑100 (MQ=255)

                                                                     cAACCACCAGACCCATCTGTTCCATCACCAAACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCACTGc                              >  1:1129958/1‑100 (MQ=255)

                                                                     cAACCACCAGACCCATCTGTTCCATCACCAAACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCACTGc                              >  2:1753893/1‑100 (MQ=255)

                                                                      aaCCACCAGACCCATCTGTTCCATCACCAGACACATTTAAGAACGAATATCCTAAGCCGAAACTACCGGAGGAACAAAAAAGTAAACACCTTTCACTAcc                             >  1:1757088/1‑100 (MQ=255)

                                                                       accaccCGACCCATCTGTTCCATCACAAGACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGATCCGGGAAGTAACCGCCTTTCACTGcc                             <  2:1798008/99‑1 (MQ=255)

                                                                       accaccAGACCCATCTGTTCCATCACCATACACATTTTAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCg                            <  1:1826915/100‑1 (MQ=255)

                                                                       accaccAGACCCATCTGTTCCATCACCAAACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGTGAAGTAACCGCCTTTCACTGCCg                            <  2:2193295/100‑1 (MQ=255)

                                                                       accaccAAACCCATCTGTTCCATCACCAAACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGATGAAACCGGGAAGTAACCGCCTTTCACTGCCg                            >  1:2053992/1‑100 (MQ=255)

                                                                        ccaccaGACCCATCTGTTCCATCACCAAACACATTTCAGAACGAATATCATGAGCAGAGTATACAGGTGGAACAGGGAAGTAACCGCCTTTAACTGCCgg                           >  1:3109968/1‑100 (MQ=255)

                                                                         caccaGACCCATCTGTTCCATCACCAGACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCGGa                          >  1:2796103/1‑100 (MQ=255)

                                                                         caccaGACCCATCTGTTCCATCACCAAACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCGGa                          >  2:190438/1‑100 (MQ=255)

                                                                           ccaGACCCATCTGTTCCATCACCAAACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCGGACg                        >  2:1259218/1‑100 (MQ=255)

                                                                             aGACCCATCTGTTCCATCACCAAACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCGGACGGt                      <  2:424363/100‑1 (MQ=255)

                                                                                cccTTCTGTTCTTTCACCTAACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCGGACgg                       <  1:725928/96‑1 (MQ=255)

                                                                                cccATCTGTTCCATCACCAGACACATTTCAGAACGAATATCCTGAGCCGAGACTACCGGTGGAACCGGGAAGTAAACGCCTTTAACTGCCGGACGGTGAc                   >  1:3521400/1‑100 (MQ=255)

                                                                                  cATCTGTTCCATCACCAGACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCGGACGGTGACCt                 >  1:3486869/1‑100 (MQ=255)

                                                                                  cATCTATTCCATCACCAGACACATTTCAGAACGAATAGTCTGAGACGCGTCTACCGGTGGAACCGGGAAGTATCCGCTTTTCACTGCTGGACTGTGa                    <  2:2949803/97‑1 (MQ=255)

                                                                                     cTGTTCCATCACTAAACACATTTCAGAACGAATATCCTGAGCCTAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCGGACGGTGACCTTTg              <  1:1344233/100‑1 (MQ=255)

                                                                                      tGTTCCATCACCAGACACATTTCAGAACGAATATCCTGAGCCCAGGCTACCGACGGAACGTGGAAGTAACCGCCTTTCAC‑GCCgggcc                        >  2:2449418/1‑85 (MQ=255)

                                                                                      tGTTCCATCACCAAACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCGGACGGTGACGTTTGt             >  2:3083994/1‑100 (MQ=255)

                                                                                      tGTTCCATCACCAAACACATTTCAGAACGAATATCCTGAGCCGAGACTACCGGTGGAACCGGGAAATAACCGCCTTTCACTGCCGGACGATGACCTTTGt             >  1:2126642/1‑100 (MQ=255)

                                                                                       gTTCCATCACCAGACGCATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCGGACGGTGACCTTTGtt            <  1:3873429/100‑1 (MQ=255)

                                                                                       gTTCCATCACCAGACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCGGACGGTGACCTTTGtt            >  2:2041039/1‑100 (MQ=255)

                                                                                       gTTCCATCACCAAACACATTTCAGAACGAATATCCTGAGCCTAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCGGACGGTGACCTTTGtt            >  2:2131926/1‑100 (MQ=255)

                                                                                       gTTCCATCACCAAACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCGGACGGTGACCTTTGtt            >  1:2919180/1‑100 (MQ=255)

                                                                                       gTTCCATCACCAAACACATTTCAGAACGAATATCCTGAGCCGAGACTACCGGTGGAAACGGGAAGTAACCGCCTTTCACTGCCGGACGGTGACCTTTGtt            >  2:3644303/1‑100 (MQ=255)

                                                                                        ttCCATCCCTAGACACTTTTCAGAACGAATATCCTGTCCCGAGTCTACCTGTGGAATCCTTAAGTATCCCCCTTTCCCTTCCTCACGGTGACCTTTGTTa           <  2:3469593/100‑1 (MQ=255)

                                                                                        ttCCATCACCTAACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTATCCGCCTTTCACTGCCGGACGGTGACCTTTGTTa           <  2:2579123/100‑1 (MQ=255)

                                                                                        ttCCATCACCAAACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCGGACGGTGACCTTTGTTa           <  2:2666956/100‑1 (MQ=255)

                                                                                        ttCCATCACCAAACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCGGACGGTGACCTTTGTTa           <  1:1526716/100‑1 (MQ=255)

                                                                                          ccATCACCAGACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCGGACGGTGACCTTTGTTacc         >  2:770439/1‑100 (MQ=255)

                                                                                            aTCACCAGACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCGGACGGTGACCTTTGTTaccac       <  2:3121417/100‑1 (MQ=255)

                                                                                            aTCACCAAACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCGGACGGTGACCTTTGTTaccac       >  2:667016/1‑100 (MQ=255)

                                                                                             tCACCAGACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCGGACGGTGACCTTTGTTaccacc      >  1:3500900/1‑100 (MQ=255)

                                                                                              cACCAAACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCGGACGGTGACCTTTGTTACCACCt     >  2:3822681/1‑100 (MQ=255)

                                                                                              cACCAAACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCGGACGGTGACCTTTGTTACCACCt     <  1:1657615/100‑1 (MQ=255)

                                                                                               aCCAAACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACAGCATTTCACTGCCGGACGGTGACCTTTGTTACCACCtt    >  1:1454628/1‑100 (MQ=255)

                                                                                                 cAGACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCGGACGGTGACCTTTGTTACCACCTTCg  <  2:2060521/100‑1 (MQ=255)

                                                                                                 cAGACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCGGACGGTGACCTTTGTTACCACCTTCg  <  2:2870232/100‑1 (MQ=255)

                                                                                                 cAGACACATTTCAGAAAGAATATCCTGAGCCGAGGCTACCGGTGGAACCGGGAAGTAACCCCCTTTCACTGCCGGACGGTGACCTTTGTTACCACCTTCg  <  2:1293492/100‑1 (MQ=255)

TGGTCATGGTATTGAAGCGGGTAGCCACTTCGTTCTGACCAGCAGTCGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAGACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCGGACGGTGACCTTTGTTACCACCTTCG  >  NC_000913/4057344‑4057540

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