BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,263,134	A→G	V179A (GTG→GCG)	tdcB ←	threonine dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,263,134	0	A	G	97.9%	602.6 / 0.1	192	V179A (GTG→GCG)	tdcB	threonine dehydratase
Reads supporting (aligned to +/- strand): ref base A (3/1); new base G (106/82); total (109/83)
Fisher's exact test for biased strand distribution p-value = 6.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 201 total aligned reads displayed.

TGTACGCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATA  >  NC_000913/3263044‑3263223

tggaCGCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCCCAATcgcg                                                                                         <  2:3376368/90‑1 (MQ=255)

 gTACGCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGt                                                                                        <  2:1317097/93‑1 (MQ=255)

    cgcgAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTAt                                                                                     <  1:4236069/93‑1 (MQ=255)

     gCCAATAACACGAATGGTGGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATc                                                                                    <  1:4496705/93‑1 (MQ=255)

     gCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATc                                                                                    <  2:2832006/93‑1 (MQ=255)

      ccAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCg                                                                                   >  2:1036662/1‑93 (MQ=255)

       cAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGa                                                                                  <  1:2146762/93‑1 (MQ=255)

       cAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGa                                                                                  >  2:3027671/1‑93 (MQ=255)

       cAATAACACGAATGGGCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTACACCGCCACCACCAATTGGCACAATCACGTTATCgc                                                                                  >  2:2331448/1‑92 (MQ=255)

        aaTAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGAc                                                                                 >  2:2886838/1‑93 (MQ=255)

           aaCACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGAcatc                                                                              <  1:3750278/93‑1 (MQ=255)

           aaCACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGAcatc                                                                              >  2:1899341/1‑93 (MQ=255)

             cacGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGAcatcat                                                                            >  2:107997/1‑93 (MQ=255)

             cacGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGAcatcat                                                                            >  1:3772863/1‑93 (MQ=255)

                gAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATaga                                                                         >  1:859354/1‑93 (MQ=255)

                gAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATaga                                                                         <  1:1997500/93‑1 (MQ=255)

                 aaTGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATagag                                                                        >  1:1208880/1‑93 (MQ=255)

                  aTGGTCGGGTTAATAGATTTAATTGCCACCTCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAgaga                                                                       <  1:2272453/93‑1 (MQ=255)

                  aTGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAgaga                                                                       <  1:1615961/93‑1 (MQ=255)

                  aTGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAgaga                                                                       <  1:3827011/93‑1 (MQ=255)

                  aTGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAgaga                                                                       >  1:2868594/1‑93 (MQ=255)

                   tGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGAt                                                                      <  1:79885/93‑1 (MQ=255)

                    ggTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATc                                                                     <  1:2544206/93‑1 (MQ=255)

                    ggTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATc                                                                     <  2:2336024/93‑1 (MQ=255)

                      tCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCtt                                                                   >  2:4586891/1‑93 (MQ=255)

                      tCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCtt                                                                   >  2:581514/1‑93 (MQ=255)

                       cGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCATAATCGCGTTATCGACATCATAGAGATCTTc                                                                  <  2:340900/93‑1 (MQ=255)

                       cGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTc                                                                  <  2:3103159/93‑1 (MQ=255)

                        gggTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACACCATAGAGATCTTcc                                                                 <  2:3117880/93‑1 (MQ=255)

                          gTTAATAGATTTAATTGCCCCCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGCCATCATATAGACCTTCCAt                                                               >  1:3598151/1‑93 (MQ=255)

                          gTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCAt                                                               >  2:2030290/1‑93 (MQ=255)

                          gTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTACAt                                                               <  1:695159/93‑1 (MQ=255)

                            tAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATaa                                                             <  1:1559588/93‑1 (MQ=255)

                              aTAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAAtt                                                           <  2:1894426/93‑1 (MQ=255)

                              aTAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAAtt                                                           >  1:3242406/1‑93 (MQ=255)

                              aTAGATTTAATTGCCACAGCAATACCAGAAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAAtt                                                           <  1:4207252/93‑1 (MQ=255)

                               tAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAAttt                                                          >  1:4378696/1‑93 (MQ=255)

                               tAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAAttt                                                          >  2:4503557/1‑93 (MQ=255)

                                 gATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTcc                                                        >  2:1985184/1‑93 (MQ=255)

                                 gATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTcc                                                        <  2:4639361/93‑1 (MQ=255)

                                   tttAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAg                                                      >  2:4214307/1‑93 (MQ=255)

                                    ttAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGa                                                     <  2:4361011/93‑1 (MQ=255)

                                    ttAATGGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGa                                                     <  2:2638111/93‑1 (MQ=255)

                                     tAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGAc                                                    >  1:1987468/1‑93 (MQ=255)

                                       aTTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATAATAGAGCTCTTCCATAATTTCCAGACCa                                                  >  1:2587847/1‑93 (MQ=255)

                                         tGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAAt                                                <  2:4555447/93‑1 (MQ=255)

                                         tGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAAt                                                <  2:115733/93‑1 (MQ=255)

                                         tGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAAt                                                <  1:2479937/93‑1 (MQ=255)

                                           ccACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCg                                              <  2:3299164/93‑1 (MQ=255)

                                            cACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGt                                             <  1:12572/93‑1 (MQ=255)

                                            cACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGt                                             >  2:1468868/1‑93 (MQ=255)

                                             cccGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGtt                                            <  2:1592458/92‑1 (MQ=255)

                                             aCCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGtt                                            >  2:2765798/1‑93 (MQ=255)

                                             aCCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGtt                                            <  2:186838/93‑1 (MQ=255)

                                               ctcAATACCAGCAATTACACCGCCACCACCAATGGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTcc                                          <  1:308712/91‑1 (MQ=255)

                                                gCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTccc                                         >  2:968992/1‑93 (MQ=255)

                                                gCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTccc                                         <  1:1350096/93‑1 (MQ=255)

                                                gCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTccc                                         <  1:1445519/93‑1 (MQ=255)

                                                   aTACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTgg                                      <  1:1496206/93‑1 (MQ=255)

                                                   aTACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTgg                                      >  1:2620756/1‑93 (MQ=255)

                                                    tACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGc                                     >  1:3890024/1‑93 (MQ=255)

                                                     aCCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGcc                                    <  2:16386/93‑1 (MQ=255)

                                                      ccAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCa                                   >  1:2354715/1‑93 (MQ=255)

                                                      ccAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCa                                   >  1:2452649/1‑93 (MQ=255)

                                                       cAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAg                                  >  2:3652183/1‑93 (MQ=255)

                                                        aGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGc                                 >  2:236673/1‑93 (MQ=255)

                                                            aTTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATc                             <  1:2409375/93‑1 (MQ=255)

                                                              tAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCAc                           >  1:1270083/1‑93 (MQ=255)

                                                                aaCCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACtt                         >  2:1998575/1‑93 (MQ=255)

                                                                aaCCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGAACTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACtt                         <  2:2797728/93‑1 (MQ=255)

                                                                 aCCGCCACCCCCAAATGGCCACATCGGGTTTTCGACATTATAGGGATTTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACttt                        <  2:2355675/93‑1 (MQ=255)

                                                                 aCCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACttt                        >  1:3806106/1‑93 (MQ=255)

                                                                  ccGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAAAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCACTCACtttt                       >  2:3974182/1‑92 (MQ=255)

                                                                   cGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCg                      <  1:2102215/93‑1 (MQ=255)

                                                                   cGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCg                      >  1:1012008/1‑93 (MQ=255)

                                                                   cGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCg                      <  2:4100111/93‑1 (MQ=255)

                                                                   cGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCg                      <  2:3935336/93‑1 (MQ=255)

                                                                       accaccaATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGatc                  >  2:3534699/1‑93 (MQ=255)

                                                                        ccaccaATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGatca                 >  1:2182041/1‑93 (MQ=255)

                                                                          accaATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGatcatc               >  2:3897248/1‑93 (MQ=255)

                                                                          accaATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGatcatc               <  1:4472741/93‑1 (MQ=255)

                                                                           ccaATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCg              <  1:1795970/93‑1 (MQ=255)

                                                                             aaTTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTa            >  1:3418583/1‑93 (MQ=255)

                                                                              aTTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTaa           <  2:3072045/93‑1 (MQ=255)

                                                                                tGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAgg         <  2:2714714/93‑1 (MQ=255)

                                                                                 ggCACAATCGCGTTATCGACATCATAGATATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATAATCGTAAGGt        <  1:312671/93‑1 (MQ=255)

                                                                                 ggCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGt        <  2:2464813/93‑1 (MQ=255)

                                                                                  gCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTg       <  1:1404205/93‑1 (MQ=255)

                                                                                  gCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTg       <  2:4040650/93‑1 (MQ=255)

                                                                                   cacaATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTgg      <  2:3246138/93‑1 (MQ=255)

                                                                                   cacaATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTgg      <  2:3771181/93‑1 (MQ=255)

                                                                                   cacaATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTgg      <  2:2188523/93‑1 (MQ=255)

                                                                                    acaATCGCGTTATCGACATCATATAGATCTTCCATCATTTCCAGACCAATCGATCCCTGGCCAGCAATCACTTTCGGGTCATCGTAAGGTggg     >  2:3212659/1‑93 (MQ=255)

                                                                                    acaATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTggg     <  1:883691/93‑1 (MQ=255)

                                                                                      aaTCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGaa   >  2:1103678/1‑92 (MQ=255)

                                                                                      aaTCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGAt   >  1:3706676/1‑93 (MQ=255)

                                                                                      aaTCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGAt   >  2:1753959/1‑93 (MQ=255)

                                                                                      aaTCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGAt   >  1:2779435/1‑93 (MQ=255)

                                                                                       aTCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATa  >  1:295799/1‑93 (MQ=255)

                                                                                       aTCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATa  >  2:754255/1‑93 (MQ=255)

TGTACGCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATA  >  NC_000913/3263044‑3263223

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