BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,980,709	A→G	M191T (ATG→ACG)	kduD ←	2‑deoxy‑D‑gluconate 3‑dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,980,709	0	A	G	85.5%	668.5 / 102.5	283	M191T (ATG→ACG)	kduD	2‑deoxy‑D‑gluconate 3‑dehydrogenase
Reads supporting (aligned to +/- strand): ref base A (20/21); new base G (120/122); total (140/143)
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.09e-03

Only 100 of 285 total aligned reads displayed.

GGCAGTCCCCAACGACCAGCTGGAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCATGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCATTCGTTCGCCATCAATCGCGTCACACCCATCACGCCGCTTTTTGA  >  NC_000913/2980619‑2980800

 gCAGTCCCCAACGACCAGCTGGAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCgtgt                                                                                          >  2:7128832/1‑93 (MQ=255)

 gCAGTCCCCAACGACCAGCTGGAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCATGt                                                                                          <  2:6191771/93‑1 (MQ=255)

  cAGTCCCCAACGACCAGCTGGAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTa                                                                                         >  1:601661/1‑93 (MQ=255)

    gTCCCCAACGACCAGCTGGAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAAc                                                                                       >  1:2976037/1‑93 (MQ=255)

     tCCCCAACGACCAGCTGGAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAAcc                                                                                      <  1:5378848/93‑1 (MQ=255)

      ccccAACGACCAGCTGGAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAAccc                                                                                     <  1:866400/93‑1 (MQ=255)

       cccAACGACCAGCTGGAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCg                                                                                    <  2:3774531/93‑1 (MQ=255)

       cccAACGACCAGCTGGAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCATGTAACCCg                                                                                    <  1:7528770/93‑1 (MQ=255)

         cAACGACCAGCTGGAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCggg                                                                                  >  2:1292146/1‑93 (MQ=255)

         cAACGACCAGCTGGAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCggg                                                                                  >  1:6931664/1‑93 (MQ=255)

         cAACGACCAGCTGGAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCggg                                                                                  <  1:5808500/93‑1 (MQ=255)

         cAACGACCAGCTGGAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCggg                                                                                  <  2:1170911/93‑1 (MQ=255)

         cAACGACCAGCTGGAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCggg                                                                                  <  2:2458844/93‑1 (MQ=255)

         cAACGACCAGCTGGAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCggg                                                                                  <  2:4713142/93‑1 (MQ=255)

          aaCGACCAGCTGGAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCgggg                                                                                 >  1:7071447/1‑93 (MQ=255)

          aaCGACCAGCTGGAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCgggg                                                                                 <  2:688273/93‑1 (MQ=255)

           aCGACCAGCTGGAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGc                                                                                >  1:1014795/1‑93 (MQ=255)

           aCGACCAGCTGGAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGc                                                                                >  1:4515612/1‑93 (MQ=255)

             gACCAGCTGGAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTa                                                                              >  1:4688756/1‑93 (MQ=255)

              aCCAGCTGGAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTAt                                                                             <  2:2483107/93‑1 (MQ=255)

               ccAGCTGGAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATc                                                                            <  1:982415/93‑1 (MQ=255)

               ccAGCTGGAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATc                                                                            <  2:77879/93‑1 (MQ=255)

               ccAGCTGGAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATc                                                                            >  1:1686931/1‑93 (MQ=255)

               ccAGCTGGAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCATGTAACCCGGGGCTATc                                                                            >  2:4348524/1‑93 (MQ=255)

                 aGCTGGAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGc                                                                          >  1:1273666/1‑93 (MQ=255)

                 aGCTGGAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGc                                                                          <  2:541679/93‑1 (MQ=255)

                 aGCTGGAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGc                                                                          >  1:7363910/1‑93 (MQ=255)

                   cTGGAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGcat                                                                        <  2:2754791/93‑1 (MQ=255)

                   cTGGAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGcat                                                                        >  2:3906000/1‑93 (MQ=255)

                    tGGAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGcatt                                                                       >  2:4494697/1‑93 (MQ=255)

                    tGGAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGcatt                                                                       <  2:2886953/93‑1 (MQ=255)

                    tGGAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGcatt                                                                       >  2:3585713/1‑93 (MQ=255)

                    tGGAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGcann                                                                       >  2:1395654/1‑91 (MQ=255)

                     ggAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGcatta                                                                      >  1:4104688/1‑93 (MQ=255)

                      gAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGcattaa                                                                     >  2:1827663/1‑93 (MQ=255)

                       aaTGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGcattaac                                                                    <  1:2454613/93‑1 (MQ=255)

                       aaTGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGcattaac                                                                    <  1:3192827/93‑1 (MQ=255)

                       aaTGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCATGTAACCCGGGGCTATCGcattaac                                                                    <  1:1904757/93‑1 (MQ=255)

                         tGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGcattaacat                                                                  <  2:3668842/93‑1 (MQ=255)

                          gCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGcattaacatt                                                                 <  1:5003052/93‑1 (MQ=255)

                           cGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGcattaacatta                                                                >  2:43376/1‑93 (MQ=255)

                               cGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAAtgt                                                            >  2:3351568/1‑93 (MQ=255)

                               cGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAAtgt                                                            <  2:1759319/93‑1 (MQ=255)

                                gagaATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAAtgtt                                                           <  2:3199733/93‑1 (MQ=255)

                                gagaATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAAtgtt                                                           >  1:7023136/1‑93 (MQ=255)

                                  gaATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAAtgttgt                                                         <  1:2273269/93‑1 (MQ=255)

                                  gaATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAAtgttgt                                                         <  1:1088919/93‑1 (MQ=255)

                                    aTTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTgtgt                                                       <  2:6552647/93‑1 (MQ=255)

                                     tttCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGtt                                                      <  1:6807952/93‑1 (MQ=255)

                                     tttCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGtt                                                      <  2:6712880/93‑1 (MQ=255)

                                     tttCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGtt                                                      <  1:6339532/93‑1 (MQ=255)

                                       tCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTa                                                    >  1:5091405/1‑93 (MQ=255)

                                       tCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTa                                                    <  2:4591859/93‑1 (MQ=255)

                                           cgcTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGccc                                                >  1:2809551/1‑93 (MQ=255)

                                            gcTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCa                                               >  1:5248336/1‑93 (MQ=255)

                                            gcTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCa                                               >  2:21537/1‑93 (MQ=255)

                                             cTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCAt                                              >  1:3875411/1‑93 (MQ=255)

                                               aCGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCAttc                                            >  1:2163644/1‑93 (MQ=255)

                                               aCGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTAGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTTTTTTTAAGCCAttc                                            >  1:6175786/1‑93 (MQ=255)

                                                cGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCAttcg                                           >  2:2716693/1‑93 (MQ=255)

                                                 gttgttCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCAttcgt                                          <  2:5707191/93‑1 (MQ=255)

                                                  ttgttCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCAttcgtt                                         <  1:4276531/93‑1 (MQ=255)

                                                  ttgttCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCAttcgtt                                         <  2:48692/93‑1 (MQ=255)

                                                   tgttCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCATGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCAttcgttc                                        >  1:4461200/1‑93 (MQ=255)

                                                    gttCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCAttcgttcg                                       >  1:379881/1‑93 (MQ=255)

                                                        aTCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCATTCGTTCGCCAt                                   >  1:2241231/1‑93 (MQ=255)

                                                          cTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCATGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCATTCGTTCGCCATCa                                 >  1:2065573/1‑93 (MQ=255)

                                                           tGCCCGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCTTTCTTTGCCCTTCaa                                >  1:1462936/1‑93 (MQ=255)

                                                             cccGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCATTCGTTCGCCATCAATc                              <  1:6624203/93‑1 (MQ=255)

                                                             cccGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCATTCGTTCGCCATCAATc                              <  1:3432669/93‑1 (MQ=255)

                                                             cccGTAGTTGTTGAGTATTGTTGGTCGCCATGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCATTCGTTCGCCATCAATc                              <  2:723564/93‑1 (MQ=255)

                                                              ccGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCATTCGTTCGCCATCAATcg                             >  1:7411128/1‑93 (MQ=255)

                                                               cGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCATTCGTTCGCCATCAATcgc                            <  2:4046336/93‑1 (MQ=255)

                                                               cGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCATTCGTTCGCCATCAATcgc                            <  1:1985092/93‑1 (MQ=255)

                                                               cGTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCATTCGTTCGCCATCAATcgc                            <  1:4027375/93‑1 (MQ=255)

                                                                gTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTATCATTAATGTTGTGTTTAGCCCATTCGTTCGCCATCAATcgcg                           <  2:5956405/93‑1 (MQ=255)

                                                                gTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCATTCGTTCGCCATCAATcgcg                           >  1:4881627/1‑93 (MQ=255)

                                                                gTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCATTCGTTCGCCATCAATcgcg                           >  2:3822567/1‑93 (MQ=255)

                                                                gTAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCATTCGTTCGCCATCAATcgcg                           >  2:5195184/1‑93 (MQ=255)

                                                                 tAGTTGTTGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCATTCGTTCGCCATCAATCGCGt                          <  1:3025819/93‑1 (MQ=255)

                                                                   gttgttGAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCATTCGTTCGCCATCAATCGCGTca                        <  1:1146797/93‑1 (MQ=255)

                                                                    ttgttgAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCATTCGTTCGCCATCAATCGCGTcac                       <  1:5857968/93‑1 (MQ=255)

                                                                    ttgttgAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCATTCGTTCGCCATCAATCGCGTcac                       >  1:5156778/1‑93 (MQ=255)

                                                                    ttgttgAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCATTCGTTCGCCATCAATCGCGTcac                       >  1:5971981/1‑93 (MQ=255)

                                                                      gttgAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCATTCGTTCGCCATCAATCGCGTCacac                     >  1:6374980/1‑93 (MQ=255)

                                                                       ttgAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCATTCGTTCGCCATCAATCGCGTCACAcc                    <  1:1849569/93‑1 (MQ=255)

                                                                       ttgAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCATTCGTTCGCCATCAATCGCGTCACAcc                    >  1:4351971/1‑93 (MQ=255)

                                                                         gAGTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCATTCGTTCGCCATCAATCGCGTCACACCCa                  <  1:54558/93‑1 (MQ=255)

                                                                           gTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCATTCGTTCGCCATCAATCGCGTCACACCCATc                <  2:232894/93‑1 (MQ=255)

                                                                           gTATTGTTGGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCATTCGTTCGCCATCAATCGCGTCACACCCATc                >  2:6169464/1‑93 (MQ=255)

                                                                            tATTGTTGGTCGCCATGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCATTCGTTCGCCATCAATCGCGTCACACCCATCa               <  1:3582647/93‑1 (MQ=255)

                                                                               tgttgGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCATTCGTTCGCCATCAATCGCGTCACACCCATCAcgc            <  2:5482077/93‑1 (MQ=255)

                                                                               tgttgGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCATTCGTTCGCCATCAATCGCGTCACACCCATCAcgc            <  2:4225411/93‑1 (MQ=255)

                                                                               tgttgGTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCATTCGTTCGCCATCAATCGCGTCACACCCATCAcgc            >  2:7106095/1‑93 (MQ=255)

                                                                                gttgGTCACCATGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCATTCGTTCGCCATCAATCGCGTCACACCCATCAcgcc           <  1:5917465/93‑1 (MQ=255)

                                                                                   ggTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCATTCGTTCGCCATCAATCGCGTCACACCCATCACGCCGCt        >  1:7384020/1‑93 (MQ=255)

                                                                                    gTCGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCATTCGTTCGCCATCAATCGCGTCACACCCATCACGCCGCtt       >  1:2506668/1‑93 (MQ=255)

                                                                                      cGCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCATTCGTTCGCCATCAATCGCGTCACACCCATCACGCCGCtttt     <  2:5018679/93‑1 (MQ=255)

                                                                                       gCCGTGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCATTCGTTCGCCATCAATCGCGTCACACCCATCACGCCGCttttt    >  2:6072473/1‑93 (MQ=255)

                                                                                         cATGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCATTCGTTCGCCATCAATCGCGTCACACCCATCACGCCGCTTTTTGa  <  1:2921870/93‑1 (MQ=255)

GGCAGTCCCCAACGACCAGCTGGAATGCGGTCGAGAATTTCCGCGCTACGTTGTTCATCTGCCCGTAGTTGTTGAGTATTGTTGGTCGCCATGTAACCCGGGGCTATCGCATTAACATTAATGTTGTGTTTAGCCCATTCGTTCGCCATCAATCGCGTCACACCCATCACGCCGCTTTTTGA  >  NC_000913/2980619‑2980800

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