BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,319,081	G→A	R340H (CGT→CAT)	atoS →	AtoS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,319,081	0	G	A	96.7%	871.8 / 8.7	240	R340H (CGT→CAT)	atoS	AtoS
Reads supporting (aligned to +/- strand): ref base G (3/5); new base A (118/114); total (121/119)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 246 total aligned reads displayed.

TTGTTCGACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCGTGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAACACGCACGGTGAAATGATAGGTGCTTTGGTGATTTTCTCTGATT  >  NC_000913/2318990‑2319170

 tgttCGACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATg                                                                                         >  1:1887794/1‑93 (MQ=255)

  gttCGACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGa                                                                                        <  2:4949238/93‑1 (MQ=255)

  gttCGACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACAGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGa                                                                                        >  1:1041400/1‑93 (MQ=255)

    tCGACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGAcc                                                                                      <  2:3615962/93‑1 (MQ=255)

     cGACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCg                                                                                     >  1:1452173/1‑93 (MQ=255)

       aCAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCa                                                                                   >  2:157945/1‑93 (MQ=255)

       aCAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCa                                                                                   >  2:570057/1‑93 (MQ=255)

        cAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCAc                                                                                  >  1:4202225/1‑93 (MQ=255)

        cAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCAc                                                                                  >  1:3742535/1‑93 (MQ=255)

        cAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCAc                                                                                  >  2:279110/1‑93 (MQ=255)

        cAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCAc                                                                                  >  2:3633336/1‑93 (MQ=255)

          aTACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCa                                                                                >  2:1271047/1‑93 (MQ=255)

          aTACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCa                                                                                >  1:1716130/1‑93 (MQ=255)

           tACTCAGGTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCAt                                                                               <  2:2100206/93‑1 (MQ=255)

            aCTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCAtt                                                                              <  1:4443702/93‑1 (MQ=255)

             cTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTg                                                                             <  2:5255045/93‑1 (MQ=255)

                 gTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACt                                                                         >  2:1791268/1‑93 (MQ=255)

                 gTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACt                                                                         <  2:201862/93‑1 (MQ=255)

                 gTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACt                                                                         >  2:3548726/1‑93 (MQ=255)

                  ttCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTc                                                                        <  1:5565774/93‑1 (MQ=255)

                  ttCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTc                                                                        >  1:5145739/1‑93 (MQ=255)

                   tCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCa                                                                       <  2:3425064/93‑1 (MQ=255)

                    cTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAg                                                                      <  2:5306764/93‑1 (MQ=255)

                    cTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAg                                                                      <  2:1033026/93‑1 (MQ=255)

                     tACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAgt                                                                     <  1:3018566/93‑1 (MQ=255)

                       cAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAgtgt                                                                   <  2:4712740/93‑1 (MQ=255)

                        aGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTc                                                                  >  1:3352612/1‑93 (MQ=255)

                         gTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCa                                                                 <  2:1802459/93‑1 (MQ=255)

                           ccAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACt                                                               >  2:5570478/1‑93 (MQ=255)

                            cAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTa                                                              <  1:4147178/93‑1 (MQ=255)

                            cAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTa                                                              <  2:4409494/93‑1 (MQ=255)

                            cAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTa                                                              >  2:1709570/1‑93 (MQ=255)

                            cAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTa                                                              <  1:1138280/93‑1 (MQ=255)

                             aGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTAc                                                             <  2:2371981/93‑1 (MQ=255)

                              gTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCCCTAcc                                                            >  2:2821444/1‑93 (MQ=255)

                              gTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTAcc                                                            >  2:991619/1‑93 (MQ=255)

                               tACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCa                                                           >  2:5778534/1‑93 (MQ=255)

                               tACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCa                                                           >  2:5364791/1‑93 (MQ=255)

                               tACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCa                                                           <  2:5677035/93‑1 (MQ=255)

                                aCTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAg                                                          <  1:4321247/93‑1 (MQ=255)

                                aCTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAg                                                          >  1:4690554/1‑93 (MQ=255)

                                 cTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGt                                                         >  1:1498610/1‑93 (MQ=255)

                                 cTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGt                                                         >  2:3149931/1‑93 (MQ=255)

                                  tGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTc                                                        >  1:5437661/1‑93 (MQ=255)

                                  gggATACGCTGGAACAGGGCACCGAACATGTGGCGCTGGAGATCAGTTTGCCAGGTCATGACCGCACCATTGAACCCAGTGTCACTACCAGTc                                                        >  1:5492564/2‑93 (MQ=255)

                                   ggATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCg                                                       >  2:3892932/1‑93 (MQ=255)

                                   ggATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCg                                                       <  2:5165768/93‑1 (MQ=255)

                                   ggATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCg                                                       >  2:3333848/1‑93 (MQ=255)

                                    gATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGt                                                      >  2:4371535/1‑93 (MQ=255)

                                     aTACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTa                                                     <  2:4009575/93‑1 (MQ=255)

                                     aTACGCTCGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTa                                                     <  2:5514226/93‑1 (MQ=255)

                                      tACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTAt                                                    <  1:1245494/93‑1 (MQ=255)

                                          cTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCAt                                                >  1:4860260/1‑93 (MQ=255)

                                           tGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATa                                               <  1:96385/93‑1 (MQ=255)

                                           tGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATa                                               <  2:894294/93‑1 (MQ=255)

                                           tGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATa                                               >  1:597399/1‑93 (MQ=255)

                                            ggAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATaa                                              <  1:228075/93‑1 (MQ=255)

                                              aaCATTGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAaca                                            <  2:1895368/93‑1 (MQ=255)

                                              aaCATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAaca                                            <  1:2419822/93‑1 (MQ=255)

                                              aaCATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAaca                                            <  2:3669458/93‑1 (MQ=255)

                                              aaCATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAaca                                            >  1:222773/1‑93 (MQ=255)

                                              aaCATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAaca                                            <  1:2931327/93‑1 (MQ=255)

                                               aCATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAacac                                           >  2:6077318/1‑93 (MQ=255)

                                                 aTGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAAcacgc                                         <  1:616223/93‑1 (MQ=255)

                                                    gCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAAcacgcacg                                      <  1:5363628/93‑1 (MQ=255)

                                                    gCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAAcacgcacg                                      <  1:4262188/93‑1 (MQ=255)

                                                      aCCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAACACGCACGGt                                    <  2:3679143/93‑1 (MQ=255)

                                                        cGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAACACGCACGGTGa                                  <  1:3643065/93‑1 (MQ=255)

                                                        cGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAACACGCACGGTGa                                  >  2:3141418/1‑93 (MQ=255)

                                                          aaCATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAACACGCACGGTGaaa                                <  2:3196979/93‑1 (MQ=255)

                                                          aaCATGGGGCGCTGGAGATCAGGTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAACACGCACGGTGaaa                                >  1:3440261/1‑93 (MQ=255)

                                                            cATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAacac                                           >  2:833146/1‑80 (MQ=255)

                                                            cATGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAacac                                           <  1:833146/80‑1 (MQ=255)

                                                             aTGTGGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAACACGCACGGTGAAATGa                             >  1:2688118/1‑93 (MQ=255)

                                                              tgtgGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCTTAACACGCACGGTGAAATGAt                            <  1:6055458/93‑1 (MQ=255)

                                                              tgtgGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAACACGCACGGTGAAATGAt                            >  2:3492920/1‑93 (MQ=255)

                                                               gtgGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAACACGCACGGTGAAATGATa                           <  2:6005945/93‑1 (MQ=255)

                                                                tgGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAACACGCACGGTGAAATGATAg                          <  1:2346522/93‑1 (MQ=255)

                                                                tgGCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAACACGCACGGTGAAATGATAg                          <  1:3931691/93‑1 (MQ=255)

                                                                 tgcgcTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAACACGCACGGTGAAATGATAgg                         <  2:5175901/92‑1 (MQ=255)

                                                                 ggCGCTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAACACGCACGGTGAAATGATAgg                         <  2:4519572/93‑1 (MQ=255)

                                                                    gcTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAACACGCACGGTGAAATGATAGGTGc                      >  2:2318501/1‑93 (MQ=255)

                                                                     cTGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAACACGCACGGTGAAATGATAGGTGCt                     <  2:11879/93‑1 (MQ=255)

                                                                      tGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAACACGCACGGTGAAATGATAGGTGCtt                    >  1:253631/1‑93 (MQ=255)

                                                                      tGGAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAACACGCACGGTGAAATGATAGGTGCtt                    >  2:4160112/1‑93 (MQ=255)

                                                                       ggAGATCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAACACGCACGGTGAAATGATAGGTGCttt                   <  1:775587/93‑1 (MQ=255)

                                                                        gagaTCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAACACGCACGGTGAAATGATAGGTGCTTTg                  >  2:5846761/1‑93 (MQ=255)

                                                                         agaTCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAACACGCACGGTGAAATGATAGGTGCTTTgg                 >  1:1013374/1‑93 (MQ=255)

                                                                          gaTCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAACACGCACGGTGAAATGATAGGTGCTTTGGt                >  2:3870515/1‑93 (MQ=255)

                                                                            tCAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAACACGCACGGTGAAATGATAGGTGCTTTGGTGa              <  2:4035503/93‑1 (MQ=255)

                                                                            tCAATTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAACACGCACGGTGAAATGATAGGTGCTTTGGTGa              <  2:653264/93‑1 (MQ=255)

                                                                             cAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAACACGCACGGTGAAATGATAGGTGCTTTGGTGAt             <  2:4336166/93‑1 (MQ=255)

                                                                             cAGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAACACGCACGGTGAAATGATAGGTGCTTTGGTGAt             <  2:3117272/93‑1 (MQ=255)

                                                                              aGTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAACACGCACGGTGAAATGATAGGTGCTTTGGTGAtt            <  1:5625919/93‑1 (MQ=255)

                                                                               gTTTTCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAACACGCACGGTGAAATGATAGGTGCTTTGGTGAttt           >  2:412696/1‑93 (MQ=255)

                                                                                   tCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAACACGCACGGTGAAATGATAGGTGCTTTGGTGATTTtctc       >  1:5287496/1‑93 (MQ=255)

                                                                                   tCCAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAACACGCACGGTGAAATGATAGGTGCTTTGGTGATTTtctc       <  1:5034458/93‑1 (MQ=255)

                                                                                     cAGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAACACGCACGGTGAAATGATAGGTGCTTTGGTGATTTTCTCTg     >  1:4282919/1‑93 (MQ=255)

                                                                                      aGGTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAACACGCACGGTGAAATGATAGGTGCTTTGGTGATTTTCTCTGa    >  2:2960864/1‑93 (MQ=255)

                                                                                       ggTCATGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAACACGCACGGTGAAATGATAGGTGCTTTGGTGATTTTCTCTGAt   >  2:5003781/1‑93 (MQ=255)

TTGTTCGACAATACTCAGTTCTACAGTCCAGTACTGGATACGCTGGAACATGGCACCGAACATGTGGCGCTGGAGATCAGTTTTCCAGGTCGTGACCGCACCATTGAACTCAGTGTCACTACCAGTCGTATTCATAACACGCACGGTGAAATGATAGGTGCTTTGGTGATTTTCTCTGATT  >  NC_000913/2318990‑2319170

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