BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,677,270	T→C	T147T (ACT→ACC)	kdgK →	ketodeoxygluconokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,677,270	0	T	C	99.5%	1378.9 / ‑3.6	406	T147T (ACT→ACC)	kdgK	ketodeoxygluconokinase
Reads supporting (aligned to +/- strand): ref base T (2/0); new base C (188/215); total (191/215)
Fisher's exact test for biased strand distribution p-value = 2.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.19e-01

Only 100 of 420 total aligned reads displayed.

GGCATTAACCGTTCATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACTTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTACGTTCTACTACTGGCG  >  NC_000913/3677180‑3677359

ggCATTAACCGTTCATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTc                                                                                         <  1:3963600/93‑1 (MQ=255)

ggCATTAACCGTTCATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTc                                                                                         >  1:822356/1‑93 (MQ=255)

ggCATTAACCGTTCATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTc                                                                                         >  2:1970999/1‑93 (MQ=255)

 gCATTAACCGTTCATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTcc                                                                                        <  2:3882691/93‑1 (MQ=255)

  cATTAACCGTTCATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTccc                                                                                       <  1:3733498/93‑1 (MQ=255)

   aTTAACCGTTCATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCt                                                                                      >  1:7673880/1‑93 (MQ=255)

     tAACCGTTCATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGa                                                                                    <  2:7720960/93‑1 (MQ=255)

        ccGTTCATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGAccc                                                                                 >  1:1695805/1‑93 (MQ=255)

        ccGTTCATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGAccc                                                                                 <  2:6732528/93‑1 (MQ=255)

        ccGTTCATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGAccc                                                                                 >  2:158758/1‑93 (MQ=255)

        ccGTTCATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGACACCTCCCTGAccc                                                                                 >  1:6500008/1‑93 (MQ=255)

          gTTCATTGCGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCaa                                                                               <  2:2773341/93‑1 (MQ=255)

          gTTCATTACGTAACGGCGCTGGGAGCGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCaa                                                                               <  1:6260360/93‑1 (MQ=255)

          gTTCATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCaa                                                                               <  2:5987340/93‑1 (MQ=255)

           ttCATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAAc                                                                              <  1:5318522/93‑1 (MQ=255)

           ttCATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAAc                                                                              >  1:1222274/1‑93 (MQ=255)

             taTTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACgg                                                                            <  1:1076251/92‑1 (MQ=255)

             cATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACgg                                                                            <  2:880216/93‑1 (MQ=255)

             cATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACgg                                                                            >  2:5780531/1‑93 (MQ=255)

             cATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACgg                                                                            >  2:8434887/1‑93 (MQ=255)

              aTTACGTAACGGCGCTGGGAACGGATAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGa                                                                           <  1:6046469/93‑1 (MQ=255)

              aTTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGAACGGCGAGAACGTTGATACCTCCCTGACCCAACGGa                                                                           <  1:4881015/93‑1 (MQ=255)

               ttACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGAt                                                                          >  1:3358044/1‑93 (MQ=255)

                tACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATg                                                                         <  2:8056240/93‑1 (MQ=255)

                 aCGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATgg                                                                        <  1:7011238/93‑1 (MQ=255)

                 aCGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATgg                                                                        <  1:3247232/93‑1 (MQ=255)

                 aCGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATgg                                                                        <  1:8741624/93‑1 (MQ=255)

                  cGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGa                                                                       <  1:5463745/93‑1 (MQ=255)

                  cGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGa                                                                       <  2:5173065/93‑1 (MQ=255)

                    tAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCCACGGATGGaaa                                                                     >  2:3588157/1‑93 (MQ=255)

                    tAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGaaa                                                                     >  1:6309080/1‑93 (MQ=255)

                         gcgcTGGGAACGGACAGTTGTAGCCAGCTGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGt                                                                <  2:4629522/93‑1 (MQ=255)

                          cgcTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTc                                                               <  1:844116/93‑1 (MQ=255)

                           gcTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCt                                                              >  2:5070305/1‑93 (MQ=255)

                            cTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTg                                                             >  1:3179147/1‑93 (MQ=255)

                             tGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGc                                                            <  1:1852197/93‑1 (MQ=255)

                             tGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGc                                                            <  1:2238676/93‑1 (MQ=255)

                              gggAACGCACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCCCCCTGACCCAACGGATGGAAAACCGTCTGcc                                                           <  2:1395587/93‑1 (MQ=255)

                               ggAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCg                                                          >  2:641038/1‑93 (MQ=255)

                                gAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCgg                                                         <  2:3274697/93‑1 (MQ=255)

                                 aaCGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCggg                                                        <  2:1477043/93‑1 (MQ=255)

                                   cGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGcc                                                      <  2:1732630/93‑1 (MQ=255)

                                   cGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGcc                                                      <  2:1267903/93‑1 (MQ=255)

                                    ggACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCt                                                     >  1:7620375/1‑93 (MQ=255)

                                    ggACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCt                                                     <  2:2970584/93‑1 (MQ=255)

                                       cAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTta                                                  <  2:4965786/93‑1 (MQ=255)

                                       cAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTta                                                  >  2:5847671/1‑93 (MQ=255)

                                       cAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTta                                                  >  2:8344828/1‑93 (MQ=255)

                                          ttttAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTtacta                                               <  2:5360651/93‑1 (MQ=255)

                                          ttttAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTtacta                                               >  2:6972128/1‑93 (MQ=255)

                                           tttAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTtactac                                              >  2:2289392/1‑93 (MQ=255)

                                            ttAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGCAAACCGTCTGGCGCGCCTTTACTACa                                             >  2:815763/1‑93 (MQ=255)

                                            ttAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACa                                             <  2:4058650/93‑1 (MQ=255)

                                            ttAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACa                                             <  2:7188866/93‑1 (MQ=255)

                                             tAGCCAGCAGATGCGGGACGCCTGTCACGGCGAGAACGTTGATACCTCCCTGCCCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACAt                                            <  2:2893203/93‑1 (MQ=255)

                                              aGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACAtt                                           <  2:413821/93‑1 (MQ=255)

                                               gCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCGGGGCCTTTACTACAtta                                          >  1:2397967/1‑92 (MQ=255)

                                               gCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTg                                          >  1:2832728/1‑93 (MQ=255)

                                               gCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTg                                          <  1:43223/93‑1 (MQ=255)

                                                ccAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGa                                         >  1:8631919/1‑93 (MQ=255)

                                                 cagcagATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGaa                                        <  2:1195000/93‑1 (MQ=255)

                                                 cagcagATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGaa                                        <  2:856138/93‑1 (MQ=255)

                                                  agcagATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGaaa                                       <  2:566318/93‑1 (MQ=255)

                                                  agcagATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGaaa                                       <  1:6820851/93‑1 (MQ=255)

                                                  agcagATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGaaa                                       >  1:4164744/1‑93 (MQ=255)

                                                  agcagATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGaaa                                       >  1:1570482/1‑93 (MQ=255)

                                                   gcagATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAAc                                      >  1:970140/1‑93 (MQ=255)

                                                     agATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCg                                    <  2:6827805/93‑1 (MQ=255)

                                                       aTGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGAc                                  >  1:2616401/1‑93 (MQ=255)

                                                        tGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACa                                 <  2:4798128/93‑1 (MQ=255)

                                                        tGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACa                                 <  1:3624052/93‑1 (MQ=255)

                                                         gCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAg                                <  2:31729/93‑1 (MQ=255)

                                                          cTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGc                               <  1:20790/93‑1 (MQ=255)

                                                           tGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCa                              >  1:5233173/1‑93 (MQ=255)

                                                           tGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCa                              <  2:616453/93‑1 (MQ=255)

                                                            ggACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACCGCAc                             <  1:7076857/93‑1 (MQ=255)

                                                             gACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCAcc                            <  2:4311429/93‑1 (MQ=255)

                                                              aCGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCg                           <  2:4373292/93‑1 (MQ=255)

                                                               cGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTCCCGGGCCTTTACTACATTGAAACCGACAGCACCgg                          <  2:3100667/93‑1 (MQ=255)

                                                                 ccTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCg                        >  2:7795042/1‑93 (MQ=255)

                                                                 ccTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCg                        >  2:1906958/1‑93 (MQ=255)

                                                                 ccTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCg                        >  2:2148844/1‑93 (MQ=255)

                                                                  cTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGa                       >  1:6207983/1‑93 (MQ=255)

                                                                    ggCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGc                     >  1:1994904/1‑93 (MQ=255)

                                                                    ggCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGc                     <  2:7606700/93‑1 (MQ=255)

                                                                       aCGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTa                  <  1:7896775/93‑1 (MQ=255)

                                                                       aCGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTa                  <  1:3612482/93‑1 (MQ=255)

                                                                       aCGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTa                  >  2:5292074/1‑93 (MQ=255)

                                                                        cGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTAc                 <  1:383791/93‑1 (MQ=255)

                                                                            gagaACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTACGTTc             >  1:5686582/1‑93 (MQ=255)

                                                                             agaACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTACGTTct            >  1:73395/1‑93 (MQ=255)

                                                                             agaACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTACGTTct            >  1:8580056/1‑93 (MQ=255)

                                                                              gaACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTACGTTcta           >  1:2299788/1‑93 (MQ=255)

                                                                              gaACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGGACGCTcta           >  1:4239635/1‑93 (MQ=255)

                                                                                 cGTTGATGCCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTACGTTctacta        <  1:6779183/93‑1 (MQ=255)

                                                                                  gTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTACGTTCtactac       <  1:5909779/93‑1 (MQ=255)

                                                                                  gTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGCCGAGCGTACGTTCtactac       >  2:7456099/1‑93 (MQ=255)

                                                                                   ttGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTACGTTCTactact      <  1:8076185/93‑1 (MQ=255)

                                                                                    tGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTACGTTCTACTACTg     <  2:8677714/93‑1 (MQ=255)

                                                                                      aTACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTACGTTCTACTACTGGc   <  1:4627473/93‑1 (MQ=255)

GGCATTAACCGTTCATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACTTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTACGTTCTACTACTGGCG  >  NC_000913/3677180‑3677359

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