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	87.1%	703.3 / 88.7	287	T147T (ACT→ACC)	kdgK	ketodeoxygluconokinase
Reads supporting (aligned to +/- strand): ref base T (21/16); new base C (118/132); total (139/148)
Fisher's exact test for biased strand distribution p-value = 2.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.05e-01

Only 100 of 288 total aligned reads displayed.

CGGCATTAACCGTTCATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACTTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTACGTTCTACTACTGGCG  >  NC_000913/3677179‑3677359

 ggCATTAACCGTTCATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACTTc                                                                                         <  1:606463/93‑1 (MQ=255)

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

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

  gCATTAACCGTTCATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTcc                                                                                        >  1:2676106/1‑93 (MQ=255)

   cATTAACCGTTCATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACTTccc                                                                                       <  1:1733945/93‑1 (MQ=255)

    aTTAACCGTTCATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCt                                                                                      <  1:2471913/93‑1 (MQ=255)

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

     ttAACCGTTCATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTg                                                                                     >  1:1700151/1‑93 (MQ=255)

        aCCGTTCATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACTTCCCTGAcc                                                                                  >  1:1003385/1‑93 (MQ=255)

        aCCGTTCATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGAcc                                                                                  >  1:5265751/1‑93 (MQ=255)

        aCCGTTCATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGAcc                                                                                  <  2:4250061/93‑1 (MQ=255)

        aCCGTTCATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGAcc                                                                                  <  2:1501513/93‑1 (MQ=255)

         ccGTTCATTACGAAACGGCGCTGGGAACGGACAGTTTTAGCCACCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGAccc                                                                                 >  1:1242506/1‑93 (MQ=255)

          cGTTCATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACTTCCCTGACCCa                                                                                >  1:2777580/1‑93 (MQ=255)

          cGTTCATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCa                                                                                >  2:5278110/1‑93 (MQ=255)

           gTTCATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACTTCCCTGACCCaa                                                                               >  1:5331206/1‑93 (MQ=255)

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

             tCATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACg                                                                             <  1:532507/93‑1 (MQ=255)

               aTTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGa                                                                           <  2:2769049/93‑1 (MQ=255)

               aTTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGa                                                                           <  2:3902413/93‑1 (MQ=255)

               aTTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGa                                                                           >  1:2580344/1‑93 (MQ=255)

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

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

                ttACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGAt                                                                          <  2:3410100/93‑1 (MQ=255)

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

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

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

                      aaCGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGaaaa                                                                    <  2:3147282/93‑1 (MQ=255)

                          gcgcTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGt                                                                <  1:3179046/93‑1 (MQ=255)

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

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

                            gcTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCt                                                              <  1:2362456/93‑1 (MQ=255)

                             cTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTg                                                             <  2:6710542/93‑1 (MQ=255)

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

                               gggAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATCCCTCCCTGACCCAACGGATGGAAAACCGTCTGcc                                                           <  1:7453164/93‑1 (MQ=255)

                               gggAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGcc                                                           <  1:4438321/93‑1 (MQ=255)

                               gggAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGcc                                                           <  2:7191749/93‑1 (MQ=255)

                                ggAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCg                                                          <  1:6859972/93‑1 (MQ=255)

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

                                 gAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACTTCCCTGACCCAACGGATGGAAAACCGTCTGCCgg                                                         >  1:942568/1‑93 (MQ=255)

                                  aaCGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACTTCCCTGACCCAACGGATGGAAAACCGTCTGCCggg                                                        >  2:1157566/1‑93 (MQ=255)

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

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

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

                                           ttttAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGACTGCCGGGCCTTtacta                                               >  2:5248857/1‑93 (MQ=255)

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

                                              tAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACTTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACAt                                            <  1:3682745/93‑1 (MQ=255)

                                              tAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACAt                                            >  1:5133155/1‑93 (MQ=255)

                                                 ccAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGCCCCAACGGATGGAAAACCGTCTCCCGGGCCTTTCTTACAttta                                         >  1:6561937/1‑91 (MQ=255)

                                                 ccAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGa                                         <  2:469106/93‑1 (MQ=255)

                                                  cagcagATGCTGGACGCCTGGCACGGCGAGAACGTTGATACTTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGaa                                        <  2:4435149/93‑1 (MQ=255)

                                                  cagcagATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGaa                                        >  1:5534331/1‑93 (MQ=255)

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

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

                                                    gcagATGCTGGACGCCTGGCACGGCGAGAACGTTGATACTTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAAc                                      <  1:2886808/93‑1 (MQ=255)

                                                    gcagATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAAc                                      <  1:6156762/93‑1 (MQ=255)

                                                     cagATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAAcc                                     >  1:5405944/1‑93 (MQ=255)

                                                      agATGCTGGACGCCTGGCACGGCGAGCACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCg                                    >  1:4032853/1‑93 (MQ=255)

                                                       gATGCTGGACGCCTGGCACGGCGAGAACGTTGATACTTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGa                                   <  2:2295158/93‑1 (MQ=255)

                                                       gATGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGa                                   >  2:6271038/1‑93 (MQ=255)

                                                        aTGCTGGACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGAc                                  <  2:3308211/93‑1 (MQ=255)

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

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

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

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

                                                             ggACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCAc                             >  1:4727083/1‑93 (MQ=255)

                                                             ggACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCAc                             <  2:645283/93‑1 (MQ=255)

                                                             ggACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCAc                             <  1:5783783/93‑1 (MQ=255)

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

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

                                                              gACGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCAcc                            >  1:1515903/1‑93 (MQ=255)

                                                               aCGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCg                           >  2:2013178/1‑93 (MQ=255)

                                                               aCGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCg                           >  2:3033825/1‑93 (MQ=255)

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

                                                                cGCCTGGCACGGCGAGAACGTTGATACTTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCgg                          <  1:6629478/93‑1 (MQ=255)

                                                                cGCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCgg                          <  2:736438/93‑1 (MQ=255)

                                                                 gCCTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGc                         >  1:1079709/1‑93 (MQ=255)

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

                                                                   cTGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGa                       <  1:7068948/93‑1 (MQ=255)

                                                                    tGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAg                      <  1:3928478/93‑1 (MQ=255)

                                                                    tGGCACGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAg                      >  1:5676499/1‑93 (MQ=255)

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

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

                                                                         cGGCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTAc                 >  1:22652/1‑93 (MQ=255)

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

                                                                          ggCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTACg                <  1:4310657/93‑1 (MQ=255)

                                                                           gCGAGAACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTACGt               >  2:4835583/1‑93 (MQ=255)

                                                                            cGAGAACGTTGATACTTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTACGtt              <  2:1358858/93‑1 (MQ=255)

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

                                                                               gaACGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTACGTTcta           <  2:3437753/93‑1 (MQ=255)

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

                                                                                aaCGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTACGTTctac          <  1:4807541/93‑1 (MQ=255)

                                                                                 aCGTTGATACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTACGTTctact         >  2:3614559/1‑93 (MQ=255)

                                                                                   gTTGATACTTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTACGTTCtactac       <  2:4933261/93‑1 (MQ=255)

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

                                                                                   gTTGATACCGCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTACGTTCtactac       <  1:7251739/93‑1 (MQ=255)

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

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

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

                                                                                        tACCTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTACGTTCTACTACTGGCg  <  2:2937976/93‑1 (MQ=255)

CGGCATTAACCGTTCATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACTTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTACGTTCTACTACTGGCG  >  NC_000913/3677179‑3677359

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