BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,078,221	C→T	Q478Q (CAG→CAA)	tktA ←	transketolase 1 isozyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,078,221	0	C	T	99.7%	1107.0 / ‑5.5	309	Q478Q (CAG→CAA)	tktA	transketolase 1 isozyme
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (170/138); total (170/139)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

Only 100 of 322 total aligned reads displayed.

CTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGCCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGC  >  NC_000913/3078131‑3078310

 tCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGc                                                                                        >  1:3873991/1‑93 (MQ=255)

   aaCACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTc                                                                                      <  2:2127843/93‑1 (MQ=255)

   aaCACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTc                                                                                      >  2:5454718/1‑93 (MQ=255)

     cacCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCaa                                                                                    >  2:7737304/1‑93 (MQ=255)

      acCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCacc                                                                                   >  2:1615720/1‑91 (MQ=255)

      acCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAAc                                                                                   >  2:700236/1‑93 (MQ=255)

      acCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAAc                                                                                   <  2:3019285/93‑1 (MQ=255)

       ccGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGCCTTGCTCAAcg                                                                                  >  2:3712346/1‑92 (MQ=255)

        cGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCg                                                                                 <  1:4262268/93‑1 (MQ=255)

         gTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCgg                                                                                >  1:5459786/1‑93 (MQ=255)

         gTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCgg                                                                                >  1:1214707/1‑93 (MQ=255)

          tATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGc                                                                               >  1:5387293/1‑93 (MQ=255)

           aTTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCt                                                                              <  2:6200221/93‑1 (MQ=255)

           aTTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCt                                                                              <  2:7679456/93‑1 (MQ=255)

            tttCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTg                                                                             >  1:6164688/1‑93 (MQ=255)

            tttCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTg                                                                             >  2:2872340/1‑93 (MQ=255)

            ttgccACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTg                                                                             <  1:1986553/90‑1 (MQ=255)

              tCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGt                                                                           <  1:549884/93‑1 (MQ=255)

              tCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGt                                                                           <  1:5512473/93‑1 (MQ=255)

                cACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGtaa                                                                         >  1:2754967/1‑91 (MQ=255)

                cACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGa                                                                         <  2:3431885/93‑1 (MQ=255)

                cACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGa                                                                         >  1:6282273/1‑93 (MQ=255)

                 aCGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAg                                                                        >  1:4186628/1‑93 (MQ=255)

                  cgcgACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGt                                                                       <  1:3681189/93‑1 (MQ=255)

                    cgACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCg                                                                     <  1:4786350/93‑1 (MQ=255)

                    cgACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCg                                                                     <  1:2251385/93‑1 (MQ=255)

                     gACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCgg                                                                    >  2:4701155/1‑93 (MQ=255)

                      aCCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCggg                                                                   >  1:3357524/1‑92 (MQ=255)

                        cgcgGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGgcc                                                                 <  2:132601/93‑1 (MQ=255)

                         gcgGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCggggcg                                                                >  2:3575246/1‑89 (MQ=255)

                         gcgGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCg                                                                >  2:4233541/1‑93 (MQ=255)

                          cgGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGt                                                               <  2:2522443/93‑1 (MQ=255)

                           ggATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTc                                                              <  1:6944616/93‑1 (MQ=255)

                           ggATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTc                                                              >  2:6567852/1‑93 (MQ=255)

                           ggATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTc                                                              <  1:1438616/93‑1 (MQ=255)

                           ggATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTc                                                              >  1:1108607/1‑93 (MQ=255)

                           ggATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTc                                                              >  1:5728519/1‑93 (MQ=255)

                            gATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCt                                                             <  2:810871/93‑1 (MQ=255)

                             aTTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCtt                                                            >  1:2490958/1‑93 (MQ=255)

                                  aCCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGccc                                                       <  1:236775/93‑1 (MQ=255)

                                  aCCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGccc                                                       <  1:219506/93‑1 (MQ=255)

                                   ccTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCa                                                      <  1:7057306/93‑1 (MQ=255)

                                      ggTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGAc                                                   <  2:1064709/93‑1 (MQ=255)

                                         cacaCGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGa                                                <  2:6022947/93‑1 (MQ=255)

                                         cacaCGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGa                                                <  1:4697817/93‑1 (MQ=255)

                                          acacGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGGCGTCTTCGCCCAGACAGAt                                               >  1:2888376/1‑93 (MQ=255)

                                          acacGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGag                                               >  1:3615097/1‑92 (MQ=255)

                                          acacGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGAt                                               <  1:4967448/93‑1 (MQ=255)

                                           cacGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATg                                              >  1:1566853/1‑93 (MQ=255)

                                           cacGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATg                                              <  1:253827/93‑1 (MQ=255)

                                           cacGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATg                                              >  1:7267344/1‑93 (MQ=255)

                                            acGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATgg                                             <  1:2588249/93‑1 (MQ=255)

                                             cGGCCGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGCCCGATGGa                                            <  2:3320493/93‑1 (MQ=255)

                                              ggACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGGCTTCGCCCAGACCGATGGAg                                           >  1:4491460/1‑93 (MQ=255)

                                               gACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATggcgg                                          >  2:7186760/1‑90 (MQ=255)

                                               gACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGt                                          <  1:3568759/93‑1 (MQ=255)

                                                   ccATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTg                                      <  2:2956483/93‑1 (MQ=255)

                                                   ccATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTg                                      <  2:5663404/93‑1 (MQ=255)

                                                   ccATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTg                                      >  2:7421555/1‑93 (MQ=255)

                                                    cATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTgg                                     >  2:2243768/1‑93 (MQ=255)

                                                    cATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTgg                                     <  2:5183715/93‑1 (MQ=255)

                                                      tGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGgt                                   <  1:7409689/93‑1 (MQ=255)

                                                       gTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGgtg                                  >  1:3260428/1‑93 (MQ=255)

                                                        tAGACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGgtgt                                 <  1:6014126/93‑1 (MQ=255)

                                                          gACATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGCAGTGGGTGTaa                               >  1:1018941/1‑93 (MQ=255)

                                                           aCATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTaaa                              <  1:89769/93‑1 (MQ=255)

                                                            cATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAAc                             <  2:3983896/93‑1 (MQ=255)

                                                            cATGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAAc                             >  2:4310976/1‑93 (MQ=255)

                                                              tGTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCa                           >  2:6428574/1‑93 (MQ=255)

                                                               gtgcGGGGTTACTCGCAGAGAAGCGACTTGCTCATCCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCAt                          <  2:1872606/90‑1 (MQ=255)

                                                               gTTCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCAt                          <  2:2679173/93‑1 (MQ=255)

                                                                 tCGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTTGTGGGTGTAAACCATCa                        <  2:5798280/93‑1 (MQ=255)

                                                                  cGGGGTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCAc                       >  2:3681225/1‑93 (MQ=255)

                                                                   ggggTTACGCGGAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCAcc                      <  2:3905908/93‑1 (MQ=255)

                                                                   ggggTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCAcc                      <  1:3016701/93‑1 (MQ=255)

                                                                    gggTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCt                     <  1:1909594/93‑1 (MQ=255)

                                                                    gggTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCt                     >  1:6743775/1‑93 (MQ=255)

                                                                     ggTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTg                    >  2:1519733/1‑93 (MQ=255)

                                                                      gTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGa                   <  1:7771525/93‑1 (MQ=255)

                                                                      gTTACGCGCAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGa                   >  2:5941968/1‑93 (MQ=255)

                                                                              cAGAGAAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGttt           <  1:4856110/93‑1 (MQ=255)

                                                                               agagAAGCGACTTGCTCAACCGGCTGGGGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGACGTGGGTGTAAACCATCACCTGACGCTGTTtc          >  1:652465/1‑93 (MQ=255)

                                                                                gagaAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTtca         >  2:6883309/1‑93 (MQ=255)

                                                                                gagaAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTtca         <  1:5130120/93‑1 (MQ=255)

                                                                                gagaAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTtca         <  2:5563326/93‑1 (MQ=255)

                                                                                gagaAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTtca         <  2:1957146/93‑1 (MQ=255)

                                                                                 agaAGCGCCTTGCTCAACCGGGTGGTGAGTCGGGCCGTCTTCCCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTtcat        <  2:1566088/93‑1 (MQ=255)

                                                                                 agaAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTtcat        <  2:6131921/93‑1 (MQ=255)

                                                                                 agaAGCGACTTGCTCAACCGGCTGGGGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTtcat        >  1:1306811/1‑93 (MQ=255)

                                                                                  gaAGCGACTTGCTCATCCGGCTGGTGAGTCGGGCCTTCTTCGCCCAGACCGATGGAGTCGTGTGTGTAACCCATCACCTGACGCTGTTtcatc       <  2:3752118/93‑1 (MQ=255)

                                                                                  gaAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTtcatc       >  2:4658202/1‑93 (MQ=255)

                                                                                  gaAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTtcatc       <  1:4409859/93‑1 (MQ=255)

                                                                                  gaAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTtcatc       >  2:7358023/1‑93 (MQ=255)

                                                                                  gaAGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTtcatc       <  1:1052375/93‑1 (MQ=255)

                                                                                   aaGCGACTTGTTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTtcatca      <  2:7525087/93‑1 (MQ=255)

                                                                                   aaGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCATCTGACGCTGTTtcatca      >  2:4776393/1‑93 (MQ=255)

                                                                                   aaGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTtcatca      <  2:4461084/93‑1 (MQ=255)

                                                                                    aGCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTtcatctg     >  2:1381149/1‑91 (MQ=255)

                                                                                     gCGACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAgc    <  2:6643455/93‑1 (MQ=255)

                                                                                       gACTTGCTCAACCGGCTGGTGAGTCGGGCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAgcgc  >  2:7808492/1‑93 (MQ=255)

CTCAACACCGTATTTCCACGCGACCGCGGATTCAACCTGGTCACACGGACGCCATGTAGACATGTTCGGGGTTACGCGCAGAGAAGCGACCTGCTCAACCGGCTGGTGAGTCGGCCCGTCTTCGCCCAGACCGATGGAGTCGTGGGTGTAAACCATCACCTGACGCTGTTTCATCAGCGC  >  NC_000913/3078131‑3078310

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