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	99.4%	1563.2 / NA	488	V179A (GTG→GCG)	tdcB	threonine dehydratase
Reads supporting (aligned to +/- strand): ref base A (0/1); major base G (246/238); minor base T (1/2); total (247/241)
Fisher's exact test for biased strand distribution p-value = 6.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 518 total aligned reads displayed.

TGTACGCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATA  >  NC_000913/3263044‑3263223

tGTACGCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATcgcg                                                                                         >  1:7567629/1‑93 (MQ=255)

tGTACGCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATcgcg                                                                                         >  2:7626803/1‑93 (MQ=255)

tGTACGCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATcgcg                                                                                         >  1:9339443/1‑93 (MQ=255)

tGTACGCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATcgcg                                                                                         >  2:3343133/1‑93 (MQ=255)

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

    cGCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCTCCGTAATACCAGCAATTAATCCGCCACCACCAATTGGCACAATCGCGTTAt                                                                                     <  1:5038011/93‑1 (MQ=255)

    cGCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTAt                                                                                     <  1:3594162/93‑1 (MQ=255)

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

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

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

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

         aTAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGAca                                                                                >  1:9654835/1‑93 (MQ=255)

          tAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGAcat                                                                               >  2:5803464/1‑93 (MQ=255)

          tAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGAcat                                                                               >  1:4074364/1‑93 (MQ=255)

            acacGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGAcatca                                                                             >  1:8551509/1‑93 (MQ=255)

              atgAATGGTCGGTGTAATAAATTTAATCGCCTCCGCAATACCAGCAATTAAACCGCCACCACCAAGTGGCACAATCGCGTTATCGACATCATa                                                                           <  2:2330296/91‑1 (MQ=255)

              acGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATa                                                                           >  2:599102/1‑93 (MQ=255)

               cGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATag                                                                          <  1:2538283/93‑1 (MQ=255)

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

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

                 aaTGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATagag                                                                        <  1:3000495/93‑1 (MQ=255)

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

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

                    ggTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGATATCGACATCATAGAGATc                                                                     >  1:4465512/1‑93 (MQ=255)

                      tCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCATTATCGACATCATAGAGATCtt                                                                   >  1:1108774/1‑93 (MQ=255)

                        gggTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCCATTGGCACAATCGCGTTATCGACATCATAGAGATCTTcc                                                                 >  2:3642633/1‑93 (MQ=255)

                         ggTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCa                                                                <  1:8265723/93‑1 (MQ=255)

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

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

                           ttAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATa                                                              >  1:6199233/1‑93 (MQ=255)

                           ttAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATa                                                              >  1:6769109/1‑93 (MQ=255)

                           ttAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATa                                                              >  2:4235449/1‑93 (MQ=255)

                           ttAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATa                                                              <  2:9415107/93‑1 (MQ=255)

                            tAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATaa                                                             >  2:8934892/1‑93 (MQ=255)

                            tAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATaa                                                             >  1:2302775/1‑93 (MQ=255)

                            tAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATaa                                                             >  2:7862224/1‑93 (MQ=255)

                             aaTAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAAt                                                            >  2:2695419/1‑93 (MQ=255)

                             aaTAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAAt                                                            >  2:8088838/1‑93 (MQ=255)

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

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

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

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

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

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

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

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

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

                                   tttAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAg                                                      <  1:206333/93‑1 (MQ=255)

                                       aTTGCCACCGCAATACCAGCAATTAACCCGCCACCGCCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCa                                                  <  1:7763162/93‑1 (MQ=255)

                                       aTTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCa                                                  >  1:8328541/1‑93 (MQ=255)

                                       aTTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCa                                                  >  1:2681235/1‑93 (MQ=255)

                                        ttGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCaa                                                 >  2:6235974/1‑93 (MQ=255)

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

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

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

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

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

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

                                              ccGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTc                                           <  1:4126340/93‑1 (MQ=255)

                                              ccGCAATACCAGCAAATAAACCGCAACCGCCAAATGGCACATTCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGtta                                           <  1:8975835/93‑2 (MQ=255)

                                               cGCAATACCAGCAATTACACCGCCACCAACAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTcc                                          <  1:1712260/93‑1 (MQ=255)

                                               cGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTcc                                          >  2:2205326/1‑93 (MQ=255)

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

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

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

                                                 cAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCt                                        >  1:5325798/1‑93 (MQ=255)

                                                      ccAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCa                                   <  1:5512888/93‑1 (MQ=255)

                                                       cAGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAg                                  <  2:3785300/93‑1 (MQ=255)

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

                                                        aGCAATTAAGCCGCCCCCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGc                                 <  2:631991/93‑1 (MQ=255)

                                                        aGCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGc                                 <  1:5882349/93‑1 (MQ=255)

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

                                                         gCAATTAACCCGCCACCACCAA‑TGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCaa                               <  2:6464885/93‑1 (MQ=255)

                                                         gCAATTAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCa                                >  2:5262333/1‑93 (MQ=255)

                                                           aaTTAAACCGCCACCACCAATTGGCACAATGGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAAt                              <  2:9821378/93‑1 (MQ=255)

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

                                                             ttAAACCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTGCCAGGCCAATCGTTCCCTGGCCAGCAATCa                            <  2:4727308/93‑1 (MQ=255)

                                                               aaaCCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACt                          <  1:9450969/93‑1 (MQ=255)

                                                               aaaCCGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACt                          <  2:6743049/93‑1 (MQ=255)

                                                                  ccGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTc                       <  1:9121954/93‑1 (MQ=255)

                                                                  ccGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTc                       <  2:8449926/93‑1 (MQ=255)

                                                                  ccGCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTc                       >  2:7292100/1‑93 (MQ=255)

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

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

                                                                    gCCACCACCAATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCgg                     <  2:1266089/93‑1 (MQ=255)

                                                                     ccaccaccaATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGa                    <  1:3914144/93‑1 (MQ=255)

                                                                     ccaccaccaATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGa                    >  2:4072962/1‑93 (MQ=255)

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

                                                                         caccaATTGGCACAATCGCGTTATCGACGTCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGatcat                >  1:9601059/1‑93 (MQ=255)

                                                                         caccaATTGGCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGatcat                <  1:3214194/93‑1 (MQ=255)

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

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

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

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

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

                                                                               agggCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGCCCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAg          <  2:9144043/91‑1 (MQ=255)

                                                                                 ggCACAATCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAggg        >  1:8212114/1‑92 (MQ=255)

                                                                                      aaTCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGAt   <  2:1345101/93‑1 (MQ=255)

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

                                                                                      aaTCGCGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGAt   <  1:3104445/93‑1 (MQ=255)

TGTACGCCAATAACACGAATGGTCGGGTTAATAGATTTAATTGCCACCGCAATACCAGCAATTAAACCGCCACCACCAATTGGCACAATCACGTTATCGACATCATAGAGATCTTCCATAATTTCCAGACCAATCGTTCCCTGGCCAGCAATCACTTTCGGATCATCGTAAGGTGGGATA  >  NC_000913/3263044‑3263223

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